JP5722984B2 - Medical image processing apparatus and medical image processing program - Google Patents
Medical image processing apparatus and medical image processing program Download PDFInfo
- Publication number
- JP5722984B2 JP5722984B2 JP2013252829A JP2013252829A JP5722984B2 JP 5722984 B2 JP5722984 B2 JP 5722984B2 JP 2013252829 A JP2013252829 A JP 2013252829A JP 2013252829 A JP2013252829 A JP 2013252829A JP 5722984 B2 JP5722984 B2 JP 5722984B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- volume data
- medical image
- image processing
- time
- 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.)
- Active
Links
- 238000012545 processing Methods 0.000 title claims description 179
- 238000004458 analytical method Methods 0.000 claims description 104
- 238000012937 correction Methods 0.000 claims description 84
- 238000000034 method Methods 0.000 claims description 45
- 230000035772 mutation Effects 0.000 claims description 7
- 238000002591 computed tomography Methods 0.000 description 29
- 238000003384 imaging method Methods 0.000 description 19
- 101100125299 Agrobacterium rhizogenes aux2 gene Proteins 0.000 description 15
- 238000010586 diagram Methods 0.000 description 12
- 230000002861 ventricular Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 8
- 239000002872 contrast media Substances 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 7
- 210000000056 organ Anatomy 0.000 description 7
- 230000002123 temporal effect Effects 0.000 description 7
- 238000002595 magnetic resonance imaging Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000009877 rendering Methods 0.000 description 5
- 230000010412 perfusion Effects 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 206010019280 Heart failures Diseases 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 230000032683 aging Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000004217 heart function Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明は、医療用画像処理装置及び医療用画像処理プログラムに関する。 The present invention relates to a medical image processing apparatus and a medical image processing program.
Ｘ線ＣＴ装置およびＭＲＩ装置等の各種画像診断機器は、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行うものである。例えば、心臓や血管等の循環器系およびその他の動きのある臓器の場合、それらを構成する組織の動きを断層像により観察して、それら臓器等の機能を診断することが行なわれている。このような画像診断機器により、高解像度の画像が短時間に得られるようになったことに伴い、取得したボリュームデータから、臓器や腫瘍等の関心領域（Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ）を抽出し、これらを見やすく可視化し、あるいは定量化して面積や体積を量ることにより病変部の診断に役立てている。 Various diagnostic imaging apparatuses such as an X-ray CT apparatus and an MRI apparatus generate a tomographic image of an examination site from volume data of a subject and perform diagnosis by displaying it on a monitor screen. For example, in the case of a circulatory system such as the heart and blood vessels and other organs with movement, the movement of tissues constituting them is observed with a tomographic image to diagnose the function of these organs and the like. With such diagnostic imaging equipment, high-resolution images can be obtained in a short period of time, and regions of interest (Region of Interest) such as organs and tumors are extracted from the acquired volume data. It is useful for diagnosis of lesions by visualizing or quantifying for easy viewing and measuring the area and volume.
従来例１として、特許文献１に開示されている医療用画像処理装置は、各フェイズの複数のボリュームデータから新たに一の画像を生成し、各フェイズ現れる組織の変化を一の画像に表示することができる。しかし、各フェイズ間を補うボリュームデータを生成しているものではない。ここで、フェイズ（Ｐｈａｓｅ）とは、同一の対象を同一手法で短時間のうちに撮像した複数のボリュームデータのうちの１つのボリュームデータを指す。
As a conventional example 1, the medical image processing apparatus disclosed in
従来例２として、特許文献２に開示されている医療用画像処理装置は、各フェイズの抽出領域から構造情報（たとえば、血管の中心線）を求めて、各フェイズの異常個所を比較、検出している。そのため、少ない計算量で領域抽出を行うことができる。しかし、各フェイズ間を補うボリュームデータを生成し、各フェイズの領域抽出を行っていない。
As Conventional Example 2, the medical image processing apparatus disclosed in
従来例３として、非特許文献１では、いわゆる、モーフィングにより、医療画像とは異なる通常の画像で構成される動画像のフレーム間隔（例えば、１／６０秒）の中間位置で新たに画像を補間している。
As Conventional Example 3, in
しかし、従来例１では、複数のボリュームデータから高画質な一の画像を生成しているが、複数のボリュームデータそれぞれから作成する個別の画像の画質は改善しない。これは特に動画となったときに顕著である。また、従来例２は複数のボリュームデータに含まれる組織の対応関係を取得するに留まる。また、従来例１及び従来例２で開示されている技術では、各フェイズ間を補うボリュームデータを生成しているものではないので、各フェイズの位置合わせの精度を高めるには別の方法を検討する必要がある。また、従来例３では、新たに画像を補間するタイミングは、動画像（たとえば、テレビ映像）のフレーム間隔（例えば、１／６０秒）の中間位置に固定されており、この位置を変更することはまでは示唆されていない。さらに、新たに画像を補間しても画質を改善しているわけではない。 However, in Conventional Example 1, one image with high image quality is generated from a plurality of volume data, but the image quality of individual images created from each of the plurality of volume data is not improved. This is particularly noticeable when it comes to moving images. Further, Conventional Example 2 only acquires the correspondence relationship between the organizations included in the plurality of volume data. In addition, the techniques disclosed in Conventional Example 1 and Conventional Example 2 do not generate volume data that compensates for each phase, so another method is considered to improve the accuracy of positioning in each phase. There is a need to. Further, in Conventional Example 3, the timing for newly interpolating an image is fixed at an intermediate position of a frame interval (for example, 1/60 second) of a moving image (for example, a television image), and this position is changed. Until is not suggested. Furthermore, even if a new image is interpolated, the image quality is not improved.
本発明の目的は、動き解析情報及び画像の画質情報に基づき、画像の変化を正確に表示することができる医療用画像処理装置及び医療用画像処理プログラムを提供することである。 An object of the present invention is to provide a medical image processing apparatus and a medical image processing program capable of accurately displaying image changes based on motion analysis information and image quality information.
本発明は、時刻情報を有する複数のボリュームデータのそれぞれを可視化する医療用画像処理装置であって、前記複数のボリュームデータを生成し、それぞれの画質情報を生成する第１生成部と、前記複数のボリュームデータのそれぞれの動き解析情報を生成する第２生成部と、前記複数のボリュームデータの内の少なくともひとつのボリュームデータについて、当該ボリュームデータのボクセルと、当該ボリュームデータとは異なる他のボリュームデータのボクセルとを動き解析情報により対応付け、前記対応付けしたボクセルの信頼度を前記画質情報に基づいて決定し、前記対応付けしたボクセルのボクセル値を、前記対応付けしたボクセルの信頼度に応じて変わるパラメータのフィルタを用いて補正する補正処理部と、前記補正処理部により補正したボクセル値に基づき、前記複数のボリュームデータを可視化した複数の画像を表示するための処理を行う表示処理部と、を備える医療用画像処理装置を提供する。 The present invention is a medical image processing apparatus that visualizes each of a plurality of volume data having time information, the first generation unit generating the plurality of volume data and generating the respective image quality information, and the plurality of the volume data A second generation unit that generates motion analysis information for each of the volume data, and for at least one of the plurality of volume data, voxels of the volume data and other volume data different from the volume data Is associated with the motion analysis information, the reliability of the associated voxel is determined based on the image quality information, and the voxel value of the associated voxel is determined according to the reliability of the associated voxel. Correction processing unit that corrects using a parameter filter that changes, and the correction processing unit Based on a more corrected voxel values, and a display processing unit that performs processing for the plurality of volume data for displaying a plurality of images visualized, to provide a medical image processing apparatus comprising a.
上記医療用画像処理装置では、更に、前記ボリュームデータのボクセルと、当該ボリュームデータとは異なる他のボリュームデータのボクセルとの間を補間する補間ボリュームデータを生成する第３生成部を備える。 The medical image processing apparatus further includes a third generation unit that generates interpolated volume data for interpolating between the voxel of the volume data and a voxel of other volume data different from the volume data .
上記医療用画像処理装置では、前記補正処理部は、更に、前記動き解析情報から、前記複数の画像のうちの少なくとも一つの画像において、部分的な変異又は変形を取得し、前記変異又は変形から、前記補間ボリュームデータのボクセル値の変化を予測し、当該予測したボクセル値の変化に基づき、前記フィルタを補正する。 In the medical image processing apparatus, the correction processing unit further acquires a partial variation or deformation in at least one of the plurality of images from the motion analysis information, and from the variation or deformation A change in the voxel value of the interpolated volume data is predicted, and the filter is corrected based on the predicted change in the voxel value .
上記医療用画像処理装置では、前記第３生成部は、前記ボリュームデータと前記補間ボリュームデータの間隔と、前記補間ボリュームデータと前記ボリュームデータとは異なる前記他のボリュームデータの間隔とは、非等間隔となるように、前記補間ボリュームデータを生成する。 In the medical image processing apparatus, the third generation unit is configured such that an interval between the volume data and the interpolated volume data and an interval between the interpolated volume data and the other volume data different from the volume data are unequal. The interpolated volume data is generated so as to have an interval.
また、上記医療用画像処理装置では、前記第１生成部は前記画質情報を医療画像撮影装置から取得した情報を少なくとも用いて生成する。 In the medical image processing apparatus, the first generation unit generates the image quality information using at least information acquired from a medical image photographing apparatus.
また、本発明は、医療用画像処理プログラムであって、コンピュータを、上記医療用画像処理装置として機能させる医療用画像処理プログラムを提供する。 The present invention also provides a medical image processing program for causing a computer to function as the medical image processing apparatus.
本発明に係る医療用画像処理装置及び医療用画像処理プログラムによれば、動き解析情報及び画像の画質情報に基づき、複数のボリュームデータからそれぞれ画像を作成するときの画質を向上させることができる。 According to the medical image processing apparatus and the medical image processing program according to the present invention, it is possible to improve image quality when an image is created from each of a plurality of volume data based on motion analysis information and image quality information.
以下、本発明の実施形態について、図面を参照して説明する。 Hereinafter, embodiments of the present invention will be described with reference to the drawings.
図１を参照して、実施の形態１に係る医療用画像処理装置について説明する。図１に示す医療用画像処理装置１００は、ボリュームデータ生成部１０１と、動き解析情報生成部１０３と、補間ボリュームデータ生成部１０５と、画像生成部１０７と、表示処理部１０９と、を備える。医療用画像処理装置１００には、医療用画像処理装置１００で生成された画像を表示する表示部１１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部１１３と、が設けられている。
With reference to FIG. 1, the medical image processing apparatus according to the first embodiment will be described. A medical image processing apparatus 100 illustrated in FIG. 1 includes a volume
図１に示すように、医療用画像処理装置１００は、コンピュータ断層撮影装置に接続されている。コンピュータ断層撮影装置は、被検体の組織等を可視化するものである。Ｘ線源１からは同図に鎖線で示す縁部ビームを有するピラミッド状のＸ線ビーム束２が放射される。Ｘ線ビーム束２は、例えば患者３である被検体を透過しＸ線検出器４に照射される。Ｘ線源１及びＸ線検出器４は、本実施形態の場合にはリング状のガントリー５に互いに対向配置されている。リング状のガントリー５は、このガントリーの中心点を通るシステム軸線６に対して、同図に示されていない保持装置に回転可能（矢印ａ参照）に支持されている。
As shown in FIG. 1, the medical image processing apparatus 100 is connected to a computed tomography apparatus. The computer tomography apparatus visualizes the tissue of a subject. The
患者３は、Ｘ線が透過するテーブル７上に寝ている。このテーブル７は、図示されていない支持装置によりシステム軸線６に沿って移動可能（矢印ｂ参照）に支持されている。従って、Ｘ線源１及びＸ線検出器４は、システム軸線６に対して回転可能でありかつシステム軸線６に沿って患者３に対して相対的に移動可能である測定システムを構成するので、システム軸線６に関して種々の投影角及び種々の位置のもとで患者３にＸ線を投射する。その際に発生するＸ線検出器４の出力信号は、ボリュームデータ生成部１０１に供給され、ボリュームデータに変換される。
The patient 3 lies on the table 7 that transmits X-rays. The table 7 is supported by a support device (not shown) so as to be movable along the system axis 6 (see arrow b). Thus, the
シーケンス走査の場合、患者３の層毎の走査が行なわれる。その際に、Ｘ線源１及びＸ線検出器４はシステム軸線６を中心に患者３の周りを回転し、Ｘ線源１及びＸ線検出器４を含む測定システムは患者３の２次元断層を走査するために多数の投影を撮影する。その際に取得された測定値から、走査された断層を表示する断層像が再構成される。相連続する断層の走査の間に、患者３はその都度システム軸線６に沿って移動される。この過程は全ての関心断層が捕捉されるまで繰り返される。
In the case of a sequence scan, a scan for each layer of the patient 3 is performed. At that time, the
スパイラル走査中は、Ｘ線源１及びＸ線検出器４を含む測定システムはシステム軸線６を中心に回転し、テーブル７は連続的に矢印ｂの方向に移動する。すなわち、Ｘ線源１及びＸ線検出器４を含む測定システムは、患者３に対して相対的に連続的にスパイラル軌道上を、患者３の関心領域が全部捕捉されるまで移動する。本実施形態の場合、同図に示されたコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給される。
During spiral scanning, the measurement system including the
ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部１０３及び画像生成部１０７へ出力する。なお、ボリュームデータ生成部１０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータを保持することができる。
The volume
ここで、本実施の形態では、ボリュームデータ生成部１０１は、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
Here, in the present embodiment, the volume
以下、説明のため、ボリュームデータ生成部１０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。
なお、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
Hereinafter, for the sake of explanation, the volume
In the following, for explanation, the volume data of time information t k, described omitted volume data t k as appropriate.
動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成し（図２参照）、補間ボリュームデータ生成部１０５に出力する。
Of the volume data output from the volume
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。 Here, the motion analysis information refers to information on a correspondence relationship between corresponding positions or corresponding objects on a plurality of volume data, and information on a process in which the positions and objects move and change. In addition, each volume data pixel serves as an index indicating a position at an arbitrary time between time k-1 and time k .
図２は、動き解析情報生成部１０３の動作を概念的に示す図である。図２に示すように、動き解析情報生成部１０３は、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを、動き解析情報として生成している。
FIG. 2 is a diagram conceptually illustrating the operation of the motion analysis
補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、画像生成部１０７に出力する。
Interpolated volume data generation unit 105, based on the motion analysis information outputted from the motion
ここで、補間ボリュームデータとは、時刻ｔｋ−１と時刻ｔｋとの間の特定の時刻ｔａｕｘにおけるボリュームデータである。つまり、補間ボリュームデータ生成部１０５が補間ボリュームデータｔａｕｘを生成することは、実際はコンピュータ断層撮影装置で撮像されていない時刻ｔａｕｘにおいて、あたかもコンピュータ断層撮影装置により撮像して、ボリュームデータ生成部１０１がボリュームデータを生成したことと、同等となる。
Here, the interpolated volume data is volume data at a particular time t aux between time t k-1 and time t k. That is, the interpolated volume data generation unit 105 generates the interpolated volume data t aux at the time t aux that is not actually captured by the computed tomography apparatus as if it was captured by the computed tomography apparatus, and the volume
図３は、補間ボリュームデータ生成部１０５の動作を概念的に示す図である。図３に示すように、補間ボリュームデータ生成部１０５は、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。 FIG. 3 is a diagram conceptually illustrating the operation of the interpolation volume data generation unit 105. As shown in FIG. 3, the interpolation volume data generating unit 105 interpolates between the volume data t k-1 and the volume data t k including information of a three-dimensional image, the interpolation volume data t aux and its time information t aux .
ここで、本実施の形態では、図３に示すように補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、後述する。 In the present embodiment, the time information t aux with interpolation volume data t aux as shown in FIG. 3, not in the middle of the time information and the time information t k-1 and time information t k. In other words, the distance A between the time information t k-1 of the volume data t k-1 which is generated in the time information t aux and immediately before the timing of the interpolation volume data, the time information t aux and immediately after the timing of the interpolated volume data different spacing and distance B between the time information t k of the generated volume data to, i.e., in the non-equal interval. The reason for this will be described later.
画像生成部１０７は、ボリュームデータ生成部１０１から出力されるボリュームデータｔｋ−１、ｔｋ、及び時刻情報ｔｋ−１、ｔｋ、並びに補間ボリュームデータ生成部１０５から出力される補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。
The
表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部１０７で生成された３次元画像のシーケンスを生成して、表示部１１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
The
表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部１０７で生成された３次元画像のシーケンスを、表示処理部１０９での処理の後、動画像として再生することができる。
The
操作部１１３は、表示部１１１に表示された画像を加工するための各種のインターフェースを備えており、例えば、補間画像ｔａｕｘを生成するための補間パラメータ（時刻情報ｔａｕｘなど）を設定する。そして、ユーザは、表示部１１１を見ながら、操作部１１３により、医療用画像処理装置１００の各部を操作することができる。そのため、ユーザは、インタラクティブに補間パラメータを変更できる。したがって、本実施の形態の医療用画像処理装置において、画像生成部１０７は、操作部１１３から指示に基づき、インタラクティブに補間画像ｔａｕｘを修正することができる。
The
ここで、補間パラメータとして、たとえば、（１）位置補間アルゴリズム（リニア／キュービック／スプライン）、（２）ボクセル値の補間のアルゴリズム（リニア／キュービック／スプライン）、（３）補間により生成される３Ｄ画像のピッチ（時間間隔）や数、並びに（４）４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。 Here, as interpolation parameters, for example, (1) position interpolation algorithm (linear / cubic / spline), (2) voxel value interpolation algorithm (linear / cubic / spline), (3) 3D image generated by interpolation (4) 4D filter methods include smoothing, edge enhancement, maximum / minimum value projection, difference, cumulative addition, histogram matching, and combinations thereof.
図４を参照して、画像生成部１０７が生成する３次元画像について説明する。図４は、画像生成部１０７が生成した各画像を模式的に示す図である。図４に示すように、実画像ｔｋ−１と実画像ｔｋとの間で、時刻情報ｔａｕｘを有する補間画像ｔａｕｘが、画像生成部１０７によって生成される。以下、説明のためボリュームデータｔｋ−１、ｔｋから生成された、時刻ｔｋ−１、ｔｋにおける実際の３次元画像を、それぞれ実画像ｔｋ−１と実画像ｔｋと称し、補間ボリュームデータｔａｕｘから生成された、時刻ｔａｕｘにおける３次元の補間画像を、補間画像ｔａｕｘと称する。
A three-dimensional image generated by the
ここで、図４に示すように、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔となっている。本実施の形態において、ボリュームデータｔｋ−１と、時刻情報ｔｋ−１の次のタイミングで生成されるボリュームデータｔｋとの間隔が、およそ０．１秒から数年単位である。これは、画像生成部１０７で生成される医療画像（実画像ｔｋ−１、実画像ｔｋ、及び補間画像ｔａｕｘ）とは異なる通常の画像で構成される動画像（たとえば、テレビなどの映像）のフレーム間隔（例えば、１／６０秒）と比べて非常に大きいことに起因する。
Here, as shown in FIG. 4, a distance C between the time information t k-1 of the generated time information t aux and immediately before the timing of the interpolated image t aux actual image t k-1, the interpolation image t aux time information t aux between the distance D between the time information t k of the actual image t k generated in the timing immediately after that in the non-equal interval. In this embodiment, the volume data t k-1, the interval between the volume data t k generated at the next timing of the time information t k-1 is the number yearly approximately 0.1 seconds. This is a moving image (for example, a television or the like) composed of a normal image different from the medical image (actual image t k−1 , actual image t k , and interpolation image t aux ) generated by the
つまり、本実施の形態では、実画像ｔｋ−１、と実画像ｔｋとの間の間隔（Ｃ＋Ｄ）が、例えば年単位であって、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経年変化を観察するのに、その中間位置での補間画像では十分ではないからである。そのため、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔とし、実画像ｔｋ−１、と実画像ｔｋとの間で被検体である患者３の検査部位（関心領域）の経年変化を中間位置とは異なる、実画像ｔｋ−１、と実画像ｔｋとの間の様々な位置で補間画像を生成することを可能としている。また、撮影に要する放射線量を最小にする目的で、関心の最も高い時刻では撮影の時間間隔を密にし、それ以外では疎に撮影することが行われている。また、心電図など被検体のリアルタイム情報を取得する別途の機器で指定されるタイミングで撮影することもある。これらによって撮影間隔が非等間隔となりうる。リアルタイム情報には他には例えば、造影剤の注入のタイミングもある。 That is, in this embodiment, the actual image t k-1, the distance between the actual image t k (C + D) is, for example, a yearly, the actual image t k-1 and the actual image t k This is because the interpolated image at the intermediate position is not sufficient for observing the secular change of the region of interest of the patient 3 which is the subject occurring between. Therefore, the distance C between the time information t k-1 of the generated time information t aux and immediately before the timing of the interpolated image t aux actual image t k-1, the interpolation image t aux time information t aux and immediately the spacing D between the time information t k of the actual image t k generated in timing is a non-equal interval, the testing area of the patient 3 is a subject between the actual image t k-1, and the actual image t k different from the aging of the intermediate position (ROI), it is made possible to generate an interpolated image at different positions between an actual image t k-1, and the actual image t k. For the purpose of minimizing the amount of radiation required for imaging, the imaging time interval is dense at the time of highest interest, and sparse imaging is performed at other times. In addition, imaging may be performed at a timing designated by a separate device that acquires real-time information of the subject such as an electrocardiogram. As a result, the photographing intervals can be non-uniform. Other real-time information includes, for example, the timing of contrast medium injection.
上述のように、本実施の形態では、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、非等間隔であるが、これに限らない。実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経時変化を観察するのに適したタイミングであれば、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、等間隔に設定されてもよい。すわなち、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の検査部位（関心領域）の経時変化を観察するのに適したタイミングであれば、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成しても良い。言い換えると、本実施の形態において、ボリュームデータｔａｕｘの時刻ｔａｕｘは、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。これにより、本実施の形態では、Ｘ線の線量を削減し、被検体である患者３への影響を極力抑えることができる。
As described above, in this embodiment, the distance C between the time information t k-1 of the real image t k-1 which is generated in the time information t aux and immediately before the timing of the interpolated image t aux, interpolated image t the spacing D between the time information t k of the actual image t k which is generated in the time information t aux and immediately after the timing of the aux is is a non-equidistant, not limited to this. If timing suitable for observing a change with time of a region of interest of the patient 3 is a subject that is generated between the actual image t k-1 and the actual image t k, the time information of the interpolation image t aux t aux and the distance C between the time information t k-1 of the real image t k-1 which is generated in the timing immediately before the interpolation image t aux time information t aux and the actual image t k generated in timing immediately after the spacing D between the time information t k of, may be set at equal intervals. Nachi Suwa, if a timing suitable for the observation of the time course of the testing area of the patient 3 is a subject that is generated between the actual image t k-1 and the actual image t k (region of interest) The
（任意の時刻の補間画像）
図５は、時刻ｔｋ−１と時刻ｔｋ−１との任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成する一例を模式的に示す図である。図５に示すように、画像生成部１０７は、実画像ｔｋ−１を生成した時刻０．０ｓｅｃと、実画像ｔｋを生成した時刻情報１．５ｓｅｃとの間の任意の時刻情報１．０ｓｅｃの補間画像ｔａｕｘを生成する。そのため、実画像を生成した時刻の間の任意の時刻で補間画像を生成することができ、実画像ｔｋ−１、と実画像ｔｋとの間で発生している患者３の診断範囲の経時変化を観察することができる。
(Interpolated image at any time)
FIG. 5 is a diagram schematically illustrating an example in which the
（補間画像とモデルとの比較）
図６は、心室容量の大きさの経時変化を示すモデルと画像生成部１０７が作成した補間画像ｔａｕｘとを比較した図である。上述のように、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の時刻であれば、任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成することが可能なので、たとえば図６に示すように、心室容量の大きさの経時変化を示すモデルの最大容量を示す時刻で、画像生成部１０７が補間画像ｔａｕｘを作成し、心室容量の大きさの経時変化を示すモデルと補間画像ｔａｕｘにあらわれる心室の大きさとを比較することができる。
(Comparison between interpolation image and model)
FIG. 6 is a diagram comparing a model showing temporal changes in the size of the ventricular volume and the interpolation image t aux created by the
心室容量の大きさの経時変化を示すモデルは、図６に示すように線形的に変化するものではないので、２つの実画像ｔｋ−１、ｔｋ間の中間位置では、心室容量は最大とならない、したがって、本実施の形態では、モデル上、心室容量が最大となる特徴を有する時刻で、実画像ｔｋ−１、ｔｋから補間画像ｔａｕｘを生成する。図６に示すように、モデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃである。したがって、画像生成部１０７は、時刻情報ｔａｕｘをモデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃとし、補間画像ｔａｕｘを生成する。
Since the model showing the temporal change in the magnitude of the ventricular volume does not change linearly as shown in FIG. 6, the ventricular volume is the maximum at the intermediate position between the two real images t k−1 and t k. Therefore, in the present embodiment, the interpolated image t aux is generated from the real images t k−1 and t k at the time having the characteristic that the ventricular volume is maximized on the model. As shown in FIG. 6, on the model, the time when the ventricular volume becomes maximum is t = 0.78 sec. Therefore, the
上述のように、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定される時刻を補間画像の時刻情報とすることで、画像生成部１０７は、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を画像として生成し、表示部１１１が可視化することができる。
As described above, the
なお、心室容量の大きさの経時変化を示すモデルから、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定しているが、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を推定する手段は当該モデルだけに限らない。例えば、被検体である患者３に対して造影剤などが注入された外的要因が発生した時刻を、画像上に特徴が現れる時刻であると推定することも可能できる。そのため、補間画像ｔａｕｘを作成し、実画像ｔｋ−１、ｔｋ、及び補間画像ｔａｕｘの一連の画像の変化から、被検体（患者３）の検査部位（関心領域）の外的要因による変化を観察することができる。また、時刻情報ｔａｕｘモデルのみならず外部機器からの入力によって定めることが出来る。これは心電図など被検体のリアルタイム情報を実画像の取得と同時に行い、これを参照することによって求めることが出来る。例えば、心電図上の波形の特徴的な時刻や、造影剤の注入を行ったタイミングなどが考えられる。 Although it is estimated from the model showing the temporal change in the size of the ventricular volume, the characteristics of the examination site (region of interest) of the subject to be observed (patient 3) appear, the subject to be observed ( The means for estimating the characteristics of the examination region (region of interest) of the patient 3) is not limited to the model. For example, it is possible to estimate that the time when an external factor in which a contrast medium or the like is injected into the patient 3 as the subject is the time when the feature appears on the image. Therefore, an interpolation image t aux is created, and external factors of the examination region (region of interest) of the subject (patient 3) are determined from a series of image changes of the real images t k−1 , t k and the interpolation image t aux. Can observe the change. Further, it can be determined not only by the time information aux model but also by input from an external device. This can be obtained by performing real-time information on the subject such as an electrocardiogram simultaneously with the acquisition of the real image and referring to it. For example, the characteristic time of the waveform on the electrocardiogram, the timing when the contrast medium is injected, and the like can be considered.
（異なる画像セットの同期）
また、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成し、複数の３次元画像セットから正確に同期した画像のシーケンスを生成することができる。図７は、画像生成部１０７が生成した２組の時系列に並んだ３次元画像セットを示す。また、図８は、図７に示す２組みの３次元画像セットを同期させた画像のシーケンスを示す図である。
(Synchronize different image sets)
In addition, the
図７に示す３次元画像セットＢでは、第１の時刻情報ｔｋ−１（＝０．２ｓｅｃ）を有する第１実画像ｔｋ−１、及び第２の時刻情報ｔｋ（＝１．７ｓｅｃ）を有する第２実画像ｔｋ、並びに、第１実画像ｔｋ−１及び第２実画像ｔｋから生成され、時刻情報ｔａｕｘ２を有する第２補間画像ｔａｕｘ２が時系列に並んでいる。また、図７に示す３次元画像セットＡでは、第３の時刻情報ｔｋ−１（＝０．０ｓｅｃ）を有する第３実画像ｔｋ−１、及び第４の時刻情報ｔｋ（＝１．５ｓｅｃ）を有する第４実画像ｔｋ、並びに第３実画像ｔｋ−１及び第４実画像ｔｋから生成され、時刻情報ｔａｕｘ１を有する第１補間画像ｔａｕｘ１が時系列に並んでいる。 In the three-dimensional image set B shown in FIG. 7, the first real image t k−1 having the first time information t k−1 (= 0.2 sec) and the second time information t k (= 1.7 sec). second actual image t k with), and is generated from the first actual image t k-1 and the second real image t k, the second interpolation image t aux2 with time information t aux2 are arranged in time series . In the three-dimensional image set A shown in FIG. 7, the third real image t k−1 having the third time information t k−1 (= 0.0 sec) and the fourth time information t k (= 1). fourth actual image t k with .5sec), and the third is generated from the actual image t k-1 and the fourth actual image t k, the first interpolated image t aux1 with time information t aux1 is arranged in chronological Yes.
そして、図８に示すように、画像生成部１０７は、各画像の保持する時刻情報に基づき、３次元画像セットＢの第２補間画像ｔａｕｘ２を有する時刻情報ｔａｕｘ２を、３次元画像セットＡに時間的に同期させて、新たな画像のシーケンスを生成する。
Then, as illustrated in FIG. 8, the
したがって、本実施の形態では、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成することができる。そして、表示処理部１０９は、複数の３次元画像セットから正確に時間的に同期した画像のシーケンスを生成することができる。
Therefore, in the present embodiment, the
また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、観察対象部位（心臓）のマルチフェイズデータでの同一フェイズ間のずれの補正も可能となる。また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、「コンピュータ断層撮影装置から取得した検査の画像取得時間情報」と「表示処理部１０９で生成される再生フレームレート」を元に、表示部１１１で正確な実時間で動画再生することができる。
In the present embodiment, if the
さらに、本実施の形態では、画像生成部１０７は、図７に示す一方の３次元画像セットＢにだけ存在する特定の取得時刻（例えば、ｔｋ＝１．７ｓｅｃ）の実画像ｔｋに合わせて、他方の３次元画像セットＡの第３実画像ｔｋ−１及び第４実画像ｔｋから新たな補間画像を作成することにより、より正確な３次元画像間の比較を可能にする。
Further, in this embodiment, the
図９に、画像生成部１０７が時間的に等間隔に複数の補間画像を実画像に補った場合の画像のシーケンスを示す。また、画像生成部１０７は、図９に示すように、同一の被検体から同一の検査部位を対象として、撮影間隔の異なるコンピュータ断層撮影装置により供給された複数のボリュームデータｔ１（＝０．００ｓｅｃ）、ｔ２（＝０．１７ｓｅｃ）、ｔ３（＝０．３０ｓｅｃ）から生成された複数の実画像ｔ１、ｔ２、ｔ３のうち実画像ｔ１、ｔ３に対して、時間的に等間隔に補う、時刻情報ｔ４（＝０．１０ｓｅｃ）を有する補間画像ｔａｕｘ４、及び時刻情報ｔ５（＝０．２０ｓｅｃ）を有する補間画像ｔａｕｘ５を生成することができる。そのため、本実施の形態では、時間的に等間隔でない実画像に対して、補間画像を等時間間隔に生成するので、表示部１１１が表示する医療用画像の表示や解析の質・能力を高めることができる。
FIG. 9 shows a sequence of images when the
図１０に、実画像に対して一時的に複数の補間画像補い、表示部１１１に表示した例を示す。画像生成部１０７は、図１０に示すように、表示部１１１に３次元画像のシーケンスを表示している時のみ、一時的に、画像生成部１０７が複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を生成し、各フレームごとに、レンダリング後に複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を破棄することで、メモリの使用量を削減でき、本実施の形態の医療用画像処理装置１００への負荷を軽減することができる。
FIG. 10 shows an example in which a plurality of interpolated images are temporarily supplemented with respect to the actual image and displayed on the
次に、図１１を参照して、本実施の形態に係る医療用画像処理装置１００の動作処理フローを説明する。図１１は、医療用画像処理装置１００の動作処理フローである。 Next, an operation processing flow of the medical image processing apparatus 100 according to the present embodiment will be described with reference to FIG. FIG. 11 is an operation processing flow of the medical image processing apparatus 100.
まず、ステップＳＴ１で、ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータ、及び各ボリュームデータの時刻情報ｔｎを生成する。そして、ステップＳＴ２へ遷移する。
First, in step ST1, the volume
次に、ステップＳＴ２で、動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。そして、ステップＳＴ３へ遷移する。
Next, in step ST2, the motion
次に、ステップＳＴ３で、補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。そして、ステップＳＴ４へ遷移する。
Next, in step ST3, the interpolated volume data generation unit 105, based on the motion analysis information outputted from the motion
次に、ステップＳＴ４で、画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。そして、ステップＳＴ５へ遷移する。
Next, in step ST4, the
次に、ステップＳＴ５で、表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。そして、ステップＳＴ６へ遷移する。
Next, in step ST < b> 5, the
最後に、ステップＳＴ６で、表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。そして、動作処理が終了する。
Finally, in step ST6, the
以上、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置１００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。 As described above, according to the medical image processing apparatus 100 according to the present embodiment, it is possible to accurately display a change in a specific region in an image by interpolating an interpolated image between real images based on motion analysis information. it can. Therefore, according to the medical image processing apparatus 100 according to the present embodiment, it is possible to accurately observe the examination region (region of interest) of the patient 3 as the subject.
例えば、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。 For example, if the medical image processing apparatus 100 interpolates an interpolated image between real images based on motion analysis information, it is possible to satisfactorily observe a cardiac function failure or the like only by performing cine reproduction using a conventional rendering method. Further, if the medical image processing apparatus 100 interpolates interpolated images between real images based on the motion analysis information, it is possible to observe not only those with short time intervals such as cardiac multiphase data but also aging of lesions. It can also be applied to long time interval data.
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、セグメンテーションにより隣接する関心領域（臓器）を設定し、動き解析情報により複数の臓器の動きが同一の部分を推定、臓器間癒着の可能性あり、として表示部１１１に表示することができる。
Furthermore, when the medical image processing apparatus 100 interpolates interpolated images between real images based on motion analysis information, it is possible to perform reproducible fully automated quantification in strain analysis, perfusion analysis, cardiac function analysis, and the like. Anyone can get accurate results at a lower price than before. Furthermore, when the medical image processing apparatus 100 interpolates interpolated images between real images based on motion analysis information, adjacent regions of interest (organs) are set by segmentation, and movements of a plurality of organs are determined by motion analysis information. Can be displayed on the
また、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間すれば、表示部１１１のノイズ軽減、ノイズのばらつき調整が可能である。
Further, according to the medical image processing apparatus 100 according to the present embodiment, if the interpolated image is interpolated between the actual images based on the motion analysis information, the noise of the
例えば、医療用画像処理装置１００が生成した動き解析情報によって、実画像及び補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置１００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。 For example, the noise can be reduced without losing the original detailed image information by specifying the corresponding pixels of the actual image and the interpolated image based on the motion analysis information generated by the medical image processing apparatus 100. Further, even if a phase with a large amount of X-ray irradiation and a phase with a small amount of X-ray irradiation are mixed in the CT multi-phase examination by the medical image processing apparatus 100, the pixel values of the previous and subsequent phases are used for the phase with a small amount of X-ray irradiation. By doing so, a good interpolation image can be obtained.
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。 Furthermore, if the medical image processing apparatus 100 interpolates an interpolated image between real images based on motion analysis information, the X-ray irradiation amount is increased only in an important phase, and a real image is generated. By generating an interpolated image, it is possible to reduce the amount of X-ray exposure of the patient 3.
なお、本実施の形態に係る医療用画像処理装置１００では、図１に示すコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給されているが、これに限らない。ボリュームデータ生成部１０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
In the medical image processing apparatus 100 according to the present embodiment, a large number of continuous tomographic signals in the diagnosis range of the patient 3 are supplied to the volume
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。 In this embodiment, the time t k of the volume data t k is the time after a lapse from the volume data t k-1 a predetermined time, the predetermined time is arbitrary time, for example, the subject It can be set as appropriate depending on the time when a change in the examination site (region of interest) of (patient 3) is predicted to occur.
なお、本実施の形態において、表示処理部１０９が、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部１１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
In the present embodiment, the
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻ｔに基づき、生成しているが、これに限らない。 In the present embodiment, the time information at which the volume data is generated is generated based on the time t supplied by the computed tomography apparatus, but the present invention is not limited to this.
なお、本実施の形態の説明に用いた各機能ブロック（図１参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。 Each functional block (see FIG. 1) used in the description of the present embodiment is typically realized as an LSI that is an integrated circuit. In particular, calculation processing for generating an actual image and an interpolation image can be performed by a GPU (Graphic Processing Unit). The GPU is an arithmetic processing unit designed specifically for image processing as compared with a general-purpose CPU (Central Processing Unit), and is mounted on a computer separately from a normal CPU.
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。 In the present embodiment, the example used for observing the heart is shown, but the present invention is not limited to this. For example, it can be used for observation of changes in tumor tissue over time, observation of blood flow perfusion, and the like.
（実施の形態２）
上述した実施の形態１に係る医療用画像処理装置１００によれば、動き解析情報に基づき、補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示する。しかし、補間画像を実画像間に補間する場合、実画像の画質によっては、単純に補間画像を補間するだけでは、画像内の特定の領域の変化を正確に表示できない場合もある。
(Embodiment 2)
According to the medical image processing apparatus 100 according to the first embodiment described above, a change in a specific region in an image is accurately displayed by interpolating an interpolated image between real images based on the motion analysis information. However, when interpolating an interpolated image between real images, depending on the quality of the real image, it may not be possible to accurately display a change in a specific area in the image simply by interpolating the interpolated image.
そこで、実施の形態２に係る医療用画像処理装置３００では、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
Therefore, in the medical image processing apparatus 300 according to
実施の形態２に係る医療用画像処理装置３００では、必ずしも補間画像を生成する必要はない。動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質を補正して、２以上の実画像内の特定の領域の変化を正確に表示することができれば、補間画像を生成する必要はない。
In the medical image processing apparatus 300 according to
図１２を参照して、本実施の形態２に係る医療用画像処理装置３００の構成について説明する。図１２に示す医療用画像処理装置３００は、ボリュームデータ生成部３０１と、動き解析情報生成部３０３と、補間ボリュームデータ生成部３０５と、補正処理部３０６と、画像生成部３０７と、表示処理部３０９と、を備える。さらに、医療用画像処理装置３００には、医療用画像処理装置３００で生成された画像を表示する表示部３１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部３１３と、が設けられている。
With reference to FIG. 12, the structure of the medical image processing apparatus 300 according to the second embodiment will be described. A medical image processing apparatus 300 illustrated in FIG. 12 includes a volume
ボリュームデータ生成部３０１は、実施の形態１と同様、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部３０３と、補正処理部３０６と、画像生成部３０７とへ出力する。
As in the first embodiment, the volume
また、ボリュームデータ生成部３０１は、コンピュータ断層撮影装置から、装置の撮像条件（例えば、Ｘ線量、撮影時間、造影剤の注入の時刻）を取得し、画像生成部３０７へ出力する。撮像条件により、のちに生成される実画像の画質が大きく変化する。なお、コンピュータ断層撮影装置の撮像条件以外にも、コンピュータ断層撮影装置の撮影対象の条件（患者の動きや、患者の体内に含まれる金属など）により、のちに生成される実画像の画質が大きく変化する。
Further, the volume
なお、ボリュームデータ生成部３０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータ、コンピュータ断層撮影装置の撮像条件、コンピュータ断層撮影装置の撮影対象の条件を保持することができる。
The volume
ここで、本実施の形態では、ボリュームデータ生成部３０１は、実施の形態１と同様、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
Here, in the present embodiment, the volume
以下、説明のため、ボリュームデータ生成部３０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。また、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
Hereinafter, for the sake of explanation, the volume
動き解析情報生成部３０３は、ボリュームデータ生成部３０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。
Of the volume data output from the volume
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。 Here, the motion analysis information refers to information on a correspondence relationship between corresponding positions or corresponding objects on a plurality of volume data, and information on a process in which the positions and objects move and change. In addition, each volume data pixel serves as an index indicating a position at an arbitrary time between time k-1 and time k .
そして、動き解析情報生成部３０３は、生成した動き解析情報を、補間ボリュームデータ生成部３０５と画像生成部３０７とに出力する。動き解析情報生成部３０３が生成する動き解析情報は、実施の形態１の動き解析情報生成部１０３と同様（図２参照）、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを示す。
Then, the motion analysis
補間ボリュームデータ生成部３０５は、操作部３１３の入力指示に基づき、必要に応じて、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、補正処理部３０６に出力する。
Interpolated volume
ここで、本実施の形態では、実施の形態１と同様、補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、実施の形態１で上述した理由と同じであるため、詳細な説明は省略する。 In the present embodiment, as in the first embodiment, the time information t aux with interpolation volume data t aux is not in the middle of the time information and the time information t k-1 and time information t k. In other words, the distance A between the time information t k-1 of the volume data t k-1 which is generated in the time information t aux and immediately before the timing of the interpolation volume data, the time information t aux and immediately after the timing of the interpolated volume data different spacing and distance B between the time information t k of the generated volume data to, i.e., in the non-equal interval. Since this reason is the same as the reason described in the first embodiment, a detailed description is omitted.
本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像（入力画像）の画質、又は２つの実画像間を補間する補間画像の画質を補正する。
As one feature of the medical image processing apparatus 300 according to the present embodiment, the
＜補正処理１：画質の均一化＞
補正処理の一例として、コンピュータ断層撮影装置の撮影対象の条件などから、補正処理部３０６の補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一でない場合、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、ボリュームデータ生成部３０１から入力されたボリュームデータｔｋのボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋ、の画素値を求める。詳細は、図１３、図１４を参照して後述する。
<Correction process 1: Uniform image quality>
As an example of the correction process, and the like imaging target condition of the computer tomography apparatus, when the image quality of the image at time t k-1, t k before correction process of the
図１３に、補正処理前後の実画像ｔｋ−１、ｔｋの例（１）を示す。なお、補正処理前の実画像（入力画像）は、説明のため示しているが、画像生成部３０７が実際に生成する画像（出力画像）は、補正処理部３０６で補正処理が行われた後の実画像ｔｋ−１、ｔｋである。図１３に示すように、補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一ではないが、補正処理部３０６で補正処理が行われた後は、実画像ｔｋ−１、ｔｋの画質は均一となる。そのため、本実施の形態に係る医療用画像処理装置３００は、動き解析情報及び画像の画質情報に基づき、複数のボリュームデータからそれぞれ画像を作成するときの画質を向上させることができる。さらに、本実施の形態に係る医療用画像処理装置３００は、補正処理後の各画像の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の経時間変化を、正確に表示することができる。
なお、必要に応じて、補正処理部３０６は、補間ボリュームデータ生成部３０５から出力された補間ボリュームデータｔａｕｘを、実画像ｔｋ−１、ｔｋと同様に補正処理を行う。
FIG. 13 shows an example (1) of the real images t k−1 and t k before and after the correction process. Note that the actual image (input image) before the correction processing is shown for explanation, but the image (output image) actually generated by the
If necessary, the
次に、図１４を参照して、実画像ｔｋ−１、ｔｋの画質を均一とするための補正処理について説明する。図１４は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（１）である。 Next, with reference to FIG. 14, a correction process for making the image quality of the actual images t k−1 and t k uniform will be described. FIG. 14 is an example (1) of the correction process flow of the medical image processing apparatus 300 according to the present embodiment.
図１４に示すように、ステップＳＴ３１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ３２へ遷移する。 As shown in FIG. 14, in step ST31, volume data t k and t k−1 are acquired for two or more input phases corresponding to before and after the output phase. And it changes to step ST32.
ステップＳＴ３２では、ボリュームデータｔｋ、ｔｋ−１それぞれについて、ボクセル位置ごとに、信頼度（０〜１）を決定する。ここで、ボリュームデータｔｋ、ｔｋ−１についての信頼度（信頼性の高さ）は、出力フェイズとの時間的距離、ボリュームデータｔｋ、ｔｋ−１を撮影した時の撮影条件（Ｘ線照射量など）、当該ボリュームデータｔｋ、ｔｋ−１のＳ／Ｎ比、近傍の領域の移動速度や画素値の特徴などにより、決定することが考えられる。そして、ステップＳＴ３３へ遷移する。 In step ST32, the reliability (0 to 1) is determined for each voxel position for each of the volume data t k and t k−1 . Here, the reliability of the volume data t k, t k-1 (high reliability) is the temporal distance between the output phase, volume data t k, photographing conditions at the time of photographing a t k-1 ( X-ray irradiation amount, etc.), the S / N ratio of the volume data t k , t k−1 , the moving speed of the neighboring region, the feature of the pixel value, etc. And it changes to step ST33.
ステップＳＴ３３では、出力フェイズの全ボクセルについて、ボリュームデータｔｋ、ｔｋ−１の画素値を、動き解析情報及びボリュームデータｔｋ、ｔｋ−１より求める。
そして、ステップＳＴ３４へ遷移する。
ステップＳＴ３４では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報により求める。そして、ステップＳＴ３５へ遷移する。
In step ST33, the pixel values of the volume data t k and t k−1 are obtained from the motion analysis information and the volume data t k and t k−1 for all voxels in the output phase.
And it changes to step ST34.
In step ST34, for all voxels in the output phase, corresponding positions in the preceding and following input phases are obtained from the motion analysis information. And it changes to step ST35.
ステップＳＴ３５では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ３６へ遷移する。
ステップＳＴ３６では、入力フェイズに対応する入力画素値の信頼度を取得する。
そして、ステップＳＴ３７へ遷移する。
In step ST35, the pixel value at each corresponding position is sampled on the input image, and the sampled value is set as the input pixel value corresponding to the input phase. And it changes to step ST36.
In step ST36, the reliability of the input pixel value corresponding to the input phase is acquired.
And it changes to step ST37.
ステップＳＴ３７では、入力フェイズに対応する入力画素値に、取得した信頼度を乗じ、当該入力画素値について、各入力フェイズで加算したものを、出力フェイズで、補正後のボリュームデータｔｋ、ｔｋ−１の画素値とする。そして、処理が終了する。 In step ST37, the input pixel value corresponding to the input phase is multiplied by the obtained reliability, and the input pixel value is added in each input phase, and the corrected volume data t k and t k in the output phase. The pixel value is -1 . Then, the process ends.
ここで、ステップＳＴ３２における、ボリュームデータｔｋ、ｔｋ−１の信頼度（信頼性の高さ）について、画素値の信頼性の高さに応じ、４Ｄフィルタのパラメータを各実画像ｔｋ、ｔｋ−１毎、又はボクセル毎に変えることで、全体の適正なＳ／Ｎを保ちつつ、元画像の空間分解能を最大限活用することができる。 Here, with regard to the reliability (high reliability) of the volume data t k and t k−1 in step ST32, the parameters of the 4D filter are set to the respective real images t k , according to the reliability of the pixel values. By changing every t k−1 or every voxel, the spatial resolution of the original image can be utilized to the maximum while maintaining the entire appropriate S / N.
画素値の信頼性は、時刻、ボクセルの位置、画素値、収縮度、速度、加速度、近傍の画素の変位から求めた回転量、変形の度合いなどの関数である。また、近傍の領域の特徴（高画素値の領域の有無、など）が、画素値の信頼性に影響することもある。 The reliability of the pixel value is a function of time, the position of the voxel, the pixel value, the contraction degree, the speed, the acceleration, the rotation amount obtained from the displacement of the neighboring pixel, the degree of deformation, and the like. In addition, characteristics of neighboring areas (such as the presence or absence of a high pixel value area) may affect the reliability of pixel values.
４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。 Examples of 4D filter methods include smoothing, edge enhancement, maximum / minimum value projection, difference, cumulative addition, histogram matching, and combinations thereof.
＜補正処理２：画像の特性の均一＞
また、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、並びにボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求める。詳細は、図１５を参照して後述する。
<Correction process 2: Uniformity of image characteristics>
The
なお、補正処理部３０６が画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求めるのは、たとえば、コンピュータ断層撮影装置の撮影対象の条件(画質情報の元となる条件)などから、生成される画像の特性（明度、ダイナミックレンジ）が均一でない場合が考えられる。
Note that the
ここで、画素値の値的な歪み量とは、対応時刻、変位後の位置、対応位置の入力画素値、変位場の収縮度、変位の大きさ（速度）やその加速度、近傍での最大画素値や、その最大画素値のボクセルから当該ボクセルまでの距離、などの関数である。 Here, the amount of distortion of the pixel value is the corresponding time, the position after displacement, the input pixel value at the corresponding position, the contraction degree of the displacement field, the magnitude (speed) of the displacement, its acceleration, and the maximum in the vicinity. These are functions such as the pixel value and the distance from the voxel having the maximum pixel value to the voxel.
ボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルと、動き解析情報から画素値の歪の量の関数を求める方法としては、対応する位置の画素値の変化量を画像内で統計処理する、などが考えられる。 Volume data t k input from the volume data generation unit 301, and the voxel of the voxel and the volume data t k-1 of the method for determining the function of the amount of distortion of the pixel values from the motion analysis information, the corresponding position pixel For example, the amount of change in the value is statistically processed in the image.
そして、補正処理部３０６で求めた画素値の値的な歪み量に基づき、実画像ｔｋ−１、ｔｋを補正処理すると、生成される画像の特性（明度、ダイナミックレンジ）が均一な画像を作成することができる。 Then, the correction processing unit on the basis of the value distortion amount of the pixel value obtained at 306, when the correction processing of the actual image t k-1, t k, characteristics of the image to be generated (lightness, dynamic range) of uniform image Can be created.
＜補正処理フロー（２）：画像の特性（明度、ダイナミックレンジ）の均一化＞
ここで、図１５を参照して、実画像ｔｋ−１、ｔｋの画質、における特性（明度、ダイナミックレンジ）を均一とするための、補正処理部３０６で補正処理について説明する。図１５は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（２）である。
<Correction processing flow (2): Uniformity of image characteristics (brightness, dynamic range)>
Here, with reference to FIG. 15, correction processing by the
図１５に示すように、ステップＳＴ４１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ４２へ遷移する。 As shown in FIG. 15, in step ST41, volume data t k and t k−1 are acquired for two or more input phases corresponding to before and after the output phase. And it changes to step ST42.
ステップＳＴ４２では、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１と動き解析情報から、画素値の値的な歪の量の関数を算出する。そして、ステップＳＴ４３へ遷移する。
ステップＳＴ４３では、歪の量の関数から、ボリュームデータｔｋ、ｔｋ−１の各ボクセル位置について、画素値の補正係数を算出する。そして、ステップＳＴ４４へ遷移する。
In step ST42, for two or more input phases, the volume data t k, from t k-1 and the motion analysis information, calculates the function of the amount of valued distortion of pixel values. And it changes to step ST43.
In step ST43, a correction coefficient for the pixel value is calculated for each voxel position of the volume data t k and t k−1 from the function of the amount of distortion. And it changes to step ST44.
ステップＳＴ４４では、出力フェイズの全ボクセルについて、動き解析情報及びボリュームデータｔｋ、ｔｋ−１から、ボリュームデータｔｋ、ｔｋ−１の画素値を算出する。そして、ステップＳＴ４５へ遷移する。
ステップＳＴ４５では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報から算出する。そして、ステップＳＴ４６へ遷移する。
In step ST44, for all voxels in the output phase is calculated motion analysis information and volume data t k, from t k-1, volume data t k, the pixel value of t k-1. And it changes to step ST45.
In step ST45, for all voxels in the output phase, corresponding positions in the previous and subsequent input phases are calculated from the motion analysis information. And it changes to step ST46.
ステップＳＴ４６では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ４７へ遷移する。
ステップＳＴ４７では、上記入力画素値に、対応ボクセル位置の補正係数による変換をし、これをそれぞれの入力フェイズについて加算したものを、出力フェイズの画素値とする。そして、処理が終了する。
In step ST46, the pixel value at each corresponding position is sampled on the input image, and the sampled value is set as the input pixel value corresponding to the input phase. And it changes to step ST47.
In step ST47, the pixel value of the output phase is obtained by converting the input pixel value by the correction coefficient of the corresponding voxel position and adding it for each input phase. Then, the process ends.
＜補正処理３：変形又は変異による画素値の変化＞
また、本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報に基づき、動き解析情報から得られる入力画像内の特定の領域に関する変形又変位による画素値の変化を見込んだ補正処理を行うこともできる。図１６を参照して、説明する。
<Correction process 3: Change in pixel value due to deformation or mutation>
In addition, as one feature of the medical image processing apparatus 300 according to the present embodiment, the
図１６は、実画像ｔｋ−１、ｔｋの補正処理の例（２）を説明するための図である。説明のため、図１６では、入力画像ｔｋ−１、ｔｋ内の特定の領域を領域Ｒ１とする。図１６に示すように、領域Ｒ１は、実画像ｔｋ−１から実画像ｔｋへ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ１では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像ｔｋ上、領域Ｒ１の画素値が大きくなってしまう。これは、造影剤の注入の時刻情報を加味すると更に効果的である。 Figure 16 is a diagram for explaining an example of correction processing of the actual image t k-1, t k ( 2). For the sake of explanation, in FIG. 16, a specific region in the input images t k−1 and t k is defined as a region R1. As shown in FIG. 16, region R1 is reduced during the transition from the actual image t k-1 to the actual image t k. If region R is reduced, such as a contrast agent in the region R1 is concentrated If no correction process is executed, the actual image t k, the pixel value of the region R1 increases. This is more effective when the time information of contrast medium injection is taken into account.
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる動き解析情報に基づき、領域Ｒ１の変形変異を取得し、実画像ｔｋに対して変形変異を見込んだ補正処理を行う。
Therefore, the
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ１に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
As described above, the
＜補正処理４：代表値での置き換え＞
図１７は、実画像ｔｋ−１、ｔｋ、ｔｋ＋１の補正処理の例（３）を説明するための図である。説明のため、図１７では、実画像（入力画像）ｔｋ−１、ｔｋ、ｔｋ＋１内の特定の領域を領域Ｒ２とする。図１７に示すように、領域Ｒ２は、実画像ｔｋ−１から実画像ｔｋ＋１へ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ２では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像（入力画像）ｔｋ−１、よりも、実画像ｔｋ、ｔｋ＋１の領域Ｒ２の画素値が大きくなってしまう。
<Correction process 4: Replacement with representative value>
FIG. 17 is a diagram for explaining an example (3) of the correction process for the real images t k−1 , t k , and t k + 1 . For the sake of explanation, in FIG. 17, specific regions in the real images (input images) t k−1 , t k , and t k + 1 are defined as a region R2. As illustrated in FIG. 17, the region R2 is reduced during the transition from the real image tk -1 to the real image tk + 1 . When the region R is reduced, the contrast agent and the like are concentrated in the region R2. If no correction process is performed, the region R2 of the real images t k and t k + 1 rather than the real image (input image) t k−1 . The pixel value of becomes large.
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋについて実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値を比較し、もっともＳ／Ｎの高い画素の画素値を、実画像ｔｋにおける画素値とする。
Therefore, the
また、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋに対して、実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値で、Ｓ／Ｎを加味した加重平均を取り、それを新しい画素値とすることも可能である。
The
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ２に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
As described above, the
＜補正処理５：２以上実画像の画質の均一化＞
図１８は、実画像及び補間画像の補正処理の例（４）を説明するための図である。補正処理の他の例として、図１８に示すように、上述した補正処理を組み合わせて、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋ間の補間画像tａｕｘ２について、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び画質を均一化する補正処理（図１３、図１４参照）後の実画像ｔｋ−１、ｔｋの画質情報などに基づき、補間画像tａｕｘ２の補正処理を行うことができる。
<Correction Process 5: Uniformity of Image Quality of Two or More Real Images>
FIG. 18 is a diagram for explaining an example (4) of the correction process of the real image and the interpolated image. Other examples of correction processing, as shown in FIG. 18, a combination of above-described correction processing, the
補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１６を参照して説明した方法を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋに関する動き解析情報に基づき、図１８に示す領域Ｒ３の変形変異を取得し、補間画像tａｕｘ２に対して変形変異を見込んだ補正処理を行う。
When performing the correction processing of the interpolated image t aux2 , for example, when the method described with reference to FIG. 16 is used, the
また、補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１７を参照して説明した方法を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、補間画像ｔａｕｘ２間、及び補間像ｔａｕｘ２、実画像ｔｋ間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、補間実画像ｔａｕｘ２について実画像ｔｋ−１、ｔｋ、及び補間画像ｔａｕｘ２で対応する画素値を比較し、もっともＳ／Ｎの高い画素の画素値を、実画像ｔａｕｘ２における画素値とする。この場合、比較する補間画像ｔａｕｘ２の画素値をゼロとすれば良い。
Further, when performing the correction processing of the interpolated image t aux2 , for example, when the method described with reference to FIG. 17 is used, the
＜補正処理６：２番目以降の画像の画素値の差分の保存＞
図１９は、実画像の補正処理の例（５）を説明するための図である。補正処理の他の例として、図１９に示すように、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋについて、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び実画像ｔｋ−１、ｔｋの画質情報などに基づき、実画像ｔｋの補正処理を行い、実画像ｔｋ−１から数えて２番目の実画像ｔｋについては、実画像ｔｋ−１の画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。同様に、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ、ｔｋ＋１について、実画像ｔｋ、ｔｋ＋１に関する動き解析情報、及び実画像ｔｋ、ｔｋ＋１の画質情報などに基づき、実画像ｔｋ＋１の補正処理を行い、実画像ｔｋ＋１については、実画像ｔｋの画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。
<Correction Process 6: Saving Differences in Pixel Values of Second and Later Images>
FIG. 19 is a diagram for explaining an example (5) of a real image correction process. As another example of the correction process, as illustrated in FIG. 19, the
画像生成部３０７は、補正処理部３０６から出力される補正処理後のボリュームデータｔｋ−１、ｔｋ、とその時刻情報ｔｋ−１、ｔｋ、補正処理後の補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。そして、画像生成部３０７は、補正処理後の、補間ボリュームデータｔａｕｘ、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成し、時系列に並べる。
The
表示処理部３０９は、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部３０７で生成された３次元画像のシーケンスを生成して、表示部３１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
The
表示部３１１は、表示処理部３０９で適宜処理された、時系列の実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部３０７で生成された３次元画像のシーケンスを、表示処理部３０９での処理の後、動画像として再生することができる。
The
以上、実施の形態２に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置３００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。 As described above, according to the medical image processing apparatus 300 according to the second embodiment, based on the image quality information of two or more real images in addition to the motion analysis information, the image quality of two or more real images or between two real images is obtained. By correcting the image quality of the interpolated image to be interpolated, it is possible to accurately display a change in a specific region (for example, an examination site (region of interest) of the subject (patient 3)) in two or more actual images. Therefore, according to the medical image processing apparatus 300 according to the present embodiment, it is possible to accurately observe the examination site (region of interest) of the patient 3 as the subject.
例えば、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置３００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。 For example, the medical image processing apparatus 300 corrects the image quality of two or more real images or the image quality of an interpolated image that interpolates between two real images based on image quality information of two or more real images in addition to motion analysis information. If this is the case, it is possible to satisfactorily observe a cardiac function failure or the like simply by reproducing the cine using the conventional rendering method. In addition, if the medical image processing apparatus 300 interpolates interpolated images between real images based on motion analysis information, it is possible to observe not only those with short time intervals such as cardiac multiphase data but also changes over time of lesions. It can also be applied to long time interval data.
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、複数のボリュームデータに対してセグメンテーション処理を行うときに処理結果を安定させることが出来る。また、臓器間癒着の判断を自動処理で行う場合も医師等が動画を直接観察する場合も容易になる。 Further, the medical image processing apparatus 300 corrects the image quality of two or more real images or the image quality of an interpolated image that interpolates between two real images based on the image quality information of two or more real images in addition to the motion analysis information. In this way, anyone can obtain accurate and accurate results at a lower cost than in the past by fully automated quantification with reproducibility in strain analysis, perfusion analysis, cardiac function analysis, and the like. Furthermore, the medical image processing apparatus 300 uses the image quality information of two or more actual images in addition to the motion analysis information, and the image quality of the two or more actual images or the image quality of the interpolated image that interpolates between the two actual images. Is corrected, the processing result can be stabilized when segmentation processing is performed on a plurality of volume data. In addition, it is easy to determine the adhesion between organs by automatic processing and when a doctor or the like directly observes a moving image.
また、本実施の形態に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、表示部３１１のノイズ軽減、ノイズのばらつき調整が可能である。
Further, according to the medical image processing apparatus 300 according to the present embodiment, based on the image quality information of two or more real images in addition to the motion analysis information, the image quality of two or more real images or between two real images is determined. If the image quality of the interpolated image to be interpolated is corrected, the noise of the
例えば、医療用画像処理装置３００が生成した動き解析情報によって、実画像又は補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置３００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。 For example, the noise can be reduced without losing the original detailed image information by specifying the corresponding pixel of the actual image or the interpolated image based on the motion analysis information generated by the medical image processing apparatus 300. Moreover, even if a phase with a large amount of X-ray irradiation and a phase with a small amount are mixed in the CT multi-phase examination by the medical image processing apparatus 300, the pixel values of the preceding and succeeding phases are used for the phase with a small amount of X-ray irradiation. By doing so, a good interpolation image can be obtained.
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。 Further, the medical image processing apparatus 300 corrects the image quality of two or more real images or the image quality of an interpolated image that interpolates between two real images based on the image quality information of two or more real images in addition to the motion analysis information. In this case, the X-ray exposure dose of the patient 3 can be reduced by generating an actual image by increasing the X-ray irradiation amount only in an important phase, and generating an interpolated image in other phases. .
なお、本実施の形態に係る医療用画像処理装置３００では、コンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部３０１に供給されているが、これに限らない。ボリュームデータ生成部３０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
In the medical image processing apparatus 300 according to the present embodiment, a large number of continuous tomographic signals in the diagnosis range of the patient 3 are supplied to the volume
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。 In this embodiment, the time t k of the volume data t k is the time after a lapse from the volume data t k-1 a predetermined time, the predetermined time is arbitrary time, for example, the subject It can be set as appropriate depending on the time when a change in the examination site (region of interest) of (patient 3) is predicted to occur.
なお、本実施の形態において、表示処理部３０９が、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部３１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
In the present embodiment, the
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻に基づき、生成しているが、これに限らない。 In the present embodiment, the time information at which the volume data is generated is generated based on the time supplied by the computed tomography apparatus, but the present invention is not limited to this.
なお、本実施の形態の説明に用いた各機能ブロック（図１２参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。 Each functional block (see FIG. 12) used in the description of the present embodiment is typically realized as an LSI that is an integrated circuit. In particular, calculation processing for generating an actual image and an interpolation image can be performed by a GPU (Graphic Processing Unit). The GPU is an arithmetic processing unit designed specifically for image processing as compared with a general-purpose CPU (Central Processing Unit), and is mounted on a computer separately from a normal CPU.
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。 In the present embodiment, the example used for observing the heart is shown, but the present invention is not limited to this. For example, it can be used for observation of changes in tumor tissue over time, observation of blood flow perfusion, and the like.
１００、３００ 医療用画像処理装置
１０１、３０１ ボリュームデータ生成部
１０３、３０３ 動き解析情報生成部
１０５、３０５ 補間ボリュームデータ生成部
１０７、３０７ 画像生成部
１０９、３０９ 表示処理部
１１１、３１１ 表示部
１１３、３１３ 操作部
３０６ 補正処理部
100, 300 Medical
Claims (6)
前記複数のボリュームデータを生成し、それぞれの画質情報を生成する第１生成部と、
前記複数のボリュームデータのそれぞれの動き解析情報を生成する第２生成部と、
前記複数のボリュームデータの内の少なくともひとつのボリュームデータについて、当該ボリュームデータのボクセルと、当該ボリュームデータとは異なる他のボリュームデータのボクセルとを動き解析情報により対応付け、前記対応付けしたボクセルの信頼度を前記画質情報に基づいて決定し、前記対応付けしたボクセルのボクセル値を、前記対応付けしたボクセルの信頼度に応じて変わるパラメータのフィルタを用いて補正する補正処理部と、
前記補正処理部により補正したボクセル値に基づき、前記複数のボリュームデータを可視化した複数の画像を表示するための処理を行う表示処理部と、を備える医療用画像処理装置。 A medical image processing apparatus for visualizing each of a plurality of volume data having time information,
A first generation unit that generates the plurality of volume data and generates image quality information of each of the volume data;
A second generator for generating motion analysis information for each of the plurality of volume data;
For at least one volume data of the plurality of volume data, a voxel of the volume data and a voxel of other volume data different from the volume data are associated by motion analysis information, and the reliability of the associated voxel A correction processing unit that determines a degree based on the image quality information, and corrects the voxel value of the associated voxel using a parameter filter that changes according to the reliability of the associated voxel;
A medical image processing apparatus comprising: a display processing unit that performs a process for displaying a plurality of images obtained by visualizing the plurality of volume data based on the voxel values corrected by the correction processing unit.
前記ボリュームデータのボクセルと、当該ボリュームデータとは異なる他のボリュームデータのボクセルとの間を補間する補間ボリュームデータを生成する第３生成部を備える、医療用画像処理装置。 The medical image processing apparatus according to claim 1 , further comprising:
A medical image processing apparatus comprising: a third generation unit that generates interpolated volume data for interpolating between voxels of the volume data and voxels of other volume data different from the volume data.
前記補正処理部は、更に、
前記動き解析情報から、前記複数の画像のうちの少なくとも一つの画像において、部分的な変異又は変形を取得し、
前記変異又は変形から、前記補間ボリュームデータのボクセル値の変化を予測し、
当該予測したボクセル値の変化に基づき、前記フィルタを補正する医療用画像処理装置。 The medical image processing apparatus according to claim 2 ,
The correction processing unit further includes:
From the motion analysis information, in at least one of the plurality of images, obtaining a partial variation or deformation,
From the mutation or deformation, predict the change in the voxel value of the interpolated volume data ,
A medical image processing apparatus that corrects the filter based on the predicted change in the voxel value.
前記第３生成部は、
前記ボリュームデータと前記補間ボリュームデータの間隔と、前記補間ボリュームデータと前記ボリュームデータとは異なる前記他のボリュームデータの間隔とは、非等間隔となるように、前記補間ボリュームデータを生成する、医療用画像処理装置。 The medical image processing apparatus according to claim 2 or 3,
The third generator is
Generating the interpolated volume data such that an interval between the volume data and the interpolated volume data and an interval between the other volume data different from the interpolated volume data and the volume data are unequal. Image processing device.
前記第１生成部は、前記画質情報を医療画像撮影装置から取得した情報を少なくとも用いて生成する、医療画像処理装置。 The medical image processing apparatus according to claim 1,
The first generation unit is a medical image processing device that generates the image quality information using at least information acquired from a medical image capturing device.
コンピュータを、請求項１ないし５のいずれか一項の医療用画像処理装置として機能させる医療用画像処理プログラム。 A medical image processing program,
A medical image processing program for causing a computer to function as the medical image processing apparatus according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013252829A JP5722984B2 (en) | 2013-12-06 | 2013-12-06 | Medical image processing apparatus and medical image processing program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013252829A JP5722984B2 (en) | 2013-12-06 | 2013-12-06 | Medical image processing apparatus and medical image processing program |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009250859A Division JP2011092547A (en) | 2009-10-30 | 2009-10-30 | Medical image processor and medical image processing program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2014050753A JP2014050753A (en) | 2014-03-20 |
| JP5722984B2 true JP5722984B2 (en) | 2015-05-27 |
Family
ID=50609770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013252829A Active JP5722984B2 (en) | 2013-12-06 | 2013-12-06 | Medical image processing apparatus and medical image processing program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5722984B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7317639B2 (en) * | 2019-09-05 | 2023-07-31 | 富士フイルムヘルスケア株式会社 | Radiation image processing system and image processing method |
| WO2021141135A1 (en) * | 2020-01-09 | 2021-07-15 | 学校法人東京女子医科大学 | Functional ischemia detection technique by coronary artery ct4d flow imaging |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5273040A (en) * | 1991-11-14 | 1993-12-28 | Picker International, Inc. | Measurement of vetricle volumes with cardiac MRI |
| JP3350208B2 (en) * | 1994-03-23 | 2002-11-25 | 東芝医用システムエンジニアリング株式会社 | Image display device |
| US6192079B1 (en) * | 1998-05-07 | 2001-02-20 | Intel Corporation | Method and apparatus for increasing video frame rate |
| JP2003244659A (en) * | 2002-02-14 | 2003-08-29 | Toshiba Corp | Frame interpolation method |
| US6697508B2 (en) * | 2002-05-10 | 2004-02-24 | Visiongate, Inc. | Tomographic reconstruction of small objects using a priori knowledge |
| US7813473B2 (en) * | 2002-07-23 | 2010-10-12 | General Electric Company | Method and apparatus for generating temporally interpolated projections |
| US7596267B2 (en) * | 2003-02-28 | 2009-09-29 | Cedara Software Corp. | Image region segmentation system and method |
| JP4180465B2 (en) * | 2003-07-31 | 2008-11-12 | サミー株式会社 | Launch assist mechanism for ball game machines |
| JP4703193B2 (en) * | 2005-01-14 | 2011-06-15 | 株式会社東芝 | Image processing device |
| JP4327139B2 (en) * | 2005-09-05 | 2009-09-09 | ザイオソフト株式会社 | Image processing method and image processing program |
| EP1958155A2 (en) * | 2005-11-24 | 2008-08-20 | Philips Intellectual Property & Standards GmbH | Motion compensated ct reconstruction of high contrast objects |
| JP2008035895A (en) * | 2006-08-01 | 2008-02-21 | Ziosoft Inc | Image processing method and image processing program |
| JP4854546B2 (en) * | 2007-03-06 | 2012-01-18 | キヤノン株式会社 | Image processing apparatus and image processing method |
| JP2008220414A (en) * | 2007-03-08 | 2008-09-25 | Canon Inc | Apparatus and methdo for processing radiation image |
| JP4690361B2 (en) * | 2007-05-28 | 2011-06-01 | 富士フイルム株式会社 | Cardiac function analysis apparatus, method and program thereof |
| JP5514397B2 (en) * | 2007-10-03 | 2014-06-04 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Image display apparatus and X-ray tomography apparatus |
| WO2009093693A1 (en) * | 2008-01-25 | 2009-07-30 | Konica Minolta Holdings, Inc. | Image generation device, image generation method, and program |
-
2013
- 2013-12-06 JP JP2013252829A patent/JP5722984B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2014050753A (en) | 2014-03-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9251560B2 (en) | Medical image diagnosis apparatus and phase determination method using medical image diagnosis apparatus | |
| JP5643218B2 (en) | X-ray CT apparatus and image display method using X-ray CT apparatus | |
| WO2016080311A1 (en) | X-ray ct device, projection data upsampling method, and image reconstruction method | |
| KR20150085462A (en) | Medical imaging apparatus and medical image processing method thereof | |
| US9519993B2 (en) | Medical image processing apparatus | |
| US8355555B2 (en) | System and method for multi-image based virtual non-contrast image enhancement for dual source CT | |
| WO2005046478A1 (en) | Image processing method, image processing device, medical image diagnosis support system, and time-axis direction filtering method | |
| JP7292942B2 (en) | MEDICAL IMAGE PROCESSING APPARATUS, METHOD AND PROGRAM | |
| JP5156863B2 (en) | X-ray CT system | |
| JP2012000135A (en) | Multi-modality dynamic image diagnostic apparatus | |
| JP7017492B2 (en) | Tomographic image generator, method and program | |
| JP2012147934A (en) | X-ray diagnostic device, image processing device and image processing program | |
| JP5408493B2 (en) | Medical image processing apparatus and medical image processing program | |
| JP2010246862A (en) | Medical image generation apparatus and program | |
| JP5002344B2 (en) | Medical image diagnostic apparatus and medical image display apparatus | |
| JPWO2009093693A1 (en) | Image generating apparatus, image generating method, and program | |
| JP5978432B1 (en) | Medical image processing apparatus and medical image processing method | |
| US9965875B2 (en) | Virtual projection image method | |
| JP5722984B2 (en) | Medical image processing apparatus and medical image processing program | |
| JP2011067279A (en) | Medical image processor and medical image processing program | |
| JP6419551B2 (en) | X-ray diagnostic apparatus, image processing apparatus, and image processing program | |
| US20140212016A1 (en) | X-ray ct apparatus and image correction method | |
| JP2011092547A (en) | Medical image processor and medical image processing program | |
| JP6955909B2 (en) | Image processing device | |
| WO2004071301A1 (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: 20131206 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140110 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140731 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140805 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20141002 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20150209 |
|
| 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: 20150317 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150326 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5722984 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |