JP2019042497A - Medical information processing device, medical image diagnostic device and program - Google Patents

Medical information processing device, medical image diagnostic device and program Download PDF

Info

Publication number
JP2019042497A
JP2019042497A JP2018142159A JP2018142159A JP2019042497A JP 2019042497 A JP2019042497 A JP 2019042497A JP 2018142159 A JP2018142159 A JP 2018142159A JP 2018142159 A JP2018142159 A JP 2018142159A JP 2019042497 A JP2019042497 A JP 2019042497A
Authority
JP
Japan
Prior art keywords
electric signal
signal propagation
positions
myocardium
motion information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2018142159A
Other languages
Japanese (ja)
Other versions
JP7144230B2 (en
Inventor
昂彦 西岡
Takahiko Nishioka
昂彦 西岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to US16/119,527 priority Critical patent/US11484362B2/en
Publication of JP2019042497A publication Critical patent/JP2019042497A/en
Application granted granted Critical
Publication of JP7144230B2 publication Critical patent/JP7144230B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

To calculate propagation of an electric signal while reducing a patient burden.SOLUTION: According to an embodiment, a medical information processing device includes a storage part, an acquisition part, and a calculation part. The storage part stores a model which shows association between a cellular deformation and electric signal propagation. The acquisition part acquires four-dimensional image data obtained by imaging a three-dimensional structure of a heart muscle in time series, and extracts motion information of the heart muscle at a plurality of positions of the acquired four-dimensional image data. The calculation part calculates electric signal propagation at the plurality of positions by applying the model to the motion information at the plurality of positions.SELECTED DRAWING: Figure 1

Description

本発明の実施形態は、医用情報処理装置、医用画像診断装置及びプログラムに関する。   Embodiments described herein relate generally to a medical information processing apparatus, a medical image diagnostic apparatus, and a program.

心臓は、伝搬する電気信号に応じて収縮(興奮)することによって、心腔内の血液を全身に送り出す臓器である。そのため、電気信号の伝搬における異常は、頻脈等の不整脈の原因となる。そして、不整脈の治療をする場合、医師は、電気信号の伝搬に異常があるか否か、異常の原因が心臓のどの部分にあるのか、といった情報に基づいて、治療計画を立てたり手術を行なったりする。   The heart is an organ that pumps blood in the heart chamber to the whole body by contracting (exciting) in response to a propagating electrical signal. Therefore, abnormalities in the propagation of electrical signals cause arrhythmias such as tachycardia. When treating arrhythmias, doctors make treatment plans and perform surgery based on information such as whether there is an abnormality in the propagation of electrical signals and in which part of the heart the abnormality is caused. Or

ここで、従来、電気信号の伝搬の評価方法としては、心腔内に電極カテーテルを挿入し、心腔内の電位変化を記録する方法が知られている。更に、心腔内をペーシング刺激し、頻脈を誘発することで、不整脈の診断や、カテーテルアブレーションにおける対象部位の設定等が可能である。しかしながら、かかる評価方法は、心腔内に電極カテーテルを挿入したりペーシング刺激によって頻脈を誘発したりする点で、患者負担が大きかった。   Here, conventionally, as a method for evaluating propagation of an electric signal, a method is known in which an electrode catheter is inserted into a heart chamber and a potential change in the heart chamber is recorded. Further, by pacing the inside of the heart chamber and inducing tachycardia, it is possible to diagnose arrhythmia, set a target site in catheter ablation, and the like. However, such an evaluation method has a large patient burden in that an electrode catheter is inserted into the heart chamber and tachycardia is induced by pacing stimulation.

電気信号の伝搬の評価方法のうち、患者負担の小さいものとしては、例えば、心電図が挙げられる。しかしながら、心電図では、電気信号の伝搬に異常があるか否かを判断することはできるとしても、異常部位の特定には至らなかった。即ち、心電図では、不整脈を診断したりカテーテルアブレーションの対象部位を設定したりする上で、十分とはいえない場合があった。   Among the methods for evaluating the propagation of electrical signals, an example of a method with a small patient burden is an electrocardiogram. However, in the electrocardiogram, even if it can be determined whether or not there is an abnormality in the propagation of the electric signal, the abnormal part has not been specified. That is, an electrocardiogram may not be sufficient in diagnosing arrhythmia or setting a target site for catheter ablation.

特開2014−512201号公報JP 2014-512201 A 特開2011−212043号公報JP 2011-212043 A 特開平08−289877号公報Japanese Patent Laid-Open No. 08-289877

本発明が解決しようとする課題は、患者負担を低減しつつ、電気信号の伝搬を算出することである。   The problem to be solved by the present invention is to calculate the propagation of electrical signals while reducing patient burden.

実施形態の医用情報処理装置は、記憶部と、取得部と、算出部とを備える。記憶部は、細胞変形と電気信号伝搬との関連を示すモデルを記憶する。取得部は、心筋の3次元構造を時系列で撮像することにより得られた4次元画像データを取得し、取得した4次元画像データの複数位置における前記心筋の動き情報を抽出する。算出部は、前記モデルを前記複数位置における前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する。   The medical information processing apparatus according to the embodiment includes a storage unit, an acquisition unit, and a calculation unit. The storage unit stores a model indicating the relationship between cell deformation and electric signal propagation. The acquisition unit acquires four-dimensional image data obtained by imaging the three-dimensional structure of the myocardium in time series, and extracts the motion information of the myocardium at a plurality of positions of the acquired four-dimensional image data. The calculation unit calculates electric signal propagation at the plurality of positions by applying the model to the motion information at the plurality of positions.

図1は、第1の実施形態に係る医用情報処理システムの構成の一例を示すブロック図である。FIG. 1 is a block diagram illustrating an example of a configuration of a medical information processing system according to the first embodiment. 図2は、第1の実施形態に係るX線CT装置の構成の一例を示すブロック図である。FIG. 2 is a block diagram showing an example of the configuration of the X-ray CT apparatus according to the first embodiment. 図3は、第1の実施形態に係る4次元画像データの一例を示す図である。FIG. 3 is a diagram illustrating an example of four-dimensional image data according to the first embodiment. 図4は、第1の実施形態に係る電気信号の伝搬と筋細胞との動態シミュレーションの一例を示す図である。FIG. 4 is a diagram illustrating an example of electrical signal propagation and muscle cell dynamic simulation according to the first embodiment. 図5は、第1の実施形態に係る医用情報処理装置の処理の一連の流れを説明するためのフローチャートである。FIG. 5 is a flowchart for explaining a series of processing steps of the medical information processing apparatus according to the first embodiment.

以下、図面を参照して、医用情報処理装置、医用画像診断装置及びプログラムの実施形態について詳細に説明する。   Hereinafter, embodiments of a medical information processing apparatus, a medical image diagnostic apparatus, and a program will be described in detail with reference to the drawings.

(第1の実施形態)
まず、第1の実施形態について説明する。第1の実施形態では、医用情報処理装置及び医用画像診断装置を含んだ医用情報処理システムを一例として説明する。また、第1の実施形態では、被検体P1の心筋について、電気信号の伝搬を算出する場合を一例として説明する。
(First embodiment)
First, the first embodiment will be described. In the first embodiment, a medical information processing system including a medical information processing apparatus and a medical image diagnostic apparatus will be described as an example. In the first embodiment, a case where the propagation of an electric signal is calculated for the myocardium of the subject P1 will be described as an example.

図1に示すように、第1の実施形態に係る医用情報処理システム1は、医用画像診断装置10と、画像保管装置20と、医用情報処理装置30とを備える。なお、図1は、第1の実施形態に係る医用情報処理システム1の構成の一例を示すブロック図である。図1に示すように、医用画像診断装置10、画像保管装置20及び医用情報処理装置30は、ネットワークを介して相互に接続される。   As shown in FIG. 1, the medical information processing system 1 according to the first embodiment includes a medical image diagnostic device 10, an image storage device 20, and a medical information processing device 30. FIG. 1 is a block diagram illustrating an example of the configuration of the medical information processing system 1 according to the first embodiment. As shown in FIG. 1, the medical image diagnostic apparatus 10, the image storage apparatus 20, and the medical information processing apparatus 30 are connected to each other via a network.

医用画像診断装置10は、被検体P1から医用画像データを収集する装置である。例えば、医用画像診断装置10は、被検体P1の心筋について、4次元画像データ(時系列の3次元画像データ)を収集する。即ち、医用画像診断装置10は、被検体P1の心筋の3次元構造を時系列で撮像することにより、4次元画像データを収集する。例えば、医用画像診断装置10は、X線CT(Computed Tomography)装置やMRI(Magnetic Resonance Imaging)装置、超音波診断装置等である。   The medical image diagnostic apparatus 10 is an apparatus that collects medical image data from the subject P1. For example, the medical image diagnostic apparatus 10 collects four-dimensional image data (time-series three-dimensional image data) for the myocardium of the subject P1. That is, the medical image diagnostic apparatus 10 collects four-dimensional image data by imaging the three-dimensional structure of the myocardium of the subject P1 in time series. For example, the medical image diagnostic apparatus 10 is an X-ray CT (Computed Tomography) apparatus, an MRI (Magnetic Resonance Imaging) apparatus, an ultrasonic diagnostic apparatus, or the like.

画像保管装置20は、医用画像診断装置10によって収集された医用画像データを保管する装置である。例えば、画像保管装置20は、サーバ装置等のコンピュータ機器によって実現される。本実施形態では、画像保管装置20は、ネットワークを介して医用画像診断装置10から4次元画像データを取得し、取得した4次元画像データを、装置内又は装置外に設けられたメモリに記憶させる。   The image storage device 20 is a device that stores medical image data collected by the medical image diagnostic device 10. For example, the image storage device 20 is realized by a computer device such as a server device. In the present embodiment, the image storage device 20 acquires four-dimensional image data from the medical image diagnostic device 10 via a network, and stores the acquired four-dimensional image data in a memory provided inside or outside the device. .

医用情報処理装置30は、ネットワークを介して4次元画像データを取得し、取得した4次元画像データを用いた種々の処理を実行する。例えば、医用情報処理装置30は、ワークステーション等のコンピュータ機器によって実現される。本実施形態では、医用情報処理装置30は、ネットワークを介して、医用画像診断装置10又は画像保管装置20から4次元画像データを取得する。また、医用情報処理装置30は、取得した4次元画像データに基づいて被検体P1の心筋の動き情報を算出し、動き情報に対応した心筋の電気信号の伝搬(興奮伝搬)を算出する。   The medical information processing apparatus 30 acquires 4D image data via a network, and executes various processes using the acquired 4D image data. For example, the medical information processing apparatus 30 is realized by a computer device such as a workstation. In this embodiment, the medical information processing apparatus 30 acquires four-dimensional image data from the medical image diagnostic apparatus 10 or the image storage apparatus 20 via a network. The medical information processing apparatus 30 calculates myocardial motion information of the subject P1 based on the acquired four-dimensional image data, and calculates myocardial electrical signal propagation (excitation propagation) corresponding to the motion information.

図1に示すように、医用情報処理装置30は、入力インターフェース31と、ディスプレイ32と、メモリ33と、処理回路34とを有する。   As shown in FIG. 1, the medical information processing apparatus 30 includes an input interface 31, a display 32, a memory 33, and a processing circuit 34.

入力インターフェース31は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路34に出力する。例えば、入力インターフェース31は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、タッチパネル等により実現される。   The input interface 31 receives various input operations from the operator, converts the received input operations into electric signals, and outputs them to the processing circuit 34. For example, the input interface 31 is realized by a mouse, keyboard, trackball, switch, button, joystick, touch panel, or the like.

ディスプレイ32は、操作者によって参照されるモニタであり、処理回路34による制御の下、医用画像を動画像又は静止画像として表示したり、処理回路34が算出した電気信号の伝搬を表示したり、入力インターフェース31を介して操作者から各種指示や各種設定等を受け付けるためのGUI(Graphical User Interface)を表示したりする。例えば、ディスプレイ32は、液晶ディスプレイやCRT(Cathode Ray Tube)ディスプレイである。   The display 32 is a monitor that is referred to by an operator, displays a medical image as a moving image or a still image under the control of the processing circuit 34, displays the propagation of an electric signal calculated by the processing circuit 34, A GUI (Graphical User Interface) for receiving various instructions and various settings from the operator via the input interface 31 is displayed. For example, the display 32 is a liquid crystal display or a CRT (Cathode Ray Tube) display.

メモリ33は、例えば、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。例えば、メモリ33は、医用画像診断装置10又は画像保管装置20から取得した4次元画像データを記憶する。また、メモリ33は、細胞変形と電気信号伝搬との関連を示すモデルを記憶する。なお、メモリ33が記憶するモデルについては後述する。また、例えば、メモリ33は、医用情報処理装置30に含まれる回路がその機能を実現するためのプログラムを記憶する。   The memory 33 is realized by, for example, a RAM (Random Access Memory), a semiconductor memory element such as a flash memory, a hard disk, an optical disk, or the like. For example, the memory 33 stores four-dimensional image data acquired from the medical image diagnostic apparatus 10 or the image storage apparatus 20. The memory 33 stores a model indicating the relationship between cell deformation and electric signal propagation. The model stored in the memory 33 will be described later. In addition, for example, the memory 33 stores a program for a circuit included in the medical information processing apparatus 30 to realize its function.

処理回路34は、取得機能34a、算出機能34b及び表示制御機能34cを実行することで、医用情報処理装置30全体の動作を制御する。例えば、処理回路34は、取得機能34aに対応するプログラムをメモリ33から読み出して実行することにより、医用画像診断装置10又は画像保管装置20から4次元画像データを取得する。即ち、取得機能34aは、心筋の3次元構造を時系列で撮像することにより得られた4次元画像データを取得する。また、取得機能34aは、取得した4次元画像データから被検体P1の心筋の動き情報を算出する。具体的には、取得機能34aは、4次元画像データの複数位置における心筋の動き情報を抽出する。また、例えば、処理回路34は、算出機能34bに対応するプログラムをメモリ33から読み出して実行することにより、動き情報に対応した被検体P1の心筋の電気信号伝搬を算出する。具体的には、算出機能34bは、細胞変形と電気信号伝搬との関連を示すモデルを、4次元画像データの複数位置における動き情報に適用することにより、4次元画像データの複数位置における電気信号伝搬を算出する。また、例えば、処理回路34は、表示制御機能34cに対応するプログラムをメモリ33から読み出して実行することにより、被検体P1の心筋の電気信号伝搬をディスプレイ32において表示する。即ち、表示制御機能34cは、4次元画像データの複数位置における電気信号伝搬を表示する。なお、電気信号伝搬の算出及び表示については後述する。   The processing circuit 34 controls the overall operation of the medical information processing apparatus 30 by executing an acquisition function 34a, a calculation function 34b, and a display control function 34c. For example, the processing circuit 34 acquires four-dimensional image data from the medical image diagnostic apparatus 10 or the image storage apparatus 20 by reading a program corresponding to the acquisition function 34 a from the memory 33 and executing the program. That is, the acquisition function 34a acquires four-dimensional image data obtained by imaging the three-dimensional structure of the myocardium in time series. Further, the acquisition function 34a calculates myocardial motion information of the subject P1 from the acquired four-dimensional image data. Specifically, the acquisition function 34a extracts myocardial motion information at a plurality of positions of the four-dimensional image data. Further, for example, the processing circuit 34 reads the program corresponding to the calculation function 34b from the memory 33 and executes it, thereby calculating the myocardial electrical signal propagation of the subject P1 corresponding to the motion information. Specifically, the calculation function 34b applies a model indicating the relationship between cell deformation and electric signal propagation to motion information at a plurality of positions in the four-dimensional image data, thereby causing electric signals at a plurality of positions in the four-dimensional image data. Calculate the propagation. In addition, for example, the processing circuit 34 displays the electrical signal propagation of the myocardium of the subject P1 on the display 32 by reading a program corresponding to the display control function 34c from the memory 33 and executing it. That is, the display control function 34c displays electric signal propagation at a plurality of positions of the four-dimensional image data. The calculation and display of electric signal propagation will be described later.

図1に示す医用情報処理装置30においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ33へ記憶されている。処理回路34は、メモリ33からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路34は、読み出したプログラムに対応する機能を有することとなる。なお、図1においては単一の処理回路34にて、取得機能34a、算出機能34b及び表示制御機能34cが実現するものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路34を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。   In the medical information processing apparatus 30 shown in FIG. 1, each processing function is stored in the memory 33 in the form of a program that can be executed by a computer. The processing circuit 34 is a processor that realizes a function corresponding to each program by reading the program from the memory 33 and executing the program. In other words, the processing circuit 34 that has read each program has a function corresponding to the read program. In FIG. 1, it has been described that the acquisition function 34a, the calculation function 34b, and the display control function 34c are realized by a single processing circuit 34. However, the processing circuit 34 is configured by combining a plurality of independent processors. Each processor may implement a function by executing a program.

次に、4次元画像データの収集を行なう医用画像診断装置10について、図2を用いて具体的に説明する。図2においては、医用画像診断装置10の一例として、X線CT装置100について説明する。図2は、第1の実施形態に係るX線CT装置100の構成の一例を示すブロック図である。   Next, the medical image diagnostic apparatus 10 that collects four-dimensional image data will be specifically described with reference to FIG. In FIG. 2, an X-ray CT apparatus 100 will be described as an example of the medical image diagnostic apparatus 10. FIG. 2 is a block diagram showing an example of the configuration of the X-ray CT apparatus 100 according to the first embodiment.

X線CT装置100は、図2に示すように、架台装置110と、寝台装置130と、コンソール装置140とを有する。また、X線CT装置100は、図示しない心電計が接続される。心電計は、被検体P1に付された電極を介して被検体P1のECG(electro cardiogram)信号を検出して心電図を生成し、生成した心電図をX線CT装置100に送信する。なお、図2においては、非チルト状態での回転フレーム113の回転軸又は寝台装置130の天板133の長手方向をZ軸方向とする。また、Z軸方向に直交し、床面に対し水平である軸方向をX軸方向とする。また、Z軸方向に直交し、床面に対し垂直である軸方向をY軸方向とする。   As shown in FIG. 2, the X-ray CT apparatus 100 includes a gantry device 110, a couch device 130, and a console device 140. The X-ray CT apparatus 100 is connected to an electrocardiograph (not shown). The electrocardiograph detects an ECG (electro cardiogram) signal of the subject P1 via an electrode attached to the subject P1, generates an electrocardiogram, and transmits the generated electrocardiogram to the X-ray CT apparatus 100. In FIG. 2, the rotation axis of the rotating frame 113 or the longitudinal direction of the top plate 133 of the bed apparatus 130 in the non-tilt state is defined as the Z-axis direction. Further, an axial direction that is orthogonal to the Z-axis direction and is horizontal to the floor surface is taken as an X-axis direction. Further, an axial direction orthogonal to the Z-axis direction and perpendicular to the floor surface is defined as a Y-axis direction.

架台装置110は、X線管111と、X線検出器112と、回転フレーム113と、X線高電圧装置114と、制御装置115と、ウェッジ116と、コリメータ117と、データ収集回路118とを有する。   The gantry device 110 includes an X-ray tube 111, an X-ray detector 112, a rotating frame 113, an X-ray high voltage device 114, a control device 115, a wedge 116, a collimator 117, and a data acquisition circuit 118. Have.

X線管111は、熱電子を発生する陰極(フィラメント)と、熱電子の衝突を受けてX線を発生する陽極(ターゲット)とを有する真空管である。X線管111は、X線高電圧装置114から供給される高電圧を用いて、陰極から陽極に向けて熱電子を照射することにより、X線を発生する。   The X-ray tube 111 is a vacuum tube having a cathode (filament) that generates thermoelectrons and an anode (target) that generates X-rays in response to collision of the thermoelectrons. The X-ray tube 111 generates X-rays by irradiating thermionic electrons from the cathode to the anode using the high voltage supplied from the X-ray high voltage device 114.

X線検出器112は、X線管111から照射されて被検体P1を通過したX線を検出し、検出したX線量に対応した信号をデータ収集回路118へと出力する。X線検出器112は、例えば、X線管111の焦点を中心とした1つの円弧に沿ってチャネル方向に複数のX線検出素子が配列された複数のX線検出素子列を有する。X線検出器112は、例えば、チャネル方向に複数のX線検出素子が配列されたX線検出素子列がスライス方向(列方向、row方向)に複数配列された構造を有する。また、X線検出器112は、例えば、グリッドと、シンチレータアレイと、光センサアレイとを有する間接変換型の検出器である。シンチレータアレイは、複数のシンチレータを有する。シンチレータは入射X線量に応じた光子量の光を出力するシンチレータ結晶を有する。グリッドは、シンチレータアレイのX線入射側の面に配置され、散乱X線を吸収するX線遮蔽板を有する。光センサアレイは、シンチレータからの光量に応じた電気信号に変換する機能を有し、例えば、光電子増倍管(フォトマルチプライヤー:PMT)等の光センサを有する。なお、X線検出器112は、入射したX線を電気信号に変換する半導体素子を有する直接変換型の検出器であっても構わない。   The X-ray detector 112 detects X-rays irradiated from the X-ray tube 111 and passed through the subject P1, and outputs a signal corresponding to the detected X-ray dose to the data collection circuit 118. The X-ray detector 112 has, for example, a plurality of X-ray detection element arrays in which a plurality of X-ray detection elements are arranged in the channel direction along one arc centered on the focal point of the X-ray tube 111. The X-ray detector 112 has, for example, a structure in which a plurality of X-ray detection element arrays in which a plurality of X-ray detection elements are arrayed in the channel direction are arrayed in the slice direction (column direction, row direction). The X-ray detector 112 is an indirect conversion type detector having a grid, a scintillator array, and an optical sensor array, for example. The scintillator array has a plurality of scintillators. The scintillator has a scintillator crystal that outputs a photon amount of light corresponding to the incident X-ray dose. The grid is disposed on the surface on the X-ray incident side of the scintillator array and has an X-ray shielding plate that absorbs scattered X-rays. The optical sensor array has a function of converting into an electric signal corresponding to the amount of light from the scintillator, and includes an optical sensor such as a photomultiplier tube (photomultiplier: PMT). The X-ray detector 112 may be a direct conversion type detector having a semiconductor element that converts incident X-rays into electrical signals.

回転フレーム113は、X線管111とX線検出器112とを対向支持し、制御装置115によってX線管111とX線検出器112とを回転させる円環状のフレームである。例えば、回転フレーム113は、アルミニウムを材料とした鋳物である。なお、回転フレーム113は、X線管111及びX線検出器112に加えて、X線高電圧装置114やデータ収集回路118を更に支持することもできる。更に、回転フレーム113は、図2において図示しない種々の構成を更に支持することもできる。以下では、架台装置110において、回転フレーム113と共に回転移動する部分及び回転フレーム113を回転部とも記載する。   The rotating frame 113 is an annular frame that supports the X-ray tube 111 and the X-ray detector 112 so as to face each other and rotates the X-ray tube 111 and the X-ray detector 112 by the control device 115. For example, the rotating frame 113 is a casting made of aluminum. The rotating frame 113 can further support an X-ray high voltage device 114 and a data acquisition circuit 118 in addition to the X-ray tube 111 and the X-ray detector 112. Further, the rotating frame 113 can further support various configurations not shown in FIG. Hereinafter, in the gantry device 110, the part that rotates together with the rotating frame 113 and the rotating frame 113 are also referred to as a rotating unit.

なお、データ収集回路118が生成した検出データは、回転フレーム113に設けられた発光ダイオード(Light Emitting Diode: LED)を有する送信機から、光通信によって、架台装置110の非回転部分に設けられた、フォトダイオードを有する受信機に送信され、コンソール装置140へと転送される。ここで、非回転部分とは、例えば、回転フレーム113を回転可能に支持する固定フレーム等である。なお、回転フレーム113から架台装置110の非回転部分への検出データの送信方法は、光通信に限らず、非接触型のデータ伝送であれば如何なる方式を採用しても構わない。   The detection data generated by the data collection circuit 118 is provided in a non-rotating portion of the gantry device 110 by optical communication from a transmitter having a light emitting diode (LED) provided in the rotating frame 113. , Transmitted to a receiver having a photodiode, and transferred to the console device 140. Here, the non-rotating part is, for example, a fixed frame that rotatably supports the rotating frame 113. The detection data transmission method from the rotating frame 113 to the non-rotating portion of the gantry device 110 is not limited to optical communication, and any method may be adopted as long as it is non-contact type data transmission.

X線高電圧装置114は、変圧器(トランス)及び整流器等の電気回路を有し、X線管111に印加する高電圧を発生する高電圧発生装置と、X線管111が照射するX線に応じた出力電圧の制御を行うX線制御装置とを有する。高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であってもよい。なお、X線高電圧装置114は、回転フレーム113に設けられてもよいし、図示しない固定フレームに設けられても構わない。   The X-ray high-voltage device 114 has an electric circuit such as a transformer and a rectifier, and generates a high voltage to be applied to the X-ray tube 111 and X-rays irradiated by the X-ray tube 111. And an X-ray control device that controls the output voltage according to the above. The high voltage generator may be a transformer system or an inverter system. Note that the X-ray high voltage device 114 may be provided on the rotating frame 113 or may be provided on a fixed frame (not shown).

制御装置115は、モータ及びアクチュエータ等の駆動機構と、この機構を制御する回路とを含む。制御装置115は、入力インターフェース143や架台装置110に設けられた入力インターフェース等からの入力信号を受けて、架台装置110及び寝台装置130の動作制御を行う。例えば、制御装置115は、回転フレーム113の回転や架台装置110のチルト、寝台装置130及び天板133の動作等について制御を行う。一例を挙げると、制御装置115は、架台装置110をチルトさせる制御として、入力された傾斜角度(チルト角度)情報により、X軸方向に平行な軸を中心に回転フレーム113を回転させる。なお、制御装置115は架台装置110に設けられてもよいし、コンソール装置140に設けられてもよい。   The control device 115 includes a drive mechanism such as a motor and an actuator, and a circuit that controls the mechanism. The control device 115 receives an input signal from an input interface 143 or an input interface provided in the gantry device 110 and performs operation control of the gantry device 110 and the couch device 130. For example, the control device 115 controls the rotation of the rotating frame 113, the tilt of the gantry device 110, the operations of the bed device 130 and the top board 133, and the like. For example, as a control for tilting the gantry device 110, the control device 115 rotates the rotating frame 113 around an axis parallel to the X-axis direction based on the input tilt angle (tilt angle) information. The control device 115 may be provided in the gantry device 110 or may be provided in the console device 140.

ウェッジ116は、X線管111から照射されたX線量を調節するためのフィルタである。具体的には、ウェッジ116は、X線管111から被検体P1へ照射されるX線が、予め定められた分布になるように、X線管111から照射されたX線を透過して減衰するフィルタである。例えば、ウェッジ116は、ウェッジフィルタ(wedge filter)やボウタイフィルタ(bow-tie filter)であり、所定のターゲット角度や所定の厚みとなるようにアルミニウム等を加工して構成される。   The wedge 116 is a filter for adjusting the X-ray dose irradiated from the X-ray tube 111. Specifically, the wedge 116 transmits and attenuates the X-rays irradiated from the X-ray tube 111 so that the X-rays irradiated from the X-ray tube 111 to the subject P1 have a predetermined distribution. It is a filter to do. For example, the wedge 116 is a wedge filter or a bow-tie filter, and is formed by processing aluminum or the like so as to have a predetermined target angle or a predetermined thickness.

コリメータ117は、ウェッジ116を透過したX線の照射範囲を絞り込むための鉛板等であり、複数の鉛板等の組み合わせによってスリットを形成する。コリメータ117は、図示しないコリメータ調整回路によって、開口度及び位置が調整される。これにより、X線管111が発生させたX線の照射範囲が調整される。   The collimator 117 is a lead plate or the like for narrowing the irradiation range of the X-rays transmitted through the wedge 116, and forms a slit by a combination of a plurality of lead plates or the like. The collimator 117 has its opening degree and position adjusted by a collimator adjustment circuit (not shown). Thereby, the irradiation range of the X-rays generated by the X-ray tube 111 is adjusted.

データ収集回路118は、DAS(Data Acquisition System)である。データ収集回路118は、X線検出器112の各X線検出素子から出力される電気信号に対して増幅処理を行う増幅器と、電気信号をデジタル信号に変換するA/D変換器とを有し、検出データを生成する。データ収集回路118は、例えば、プロセッサにより実現される。   The data acquisition circuit 118 is a DAS (Data Acquisition System). The data acquisition circuit 118 includes an amplifier that performs amplification processing on the electric signal output from each X-ray detection element of the X-ray detector 112, and an A / D converter that converts the electric signal into a digital signal. , Generate detection data. The data collection circuit 118 is realized by a processor, for example.

寝台装置130は、スキャン対象の被検体P1を載置、移動させる装置であり、基台131と、寝台駆動装置132と、天板133と、支持フレーム134とを有する。基台131は、支持フレーム134を鉛直方向に移動可能に支持する筐体である。寝台駆動装置132は、被検体P1が載置された天板133を、天板133の長軸方向に移動する駆動機構であり、モータ及びアクチュエータ等を含む。支持フレーム134の上面に設けられた天板133は、被検体P1が載置される板である。なお、寝台駆動装置132は、天板133に加え、支持フレーム134を天板133の長軸方向に移動してもよい。   The couch device 130 is a device for placing and moving the subject P1 to be scanned, and includes a base 131, a couch driving device 132, a top board 133, and a support frame 134. The base 131 is a housing that supports the support frame 134 so as to be movable in the vertical direction. The bed driving device 132 is a driving mechanism that moves the top plate 133 on which the subject P1 is placed in the long axis direction of the top plate 133, and includes a motor, an actuator, and the like. A top plate 133 provided on the upper surface of the support frame 134 is a plate on which the subject P1 is placed. The couch driving device 132 may move the support frame 134 in the long axis direction of the top plate 133 in addition to the top plate 133.

コンソール装置140は、メモリ141と、ディスプレイ142と、入力インターフェース143と、処理回路144とを有する。   The console device 140 includes a memory 141, a display 142, an input interface 143, and a processing circuit 144.

メモリ141は、例えば、RAM、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。例えば、メモリ141は、投影データや再構成画像データ、心電計から送信された心電図等を記憶する。また、例えば、メモリ141は、X線CT装置100に含まれる回路がその機能を実現するためのプログラムを記憶する。   The memory 141 is realized by, for example, a semiconductor memory element such as a RAM or a flash memory, a hard disk, an optical disk, or the like. For example, the memory 141 stores projection data, reconstructed image data, an electrocardiogram transmitted from an electrocardiograph, and the like. Further, for example, the memory 141 stores a program for a circuit included in the X-ray CT apparatus 100 to realize its function.

ディスプレイ142は、各種の情報を表示する。例えば、ディスプレイ142は、処理回路144によって生成された4次元画像データを動画像又は静止画像として表示したり、操作者からの各種操作を受け付けるためのGUI等を表示したりする。例えば、ディスプレイ142は、液晶ディスプレイやCRTディスプレイである。   The display 142 displays various information. For example, the display 142 displays the four-dimensional image data generated by the processing circuit 144 as a moving image or a still image, or displays a GUI or the like for receiving various operations from the operator. For example, the display 142 is a liquid crystal display or a CRT display.

入力インターフェース143は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路144に出力する。例えば、入力インターフェース143は、投影データを収集する際の収集条件や、CT画像データを再構成する際の再構成条件、CT画像データから4次元画像データを生成する際の画像処理条件等を操作者から受け付ける。例えば、入力インターフェース143は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、タッチパネル等により実現される。   The input interface 143 receives various input operations from the operator, converts the received input operations into electrical signals, and outputs them to the processing circuit 144. For example, the input interface 143 operates a collection condition when collecting projection data, a reconstruction condition when reconstructing CT image data, an image processing condition when generating four-dimensional image data from CT image data, and the like. Accept from the person. For example, the input interface 143 is realized by a mouse, a keyboard, a trackball, a switch, a button, a joystick, a touch panel, or the like.

処理回路144は、X線CT装置100全体の動作を制御する。例えば、処理回路144は、収集機能144a、表示制御機能144b及び制御機能144cを有する。処理回路144は、例えば、プロセッサにより実現される。   The processing circuit 144 controls the overall operation of the X-ray CT apparatus 100. For example, the processing circuit 144 has a collection function 144a, a display control function 144b, and a control function 144c. The processing circuit 144 is realized by a processor, for example.

例えば、処理回路144は、メモリ141から収集機能144aに相当するプログラムを読み出して実行することにより、X線CT装置100を制御してスキャンを実行する。ここで、収集機能144aは、例えば、コンベンショナルスキャンやヘリカルスキャン、ステップアンドシュート方式といった種々の方式でのスキャンを実行することができる。   For example, the processing circuit 144 reads out and executes a program corresponding to the acquisition function 144a from the memory 141, thereby controlling the X-ray CT apparatus 100 and executing a scan. Here, the collection function 144a can execute scans by various methods such as a conventional scan, a helical scan, and a step-and-shoot method, for example.

具体的には、収集機能144aは、寝台駆動装置132を制御することにより、被検体P1を架台装置110の撮影口内へ移動させる。また、収集機能144aは、X線高電圧装置114を制御することにより、X線管111へ高電圧を供給させる。また、収集機能144aは、コリメータ117の開口度及び位置を調整する。また、収集機能144aは、制御装置115を制御することにより、回転フレーム113を含む回転部を回転させる。また、収集機能144aは、データ収集回路118に投影データを収集させる。   Specifically, the collection function 144a moves the subject P1 into the imaging port of the gantry device 110 by controlling the bed driving device 132. Further, the collection function 144 a controls the X-ray high voltage device 114 to supply a high voltage to the X-ray tube 111. The collection function 144a adjusts the aperture and position of the collimator 117. Further, the collection function 144 a rotates the rotating unit including the rotating frame 113 by controlling the control device 115. The collection function 144a causes the data collection circuit 118 to collect projection data.

また、例えば、収集機能144aは、データ収集回路118から出力された検出データに対して対数変換処理やオフセット補正処理、チャネル間の感度補正処理、ビームハードニング補正等の前処理を施したデータを生成する。なお、前処理を施す前のデータ(検出データ)および前処理後のデータを総称して投影データと称する場合もある。また、例えば、収集機能144aは、CT画像データを生成する。具体的には、収集機能144aは、前処理後の投影データに対して、フィルタ補正逆投影法や逐次近似再構成法等を用いた再構成処理を行ってCT画像データを生成する。また、収集機能144aは、入力インターフェース143を介して操作者から受け付けた入力操作等に基づいて、CT画像データを3次元画像データに変換する。   In addition, for example, the collection function 144a performs processing on data obtained by performing preprocessing such as logarithmic conversion processing, offset correction processing, sensitivity correction processing between channels, and beam hardening correction on the detection data output from the data collection circuit 118. Generate. Note that the data before the preprocessing (detection data) and the data after the preprocessing may be collectively referred to as projection data. For example, the collection function 144a generates CT image data. Specifically, the collection function 144a generates CT image data by performing reconstruction processing using a filter-corrected back projection method, successive approximation reconstruction method, or the like on the projection data after preprocessing. The collection function 144a converts CT image data into three-dimensional image data based on an input operation received from an operator via the input interface 143.

ここで、収集機能144aは、上述した3次元画像データの収集を時系列的に行うことにより、4次元画像データを収集する。例えば、収集機能144aは、被検体P1の心臓の収縮期における複数の心位相についてそれぞれ3次元画像データを収集することで、被検体P1の心臓の収縮を3次元的に表した4次元画像データを収集する。一例を挙げると、収集機能144aは、心電計から受信した心電図に基づいてX線の照射を制御し、収縮期における複数の心位相のそれぞれに同期した3次元画像データを収集することで、4次元画像データを収集する。そして、収集機能144aは、4次元画像データを、メモリ141に記憶させたり、画像保管装置20や医用情報処理装置30に送信したりする。   Here, the collection function 144a collects the four-dimensional image data by collecting the above-described three-dimensional image data in time series. For example, the collection function 144a collects three-dimensional image data for each of a plurality of cardiac phases in the systolic phase of the subject P1, thereby obtaining four-dimensional image data that three-dimensionally represents the heart contraction of the subject P1. To collect. For example, the collection function 144a controls the X-ray irradiation based on the electrocardiogram received from the electrocardiograph, and collects three-dimensional image data synchronized with each of a plurality of cardiac phases in the systole. Collect 4D image data. Then, the collection function 144a stores the four-dimensional image data in the memory 141 or transmits it to the image storage device 20 or the medical information processing device 30.

また、例えば、処理回路144は、メモリ141から表示制御機能144bに相当するプログラムを読み出して実行することにより、ディスプレイ142において4次元画像データを表示する。また、例えば、処理回路144は、メモリ141から制御機能144cに相当するプログラムを読み出して実行することにより、入力インターフェース143を介して操作者から受け付けた入力操作に基づいて、処理回路144の各種機能を制御する。   Further, for example, the processing circuit 144 reads out a program corresponding to the display control function 144 b from the memory 141 and executes it, thereby displaying the four-dimensional image data on the display 142. Further, for example, the processing circuit 144 reads out and executes a program corresponding to the control function 144c from the memory 141, and thereby executes various functions of the processing circuit 144 based on an input operation received from an operator via the input interface 143. To control.

図2に示すX線CT装置100においては、各処理機能がコンピュータによって実行可能なプログラムの形態でメモリ141へ記憶されている。処理回路144は、メモリ141からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路144は、読み出したプログラムに対応する機能を有することとなる。なお、図2においては、収集機能144a、表示制御機能144b及び制御機能144cの各処理機能が単一の処理回路144によって実現される場合を示したが、実施形態はこれに限られるものではない。例えば、処理回路144は、複数の独立したプロセッサを組み合わせて構成され、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。また、処理回路144が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。   In the X-ray CT apparatus 100 shown in FIG. 2, each processing function is stored in the memory 141 in the form of a program executable by a computer. The processing circuit 144 is a processor that realizes a function corresponding to each program by reading the program from the memory 141 and executing the program. In other words, the processing circuit 144 that has read each program has a function corresponding to the read program. 2 shows a case where the processing functions of the collection function 144a, the display control function 144b, and the control function 144c are realized by a single processing circuit 144, but the embodiment is not limited thereto. . For example, the processing circuit 144 may be configured by combining a plurality of independent processors, and each processor may implement each processing function by executing each program. In addition, each processing function of the processing circuit 144 may be realized by being appropriately distributed or integrated into a single or a plurality of processing circuits.

上記説明において用いた「プロセッサ」という文言は、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、あるいは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサはメモリ33又はメモリ141に保存されたプログラムを読み出し実行することで機能を実現する。なお、メモリ33又はメモリ141にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて1つのプロセッサとして構成し、その機能を実現するようにしてもよい。   The term “processor” used in the above description is, for example, a CPU (Central Processing Unit), a GPU (Graphics Processing Unit), an application specific integrated circuit (ASIC), a programmable logic device (for example, It means circuits such as a simple programmable logic device (SPLD), a complex programmable logic device (CPLD), and a field programmable gate array (FPGA). The processor implements a function by reading and executing a program stored in the memory 33 or the memory 141. Instead of storing the program in the memory 33 or the memory 141, the program may be directly incorporated in the circuit of the processor. In this case, the processor realizes the function by reading and executing the program incorporated in the circuit. Note that each processor of the present embodiment is not limited to being configured as a single circuit for each processor, but may be configured as a single processor by combining a plurality of independent circuits to realize the function. Good.

以上、医用情報処理システム1の構成の一例について説明した。かかる構成の下、医用情報処理システム1における医用情報処理装置30は、被検体P1の負担を低減しつつ、被検体P1の心筋における電気信号の伝搬を算出する。具体的には、医用情報処理装置30は、以下詳細に説明する処理回路34による処理によって、被検体P1の心筋の4次元画像データに基づく動き情報を取得し、動き情報に対応した心筋の電気信号の伝搬を算出することで、被検体P1の心筋における電気信号の伝搬を非侵襲的に算出する。以下、第1の実施形態に係る医用情報処理装置30が行う処理について詳細に説明する。   Heretofore, an example of the configuration of the medical information processing system 1 has been described. Under such a configuration, the medical information processing apparatus 30 in the medical information processing system 1 calculates the propagation of the electrical signal in the myocardium of the subject P1 while reducing the burden on the subject P1. Specifically, the medical information processing apparatus 30 acquires motion information based on the four-dimensional image data of the myocardium of the subject P1 through processing by the processing circuit 34 described in detail below, and outputs the myocardial electricity corresponding to the motion information. By calculating the propagation of the signal, the propagation of the electrical signal in the myocardium of the subject P1 is calculated non-invasively. Hereinafter, processing performed by the medical information processing apparatus 30 according to the first embodiment will be described in detail.

まず、取得機能34aは、X線CT装置100又は画像保管装置20から、被検体P1の心筋の4次元画像データを取得してメモリ33に格納する。例えば、取得機能34aは、被検体P1の心筋の4次元画像データとして、図3に示す4次元画像データI1を取得して、メモリ33に格納する。なお、図3は、第1の実施形態に係る4次元画像データI1の一例を示す図である。ここで、4次元画像データI1は、被検体P1の心筋について時系列的に収集された3次元画像データの集合であり、図3は、4次元画像データI1に含まれる3次元画像データの一つを示すものである。   First, the acquisition function 34 a acquires four-dimensional image data of the myocardium of the subject P <b> 1 from the X-ray CT apparatus 100 or the image storage apparatus 20 and stores it in the memory 33. For example, the acquisition function 34a acquires the four-dimensional image data I1 shown in FIG. 3 as the four-dimensional image data of the myocardium of the subject P1 and stores it in the memory 33. FIG. 3 is a diagram illustrating an example of the four-dimensional image data I1 according to the first embodiment. Here, the four-dimensional image data I1 is a set of three-dimensional image data collected in time series for the myocardium of the subject P1, and FIG. 3 shows one of the three-dimensional image data included in the four-dimensional image data I1. It shows one.

次に、取得機能34aは、メモリ33から4次元画像データI1を読み込んで、4次元画像データI1の複数位置における心筋の動き情報B1を抽出する。ここで、動き情報B1とは、例えば、心筋の各位置を示す座標と、心筋の各位置における移動量及び移動方向とを対応付けた数値データである。また、動き情報B1とは、例えば、心筋の各位置における移動量及び移動方向をマッピングした画像データである。   Next, the acquisition function 34a reads the four-dimensional image data I1 from the memory 33, and extracts myocardial motion information B1 at a plurality of positions of the four-dimensional image data I1. Here, the motion information B1 is, for example, numerical data in which coordinates indicating each position of the myocardium are associated with a movement amount and a movement direction at each position of the myocardium. In addition, the motion information B1 is image data in which a movement amount and a movement direction at each position of the myocardium are mapped, for example.

例えば、取得機能34aは、まず、4次元画像データI1に含まれる複数の3次元画像データの各々から、心筋に対応する領域(心筋領域)を抽出する。一例を挙げると、取得機能34aは、3次元画像データの各々に対して、CT値が空間的に連続する領域を抽出する領域拡張(region growing)法や形状テンプレートを用いたパターンマッチング法などを用いてセグメンテーション処理を行うことにより、心筋領域を抽出する。   For example, the acquisition function 34a first extracts a region (myocardial region) corresponding to the myocardium from each of a plurality of three-dimensional image data included in the four-dimensional image data I1. For example, the acquisition function 34a uses a region growing method for extracting regions where CT values are spatially continuous, a pattern matching method using a shape template, and the like for each of the three-dimensional image data. A myocardial region is extracted by performing a segmentation process using it.

次に、取得機能34aは、3次元画像データの各々から抽出した心筋領域を相互に位置合わせし、解剖学的な位置の対応を算出する。これにより、取得機能34aは、各3次元画像データが収集される間に解剖学的な位置が移動した方向及び移動量を、動き情報B1として抽出することができる。   Next, the acquisition function 34a aligns the myocardial regions extracted from each of the three-dimensional image data with each other, and calculates the correspondence between the anatomical positions. Thereby, the acquisition function 34a can extract the direction and amount of movement of the anatomical position while each piece of 3D image data is collected as the movement information B1.

次に、算出機能34bは、モデルM11を使用し、電気信号の伝搬と筋細胞との動態シミュレーションを実行して、動き情報B1に対応した被検体P1の心筋の電気信号の伝搬A1を算出する。ここで、算出機能34bによる電気信号の伝搬A1の算出について、図4を用いて具体的に説明する。図4は、第1の実施形態に係る電気信号の伝搬と筋細胞との動態シミュレーションの一例を示す図である。   Next, using the model M11, the calculation function 34b executes the electrical signal propagation and the dynamic simulation of muscle cells, and calculates the myocardial electrical signal propagation A1 of the subject P1 corresponding to the motion information B1. . Here, the calculation of the electric signal propagation A1 by the calculation function 34b will be specifically described with reference to FIG. FIG. 4 is a diagram illustrating an example of electrical signal propagation and muscle cell dynamic simulation according to the first embodiment.

例えば、算出機能34bは、まず、メモリ33からモデルM11を読み出す。ここで、モデルM11は、細胞変形と電気信号伝搬との関連を示すモデルである。例えば、モデルM11は、図4に示すように、細胞の動き(細胞が生む力や細胞の長さ等の時間変化)と、電気信号(Na+やCa2+といったイオンの濃度等)との電気生理学的な関連を示した筋細胞モデルである。なお、メモリ33は、細胞の種類(例えば、固有心筋、特殊心筋等)に応じたモデルM11を記憶する場合であってもよい。この場合、算出機能34bは、対象部位に含まれる細胞の種類に応じて、1又は複数のモデルM11をメモリ33から読み出す。   For example, the calculation function 34b first reads the model M11 from the memory 33. Here, the model M11 is a model showing the relationship between cell deformation and electric signal propagation. For example, as shown in FIG. 4, the model M11 has an electrophysiological relationship between a cell movement (a time change such as a force generated by a cell and a cell length) and an electrical signal (concentration of ions such as Na + and Ca2 +). It is a muscle cell model that showed a strong relationship. Note that the memory 33 may store a model M11 corresponding to a cell type (for example, a specific myocardium, a special myocardium, etc.). In this case, the calculation function 34b reads one or a plurality of models M11 from the memory 33 in accordance with the type of cell included in the target region.

次に、算出機能34bは、図4に示すように、4次元画像データI1に基づく動き情報B1を、モデルM11に入力して、心筋の電気信号の伝搬A1を算出する。即ち、算出機能34bは、モデルM11を、4次元画像データI1の複数位置における動き情報B1に適用することにより、4次元画像データI1の複数位置における電気信号の伝搬A1を算出する。具体的には、まず、算出機能34bは、図4に示すように、4次元画像データI1に基づいて心筋モデルM21を生成する。例えば、算出機能34bは、4次元画像データI1から抽出した心筋領域に基づいて、被検体P1の心筋の形状を示す形状データを生成し、生成した形状データを筋細胞の大きさ及び形状に応じて分割することで、心筋モデルM21を生成する。即ち、算出機能34bは、心筋モデルM21を、筋細胞の集合として生成する。   Next, as shown in FIG. 4, the calculation function 34b inputs the motion information B1 based on the four-dimensional image data I1 to the model M11 and calculates the propagation A1 of the myocardial electrical signal. That is, the calculation function 34b calculates the electric signal propagation A1 at the plurality of positions of the four-dimensional image data I1 by applying the model M11 to the motion information B1 at the plurality of positions of the four-dimensional image data I1. Specifically, first, the calculation function 34b generates a myocardial model M21 based on the four-dimensional image data I1, as shown in FIG. For example, the calculation function 34b generates shape data indicating the shape of the myocardium of the subject P1 based on the myocardial region extracted from the four-dimensional image data I1, and uses the generated shape data according to the size and shape of the myocyte. The myocardial model M21 is generated. That is, the calculation function 34b generates the myocardial model M21 as a set of myocytes.

次に、算出機能34bは、心筋モデルM21における個々の筋細胞において、その筋細胞の位置に対応した動き情報B1を、力や細胞の長さ等の時間変化としてモデルM11に入力する。これにより、算出機能34bは、心筋モデルM21における個々の筋細胞において、入力した力や細胞の長さ等の時間変化に応じた、イオン濃度の時間変化を算出する。更に、算出機能34bは、心筋モデルM21における個々の筋細胞において、イオン濃度の時間変化に応じた電気信号の時間変化を算出することで、心筋の電気信号の伝搬A1を算出する。   Next, the calculation function 34b inputs motion information B1 corresponding to the position of each muscle cell in the myocardial model M21 to the model M11 as a change over time such as force and cell length. Thereby, the calculation function 34b calculates the time change of the ion concentration according to the time change of the input force, the length of the cell, etc. in each myocyte in the myocardial model M21. Further, the calculation function 34b calculates the electrical signal propagation A1 of the myocardium by calculating the temporal change of the electrical signal according to the temporal change of the ion concentration in each myocyte in the myocardial model M21.

ここで、算出機能34bは、筋細胞間での電気信号の連続性を考慮して、電気信号の伝搬A1を算出するものであってもよい。即ち、電気信号は、筋細胞を伝搬する間に途絶等することはあるとしても、ある程度の連続性を持って流れるものであるから、算出機能34bは、筋細胞間での電気信号の連続性を実現するように、心筋モデルM21における個々の筋細胞における電気信号を調整しながら、心筋の電気信号の伝搬A1を算出してもよい。   Here, the calculation function 34b may calculate the propagation A1 of the electric signal in consideration of the continuity of the electric signal between muscle cells. In other words, even if the electrical signal is interrupted while propagating through the muscle cells, the calculation function 34b has the continuity of the electrical signal between the muscle cells because it flows with a certain degree of continuity. The myocardial electrical signal propagation A1 may be calculated while adjusting the electrical signal in each myocyte in the myocardial model M21.

なお、算出機能34bが電気信号の伝搬A1を算出する範囲は、心筋の一部でもよいし、心筋の全部でもよい。例えば、算出機能34bは、被検体P1の症例に応じた範囲において、電気信号の伝搬A1を算出する。一例を挙げると、被検体P1の心筋が興奮伝搬障害を有している場合、算出機能34bは、電気信号の伝搬A1として、心筋内部の各位置、及び、心筋の表面(心腔側の面及び外側の面)の各位置について、電気信号の時間変化をそれぞれ算出する。これにより、算出機能34bは、心筋内部において電気信号の伝搬が阻害されている場合でも、異常個所を算出することができる。一方で、心筋内部の電気信号について観察を要しない症例である場合、算出機能34bは、電気信号の伝搬A1として、心筋の心腔側の面における電気信号の時間変化のみを算出してもよい。   The range in which the calculation function 34b calculates the electric signal propagation A1 may be a part of the myocardium or the entire myocardium. For example, the calculation function 34b calculates the electric signal propagation A1 in a range corresponding to the case of the subject P1. As an example, when the myocardium of the subject P1 has an excitatory propagation disorder, the calculation function 34b uses each position inside the myocardium and the surface of the myocardium (surface on the heart chamber side) as the electric signal propagation A1. And the temporal change of the electric signal for each position on the outer surface). Thereby, the calculation function 34b can calculate an abnormal part even when propagation of an electrical signal is inhibited inside the myocardium. On the other hand, in the case where it is not necessary to observe the electrical signal inside the myocardium, the calculation function 34b may calculate only the time change of the electrical signal on the heart chamber side surface of the myocardium as the propagation A1 of the electrical signal. .

上述したように、算出機能34bは、電気信号の伝搬と筋細胞との動態シミュレーションを実行して、電気信号の伝搬A1を算出する。ここで、算出機能34bは、心筋の各位置における細胞が周辺細胞から受ける力を考慮して、電気信号の伝搬A1を算出してもよい。これは、電気信号を適切に算出するためには、各細胞が電気信号に起因して自ら行う動きをモデルM11に入力する必要があるのに対して、4次元画像データI1に基づく動き情報B1は、各細胞が電気信号に起因して自ら行う動きと、周辺細胞に引っ張られたことによる動きとの双方を含んでいるためである。   As described above, the calculation function 34b executes the electric signal propagation and the dynamic simulation of muscle cells to calculate the electric signal propagation A1. Here, the calculation function 34b may calculate the propagation A1 of the electric signal in consideration of the force that the cell at each position of the myocardium receives from the surrounding cells. This is because, in order to appropriately calculate the electric signal, each cell needs to input the movement that the cell itself performs due to the electric signal to the model M11, whereas the movement information B1 based on the four-dimensional image data I1. This is because each cell includes both a movement caused by an electric signal and a movement caused by being pulled by a surrounding cell.

例えば、算出機能34bは、4次元画像データI1に基づく心筋モデルM21を用いて、心筋の各位置における細胞が周辺細胞から受ける力を算出する。一例を挙げると、算出機能34bは、心筋モデルM21において、心筋の各位置における細胞にラプラスモデル(以下の式(1))を適用して、心筋の各位置における細胞が周辺細胞から受ける力を算出する。   For example, the calculation function 34b uses the myocardial model M21 based on the four-dimensional image data I1 to calculate the force that the cells at each position of the myocardium receive from surrounding cells. For example, in the myocardial model M21, the calculation function 34b applies the Laplace model (the following equation (1)) to the cells at each position of the myocardium, thereby calculating the force that the cells at each position of the myocardium receive from surrounding cells. calculate.

Figure 2019042497
Figure 2019042497

式(1)において、「ΔP」は心筋内外の圧力差であり、「T」は心筋の各位置における細胞が周辺細胞から受ける力(張力)であり、「R」は曲率半径である。例えば、算出機能34bは、心筋モデルM21に基づいて、心筋の形状が示す曲率半径Rを、心筋の位置ごとに取得する。また、算出機能34bは、心筋にかかる力に基づいて、圧力差ΔPを取得する。例えば、算出機能34bは、被検体P1の血圧に基づく解析を行なうことにより、圧力差ΔPを取得する。一例を挙げると、算出機能34bは、まず、心筋の各位置が、心臓におけるいずれの部屋(左心房、左心室、右心房及び右心室)に位置しているかを判定する。次に、算出機能34bは、判定した部屋における圧力の基準値(例えば、予め複数人について計測した圧力の平均値等)をメモリ33から取得する。次に、算出機能34bは、取得した基準値を被検体P1の血圧に応じて補正することで、圧力差ΔPを取得する。更に、算出機能34bは、曲率半径R及び圧力差ΔPを式(1)に代入することにより、張力Tを算出する。   In Expression (1), “ΔP” is a pressure difference inside and outside the myocardium, “T” is a force (tension) that a cell at each position of the myocardium receives from surrounding cells, and “R” is a radius of curvature. For example, the calculation function 34b acquires, for each position of the myocardium, the radius of curvature R indicated by the shape of the myocardium based on the myocardial model M21. Moreover, the calculation function 34b acquires the pressure difference ΔP based on the force applied to the myocardium. For example, the calculation function 34b acquires the pressure difference ΔP by performing an analysis based on the blood pressure of the subject P1. As an example, the calculation function 34b first determines in which room (left atrium, left ventricle, right atrium, and right ventricle) in the heart each position of the myocardium is located. Next, the calculation function 34b acquires from the memory 33 a reference value of pressure in the determined room (for example, an average value of pressures measured in advance for a plurality of persons). Next, the calculation function 34b acquires the pressure difference ΔP by correcting the acquired reference value according to the blood pressure of the subject P1. Furthermore, the calculation function 34b calculates the tension T by substituting the curvature radius R and the pressure difference ΔP into the equation (1).

次に、算出機能34bは、算出した張力Tの影響を除外するように、4次元画像データI1に基づく動き情報B1を補正する。言い換えると、算出機能34bは、4次元画像データI1に基づく動き情報B1から、周辺細胞に引っ張られたことによる動きを差分し、各細胞が電気信号に起因して自ら行う動きを算出する。そして、算出機能34bは、モデルM11を使用し、電気信号の伝搬と筋細胞との動態シミュレーションを実行して、補正後の動き情報B1(各細胞が電気信号に起因して自ら行う動き)に対応した心筋の電気信号の伝搬A1を算出する。   Next, the calculation function 34b corrects the motion information B1 based on the four-dimensional image data I1 so as to exclude the influence of the calculated tension T. In other words, the calculation function 34b subtracts the movement caused by being pulled by surrounding cells from the movement information B1 based on the four-dimensional image data I1, and calculates the movement that each cell performs due to the electrical signal. Then, the calculation function 34b uses the model M11 to execute electrical signal propagation and muscle cell dynamics simulation, and use the corrected motion information B1 (the motion that each cell performs itself due to the electrical signal). The propagation A1 of the corresponding myocardial electrical signal is calculated.

なお、これまで、心筋の各位置における細胞が周辺細胞から受ける力を、ラプラスモデルを用いて算出する場合について説明したが、実施形態はこれに限定されるものではない。例えば、算出機能34bは、有限要素法(FEM:Finite Element Method)や境界要素法等の数値解析法を用いた構造解析を行なうことにより、心筋の各位置における細胞が周辺細胞から受ける力を算出してもよい。   Heretofore, the case has been described in which the force received by the cells at each position of the myocardium from the surrounding cells is calculated using the Laplace model, but the embodiment is not limited to this. For example, the calculation function 34b calculates the force that a cell at each position of the myocardium receives from surrounding cells by performing a structural analysis using a numerical analysis method such as a finite element method (FEM) or a boundary element method. May be.

また、算出機能34bは、動き情報B1から、心筋の各細胞が電気信号に起因して自ら行なう動きの成分を抽出して、電気信号伝搬を算出してもよい。これは、4次元画像データI1に基づく動き情報B1は、各細胞が自ら行う動きの他、心臓全体が被検体P1に対して揺れ動くことによる動き成分や、血液の流出入による心臓の膨張及び収縮に伴う動き成分を含むためである。   In addition, the calculation function 34b may extract an electrical signal propagation by extracting a component of motion that each cell of the myocardium itself causes from the electrical signal from the motion information B1. This is because the movement information B1 based on the four-dimensional image data I1 includes the movement component of each cell, the movement component due to the whole heart swinging with respect to the subject P1, and the expansion and contraction of the heart due to the inflow and outflow of blood. This is because it includes a motion component associated with.

具体的には、算出機能34bは、まず、動き情報B1から、心筋の各細胞が自ら行なう動きの成分を抽出する。換言すると、算出機能34bは、動き情報B1から、心筋の各位置における細胞が発生させる力の方向の成分に対応する動き情報(以下、動き情報B1’とする)を抽出する。動き情報B1’は、部分動き情報の一例である。   Specifically, the calculation function 34b first extracts, from the movement information B1, a component of movement performed by each cell of the myocardium. In other words, the calculation function 34b extracts motion information (hereinafter referred to as motion information B1 ') corresponding to a component in the direction of force generated by cells at each position of the myocardium from the motion information B1. The motion information B1 'is an example of partial motion information.

一例を挙げると、算出機能34bは、まず、心筋の各位置における細胞の繊維方向を取得する。ここで、算出機能34bは、心筋の各位置における細胞の繊維方向を、知見に基づいて取得することができる。例えば、算出機能34bは、電気信号伝搬の算出処理に先立って、心筋の各位置における細胞の繊維方向を示す3次元モデルMfを取得し、メモリ33に格納する。一例を挙げると、算出機能34bは、標準的な心臓の形状を持つ3次元モデルに対して、標準的な心臓上の繊維方向をマッピングすることによって3次元モデルMfを生成し、メモリ33に格納する。なお、算出機能34bは、医用情報処理装置30以外の他の装置において生成された3次元モデルMfを取得して、メモリ33に格納してもよい。次に、算出機能34bは、被検体P1の心臓の形状に応じて3次元モデルMfを変形させる。例えば、算出機能34bは、4次元画像データI1に基づく心筋モデルM21と形状が一致するように、3次元モデルMfを変形させる。そして、算出機能34bは、変形後の3次元モデルMfが示す繊維方向を被検体P1の心臓に割り当てることで、心筋の各位置における細胞の繊維方向を取得する。例えば、算出機能34bは、変形後の3次元モデルMfが示す繊維方向を心筋モデルM21の各位置に割り当てることで、心筋モデルM21の各位置における細胞の繊維方向を取得する。ここで、心筋の細胞は、その繊維の方向に沿って力を発生させる性質がある。従って、算出機能34bは、細胞の繊維方向に基づいて、動き情報B1から、心筋の各位置における細胞が発生させる力の方向の成分に対応する動き情報B1’を抽出することができる。例えば、算出機能34bは、動き情報B1を、細胞の繊維方向の成分と他の方向の成分とに分解し、細胞の繊維方向の成分を動き情報B1’として抽出する。   For example, the calculation function 34b first acquires the fiber direction of the cell at each position of the myocardium. Here, the calculation function 34b can acquire the fiber direction of the cell at each position of the myocardium based on the knowledge. For example, the calculation function 34 b acquires a three-dimensional model Mf indicating the fiber direction of the cell at each position of the myocardium and stores it in the memory 33 prior to the calculation process of electrical signal propagation. For example, the calculation function 34b generates a three-dimensional model Mf by mapping a fiber direction on a standard heart with respect to a three-dimensional model having a standard heart shape, and stores it in the memory 33. To do. Note that the calculation function 34 b may acquire a three-dimensional model Mf generated by a device other than the medical information processing device 30 and store it in the memory 33. Next, the calculation function 34b deforms the three-dimensional model Mf according to the shape of the heart of the subject P1. For example, the calculation function 34b deforms the three-dimensional model Mf so that the shape matches the myocardial model M21 based on the four-dimensional image data I1. Then, the calculation function 34b acquires the fiber direction of the cell at each position of the myocardium by assigning the fiber direction indicated by the deformed three-dimensional model Mf to the heart of the subject P1. For example, the calculation function 34b acquires the fiber direction of the cell at each position of the myocardial model M21 by assigning the fiber direction indicated by the deformed three-dimensional model Mf to each position of the myocardial model M21. Here, myocardial cells have the property of generating force along the direction of their fibers. Therefore, the calculation function 34b can extract the motion information B1 'corresponding to the component in the direction of the force generated by the cell at each position of the myocardium from the motion information B1 based on the fiber direction of the cell. For example, the calculation function 34b decomposes the motion information B1 into a component in the cell fiber direction and a component in another direction, and extracts the component in the cell fiber direction as the motion information B1 '.

別の例を挙げると、算出機能34bは、動き情報B1から、心筋の曲面に垂直な成分を除外することにより、動き情報B1’を抽出する。例えば、算出機能34bは、まず、4次元画像データI1に基づく心筋モデルM21から、心筋の形状を示す曲面を抽出する。次に、算出機能34bは、抽出した曲面の各位置において法線方向を算出する。ここで、心筋の細胞は、通常、心筋の曲面に垂直な方向の力を発生することはない。即ち、心室又は心房の内部から心臓壁に対して垂直に発生する力は、通常、血液の流出入による心臓の膨張及び収縮によるものであって、心筋の細胞が発生させたものではない。従って、算出機能34bは、動き情報B1から心筋の曲面に垂直な成分を除外することにより、心筋の各位置における細胞が発生させる力の方向の成分に対応する動き情報B1’を抽出することができる。例えば、算出機能34bは、動き情報B1を、心筋の曲面に垂直な成分と他の方向の成分とに分解し、他の方向の成分を動き情報B1’として抽出する。   As another example, the calculation function 34b extracts motion information B1 'by excluding a component perpendicular to the curved surface of the myocardium from the motion information B1. For example, the calculation function 34b first extracts a curved surface indicating the shape of the myocardium from the myocardial model M21 based on the four-dimensional image data I1. Next, the calculation function 34b calculates a normal direction at each position of the extracted curved surface. Here, myocardial cells usually do not generate a force in a direction perpendicular to the curved surface of the myocardium. That is, the force generated perpendicularly to the heart wall from the inside of the ventricle or the atrium is usually due to the expansion and contraction of the heart due to the inflow and outflow of blood, and is not generated by myocardial cells. Accordingly, the calculation function 34b can extract the motion information B1 ′ corresponding to the component in the direction of the force generated by the cells at each position of the myocardium by excluding the component perpendicular to the curved surface of the myocardium from the motion information B1. it can. For example, the calculation function 34b decomposes the motion information B1 into a component perpendicular to the curved surface of the myocardium and a component in another direction, and extracts the component in the other direction as the motion information B1 '.

そして、算出機能34bは、モデルM11を使用し、電気信号の伝搬と筋細胞との動態シミュレーションを実行して、動き情報B1’に対応した心筋の電気信号の伝搬A1を算出する。即ち、算出機能34bは、モデルM11を、4次元画像データI1の複数位置における動き情報B1’に適用することにより、4次元画像データI1の複数位置における電気信号の伝搬A1を算出する。ここで、動き情報B1’を用いることにより、算出機能34bは、電気信号の伝搬A1の精度を向上させることができる。即ち、算出機能34bは、心筋の細胞が自ら行なう動きの成分を抽出して動態シミュレーションを実行することにより、動態シミュレーションの精度を向上させることができる。なお、動き情報B1’は、心筋の各位置における細胞が周辺細胞から受ける力の影響を除外するように補正したものであってもよい。   Then, the calculation function 34b calculates the propagation A1 of the myocardial electrical signal corresponding to the motion information B1 'by using the model M11 and executing the electrical signal propagation and the dynamic simulation of the muscle cell. That is, the calculation function 34b calculates the electric signal propagation A1 at the plurality of positions of the four-dimensional image data I1 by applying the model M11 to the motion information B1 'at the plurality of positions of the four-dimensional image data I1. Here, by using the motion information B1 ', the calculation function 34b can improve the accuracy of the propagation A1 of the electric signal. That is, the calculation function 34b can improve the accuracy of the dynamic simulation by extracting the component of the movement performed by the myocardial cell itself and executing the dynamic simulation. The motion information B1 'may be corrected so as to exclude the influence of the force that the cells at each position of the myocardium receive from the surrounding cells.

算出機能34bが電気信号の伝搬A1を算出した後、表示制御機能34cは、ディスプレイ32において電気信号の伝搬A1を表示する。例えば、表示制御機能34cは、ディスプレイ32において、4次元画像データI1を動画像又は静止画像として表示するとともに、電気信号の伝搬A1を表示する。なお、表示制御機能34cは、4次元画像データI1に含まれる心筋の全体を表示してもよいし、電気信号の伝搬A1に異常のある領域のみを表示してもよい。   After the calculation function 34b calculates the electric signal propagation A1, the display control function 34c displays the electric signal propagation A1 on the display 32. For example, the display control function 34c displays, on the display 32, the four-dimensional image data I1 as a moving image or a still image, and also displays electric signal propagation A1. Note that the display control function 34c may display the entire myocardium included in the four-dimensional image data I1, or may display only a region where the electric signal propagation A1 is abnormal.

一例を挙げると、まず、表示制御機能34cは、心筋の電気信号の伝搬A1に基づいて、心筋の各位置に電気信号が到達した時刻(興奮開始時刻)や、電気信号が最大となった最大興奮時刻、興奮の時間幅等を取得する。そして、表示制御機能34cは、ディスプレイ32において、取得した興奮開始時刻、最大興奮時刻、時間幅等に応じたカラーを付した4次元画像データI1を表示する。   For example, first, the display control function 34c, based on the propagation A1 of the myocardial electrical signal, the time when the electrical signal reaches each position of the myocardium (excitation start time) and the maximum when the electrical signal is maximized. Get excitement time, excitement duration, etc. Then, the display control function 34c displays, on the display 32, the four-dimensional image data I1 with a color corresponding to the acquired excitement start time, maximum excitement time, time width, and the like.

別の一例を挙げると、まず、表示制御機能34cは、心筋の電気信号の伝搬A1に基づいて、心筋の各位置における電圧の大きさ(最大値や平均値等)を取得する。そして、表示制御機能34cは、ディスプレイ32において、取得した電圧の大きさに応じたカラーを付して、4次元画像データI1を表示する。   As another example, first, the display control function 34c acquires the magnitude (maximum value, average value, etc.) of the voltage at each position of the myocardium based on the propagation A1 of the myocardial electrical signal. Then, the display control function 34c displays the four-dimensional image data I1 on the display 32 with a color according to the magnitude of the acquired voltage.

別の一例を挙げると、まず、表示制御機能34cは、心筋の電気信号の伝搬A1に基づいて、心筋の各位置における電圧の時間変化を取得する。そして、表示制御機能34cは、ディスプレイ32において、電圧の変化量をカラーで表現しながら、4次元画像データI1を動画像として表示する。また、例えば、表示制御機能34cは、ディスプレイ32において、電圧の時間変化を表したグラフを表示する。   As another example, first, the display control function 34c acquires a temporal change in voltage at each position of the myocardium based on the propagation A1 of the myocardial electrical signal. The display control function 34c displays the four-dimensional image data I1 as a moving image on the display 32 while expressing the amount of voltage change in color. Further, for example, the display control function 34 c displays a graph representing a change in voltage over time on the display 32.

別の一例を挙げると、まず、表示制御機能34cは、心筋の電気信号の伝搬A1に基づいて、心筋の各位置における電気信号の伝搬方向を取得する。そして、表示制御機能34cは、ディスプレイ32において、被検体P1の心筋の各位置における電気信号の伝搬方向を表す矢印(ベクトル)を付した4次元画像データI1を表示する。   As another example, first, the display control function 34c acquires the propagation direction of the electrical signal at each position of the myocardium based on the propagation A1 of the myocardial electrical signal. The display control function 34c displays on the display 32 the four-dimensional image data I1 with an arrow (vector) indicating the propagation direction of the electrical signal at each position of the myocardium of the subject P1.

また、表示制御機能34cは、種々の情報と併せて、電気信号の伝搬A1を表示することとしてもよい。例えば、取得機能34aは、被検体P1の心臓について収集された時系列の3次元超音波画像データを更に取得する。次に、算出機能34bは、時系列の3次元超音波画像データから、心筋の各位置における変位(Displacement)を算出する。次に、算出機能34bは、心筋の各位置について、変位を微分することにより速度を算出する。更に、算出機能34bは、心筋の各位置における速度の時間変化から、心筋の局所的な歪み(Strain)の時間変化曲線であるStrainカーブを算出する。   The display control function 34c may display the electric signal propagation A1 together with various information. For example, the acquisition function 34a further acquires time-series three-dimensional ultrasound image data collected for the heart of the subject P1. Next, the calculation function 34b calculates a displacement (Displacement) at each position of the myocardium from the time-series three-dimensional ultrasonic image data. Next, the calculation function 34b calculates the velocity by differentiating the displacement for each position of the myocardium. Further, the calculation function 34b calculates a Strain curve that is a time change curve of local myocardial strain (Strain) from the time change in velocity at each position of the myocardium.

そして、表示制御機能34cは、電気信号の伝搬A1と併せて、Strainカーブに基づく3次元超音波画像データを表示する。例えば、表示制御機能34cは、電気信号の伝搬A1に基づく興奮開始時刻に応じてカラーを付した4次元画像データI1と、Strainが最大となる時刻に応じたカラーを付した3次元超音波画像データとを、ディスプレイ32において表示する。これにより、表示制御機能34cは、被検体P1の心臓の動きについて、より多くの情報を操作者に提供することができる。   The display control function 34c displays the three-dimensional ultrasonic image data based on the strain curve together with the electric signal propagation A1. For example, the display control function 34c includes the four-dimensional image data I1 colored according to the excitement start time based on the propagation A1 of the electric signal, and the three-dimensional ultrasonic image colored according to the time when the strain becomes maximum. The data is displayed on the display 32. Thereby, the display control function 34c can provide more information about the movement of the heart of the subject P1 to the operator.

なお、3次元超音波画像データに基づくStrainカーブは、周辺細胞の動きに引きずられる影響があること、及び、変位を微分して求める事情から、ノイズが多い。これに対して、算出機能34bは、4次元画像データI1に基づく動き情報B1を電気信号に分解することにより、局所の細胞が動いた真の時刻を精度よく推定することができる。そして、表示制御機能34cは、電気信号の伝搬A1を表示することで、心筋の各細胞が動いた時刻を高精度で操作者に提示することができる。   Note that the strain curve based on the three-dimensional ultrasonic image data has a lot of noise due to the influence of being dragged by the movement of surrounding cells and the fact that the displacement is differentiated. In contrast, the calculation function 34b can accurately estimate the true time when the local cell has moved by decomposing the motion information B1 based on the four-dimensional image data I1 into an electrical signal. The display control function 34c can display the time when each cell of the myocardium moves to the operator with high accuracy by displaying the propagation A1 of the electric signal.

また、算出機能34bは、電気信号の伝搬A1を算出した後、電気信号の伝搬A1に基づいて、心筋の動き情報を算出してもよい。以下、電気信号の伝搬A1に基づいて算出される心筋の動き情報を、動き情報B2とする。   The calculation function 34b may calculate myocardial motion information based on the electric signal propagation A1 after calculating the electric signal propagation A1. Hereinafter, the motion information of the myocardium calculated based on the propagation A1 of the electrical signal is referred to as motion information B2.

例えば、算出機能34bは、動き情報B1に対応した電気信号の伝搬A1を算出する処理の逆の処理により、電気信号の伝搬A1に対応した動き情報B2を算出する。即ち、算出機能34bは、モデルM11を動き情報B1に適用することによって電気信号の伝搬A1を算出した後、モデルM11を電気信号の伝搬A1に適用することによって動き情報B2を算出する。   For example, the calculation function 34b calculates the motion information B2 corresponding to the electrical signal propagation A1 by the reverse process of the process of calculating the electrical signal propagation A1 corresponding to the motion information B1. That is, the calculation function 34b calculates the motion information B2 by applying the model M11 to the electrical signal propagation A1 after calculating the electrical signal propagation A1 by applying the model M11 to the motion information B1.

算出機能34bが動き情報B2を算出した後、表示制御機能34cは、ディスプレイ32において、動き情報B2と4次元画像データI1とを比較可能に表示する。例えば、表示制御機能34cは、4次元画像データI1に含まれる1つの3次元画像データを動き情報B2に基づいて順次変形させることで、動き情報B2に基づく4次元画像データI2を生成する。そして、表示制御機能34cは、4次元画像データI1と4次元画像データI2とを並べて表示する。   After the calculation function 34b calculates the motion information B2, the display control function 34c displays the motion information B2 and the four-dimensional image data I1 on the display 32 so that they can be compared. For example, the display control function 34c generates four-dimensional image data I2 based on the motion information B2 by sequentially deforming one three-dimensional image data included in the four-dimensional image data I1 based on the motion information B2. Then, the display control function 34c displays the 4D image data I1 and the 4D image data I2 side by side.

動き情報B2と4次元画像データI1とを比較可能に表示することにより、表示制御機能34cは、算出された電気信号の伝搬A1の評価を可能とする。即ち、電気信号の伝搬A1に基づく動き情報B2と、4次元画像データI1とが整合している場合、操作者は、電気信号の伝搬A1の精度を信頼できると判断することができる。   By displaying the motion information B2 and the four-dimensional image data I1 so as to be comparable, the display control function 34c enables evaluation of the propagation A1 of the calculated electric signal. That is, when the motion information B2 based on the electric signal propagation A1 and the four-dimensional image data I1 match, the operator can determine that the accuracy of the electric signal propagation A1 is reliable.

また、4次元画像データI1は、心筋の各細胞が電気信号に起因して自ら行なう動き以外の他の成分(心臓全体が揺れ動くことによる動き成分や、心臓の膨張及び収縮に伴う動き成分)を含んでいる。これに対して、電気信号の伝搬A1に基づく動き情報B2は、心筋の各細胞が自ら行なう動きを示すものである。従って、表示制御機能34cは、動き情報B2と4次元画像データI1とを比較可能に表示することにより、被検体P1の心筋の動きについて、他の成分の有無による違いを操作者に提示することができる。   In addition, the four-dimensional image data I1 includes other components (motion components due to the whole heart shaking and motion components associated with the expansion and contraction of the heart) other than the motion of each cell of the myocardium itself due to the electrical signal. Contains. On the other hand, the motion information B2 based on the propagation A1 of the electric signal indicates the motion that each cell of the myocardium performs by itself. Therefore, the display control function 34c displays the motion information B2 and the four-dimensional image data I1 so as to be comparable, thereby presenting to the operator the difference in the myocardial motion of the subject P1 due to the presence or absence of other components. Can do.

次に、医用情報処理装置30による処理の手順の一例を、図5を用いて説明する。図5は、第1の実施形態に係る医用情報処理装置30の処理の一連の流れを説明するためのフローチャートである。ステップS101、ステップS102及びステップS106は、取得機能34aに対応するステップである。ステップS103及びステップS104は、算出機能34bに対応するステップである。ステップS105は、表示制御機能34cに対応するステップである。   Next, an example of a processing procedure performed by the medical information processing apparatus 30 will be described with reference to FIG. FIG. 5 is a flowchart for explaining a series of processing steps of the medical information processing apparatus 30 according to the first embodiment. Step S101, step S102, and step S106 are steps corresponding to the acquisition function 34a. Steps S103 and S104 are steps corresponding to the calculation function 34b. Step S105 is a step corresponding to the display control function 34c.

まず、処理回路34は、X線CT装置100又は画像保管装置20から、被検体P1の心筋の4次元画像データI1を取得してメモリ33に格納する(ステップS101)。次に、処理回路34は、メモリ33から4次元画像データI1を読み出して、4次元画像データI1から被検体P1の心筋の動き情報B1を取得する(ステップS102)。次に、処理回路34は、モデルM11を使用し、電気信号の伝搬と筋細胞との動態シミュレーションを実行して(ステップS103)、動き情報B1に対応した電気信号の伝搬A1を算出する(ステップS104)。そして、処理回路34は、算出した電気信号の伝搬A1を、ディスプレイ32において表示する(ステップS105)。   First, the processing circuit 34 acquires the four-dimensional image data I1 of the myocardium of the subject P1 from the X-ray CT apparatus 100 or the image storage apparatus 20 and stores it in the memory 33 (step S101). Next, the processing circuit 34 reads the four-dimensional image data I1 from the memory 33 and acquires the myocardial motion information B1 of the subject P1 from the four-dimensional image data I1 (step S102). Next, the processing circuit 34 uses the model M11 to execute electric signal propagation and dynamic simulation of muscle cells (step S103), and calculates the electric signal propagation A1 corresponding to the motion information B1 (step S103). S104). Then, the processing circuit 34 displays the calculated propagation A1 of the electric signal on the display 32 (step S105).

ここで、処理回路34は、更に動き情報B1を取得したか否かを判定する(ステップS106)。動き情報B1を取得した場合(ステップS106肯定)、処理回路34は、再度ステップS103に移行する。一方で、動き情報B1を取得しなかった場合(ステップS106否定)、処理回路34は、処理を終了する。   Here, the processing circuit 34 further determines whether or not the motion information B1 has been acquired (step S106). When the motion information B1 is acquired (Yes at Step S106), the processing circuit 34 proceeds to Step S103 again. On the other hand, when the motion information B1 is not acquired (No at Step S106), the processing circuit 34 ends the process.

上述したように、第1の実施形態によれば、取得機能34aは、4次元画像データI1に基づく被検体P1の心筋の動き情報B1を取得する。算出機能34bは、細胞の動きと電気信号との関連を示すモデルM11を使用し、動き情報B1に対応した心筋の電気信号の伝搬A1を算出する。従って、第1の実施形態に係る医用情報処理装置30は、心筋における電気信号の伝搬を非侵襲的に算出することができる。例えば、医用情報処理装置30は、被検体P1の心臓に直接刺激を与えて頻脈を誘発するといった患者負担を要することなく、電気信号の伝搬を算出することができる。   As described above, according to the first embodiment, the acquisition function 34a acquires the myocardial motion information B1 of the subject P1 based on the four-dimensional image data I1. The calculation function 34b calculates the propagation A1 of the myocardial electrical signal corresponding to the motion information B1, using the model M11 indicating the relationship between the cell motion and the electrical signal. Therefore, the medical information processing apparatus 30 according to the first embodiment can non-invasively calculate the propagation of the electrical signal in the myocardium. For example, the medical information processing apparatus 30 can calculate the propagation of the electric signal without requiring a patient burden of directly stimulating the heart of the subject P1 to induce tachycardia.

また、上述したように、第1の実施形態によれば、算出機能34bは、細胞が周辺細胞から受ける張力Tを更に算出し、張力Tの影響を除外するように動き情報B1を補正して、補正後の動き情報B1に対応した電気信号の伝搬A1を算出する。従って、第1の実施形態に係る医用情報処理装置30は、心筋状態を考慮して、電気信号の伝搬をより適切に算出することができる。即ち、心筋の変形のみに着目して電気信号の伝搬を評価する場合と比較し、医用情報処理装置30は、電気信号の伝搬を精度よく算出することができる。   As described above, according to the first embodiment, the calculation function 34b further calculates the tension T that the cell receives from the surrounding cells, and corrects the motion information B1 so as to exclude the influence of the tension T. Then, the propagation A1 of the electrical signal corresponding to the corrected motion information B1 is calculated. Therefore, the medical information processing apparatus 30 according to the first embodiment can more appropriately calculate the propagation of the electric signal in consideration of the myocardial state. That is, the medical information processing apparatus 30 can calculate the propagation of the electrical signal with higher accuracy than the case where the propagation of the electrical signal is evaluated by focusing only on the deformation of the myocardium.

なお、これまで、1つの心拍について動き情報B1を取得する場合について説明したが、取得機能34aは、複数の心拍について動き情報B1を取得する場合であってもよい。例えば、まず、X線CT装置100における収集機能144aは、心電計から受信した心電図に基づいてX線の照射を制御し、複数の心位相のそれぞれに同期した投影データを収集する。次に、収集機能144aは、投影データを、正常な心拍(不整脈)について収集された投影データと、異常な心拍(整脈)について収集された投影データとに分類する。ここで、収集機能144aは、各投影データが正常な心拍について収集されたものか異常な心拍について収集されたものかを、心電図に基づいて分類してもよいし、投影データ自体に基づいて分類してもよい。   Although the case where the motion information B1 is acquired for one heartbeat has been described so far, the acquisition function 34a may be a case where the motion information B1 is acquired for a plurality of heartbeats. For example, first, the collection function 144a in the X-ray CT apparatus 100 controls X-ray irradiation based on an electrocardiogram received from an electrocardiograph, and collects projection data synchronized with each of a plurality of cardiac phases. Next, the collection function 144a classifies the projection data into projection data collected for a normal heartbeat (arrhythmia) and projection data collected for an abnormal heartbeat (arrhythmia). Here, the collection function 144a may classify whether each projection data is collected for a normal heartbeat or an abnormal heartbeat based on the electrocardiogram or based on the projection data itself. May be.

次に、収集機能144aは、正常な心拍について収集された投影データから、複数の心位相のそれぞれについて3次元画像データを生成することで、4次元画像データI11を収集する。また、収集機能144aは、異常な心拍について収集された投影データから、複数の心位相のそれぞれについて3次元画像データを生成することで、4次元画像データI12を収集する。そして、収集機能144aは、正常な心拍について収集された4次元画像データI11と、異常な心拍について収集された4次元画像データI12とを、メモリ141に記憶させたり、画像保管装置20や医用情報処理装置30に送信したりする。なお、4次元画像データI11及び4次元画像データI12は、4次元画像データI1の一例である。   Next, the collection function 144a collects four-dimensional image data I11 by generating three-dimensional image data for each of a plurality of cardiac phases from projection data collected for a normal heartbeat. The collection function 144a collects the four-dimensional image data I12 by generating three-dimensional image data for each of a plurality of cardiac phases from the projection data collected for an abnormal heartbeat. The collection function 144a stores the four-dimensional image data I11 collected for the normal heartbeat and the four-dimensional image data I12 collected for the abnormal heartbeat in the memory 141, or stores the image storage device 20 and the medical information. Or transmitted to the processing device 30. The 4D image data I11 and the 4D image data I12 are examples of the 4D image data I1.

次に、取得機能34aは、正常な心拍について収集された4次元画像データI11から動き情報B11を取得し、異常な心拍について収集された4次元画像データI12から動き情報B12を取得する。即ち、取得機能34aは、複数の心拍についてそれぞれ動き情報を取得する。次に、算出機能34bは、モデルM11を使用し、動き情報B11に対応した電気信号の伝搬A11と、動き情報B12に対応した電気信号の伝搬A12とをそれぞれ算出する。即ち、算出機能34bは、心拍ごとに、電気信号の伝搬を算出する。なお、動き情報B11及び動き情報B12は動き情報B1の一例であり、電気信号の伝搬A11及び電気信号の伝搬A12は電気信号の伝搬A1の一例である。   Next, the acquisition function 34a acquires the motion information B11 from the 4D image data I11 collected for the normal heartbeat, and acquires the motion information B12 from the 4D image data I12 collected for the abnormal heartbeat. That is, the acquisition function 34a acquires motion information for each of a plurality of heartbeats. Next, the calculation function 34b uses the model M11 to calculate the propagation A11 of the electrical signal corresponding to the motion information B11 and the propagation A12 of the electrical signal corresponding to the motion information B12. That is, the calculation function 34b calculates the propagation of the electrical signal for each heartbeat. The motion information B11 and the motion information B12 are examples of the motion information B1, and the electrical signal propagation A11 and the electrical signal propagation A12 are examples of the electrical signal propagation A1.

そして、表示制御機能34cは、電気信号の伝搬A11と電気信号の伝搬A12とを、ディスプレイ32において比較可能に表示する。例えば、表示制御機能34cは、電気信号の伝搬A11と電気信号の伝搬A12とを、ディスプレイ32において、並べて表示したり、切り替えて表示したりする。これにより、医用情報処理装置30は、心拍間の差異と電気信号の伝搬との関係を表示し、診断を支援することができる。例えば、医用情報処理装置30は、異常な心拍における電気信号の伝搬A12において、正常な心拍における電気信号の伝搬A11と異なる部分を、心拍の異常の原因として提示することができる。   The display control function 34c displays the electric signal propagation A11 and the electric signal propagation A12 on the display 32 so that they can be compared. For example, the display control function 34c displays the electric signal propagation A11 and the electric signal propagation A12 side by side on the display 32 or by switching them. Thereby, the medical information processing apparatus 30 can display the relationship between the difference between heartbeats and the propagation of an electric signal, and can support diagnosis. For example, the medical information processing apparatus 30 can present a portion of the electrical signal propagation A12 in an abnormal heartbeat that is different from the electrical signal propagation A11 in a normal heartbeat as the cause of the heartbeat abnormality.

また、これまで、1の被検体P1から動き情報B1を取得する場合について説明したが、取得機能34aは、複数の被検体P1から動き情報B1を取得する場合であってもよい。例えば、まず、X線CT装置100における収集機能144aは、心臓に疾患を有しない被検体P11の心筋の4次元画像データI13と、心臓に疾患を有する被検体P12の心筋の4次元画像データI14とを収集する。次に、収集機能144aは、4次元画像データI13と4次元画像データI14とを、メモリ141に記憶させたり、画像保管装置20や医用情報処理装置30に送信したりする。なお、被検体P11及び被検体P12は被検体P1の一例であり、4次元画像データI13及び4次元画像データI14は4次元画像データI1の一例である。   Although the case where the motion information B1 is acquired from one subject P1 has been described so far, the acquisition function 34a may be a case where the motion information B1 is acquired from a plurality of subjects P1. For example, first, the acquisition function 144a in the X-ray CT apparatus 100 performs the four-dimensional image data I13 of the myocardium of the subject P11 having no heart disease and the four-dimensional image data I14 of the myocardium of the subject P12 having a heart disease. And collect. Next, the collection function 144a stores the 4D image data I13 and the 4D image data I14 in the memory 141 or transmits them to the image storage device 20 or the medical information processing device 30. The subject P11 and the subject P12 are examples of the subject P1, and the four-dimensional image data I13 and the four-dimensional image data I14 are examples of the four-dimensional image data I1.

次に、取得機能34aは、4次元画像データI13から動き情報B13を取得し、4次元画像データI14から動き情報B14を取得する。即ち、取得機能34aは、複数の被検体P1から動き情報を取得する。次に、算出機能34bは、モデルM11を使用し、動き情報B13に対応した電気信号の伝搬A13と、動き情報B14に対応した電気信号の伝搬A14とをそれぞれ算出する。即ち、算出機能34bは、被検体ごとに、電気信号の伝搬を算出する。なお、動き情報B13及び動き情報B14は動き情報B1の一例であり、電気信号の伝搬A13及び電気信号の伝搬A14は電気信号の伝搬A1の一例である。   Next, the acquisition function 34a acquires the motion information B13 from the 4D image data I13, and acquires the motion information B14 from the 4D image data I14. That is, the acquisition function 34a acquires motion information from a plurality of subjects P1. Next, using the model M11, the calculation function 34b calculates an electric signal propagation A13 corresponding to the motion information B13 and an electric signal propagation A14 corresponding to the motion information B14. That is, the calculation function 34b calculates the propagation of the electrical signal for each subject. The motion information B13 and the motion information B14 are examples of the motion information B1, and the electrical signal propagation A13 and the electrical signal propagation A14 are examples of the electrical signal propagation A1.

そして、表示制御機能34cは、被検体ごとに算出された電気信号の伝搬A13と電気信号の伝搬A14とを、ディスプレイ32において、比較可能に表示する。例えば、表示制御機能34cは、電気信号の伝搬A13と電気信号の伝搬A14とを、ディスプレイ32において、並べて表示したり、切り替えて表示したりする。これにより、医用情報処理装置30は、被検体間の差異と電気信号の伝搬との関係を表示し、診断を支援することができる。例えば、医用情報処理装置30は、心臓に疾患を有する被検体P12の電気信号の伝搬A14において、心臓に疾患を有しない被検体P11の電気信号の伝搬A13と異なる部分を、心拍の異常の原因として提示することができる。   Then, the display control function 34c displays the electric signal propagation A13 and the electric signal propagation A14 calculated for each subject on the display 32 so that they can be compared. For example, the display control function 34c displays the electric signal propagation A13 and the electric signal propagation A14 side by side on the display 32 or by switching them. Thereby, the medical information processing apparatus 30 can display the relationship between the difference between the subjects and the propagation of the electric signal, and can support the diagnosis. For example, the medical information processing apparatus 30 determines the part of the electrical signal propagation A14 of the subject P12 having a disease in the heart that is different from the electrical signal propagation A13 of the subject P11 having no disease in the heart, Can be presented as

また、これまで、心臓の収縮期について4次元画像データI1を収集する場合について説明したが、収集機能144aは、拡張期を含めた1心拍分の4次元画像データI1を収集する場合であってもよい。例えば、被検体P1が心臓に有する疾患等に起因して、収縮期以外にも電気信号の伝搬が生じ得る場合、収集機能144aは、1心拍分の4次元画像データI1を収集する。そして、取得機能34aは、1心拍分の4次元画像データI1に基づく動き情報B1を取得し、算出機能34bは、1心拍分の動き情報B1に対応した電気信号の伝搬A1を算出する。   The case where the four-dimensional image data I1 is collected for the systole of the heart has been described so far, but the collection function 144a is a case where the four-dimensional image data I1 for one heartbeat including the diastole is collected. Also good. For example, when an electric signal can be propagated in addition to the systole due to a disease or the like that the subject P1 has in the heart, the collection function 144a collects four-dimensional image data I1 for one heartbeat. The acquisition function 34a acquires the motion information B1 based on the four-dimensional image data I1 for one heartbeat, and the calculation function 34b calculates the propagation A1 of the electrical signal corresponding to the motion information B1 for one heartbeat.

(第2の実施形態)
上述した第1の実施形態では、メモリ33に記憶されたモデルM11をそのまま使用し、心筋の電気信号の伝搬A1を算出する場合について説明した。これに対し、第2の実施形態では、被検体P1に応じたモデルM12を使用して、心筋の電気信号の伝搬A1を算出する場合について説明する。
(Second Embodiment)
In the first embodiment described above, the case where the model M11 stored in the memory 33 is used as it is and the propagation A1 of the myocardial electrical signal is calculated has been described. In contrast, in the second embodiment, a case where the propagation A1 of the myocardial electrical signal is calculated using the model M12 corresponding to the subject P1 will be described.

第2の実施形態に係る医用情報処理装置30は、図1に示した第1の実施形態に係る医用情報処理装置30と同様の構成を有し、算出機能34bによる処理の一部が相違する。そこで、第1の実施形態において説明した構成と同様の構成を有する点については、図1と同一の符号を付し、説明を省略する。   The medical information processing apparatus 30 according to the second embodiment has the same configuration as that of the medical information processing apparatus 30 according to the first embodiment shown in FIG. 1, and a part of the processing by the calculation function 34b is different. . Therefore, the same reference numerals as those in FIG. 1 are given to the components having the same configurations as those described in the first embodiment, and the description thereof is omitted.

例えば、算出機能34bは、まず、被検体P1について、イオン濃度、血圧及び病変分布の少なくとも一つを取得する。一例を挙げると、まず、操作者は、電解質測定機器を用いて、被検体P1の血液のイオン濃度を測定する。かかる電解質測定機器は、例えば、被検体P1から採取された血液を検査してイオン濃度を測定するものであってもよいし、被検体P1の指先に照射した光に基づいてイオン濃度を測定するものであってもよい。そして、算出機能34bは、入力インターフェース31を介して、操作者からイオン濃度の入力操作を受け付けることにより、被検体P1のイオン濃度を取得する。   For example, the calculation function 34b first acquires at least one of ion concentration, blood pressure, and lesion distribution for the subject P1. As an example, first, the operator measures the ion concentration of blood in the subject P1 using an electrolyte measurement device. Such an electrolyte measurement device may be, for example, a device that examines blood collected from the subject P1 and measures the ion concentration, or measures the ion concentration based on light irradiated to the fingertip of the subject P1. It may be a thing. Then, the calculation function 34b receives the ion concentration input operation from the operator via the input interface 31, thereby acquiring the ion concentration of the subject P1.

また、一例を挙げると、操作者は、血圧測定機器を用いて、被検体P1の血圧を測定する。そして、算出機能34bは、入力インターフェース31を介して、操作者から血圧の入力操作を受け付けることにより、被検体P1の血圧を取得する。また、一例を挙げると、算出機能34bは、4次元画像データI1について、病変部位を検出するコンピュータ支援診断(Computer Aided Diagnosis:CAD)処理を実行することにより、被検体P1の心筋における病変分布を取得する。   For example, the operator measures the blood pressure of the subject P1 using a blood pressure measurement device. Then, the calculation function 34b acquires the blood pressure of the subject P1 by receiving a blood pressure input operation from the operator via the input interface 31. For example, the calculation function 34b performs a computer-aided diagnosis (CAD) process for detecting a lesion site on the four-dimensional image data I1, thereby calculating a lesion distribution in the myocardium of the subject P1. get.

次に、算出機能34bは、メモリ33からモデルM11を読み出して、被検体P1のイオン濃度、血圧及び病変分布の少なくとも一つに応じてモデルM11を変形することにより、被検体P1に応じたモデルM12を生成する。例えば、心筋梗塞が生じている被検体P1の心臓において、心筋梗塞部(壊死している細胞の分布)が明らかとなっている場合、算出機能34bは、心筋梗塞部については心筋梗塞の筋細胞モデルを割り当て、心筋梗塞部以外についてはモデルM11と同様の筋細胞モデルを割り当てたモデルM12を生成する。そして、算出機能34bは、被検体P1に応じたモデルM12を使用して、心筋の電気信号の伝搬A1を算出する。   Next, the calculation function 34b reads the model M11 from the memory 33, and deforms the model M11 according to at least one of the ion concentration, blood pressure, and lesion distribution of the subject P1, and thereby the model according to the subject P1. M12 is generated. For example, when the myocardial infarction part (distribution of necrotic cells) is clear in the heart of the subject P1 in which myocardial infarction has occurred, the calculation function 34b performs myocardial infarction myocytes for the myocardial infarction part. A model is assigned, and a model M12 to which a muscle cell model similar to the model M11 is assigned except for the myocardial infarction is generated. Then, the calculation function 34b calculates the propagation A1 of the myocardial electrical signal using the model M12 corresponding to the subject P1.

これまで、被検体P1から取得される条件(イオン濃度の測定値、血圧の測定値、CAD処理による病変分布等)に応じて、モデルM12を生成する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、算出機能34bは、仮の条件に応じてモデルM12を生成する場合であってもよい。   So far, the case where the model M12 is generated according to the conditions (measured values of ion concentration, measured values of blood pressure, lesion distribution by CAD processing, etc.) acquired from the subject P1 has been described. However, the embodiment is not limited to this, and the calculation function 34b may be a case where the model M12 is generated according to provisional conditions.

例えば、まず、取得機能34aは、入力インターフェース31を介して、操作者から、被検体P1の心筋に対して病変分布を設定する操作を受け付ける。ここで、操作者は、例えば、被検体P1の病変が進行した場合や、手術により病変を除去した場合等を仮定して、病変分布を設定する。次に、算出機能34bは、設定された病変分布に応じたモデルM12を使用して、心筋の電気信号の伝搬A1を算出する。そして、表示制御機能34cは、算出された電気信号の伝搬A1を、ディスプレイ32において表示する。例えば、表示制御機能34cは、CAD処理による病変分布に応じたモデルM12を使用して算出された電気信号の伝搬A1と、操作者により設定された病変分布に応じたモデルM12を使用して算出された電気信号の伝搬A1とを、ディスプレイ32において、比較可能に表示する。   For example, first, the acquisition function 34a receives an operation for setting a lesion distribution for the myocardium of the subject P1 from the operator via the input interface 31. Here, for example, the operator sets the lesion distribution on the assumption that the lesion of the subject P1 has progressed or the lesion has been removed by surgery. Next, the calculation function 34b calculates the propagation A1 of the myocardial electrical signal using the model M12 corresponding to the set lesion distribution. Then, the display control function 34 c displays the calculated propagation A1 of the electric signal on the display 32. For example, the display control function 34c calculates using the propagation A1 of the electrical signal calculated using the model M12 corresponding to the lesion distribution by CAD processing and the model M12 corresponding to the lesion distribution set by the operator. The transmitted electric signal A1 is displayed on the display 32 so as to be comparable.

上述したように、第2の実施形態によれば、算出機能34bは、被検体P1に応じたモデルM12を使用して、心筋の電気信号の伝搬A1を算出する。従って、第2の実施形態に係る医用情報処理装置30は、被検体P1の個別の情報を用いて、心筋における電気信号の伝搬をより適切に算出することができる。   As described above, according to the second embodiment, the calculation function 34b calculates the propagation A1 of the myocardial electrical signal using the model M12 corresponding to the subject P1. Therefore, the medical information processing apparatus 30 according to the second embodiment can more appropriately calculate the propagation of the electrical signal in the myocardium using the individual information of the subject P1.

また、上述したように、第2の実施形態によれば、算出機能34bは、設定された病変分布に応じたモデルM12を使用して、心筋の電気信号の伝搬A1を算出する。従って、第2の実施形態に係る医用情報処理装置30は、現在の被検体P1における電気信号の伝搬のみならず、カテーテルアブレーションを実行した後の電気信号の伝搬や、カテーテルアブレーションを実行せずに病変が進行した場合の電気信号の伝搬等をも算出し、診断を支援することができる。   As described above, according to the second embodiment, the calculation function 34b calculates the propagation A1 of the myocardial electrical signal using the model M12 corresponding to the set lesion distribution. Therefore, the medical information processing apparatus 30 according to the second embodiment performs not only the propagation of the electrical signal in the current subject P1, but also the propagation of the electrical signal after performing the catheter ablation and the catheter ablation. Diagnosis can be supported by calculating the propagation of electrical signals when the lesion progresses.

(第3の実施形態)
さて、これまで第1〜第2の実施形態について説明したが、上述した実施形態以外にも、種々の異なる形態にて実施されてよいものである。
(Third embodiment)
The first and second embodiments have been described so far, but the present invention may be implemented in various different forms other than the above-described embodiments.

例えば、算出機能34bは、更に、心電図を用いて、電気信号の伝搬A1を算出してもよい。一例を挙げると、算出機能34bは、被検体P1の心筋の電気信号の伝搬A1を、被検体P1の心電図が示す電気信号の伝搬A2に対応するように算出する。   For example, the calculation function 34b may further calculate the electric signal propagation A1 using an electrocardiogram. For example, the calculation function 34b calculates the myocardial electrical signal propagation A1 of the subject P1 so as to correspond to the electrical signal propagation A2 indicated by the electrocardiogram of the subject P1.

具体的には、まず、取得機能34aは、被検体P1の心電図を取得する。ここで、取得機能34aは、4次元画像データI1を収集するスキャンの際中に測定された心電図を取得する場合であってもよいし、スキャンの前又は後に測定された心電図を取得する場合であってもよい。   Specifically, first, the acquisition function 34a acquires an electrocardiogram of the subject P1. Here, the acquisition function 34a may be a case of acquiring an electrocardiogram measured during a scan for collecting the four-dimensional image data I1, or a case of acquiring an electrocardiogram measured before or after the scan. There may be.

次に、算出機能34bは、心電図に基づいて、被検体P1の心筋における電気信号の伝搬A2を算出する。例えば、算出機能34bは、心電図の波形から、P波及びQRS波のそれぞれに相当する時間区間を識別する。ここで、算出機能34bは、洞結節において電気信号が発生した時間として、P波の開始時間を取得する。また、算出機能34bは、右心房を電気信号が伝搬した時間として、P波に相当する時間区間の開始後2/3を取得する。また、算出機能34bは、左心房を電気信号が伝搬した時間として、P波に相当する時間区間の終了前2/3を取得する。また、算出機能34bは、ヒス束、左脚/右脚を電気信号が伝搬した時間として、QRS波に相当する時間区間を取得する。   Next, the calculation function 34b calculates an electric signal propagation A2 in the myocardium of the subject P1 based on the electrocardiogram. For example, the calculation function 34b identifies time intervals corresponding to the P wave and the QRS wave from the waveform of the electrocardiogram. Here, the calculation function 34b acquires the start time of the P wave as the time when the electrical signal is generated in the sinus node. In addition, the calculation function 34b acquires 2/3 after the start of the time interval corresponding to the P wave as the time that the electrical signal propagates through the right atrium. In addition, the calculation function 34b acquires 2/3 before the end of the time interval corresponding to the P wave as the time that the electrical signal propagates through the left atrium. Further, the calculation function 34b acquires a time interval corresponding to the QRS wave as the time when the electrical signal propagates through the His bundle and the left leg / right leg.

また、算出機能34bは、モデルM11又はモデルM12を使用して動き情報B1に対応した心筋の電気信号の伝搬A1を算出する。また、算出機能34bは、算出した電気信号の伝搬A1と、心電図が示す電気信号の伝搬A2とを比較する。ここで、電気信号の伝搬A1と電気信号の伝搬A2との差異が小さい場合、算出機能34bは処理を終了する。   The calculation function 34b calculates the propagation A1 of the myocardial electrical signal corresponding to the motion information B1 using the model M11 or the model M12. The calculation function 34b compares the calculated electric signal propagation A1 with the electric signal propagation A2 indicated by the electrocardiogram. If the difference between the electrical signal propagation A1 and the electrical signal propagation A2 is small, the calculation function 34b ends the process.

一方で、電気信号の伝搬A1と電気信号の伝搬A2との差異が閾値より大きい場合、算出機能34bは、電気信号の伝搬A1を補正する。例えば、洞結節にて電気信号が発生してから右心房における電気信号の伝搬が開始するまでの時間について、電気信号の伝搬A2が示す値に対して電気信号の伝搬A1が示す値が著しく大きい場合、算出機能34bは、洞結節にて電気信号が発生してから右心房における電気信号の伝搬が開始するまでの時間を短縮するように、電気信号の伝搬A1を補正する。   On the other hand, when the difference between the electric signal propagation A1 and the electric signal propagation A2 is larger than the threshold value, the calculation function 34b corrects the electric signal propagation A1. For example, with respect to the time from when an electrical signal is generated at the sinus node to when the propagation of the electrical signal in the right atrium starts, the value indicated by the electrical signal propagation A1 is significantly larger than the value indicated by the electrical signal propagation A2. In this case, the calculation function 34b corrects the propagation A1 of the electrical signal so as to shorten the time from when the electrical signal is generated in the sinus node until the propagation of the electrical signal in the right atrium starts.

次に、算出機能34bは、モデルM11又はモデルM12を使用し、補正後の電気信号の伝搬A1に対応した心筋の動き情報B3を算出する。即ち、算出機能34bは、動き情報B1に対応した電気信号の伝搬A1を算出する処理の逆の処理により、電気信号の伝搬A1に対応した動き情報B3を算出する。次に、算出機能34bは、動き情報B1と動き情報B3とを比較する。   Next, the calculation function 34b uses the model M11 or the model M12 to calculate myocardial motion information B3 corresponding to the corrected electric signal propagation A1. That is, the calculation function 34b calculates the motion information B3 corresponding to the electric signal propagation A1 by the reverse process of the process of calculating the electric signal propagation A1 corresponding to the motion information B1. Next, the calculation function 34b compares the motion information B1 with the motion information B3.

ここで、動き情報B1と動き情報B3との差異が小さい場合、算出機能34bは処理を終了する。一方で、動き情報B1と動き情報B3との差異が閾値より大きい場合、算出機能34bは、動き情報B1を補正する。次に、算出機能34bは、補正後の動き情報B1に対応した心筋の電気信号の伝搬A1を算出する。更に、算出機能34bは、上述した電気信号の伝搬A1と電気信号の伝搬A2との比較、及び、動き情報B1と動き情報B3との比較を、所定の回数或いは差異が閾値を下回るまで、繰り返し実行する。これにより、算出機能34bは、動き情報B1に対応した電気信号の伝搬A1の精度を、心電図に基づいて向上させることができる。   Here, when the difference between the motion information B1 and the motion information B3 is small, the calculation function 34b ends the process. On the other hand, when the difference between the motion information B1 and the motion information B3 is larger than the threshold value, the calculation function 34b corrects the motion information B1. Next, the calculation function 34b calculates myocardial electrical signal propagation A1 corresponding to the corrected motion information B1. Further, the calculation function 34b repeats the comparison between the electric signal propagation A1 and the electric signal propagation A2 and the comparison between the motion information B1 and the motion information B3 until a predetermined number of times or a difference falls below a threshold value. Run. Thereby, the calculation function 34b can improve the accuracy of the propagation A1 of the electric signal corresponding to the motion information B1 based on the electrocardiogram.

また、算出機能34bは、電気信号の伝搬A1を算出した後、電気信号の伝搬A1が制約条件を満たすか否かを判定してもよい。即ち、算出機能34bは、算出した電気信号の伝搬A1が、実際に起こり得るものであるか否かを判定してもよい。   The calculation function 34b may determine whether or not the electric signal propagation A1 satisfies the constraint condition after calculating the electric signal propagation A1. That is, the calculation function 34b may determine whether or not the calculated propagation A1 of the electric signal can actually occur.

例えば、被検体P1の心筋を電気信号が伝搬する速度には限界があり、この限界を超えた速度で電気信号が伝搬することはない。また、例えば、電気信号の伝搬はつながっているものであり、不連続に電気信号が伝搬することはない。そこで、算出機能34bは、4次元画像データI1の複数位置における電気信号の伝搬A1について、心筋における電気信号伝搬の速度及び連続性の少なくとも一方に関する制約条件を満たすか否かを判定する。   For example, there is a limit on the speed at which the electric signal propagates through the myocardium of the subject P1, and the electric signal does not propagate at a speed exceeding this limit. Further, for example, the propagation of the electrical signal is continuous, and the electrical signal does not propagate discontinuously. Therefore, the calculation function 34b determines whether or not a constraint condition regarding at least one of the speed and continuity of electrical signal propagation in the myocardium is satisfied for the electrical signal propagation A1 at a plurality of positions of the four-dimensional image data I1.

ここで、電気信号の伝搬A1が制約条件を満たさない場合、算出機能34bは、例えば、4次元画像データI1の複数位置における電気信号の伝搬A1を再算出する。なお、算出機能34bは、算出条件を変更して、電気信号の伝搬A1を再算出することとしてもよい。例えば、算出機能34bは、使用するモデルや動き情報を変更して、電気信号の伝搬A1を再算出する。一例を挙げると、モデルM11を用いて電気信号の伝搬A1を算出していた場合において、算出機能34bはモデルM12を用いて電気信号の伝搬A1を再算出する。別の例を挙げると、動き情報B1を用いて電気信号の伝搬A1を算出していた場合において、算出機能34bは動き情報B1’を用いて電気信号の伝搬A1を再算出する。   When the electric signal propagation A1 does not satisfy the constraint condition, the calculation function 34b recalculates the electric signal propagation A1 at a plurality of positions of the four-dimensional image data I1, for example. The calculation function 34b may recalculate the electric signal propagation A1 by changing the calculation condition. For example, the calculation function 34b changes the model to be used and the motion information, and recalculates the electric signal propagation A1. For example, in the case where the electric signal propagation A1 is calculated using the model M11, the calculation function 34b recalculates the electric signal propagation A1 using the model M12. As another example, when the electric signal propagation A1 is calculated using the motion information B1, the calculation function 34b recalculates the electric signal propagation A1 using the motion information B1 '.

或いは、算出機能34bは、電気信号の伝搬A1のうち、制約条件を満たさない部分を除外することとしてもよい。即ち、算出機能34bは、算出した電気信号の伝搬A1のうち、制約条件を満たさない部分を取り除いて、最終結果としてもよい。この場合、表示制御機能34cは、電気信号の伝搬A1のうち、制約条件を満たさない部分を除く他の部分を表示する。或いは、表示制御機能34cは、電気信号の伝搬A1に制約条件を満たさない部分が含まれることが認識できる態様で、電気信号の伝搬A1を表示する。例えば、表示制御機能34cは、電気信号の伝搬A1のうち制約条件を満たさない部分を黒色として、電気信号の伝搬A1を表示する。   Alternatively, the calculation function 34b may exclude a portion of the electric signal propagation A1 that does not satisfy the constraint condition. That is, the calculation function 34b may remove a portion that does not satisfy the constraint condition from the calculated propagation A1 of the electric signal and obtain the final result. In this case, the display control function 34c displays other portions of the electric signal propagation A1 except for portions that do not satisfy the constraint conditions. Alternatively, the display control function 34c displays the electric signal propagation A1 in such a manner that it can be recognized that the electric signal propagation A1 includes a portion that does not satisfy the constraint condition. For example, the display control function 34c displays the electric signal propagation A1 by setting the portion of the electric signal propagation A1 that does not satisfy the constraint condition to black.

また、上述した実施形態では、医用情報処理装置30の処理回路34が、取得機能34a及び算出機能34bを有する場合について説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、X線CT装置100の処理回路144が、取得機能34a及び算出機能34bに相当する機能(以下、取得機能144d及び算出機能144eと記載する)を有する場合であってもよい。   In the above-described embodiment, the case where the processing circuit 34 of the medical information processing apparatus 30 has the acquisition function 34a and the calculation function 34b has been described. However, the embodiment is not limited to this. For example, the processing circuit 144 of the X-ray CT apparatus 100 may have a function corresponding to the acquisition function 34a and the calculation function 34b (hereinafter referred to as the acquisition function 144d and the calculation function 144e).

この場合、まず、収集機能144aが、被検体P1の心筋の4次元画像データI1を収集する。次に、取得機能144dは、4次元画像データI1から、被検体P1の心筋の動き情報B1を取得する。次に、算出機能144eは、モデルM11又はモデルM12を使用し、動き情報B1に対応した心筋の電気信号の伝搬A1を算出する。そして、表示制御機能144bは、算出された電気信号の伝搬A1を、ディスプレイ142において表示する。   In this case, first, the collection function 144a collects the four-dimensional image data I1 of the myocardium of the subject P1. Next, the acquisition function 144d acquires the myocardial motion information B1 of the subject P1 from the four-dimensional image data I1. Next, the calculation function 144e calculates the propagation A1 of the myocardial electrical signal corresponding to the motion information B1 using the model M11 or the model M12. Then, the display control function 144b displays the calculated propagation A1 of the electric signal on the display 142.

また、上述した実施形態では、医用情報処理装置30における取得機能34aが、4次元画像データI1から動き情報を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、X線CT装置100における取得機能144dが4次元画像データI1から被検体P1の心筋の動き情報B1を算出し、取得機能34aは、X線CT装置100から動き情報B1を取得する場合であってもよい。   In the above-described embodiment, the case where the acquisition function 34a in the medical information processing apparatus 30 calculates motion information from the four-dimensional image data I1 has been described. However, the embodiment is not limited to this. For example, the acquisition function 144d in the X-ray CT apparatus 100 calculates the myocardial motion information B1 of the subject P1 from the four-dimensional image data I1, and the acquisition function 34a acquires the motion information B1 from the X-ray CT apparatus 100. There may be.

また、上述した実施形態では、被検体P1の心筋について、電気信号の伝搬A1を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、取得機能34aは、心筋以外の筋肉を対象部位として、対象部位の4次元画像データに基づく動き情報を取得することができる。更に、算出機能34bは、モデルM11を使用し、動き情報に対応した対象部位の電気信号の伝搬を算出することができる。   In the above-described embodiment, the case where the propagation A1 of the electric signal is calculated for the myocardium of the subject P1 has been described. However, the embodiment is not limited to this. For example, the acquisition function 34a can acquire motion information based on the four-dimensional image data of the target part using a muscle other than the myocardium as the target part. Furthermore, the calculation function 34b can use the model M11 to calculate the propagation of the electrical signal of the target part corresponding to the motion information.

また、上述した実施形態では、モデルM11又はモデルM12を使用して、電気信号の伝搬A1を算出する場合について説明した。しかしながら、実施形態はこれに限定されるものではない。例えば、算出機能34bは、動き情報B1と電気信号の伝搬A1との対応関係を定めたデータベースをメモリ33から読み出し、取得機能34aが取得した動き情報B1をデータベースと比較することにより、電気信号の伝搬A1を算出する場合であってもよい。   In the above-described embodiment, the case where the propagation A1 of the electric signal is calculated using the model M11 or the model M12 has been described. However, the embodiment is not limited to this. For example, the calculation function 34b reads a database that defines the correspondence between the motion information B1 and the propagation A1 of the electrical signal from the memory 33, and compares the motion information B1 acquired by the acquisition function 34a with the database, thereby It may be a case where the propagation A1 is calculated.

上述した実施形態に係る各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。即ち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行われる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現されうる。   Each component of each device according to the above-described embodiment is functionally conceptual and does not necessarily need to be physically configured as illustrated. In other words, the specific form of distribution / integration of each device is not limited to that shown in the figure, and all or a part thereof may be functionally or physically distributed or arbitrarily distributed in arbitrary units according to various loads or usage conditions. Can be integrated and configured. Further, all or a part of each processing function performed in each device may be realized by a CPU and a program that is analyzed and executed by the CPU, or may be realized as hardware by wired logic.

また、上述した実施形態で説明した医用情報処理方法は、予め用意されたプログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。このプログラムは、インターネット等のネットワークを介して配布することができる。また、このプログラムは、ハードディスク、フレキシブルディスク(FD)、CD−ROM、MO、DVD等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。   Further, the medical information processing method described in the above-described embodiment can be realized by executing a program prepared in advance on a computer such as a personal computer or a workstation. This program can be distributed via a network such as the Internet. The program can also be executed by being recorded on a computer-readable recording medium such as a hard disk, a flexible disk (FD), a CD-ROM, an MO, and a DVD and being read from the recording medium by the computer.

以上説明した少なくとも1つの実施形態によれば、患者負担を低減しつつ、電気信号の伝搬を算出することができる。   According to at least one embodiment described above, propagation of an electric signal can be calculated while reducing patient burden.

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。   Although several embodiments of the present invention have been described, these embodiments are presented by way of example and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the spirit of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and the equivalents thereof.

1 医用情報処理システム
30 医用情報処理装置
34 処理回路
34a 取得機能
34b 算出機能
34c 表示制御機能
DESCRIPTION OF SYMBOLS 1 Medical information processing system 30 Medical information processing apparatus 34 Processing circuit 34a Acquisition function 34b Calculation function 34c Display control function

Claims (21)

細胞変形と電気信号伝搬との関連を示すモデルを記憶する記憶部と、
心筋の3次元構造を時系列で撮像することにより得られた4次元画像データを取得し、取得した4次元画像データの複数位置における前記心筋の動き情報を抽出する取得部と、
前記モデルを前記複数位置における前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する算出部と、
を備える、医用情報処理装置。
A storage unit for storing a model indicating a relationship between cell deformation and electric signal propagation;
An acquisition unit that acquires four-dimensional image data obtained by imaging a three-dimensional structure of the myocardium in time series, and extracts motion information of the myocardium at a plurality of positions of the acquired four-dimensional image data;
A calculation unit that calculates electric signal propagation at the plurality of positions by applying the model to the motion information at the plurality of positions;
A medical information processing apparatus comprising:
前記算出部は、前記心筋の各位置における細胞が周辺細胞から受ける力を更に算出し、当該力の影響を除外するように前記動き情報を補正して、前記モデルを前記複数位置における補正後の前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する、請求項1に記載の医用情報処理装置。   The calculation unit further calculates a force that a cell at each position of the myocardium receives from surrounding cells, corrects the motion information so as to exclude the influence of the force, and corrects the model after the correction at the plurality of positions. The medical information processing apparatus according to claim 1, wherein electrical signal propagation at the plurality of positions is calculated by applying the motion information. 前記算出部は、前記心筋の形状、及び、前記心筋にかかる力に基づいて、前記心筋の各位置における細胞が周辺細胞から受ける力を算出する、請求項2に記載の医用情報処理装置。   The medical information processing apparatus according to claim 2, wherein the calculation unit calculates a force that a cell at each position of the myocardium receives from surrounding cells based on the shape of the myocardium and the force applied to the myocardium. 前記算出部は、前記動き情報から、前記心筋の各位置における細胞が発生させる力の方向の成分に対応する部分動き情報を抽出し、前記モデルを前記複数位置における前記部分動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する、請求項1乃至3のいずれか一項に記載の医用情報処理装置。   The calculation unit extracts partial motion information corresponding to a component of a force direction generated by a cell at each position of the myocardium from the motion information, and applies the model to the partial motion information at the plurality of positions. The medical information processing apparatus according to any one of claims 1 to 3, wherein electric signal propagation at the plurality of positions is calculated by: 前記算出部は、前記心筋の各位置における細胞の繊維方向に基づいて、前記部分動き情報を抽出する、請求項4に記載の医用情報処理装置。   The medical information processing apparatus according to claim 4, wherein the calculation unit extracts the partial motion information based on a fiber direction of a cell at each position of the myocardium. 前記算出部は、前記動き情報から前記心筋の曲面に垂直な成分を除外することにより、前記部分動き情報を抽出する、請求項4に記載の医用情報処理装置。   The medical information processing apparatus according to claim 4, wherein the calculation unit extracts the partial motion information by excluding a component perpendicular to the curved surface of the myocardium from the motion information. 前記取得部は、更に、被検体の心電図を取得し、
前記算出部は、前記複数位置における電気信号伝搬を、前記心電図が示す電気信号の伝搬に対応するように算出する、請求項1乃至6のいずれか一項に記載の医用情報処理装置。
The acquisition unit further acquires an electrocardiogram of the subject,
The medical information processing apparatus according to any one of claims 1 to 6, wherein the calculation unit calculates electric signal propagation at the plurality of positions so as to correspond to electric signal propagation indicated by the electrocardiogram.
前記算出部は、算出した前記複数位置における電気信号伝搬について、前記心筋における電気信号伝搬の速度及び連続性の少なくとも一方に関する制約条件を満たすか否かを判定し、当該制約条件を満たさない場合には前記複数位置における電気信号伝搬を再算出する、請求項1乃至7のいずれか一項に記載の医用情報処理装置。   The calculation unit determines whether or not the calculated electric signal propagation at the plurality of positions satisfies a constraint condition regarding at least one of speed and continuity of the electric signal propagation in the myocardium, and when the constraint condition is not satisfied The medical information processing apparatus according to claim 1, which recalculates electric signal propagation at the plurality of positions. 前記算出部は、算出した前記複数位置における電気信号伝搬のうち、前記心筋における電気信号伝搬の速度及び連続性の少なくとも一方に関する制約条件を満たさない部分を除外する、請求項1乃至7のいずれか一項に記載の医用情報処理装置。   8. The calculation unit according to claim 1, wherein the calculation unit excludes a portion of the calculated electric signal propagation at the plurality of positions that does not satisfy a restriction condition regarding at least one of speed and continuity of electric signal propagation in the myocardium. The medical information processing apparatus according to one item. 前記算出部は、被検体に応じた前記モデルを使用して、前記複数位置における電気信号伝搬を算出する、請求項1乃至9のいずれか一項に記載の医用情報処理装置。   The medical information processing apparatus according to claim 1, wherein the calculation unit calculates electric signal propagation at the plurality of positions using the model corresponding to a subject. 前記算出部は、被検体のイオン濃度に応じた前記モデルを使用して、前記複数位置における電気信号伝搬を算出する、請求項10に記載の医用情報処理装置。   The medical information processing apparatus according to claim 10, wherein the calculation unit calculates electric signal propagation at the plurality of positions using the model according to an ion concentration of a subject. 前記算出部は、被検体の血圧に応じた前記モデルを使用して、前記複数位置における電気信号伝搬を算出する、請求項10又は11に記載の医用情報処理装置。   The medical information processing apparatus according to claim 10, wherein the calculation unit calculates electric signal propagation at the plurality of positions using the model according to a blood pressure of a subject. 前記算出部は、前記心筋の病変分布に応じた前記モデルを使用して、前記複数位置における電気信号伝搬を算出する、請求項10乃至12のいずれか一項に記載の医用情報処理装置。   The medical information processing apparatus according to any one of claims 10 to 12, wherein the calculation unit calculates electric signal propagation at the plurality of positions using the model according to a lesion distribution of the myocardium. 前記取得部は、更に、前記心筋に対して病変分布を設定する操作を受け付け、
前記算出部は、前記病変分布に応じた前記モデルを使用して、前記複数位置における電気信号伝搬を算出する、請求項1乃至12のいずれか一項に記載の医用情報処理装置。
The acquisition unit further accepts an operation for setting a lesion distribution for the myocardium,
The medical information processing apparatus according to any one of Claims 1 to 12, wherein the calculation unit calculates electric signal propagation at the plurality of positions using the model according to the lesion distribution.
前記取得部は、複数の心拍について前記動き情報を抽出し、
前記算出部は、心拍ごとに、前記複数位置における電気信号伝搬を算出する、請求項1乃至14のいずれか一項に記載の医用情報処理装置。
The acquisition unit extracts the motion information for a plurality of heartbeats,
The medical information processing apparatus according to any one of claims 1 to 14, wherein the calculation unit calculates electric signal propagation at the plurality of positions for each heartbeat.
前記算出部が算出した前記複数位置における電気信号伝搬を表示する表示制御部を更に備える、請求項1乃至15のいずれか一項に記載の医用情報処理装置。   The medical information processing apparatus according to claim 1, further comprising a display control unit that displays electric signal propagation at the plurality of positions calculated by the calculation unit. 前記取得部は、複数の被検体から前記動き情報を抽出し、
前記算出部は、被検体ごとに、前記複数位置における電気信号伝搬を算出し、
前記表示制御部は、被検体ごとに算出された前記複数位置における電気信号伝搬を比較可能に表示する、請求項16に記載の医用情報処理装置。
The acquisition unit extracts the motion information from a plurality of subjects,
The calculation unit calculates electrical signal propagation at the plurality of positions for each subject,
The medical information processing apparatus according to claim 16, wherein the display control unit displays the electric signal propagation at the plurality of positions calculated for each subject in a comparable manner.
表示制御部を更に備え、
前記算出部は、算出した前記複数位置における電気信号伝搬に基づいて前記心筋の動き情報を算出し、
前記表示制御部は、前記算出部が算出した動き情報と、前記4次元画像データとを比較可能に表示する、請求項1乃至17のいずれか一項に記載の医用情報処理装置。
A display control unit;
The calculation unit calculates the motion information of the myocardium based on the calculated electrical signal propagation at the plurality of positions,
The medical information processing apparatus according to any one of claims 1 to 17, wherein the display control unit displays the motion information calculated by the calculation unit and the four-dimensional image data so as to be comparable.
対象部位の4次元画像データの複数位置における前記対象部位の動き情報を取得する取得部と、
細胞変形と電気信号伝搬との関連を示すモデルを、前記複数位置における前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する算出部と、
を備える、医用情報処理装置。
An acquisition unit for acquiring movement information of the target part at a plurality of positions of the four-dimensional image data of the target part;
A calculation unit that calculates electric signal propagation at the plurality of positions by applying a model indicating the relationship between cell deformation and electric signal propagation to the motion information at the plurality of positions;
A medical information processing apparatus comprising:
細胞変形と電気信号伝搬との関連を示すモデルを記憶する記憶部と、
心筋の3次元構造を時系列で撮像することにより4次元画像データを収集する収集部と、
前記4次元画像データの複数位置における前記心筋の動き情報を抽出する取得部と、
前記モデルを前記複数位置における前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する算出部と、
を備える、医用画像診断装置。
A storage unit for storing a model indicating a relationship between cell deformation and electric signal propagation;
A collection unit that collects four-dimensional image data by imaging the three-dimensional structure of the myocardium in time series;
An acquisition unit for extracting movement information of the myocardium at a plurality of positions of the four-dimensional image data;
A calculation unit that calculates electric signal propagation at the plurality of positions by applying the model to the motion information at the plurality of positions;
A medical image diagnostic apparatus comprising:
心筋の3次元構造を時系列で撮像することにより得られた4次元画像データを取得し、取得した4次元画像データの複数位置における前記心筋の動き情報を抽出し、
細胞変形と電気信号伝搬との関連を示すモデルを前記複数位置における前記動き情報に適用することにより、前記複数位置における電気信号伝搬を算出する、
各処理をコンピュータに実行させる、プログラム。
Acquiring 4D image data obtained by imaging the 3D structure of the myocardium in time series, extracting motion information of the myocardium at a plurality of positions of the acquired 4D image data;
Calculating an electrical signal propagation at the plurality of positions by applying a model indicating the relationship between cell deformation and electrical signal propagation to the motion information at the plurality of positions;
A program that causes a computer to execute each process.
JP2018142159A 2017-08-31 2018-07-30 Medical information processing device, medical image diagnostic device and program Active JP7144230B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US16/119,527 US11484362B2 (en) 2017-08-31 2018-08-31 Medical information processing apparatus and medical information processing method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017167657 2017-08-31
JP2017167657 2017-08-31

Publications (2)

Publication Number Publication Date
JP2019042497A true JP2019042497A (en) 2019-03-22
JP7144230B2 JP7144230B2 (en) 2022-09-29

Family

ID=65813335

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018142159A Active JP7144230B2 (en) 2017-08-31 2018-07-30 Medical information processing device, medical image diagnostic device and program

Country Status (1)

Country Link
JP (1) JP7144230B2 (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0663026A (en) * 1992-08-21 1994-03-08 Nippon Koden Corp Simulator of electric phenomenon of heart
JP2000509618A (en) * 1996-05-06 2000-08-02 アリズミア リサーチ テクノロジー,インコーポレイテッド Prediction of cardiac arrhythmia by detecting weak potential
US20080177192A1 (en) * 2007-01-18 2008-07-24 General Electric Company Determination of cellular electrical potentials
JP2012183200A (en) * 2011-03-04 2012-09-27 Toshiba Corp Medical image processor and medical image processing program
US20140122048A1 (en) * 2012-10-30 2014-05-01 The Johns Hopkins University System and method for personalized cardiac arrhythmia risk assessment by simulating arrhythmia inducibility
JP2014512201A (en) * 2011-02-11 2014-05-22 ザ・ジョンズ・ホプキンス・ユニバーシティー System and method for planning cardiac surgery by patient
US20150302604A1 (en) * 2014-04-17 2015-10-22 Samsung Medison Co., Ltd. Medical imaging apparatus and method of operating the same
JP2015226711A (en) * 2014-06-02 2015-12-17 株式会社東芝 Medical image processor

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0663026A (en) * 1992-08-21 1994-03-08 Nippon Koden Corp Simulator of electric phenomenon of heart
JP2000509618A (en) * 1996-05-06 2000-08-02 アリズミア リサーチ テクノロジー,インコーポレイテッド Prediction of cardiac arrhythmia by detecting weak potential
US20080177192A1 (en) * 2007-01-18 2008-07-24 General Electric Company Determination of cellular electrical potentials
JP2014512201A (en) * 2011-02-11 2014-05-22 ザ・ジョンズ・ホプキンス・ユニバーシティー System and method for planning cardiac surgery by patient
JP2012183200A (en) * 2011-03-04 2012-09-27 Toshiba Corp Medical image processor and medical image processing program
US20140122048A1 (en) * 2012-10-30 2014-05-01 The Johns Hopkins University System and method for personalized cardiac arrhythmia risk assessment by simulating arrhythmia inducibility
US20150302604A1 (en) * 2014-04-17 2015-10-22 Samsung Medison Co., Ltd. Medical imaging apparatus and method of operating the same
JP2015226711A (en) * 2014-06-02 2015-12-17 株式会社東芝 Medical image processor

Also Published As

Publication number Publication date
JP7144230B2 (en) 2022-09-29

Similar Documents

Publication Publication Date Title
JP5562533B2 (en) X-ray computed tomography apparatus and tomography method
JP6953187B2 (en) Medical image diagnostic equipment and medical image processing equipment
JP2018083056A (en) Medical image processing apparatus, medical image processing method, and x-ray ct apparatus
US8798714B2 (en) Medical apparatus installation, and method for controlling a medical apparatus
JP2014087635A (en) Image processor and x-ray ct device
KR20170060698A (en) Computed tomography apparatus and control method for the same
JP2012130648A (en) Image processing apparatus and image processing method
US10561390B2 (en) Dose-reduced CT scan using dynamic collimation
US10610184B2 (en) Medical-information processing apparatus and X-ray CT apparatus
JP6933498B2 (en) Medical information processing equipment, X-ray CT equipment and medical information processing program
JP7032111B2 (en) Medical image processing equipment, X-ray CT equipment and medical image processing program
JP6981807B2 (en) Medical information processing device, X-ray CT device, medical information processing program, medical information processing method and medical information processing system
JP6466057B2 (en) Medical diagnostic imaging equipment
JP7313818B2 (en) MEDICAL IMAGE PROCESSING APPARATUS, MEDICAL IMAGE DIAGNOSTIC APPARATUS, AND MEDICAL IMAGE PROCESSING METHOD
JP4303947B2 (en) X-ray CT apparatus and image processing apparatus
JP7024954B2 (en) Medical image processing equipment and X-ray CT equipment
JP7144230B2 (en) Medical information processing device, medical image diagnostic device and program
US20210196221A1 (en) Medical image processing apparatus, recording medium, and system
US11484362B2 (en) Medical information processing apparatus and medical information processing method
JP5159299B2 (en) X-ray CT apparatus and medical image processing apparatus
US10842446B2 (en) Medical information processing apparatus, X-ray CT apparatus, and medical information processing method
JP2020096757A (en) Medical image processing device, medical image diagnostic device and image processing program
WO2013187461A1 (en) X-ray ct apparatus and image reconstruction method
JP6068177B2 (en) Medical image diagnostic apparatus, medical image processing apparatus, and medical image processing method
JP7139156B2 (en) X-ray CT device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210531

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220502

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220517

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220816

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220915

R150 Certificate of patent or registration of utility model

Ref document number: 7144230

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150