JP3725442B2 - 医用画像診断装置及びその方法 - Google Patents
医用画像診断装置及びその方法 Download PDFInfo
- Publication number
- JP3725442B2 JP3725442B2 JP2001120183A JP2001120183A JP3725442B2 JP 3725442 B2 JP3725442 B2 JP 3725442B2 JP 2001120183 A JP2001120183 A JP 2001120183A JP 2001120183 A JP2001120183 A JP 2001120183A JP 3725442 B2 JP3725442 B2 JP 3725442B2
- Authority
- JP
- Japan
- Prior art keywords
- organ
- contour
- dimensional
- calculating
- vertex
- 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.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title description 17
- 210000000056 organ Anatomy 0.000 claims description 64
- 238000006073 displacement reaction Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000002123 temporal effect Effects 0.000 claims description 15
- 238000003745 diagnosis Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims 1
- 230000003287 optical effect Effects 0.000 description 24
- 238000000605 extraction Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、超音波診断装置、ＭＲＩ、Ｘ線ＣＴ等で撮られた画像データから、臓器の運動状態を解析する医用画像診断装置及びその方法に関する。
【０００２】
【従来の技術】
従来、１つの断面画像から臓器輪郭を抽出し、２次元的な臓器運動を解析する手段が考案されているが、少数の複数断面画像から簡便に臓器輪郭を抽出し、３次元的な臓器運動を解析する手段はない。
【０００３】
また、従来、３次元データを基に輪郭抽出する手段や、ＭＲ画像に対し磁気標識を付加し臓器運動を解析する手法が考案されているが、磁気標識等を用いず画像データのみから臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析する手段はない。
【０００４】
【発明が解決しようとする課題】
そこで、本発明は、臓器の３次元的輪郭を生成し、臓器運動の解析を可能にし、また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することを可能にすることを目的とする。
【０００５】
さらに、少数の複数断面から簡便に３次元的な臓器輪郭を生成し、臓器運動の解析を可能にすることを目的とする。
【０００７】
【課題を解決するための手段】
請求項１の発明は、生体器官である臓器を撮影した複数時相の画像から前記臓器に関する複数時相の３次元輪郭を算出する輪郭算出手段と、算出した前記臓器に関する複数時相の３次元輪郭に基づいて、前記臓器の３次元輪郭上の各頂点の時間的変位量を算出する変位量算出手段と、算出した前記各頂点の時間的変位量から前記臓器の輪郭の運動情報を算出する運動算出手段と、を有し、前記変位量算出手段は、前記複数時相における３次元輪郭の各頂点の法線方向の移動距離を算出し、前記運動算出手段は、前記移動距離から前記各頂点の移動距離、または、移動速度を運動情報として算出することを特徴とする医用画像診断装置である。
【００１４】
本発明によれば、臓器の３次元的輪郭を生成し、臓器運動の解析が可能となる。
【００１５】
また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することが可能となる。
【００１６】
【発明の実施の形態】
＜第１実施例＞
以下、本発明の第１実施例について、図１から図１１を参照して説明する。
【００１７】
図１は第１実施例の医用画像診断装置１０のブロック図であり、図２はその動作及び操作の流れを示すフローチャートである。
【００１８】
第１実施例は、複数の断面画像から臓器輪郭を自動的に抽出し、３次元輪郭を補間生成し、その輪郭の運動情報の時間的変化を表示する例である。
【００１９】
図１に示すように、コンピュータによって構成される医用画像診断装置１０は、画像データを入力する画像入力部１２と、入力された画像データを記憶するメモリ１４と、画像上の臓器輪郭を抽出する輪郭抽出部１６と、抽出された複数輪郭を基に３次元輪郭を生成する３次元輪郭生成部１８と、３次元輪郭の変形から運動情報を算出する運動情報算出部２０と、算出された運動情報を表示する表示部２２とで構成される。なお、下記で説明する各機能はプログラムによって実現される。
【００２０】
以下、図２の処理の流れに沿って説明する。
【００２１】
まず、例えば超音波診断装置、ＭＲＩ、Ｘ線ＣＴ装置等で得られる複数断面画像データが画像入力部１２を介して入力され、メモリ１４に記憶される（ステップＳ１）。
【００２２】
この入力データは、ディジタル化された輝度値データでも良いし、アナログ信号による画像でも良い。アナログ信号の場合は、Ａ／Ｄ変換器によりディジタル化され入力される。
【００２３】
また、入力データは輝度値データのみに限らず、例えば超音波診断装置におけるエコーの生信号データやドップラ速度データや、Ｘ線ＣＴ装置での生のＣＴ値データ等であっても良い。
【００２４】
さらに、複数断面画像は３次元データから生成されたものであっても良い。
【００２５】
次に、輪郭抽出部１６において半自動的に臓器輪郭が抽出され、輪郭情報がメモリ１４に記憶される（ステップＳ２）。
【００２６】
輪郭抽出は、例えば手動で設定された初期輪郭を基に画像輝度値のエッジ探索を行う動的輪郭モデルを用いた手法や、画像輝度値の閾値処理による領域抽出後に輪郭を検出するような手法を用いると良い。
【００２７】
また、輪郭を手動でトレースしても良い。
【００２８】
この輪郭抽出処理は、複数の断面画像に対して行われる。輪郭は、例えば図３に示すように複数の頂点により表現され、それらの複数頂点の座標値が輪郭情報としてメモリ１４に記憶される。
【００２９】
次に、３次元輪郭生成部１８が複数の輪郭情報を基に、３次元輪郭を生成し、３次元輪郭情報をメモリ１４に記憶する（ステップＳ３）。
【００３０】
３次元輪郭の生成は例えば次の処理により行われる。
【００３１】
まず、入力された複数断面の位置関係により、ステップＳ２でメモリ１４に記憶された輪郭座標を３次元座標に変換する。
【００３２】
次に、複数輪郭間の輪郭情報を補間処理し３次元輪郭全体を生成する。
【００３３】
複数断面が１つの軸の周りに回転する場合について、具体的計算例を以下に示す。
【００３４】
ステップＳ２で記憶された輪郭の頂点座標を（ｘｉ，ｙｉ）と表す。ここで、複数断面で共通の回転軸がｙ軸と等しくなるように座標軸を設定する（図９参照）。断面１を基準とし、ｎ番目の断面の３次元空間内での角度をθｎとすると、次式により、各断面内の２次元輪郭頂点座標（ｘｉ，ｙｉ）は、３次元空間内の座標（ｘ’ｉ，ｙ’ｉ，ｚ’ｉ）に変換される。
【００３５】
【数１】
次に、複数断面の間の仮想的な断面での輪郭頂点座標を補間生成し３次元輪郭を生成する。図１０は、断面ｎと断面ｎ＋１間の頂点の３次元的な補間生成を示したものである。
【００３６】
図４は、断面画像から輪郭抽出を経て３次元輪郭が生成されるまでの処理を、心室を例にして模式的に示したものである。
【００３７】
また、図５は、管状の臓器の３次元輪郭を模式的に例示したものである。３次元輪郭は、複数頂点の３次元座標で表現され、各点の３次元座標が３次元輪郭情報としてメモリ１４に記憶される。
【００３８】
次に、運動情報算出部２０がステップＳ３でメモリ１４に記憶された３次元輪郭情報から、臓器輪郭の各局所部分の運動情報を算出し、メモリ１４に記憶する（ステップＳ４）。
【００３９】
「運動情報」とは、移動距離や移動速度等である。移動距離の算出は、例えば３次元輪郭の対応する頂点番号の頂点間の時間的移動距離によって算出する。
【００４０】
移動距離の他の算出方法として、３次元輪郭面の法線方向の移動距離により算出する方法や、臓器の運動モデルに基づき算出する方法を用いても良い。
【００４１】
本実施例は、運動情報の算出方法には制限されない。また、運動情報として、移動距離や移動速度以外の数値を算出しても良い。
【００４２】
最後に、算出された運動情報を表示部２２で表示する（ステップＳ５）。
【００４３】
この運動情報の表示は、３次元輪郭上に色の違いや数値等で表示しても良いし（図６参照）、例えばブルズアイマップのように模式的に示しても良い（図７参照）。また、部位毎の運動情報の時間的変化を心電図等と共にグラフ化しても良い（図８参照）。
【００４４】
以上の処理内容を行う場合の、ユーザの操作手順について図１１に基づいて以下に説明する。
【００４５】
ユーザは、まず１つの断面動画像または複数時相の画像を入力する（ステップＳ１１）。
【００４６】
次に、輪郭抽出処理を行う（ステップＳ１２）。輪郭抽出処理を自動で行う場合は、ユーザは自動輪郭抽出のための初期位置の指定のみ行えば良い。
【００４７】
次に、各断面動画像において３次元輪郭生成に用いる動画区間または複数時相を選択する（ステップＳ１３）。ここで、時相とは、時系列で入力された動画像の各画像における時間をいう。
【００４８】
ステップＳ１３は、例えば心臓の拡張末期から収縮末期までの区間の３次元輪郭を生成したいといった場合に、各断面動画像においてはフレーム番号と心臓の拡張収縮の時相との関係がばらばらであるため、使用する動画区間を指定するものである。
【００４９】
ステップＳ１１からステップＳ１３までを断面数分だけ繰り返す。
【００５０】
最後に、ユーザは表示実行を指示すると（ステップＳ１４）、３次元輪郭生成、運動情報算出が自動的に行われ、結果が表示部２２に表示される。
【００５１】
以上のように、本実施例によれば簡便な操作で３次元的な輪郭及びその運動情報を定量的に得ることができる。
【００５２】
第１の実施例では複数断面を利用する例を示したが、本実施例を用いて１断面から輪郭を抽出し３次元輪郭を補間生成しても良い。
【００５３】
また、複数断面として例えば心臓を対象とした超音波診断画像における、心尖部四腔断面像、心尖部二腔断面像、心尖部左室長軸断面像を用いると、３断面のみから左心室の代表的部位を網羅した壁運動情報を得ることができる。
【００５４】
＜第２実施例＞
次に、第２実施例について図１２から図１７を参照して説明する。
【００５５】
第２実施例は、３次元画像データからオプティカルフローを算出し、算出された速度ベクトルをもとに３次元輪郭の時間的変位量や時間的変位速度といった運動情報を算出し、その運動情報を反映するように輪郭表示をすることで、より実際に近い運動の様子を表示する例である。
【００５６】
図１２は第２実施例の医用画像診断装置１１０のブロック図であり、図１３は動作の流れを示す。
【００５７】
図１２に示すように、第２実施例のコンピュータによって構成される医用画像診断装置１１０は、画像データを入力する画像入力部１１２と、入力された画像データを記憶するメモリ１１４と、画像上の臓器輪郭を抽出する輪郭抽出部１１６と、抽出された複数輪郭を基に３次元輪郭を生成する３次元輪郭生成部１１８と、画像データからオプティカルフローを算出するオプティカルフロー算出部１１９と、３次元輪郭の変形から運動情報を算出する運動情報算出部１２０と、算出された運動情報を表示する表示部１２２とで構成される。なお、下記で説明する各機能はプログラムによって実現される。
【００５８】
以下、図１３の処理の流れに沿って説明する。
【００５９】
まず、３次元画像データが画像入力部１１２で入力されメモリ１１４に記憶される（ステップＳ２１）。
【００６０】
次に、３次元輪郭抽出部１１６で３次元画像データから３次元輪郭が抽出され、３次元輪郭情報がメモリ１１４に記憶される（ステップＳ２２）。
【００６１】
次に、オプティカルフロー算出部１１９で画像データからオプティカルフローを算出する（ステップＳ２３）。
【００６２】
オプティカルフローは、画像上の対象物のある点が時間と共に画像上のどこへ移動したかを示すものである。オプティカルフローの算出は、勾配法やブロックマッチングによる方法等の方法を用いれば良い。
【００６３】
オプティカルフローは、連続するフレーム間で算出しても良いし、基準フレームと現フレームとの間で算出しても良い。
【００６４】
次に、運動情報算出部１２０においてオプティカルフロー情報に基づき、３次元輪郭上の各部分における時間的変位量や時間的変位速度等の運動情報を算出し、メモリ１４に記憶する（ステップＳ２４）。
【００６５】
オプティカルフローは画像上の全ての部分において算出可能ではないので、オプティカルフローのない部分は、近隣のオプティカルフロー値を補間して生成すると良い。
【００６６】
例えば、図１４のように、頂点近傍のオプティカルフローの線形和や、その補間により頂点に対応したオプティカルフローを生成することができる。
【００６７】
次に、運動情報をもとに３次元輪郭情報を修正し、表示用３次元輪郭情報を生成する（ステップＳ２５）。
【００６８】
表示用３次元輪郭情報の生成の方法を以下に説明する。
【００６９】
表示用３次元輪郭の形状は３次元輪郭抽出部１１６で抽出された３次元輪郭の形状と同じであるが、２つの時相の表示用３次元輪郭を構成する各頂点の変位量が、オプティカルフローから算出された変位量に最も近くなるように各点を配置し、表示用３次元輪郭を生成する（図１５参照）。つまり、表示用３次元輪郭を構成する各点の運動が、臓器の局所運動を近似するように、表示用３次元輪郭を構成する各頂点を配置する。
【００７０】
このように、オプティカルフロー情報に基づき輪郭の頂点座標を決定すれば、第１実施例に示したような、輪郭決定後に頂点番号のみを基準に頂点を対応付ける方法に比べ、より正確な運動状態を表すことができる。また、後述する表示方法を用いれば、臓器の実際の運動の様子をより正確に知ることができる。
【００７１】
上記実施例は、３次元輪郭形状は画像データから求めた輪郭形状を保持し、頂点座標を輪郭表面上で修正する例であるが、画像データに基づく３次元輪郭よりもオプティカルフロー情報を重視して、基準フレーム以降の３次元輪郭形状をオプティカルフロー情報から生成しても良い。あるいは、画像データに基づく3次元輪郭形状とオプティカルフローに基づく３次元輪郭形状の双方の線形和により輪郭形状を決定しても良い。
【００７２】
最後に、表示部１２２がステップＳ２５で生成された表示用３次元輪郭情報と運動情報を表示する（ステップＳ２６）。
【００７３】
３次元輪郭表示は、例えば運動情報により色付けされた輪郭上に、輪郭を構成する点間を直線で結んだメッシュ表示を重ね合わせると良い。
【００７４】
この様に、輪郭を構成する各点を表示することで、臓器局所の運動情報をより正確に得ることができ、臓器のねじれ運動等を反映した表示ができる（図１６，図１７）。
【００７５】
第２実施例では３次元画像データを入力し直接３次元輪郭を抽出する例を示したが、複数の断面画像を利用し、断面内でのオプティカルフローを算出しても良い。ただし、その場合は断面内の運動は反映されるが、断面以外の方向の運動は反映されない。
【００７６】
【発明の効果】
本発明によれば、臓器の３次元的輪郭を生成し、臓器運動の解析が可能となる。また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することが可能となる。これにより、臓器の診断が大幅に効率良く適切に行え、その効果は多大である。
【図面の簡単な説明】
【図１】第１実施例の構成図である。
【図２】第１実施例の処理の流れのフローチャートである。
【図３】第１実施例の輪郭表現の一例である。
【図４】第１実施例における３次元輪郭生成までの処理の模式図である。
【図５】第１実施例における管状の３次元輪郭の模式図である。
【図６】本発明の運動情報の表示例である。
【図７】本発明の運動情報の表示例である。
【図８】本発明の運動情報の表示例である。
【図９】複数断面が１つの軸の周りに回転する場合の断面図である。
【図１０】断面ｎと断面ｎ＋１間の頂点の３次元的な補間生成を示した図である。
【図１１】ユーザの操作の流れを示すフローチャートである。
【図１２】第２実施例の構成図である。
【図１３】第２実施例の処理の流れのフローチャートである。
【図１４】頂点近傍のオプティカルフローの線形和や、その補間により頂点に対応したオプティカルフローを生成する図である。
【図１５】第２実施例の表示用３次元輪郭の生成例である。
【図１６】オプティカルフローを反映しない場合の３次元輪郭の例である。
【図１７】オプティカルフローを反映した場合の３次元輪郭の例である。
【符号の説明】
１０ 医用画像診断装置
１２ 画像入力部
１４ メモリ
１６ 輪郭抽出部
１８ ３次元輪郭生成部
２０ 運動情報算出部
２２ 表示部
【発明の属する技術分野】
本発明は、超音波診断装置、ＭＲＩ、Ｘ線ＣＴ等で撮られた画像データから、臓器の運動状態を解析する医用画像診断装置及びその方法に関する。
【０００２】
【従来の技術】
従来、１つの断面画像から臓器輪郭を抽出し、２次元的な臓器運動を解析する手段が考案されているが、少数の複数断面画像から簡便に臓器輪郭を抽出し、３次元的な臓器運動を解析する手段はない。
【０００３】
また、従来、３次元データを基に輪郭抽出する手段や、ＭＲ画像に対し磁気標識を付加し臓器運動を解析する手法が考案されているが、磁気標識等を用いず画像データのみから臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析する手段はない。
【０００４】
【発明が解決しようとする課題】
そこで、本発明は、臓器の３次元的輪郭を生成し、臓器運動の解析を可能にし、また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することを可能にすることを目的とする。
【０００５】
さらに、少数の複数断面から簡便に３次元的な臓器輪郭を生成し、臓器運動の解析を可能にすることを目的とする。
【０００７】
【課題を解決するための手段】
請求項１の発明は、生体器官である臓器を撮影した複数時相の画像から前記臓器に関する複数時相の３次元輪郭を算出する輪郭算出手段と、算出した前記臓器に関する複数時相の３次元輪郭に基づいて、前記臓器の３次元輪郭上の各頂点の時間的変位量を算出する変位量算出手段と、算出した前記各頂点の時間的変位量から前記臓器の輪郭の運動情報を算出する運動算出手段と、を有し、前記変位量算出手段は、前記複数時相における３次元輪郭の各頂点の法線方向の移動距離を算出し、前記運動算出手段は、前記移動距離から前記各頂点の移動距離、または、移動速度を運動情報として算出することを特徴とする医用画像診断装置である。
【００１４】
本発明によれば、臓器の３次元的輪郭を生成し、臓器運動の解析が可能となる。
【００１５】
また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することが可能となる。
【００１６】
【発明の実施の形態】
＜第１実施例＞
以下、本発明の第１実施例について、図１から図１１を参照して説明する。
【００１７】
図１は第１実施例の医用画像診断装置１０のブロック図であり、図２はその動作及び操作の流れを示すフローチャートである。
【００１８】
第１実施例は、複数の断面画像から臓器輪郭を自動的に抽出し、３次元輪郭を補間生成し、その輪郭の運動情報の時間的変化を表示する例である。
【００１９】
図１に示すように、コンピュータによって構成される医用画像診断装置１０は、画像データを入力する画像入力部１２と、入力された画像データを記憶するメモリ１４と、画像上の臓器輪郭を抽出する輪郭抽出部１６と、抽出された複数輪郭を基に３次元輪郭を生成する３次元輪郭生成部１８と、３次元輪郭の変形から運動情報を算出する運動情報算出部２０と、算出された運動情報を表示する表示部２２とで構成される。なお、下記で説明する各機能はプログラムによって実現される。
【００２０】
以下、図２の処理の流れに沿って説明する。
【００２１】
まず、例えば超音波診断装置、ＭＲＩ、Ｘ線ＣＴ装置等で得られる複数断面画像データが画像入力部１２を介して入力され、メモリ１４に記憶される（ステップＳ１）。
【００２２】
この入力データは、ディジタル化された輝度値データでも良いし、アナログ信号による画像でも良い。アナログ信号の場合は、Ａ／Ｄ変換器によりディジタル化され入力される。
【００２３】
また、入力データは輝度値データのみに限らず、例えば超音波診断装置におけるエコーの生信号データやドップラ速度データや、Ｘ線ＣＴ装置での生のＣＴ値データ等であっても良い。
【００２４】
さらに、複数断面画像は３次元データから生成されたものであっても良い。
【００２５】
次に、輪郭抽出部１６において半自動的に臓器輪郭が抽出され、輪郭情報がメモリ１４に記憶される（ステップＳ２）。
【００２６】
輪郭抽出は、例えば手動で設定された初期輪郭を基に画像輝度値のエッジ探索を行う動的輪郭モデルを用いた手法や、画像輝度値の閾値処理による領域抽出後に輪郭を検出するような手法を用いると良い。
【００２７】
また、輪郭を手動でトレースしても良い。
【００２８】
この輪郭抽出処理は、複数の断面画像に対して行われる。輪郭は、例えば図３に示すように複数の頂点により表現され、それらの複数頂点の座標値が輪郭情報としてメモリ１４に記憶される。
【００２９】
次に、３次元輪郭生成部１８が複数の輪郭情報を基に、３次元輪郭を生成し、３次元輪郭情報をメモリ１４に記憶する（ステップＳ３）。
【００３０】
３次元輪郭の生成は例えば次の処理により行われる。
【００３１】
まず、入力された複数断面の位置関係により、ステップＳ２でメモリ１４に記憶された輪郭座標を３次元座標に変換する。
【００３２】
次に、複数輪郭間の輪郭情報を補間処理し３次元輪郭全体を生成する。
【００３３】
複数断面が１つの軸の周りに回転する場合について、具体的計算例を以下に示す。
【００３４】
ステップＳ２で記憶された輪郭の頂点座標を（ｘｉ，ｙｉ）と表す。ここで、複数断面で共通の回転軸がｙ軸と等しくなるように座標軸を設定する（図９参照）。断面１を基準とし、ｎ番目の断面の３次元空間内での角度をθｎとすると、次式により、各断面内の２次元輪郭頂点座標（ｘｉ，ｙｉ）は、３次元空間内の座標（ｘ’ｉ，ｙ’ｉ，ｚ’ｉ）に変換される。
【００３５】
【数１】
【００３６】
図４は、断面画像から輪郭抽出を経て３次元輪郭が生成されるまでの処理を、心室を例にして模式的に示したものである。
【００３７】
また、図５は、管状の臓器の３次元輪郭を模式的に例示したものである。３次元輪郭は、複数頂点の３次元座標で表現され、各点の３次元座標が３次元輪郭情報としてメモリ１４に記憶される。
【００３８】
次に、運動情報算出部２０がステップＳ３でメモリ１４に記憶された３次元輪郭情報から、臓器輪郭の各局所部分の運動情報を算出し、メモリ１４に記憶する（ステップＳ４）。
【００３９】
「運動情報」とは、移動距離や移動速度等である。移動距離の算出は、例えば３次元輪郭の対応する頂点番号の頂点間の時間的移動距離によって算出する。
【００４０】
移動距離の他の算出方法として、３次元輪郭面の法線方向の移動距離により算出する方法や、臓器の運動モデルに基づき算出する方法を用いても良い。
【００４１】
本実施例は、運動情報の算出方法には制限されない。また、運動情報として、移動距離や移動速度以外の数値を算出しても良い。
【００４２】
最後に、算出された運動情報を表示部２２で表示する（ステップＳ５）。
【００４３】
この運動情報の表示は、３次元輪郭上に色の違いや数値等で表示しても良いし（図６参照）、例えばブルズアイマップのように模式的に示しても良い（図７参照）。また、部位毎の運動情報の時間的変化を心電図等と共にグラフ化しても良い（図８参照）。
【００４４】
以上の処理内容を行う場合の、ユーザの操作手順について図１１に基づいて以下に説明する。
【００４５】
ユーザは、まず１つの断面動画像または複数時相の画像を入力する（ステップＳ１１）。
【００４６】
次に、輪郭抽出処理を行う（ステップＳ１２）。輪郭抽出処理を自動で行う場合は、ユーザは自動輪郭抽出のための初期位置の指定のみ行えば良い。
【００４７】
次に、各断面動画像において３次元輪郭生成に用いる動画区間または複数時相を選択する（ステップＳ１３）。ここで、時相とは、時系列で入力された動画像の各画像における時間をいう。
【００４８】
ステップＳ１３は、例えば心臓の拡張末期から収縮末期までの区間の３次元輪郭を生成したいといった場合に、各断面動画像においてはフレーム番号と心臓の拡張収縮の時相との関係がばらばらであるため、使用する動画区間を指定するものである。
【００４９】
ステップＳ１１からステップＳ１３までを断面数分だけ繰り返す。
【００５０】
最後に、ユーザは表示実行を指示すると（ステップＳ１４）、３次元輪郭生成、運動情報算出が自動的に行われ、結果が表示部２２に表示される。
【００５１】
以上のように、本実施例によれば簡便な操作で３次元的な輪郭及びその運動情報を定量的に得ることができる。
【００５２】
第１の実施例では複数断面を利用する例を示したが、本実施例を用いて１断面から輪郭を抽出し３次元輪郭を補間生成しても良い。
【００５３】
また、複数断面として例えば心臓を対象とした超音波診断画像における、心尖部四腔断面像、心尖部二腔断面像、心尖部左室長軸断面像を用いると、３断面のみから左心室の代表的部位を網羅した壁運動情報を得ることができる。
【００５４】
＜第２実施例＞
次に、第２実施例について図１２から図１７を参照して説明する。
【００５５】
第２実施例は、３次元画像データからオプティカルフローを算出し、算出された速度ベクトルをもとに３次元輪郭の時間的変位量や時間的変位速度といった運動情報を算出し、その運動情報を反映するように輪郭表示をすることで、より実際に近い運動の様子を表示する例である。
【００５６】
図１２は第２実施例の医用画像診断装置１１０のブロック図であり、図１３は動作の流れを示す。
【００５７】
図１２に示すように、第２実施例のコンピュータによって構成される医用画像診断装置１１０は、画像データを入力する画像入力部１１２と、入力された画像データを記憶するメモリ１１４と、画像上の臓器輪郭を抽出する輪郭抽出部１１６と、抽出された複数輪郭を基に３次元輪郭を生成する３次元輪郭生成部１１８と、画像データからオプティカルフローを算出するオプティカルフロー算出部１１９と、３次元輪郭の変形から運動情報を算出する運動情報算出部１２０と、算出された運動情報を表示する表示部１２２とで構成される。なお、下記で説明する各機能はプログラムによって実現される。
【００５８】
以下、図１３の処理の流れに沿って説明する。
【００５９】
まず、３次元画像データが画像入力部１１２で入力されメモリ１１４に記憶される（ステップＳ２１）。
【００６０】
次に、３次元輪郭抽出部１１６で３次元画像データから３次元輪郭が抽出され、３次元輪郭情報がメモリ１１４に記憶される（ステップＳ２２）。
【００６１】
次に、オプティカルフロー算出部１１９で画像データからオプティカルフローを算出する（ステップＳ２３）。
【００６２】
オプティカルフローは、画像上の対象物のある点が時間と共に画像上のどこへ移動したかを示すものである。オプティカルフローの算出は、勾配法やブロックマッチングによる方法等の方法を用いれば良い。
【００６３】
オプティカルフローは、連続するフレーム間で算出しても良いし、基準フレームと現フレームとの間で算出しても良い。
【００６４】
次に、運動情報算出部１２０においてオプティカルフロー情報に基づき、３次元輪郭上の各部分における時間的変位量や時間的変位速度等の運動情報を算出し、メモリ１４に記憶する（ステップＳ２４）。
【００６５】
オプティカルフローは画像上の全ての部分において算出可能ではないので、オプティカルフローのない部分は、近隣のオプティカルフロー値を補間して生成すると良い。
【００６６】
例えば、図１４のように、頂点近傍のオプティカルフローの線形和や、その補間により頂点に対応したオプティカルフローを生成することができる。
【００６７】
次に、運動情報をもとに３次元輪郭情報を修正し、表示用３次元輪郭情報を生成する（ステップＳ２５）。
【００６８】
表示用３次元輪郭情報の生成の方法を以下に説明する。
【００６９】
表示用３次元輪郭の形状は３次元輪郭抽出部１１６で抽出された３次元輪郭の形状と同じであるが、２つの時相の表示用３次元輪郭を構成する各頂点の変位量が、オプティカルフローから算出された変位量に最も近くなるように各点を配置し、表示用３次元輪郭を生成する（図１５参照）。つまり、表示用３次元輪郭を構成する各点の運動が、臓器の局所運動を近似するように、表示用３次元輪郭を構成する各頂点を配置する。
【００７０】
このように、オプティカルフロー情報に基づき輪郭の頂点座標を決定すれば、第１実施例に示したような、輪郭決定後に頂点番号のみを基準に頂点を対応付ける方法に比べ、より正確な運動状態を表すことができる。また、後述する表示方法を用いれば、臓器の実際の運動の様子をより正確に知ることができる。
【００７１】
上記実施例は、３次元輪郭形状は画像データから求めた輪郭形状を保持し、頂点座標を輪郭表面上で修正する例であるが、画像データに基づく３次元輪郭よりもオプティカルフロー情報を重視して、基準フレーム以降の３次元輪郭形状をオプティカルフロー情報から生成しても良い。あるいは、画像データに基づく3次元輪郭形状とオプティカルフローに基づく３次元輪郭形状の双方の線形和により輪郭形状を決定しても良い。
【００７２】
最後に、表示部１２２がステップＳ２５で生成された表示用３次元輪郭情報と運動情報を表示する（ステップＳ２６）。
【００７３】
３次元輪郭表示は、例えば運動情報により色付けされた輪郭上に、輪郭を構成する点間を直線で結んだメッシュ表示を重ね合わせると良い。
【００７４】
この様に、輪郭を構成する各点を表示することで、臓器局所の運動情報をより正確に得ることができ、臓器のねじれ運動等を反映した表示ができる（図１６，図１７）。
【００７５】
第２実施例では３次元画像データを入力し直接３次元輪郭を抽出する例を示したが、複数の断面画像を利用し、断面内でのオプティカルフローを算出しても良い。ただし、その場合は断面内の運動は反映されるが、断面以外の方向の運動は反映されない。
【００７６】
【発明の効果】
本発明によれば、臓器の３次元的輪郭を生成し、臓器運動の解析が可能となる。また、臓器局所部分の実際の運動を反映した臓器運動を抽出、表示、解析することが可能となる。これにより、臓器の診断が大幅に効率良く適切に行え、その効果は多大である。
【図面の簡単な説明】
【図１】第１実施例の構成図である。
【図２】第１実施例の処理の流れのフローチャートである。
【図３】第１実施例の輪郭表現の一例である。
【図４】第１実施例における３次元輪郭生成までの処理の模式図である。
【図５】第１実施例における管状の３次元輪郭の模式図である。
【図６】本発明の運動情報の表示例である。
【図７】本発明の運動情報の表示例である。
【図８】本発明の運動情報の表示例である。
【図９】複数断面が１つの軸の周りに回転する場合の断面図である。
【図１０】断面ｎと断面ｎ＋１間の頂点の３次元的な補間生成を示した図である。
【図１１】ユーザの操作の流れを示すフローチャートである。
【図１２】第２実施例の構成図である。
【図１３】第２実施例の処理の流れのフローチャートである。
【図１４】頂点近傍のオプティカルフローの線形和や、その補間により頂点に対応したオプティカルフローを生成する図である。
【図１５】第２実施例の表示用３次元輪郭の生成例である。
【図１６】オプティカルフローを反映しない場合の３次元輪郭の例である。
【図１７】オプティカルフローを反映した場合の３次元輪郭の例である。
【符号の説明】
１０ 医用画像診断装置
１２ 画像入力部
１４ メモリ
１６ 輪郭抽出部
１８ ３次元輪郭生成部
２０ 運動情報算出部
２２ 表示部
Claims (6)
- 生体器官である臓器を撮影した複数時相の画像から前記臓器に関する複数時相の３次元輪郭を算出する輪郭算出手段と、
算出した前記臓器に関する複数時相の３次元輪郭に基づいて、前記臓器の３次元輪郭上の各頂点の時間的変位量を算出する変位量算出手段と、
算出した前記各頂点の時間的変位量から前記臓器の輪郭の運動情報を算出する運動算出手段と、
を有し、
前記変位量算出手段は、前記複数時相における３次元輪郭の各頂点の法線方向の移動距離を算出し、
前記運動算出手段は、前記移動距離から前記各頂点の移動距離、または、移動速度を運動情報として算出する
ことを特徴とする医用画像診断装置。 - 前記輪郭算出手段は、
前記臓器の複数断面が撮影された複数時相の２次元画像に基づいて、前記各断面における前記臓器の輪郭を抽出する輪郭抽出手段と、
抽出された複数時相の前記各断面における前記臓器における輪郭に基づいて、前記臓器に関する複数時相の３次元輪郭を生成する輪郭生成手段と、
を有する
ことを特徴とする請求項１に記載の医用画像診断装置。 - 前記輪郭算出手段は、
前記臓器が撮影された複数時相の３次元画像に基づいて、前記臓器に関する複数時相の３次元輪郭を算出する
ことを特徴とする請求項１に記載の医用画像診断装置。 - 前記算出された前記臓器の３次元輪郭上の各点の時間的変位量に基づき、時間的に変位する頂点によって表現された前記臓器の３次元輪郭を表示する表示手段を有した
ことを特徴とする請求項１に記載の医用画像診断装置。 - 生体器官である臓器を撮影した複数時相の画像から前記臓器に関する複数時相の３次元輪郭を算出する輪郭算出ステップと、
算出した前記臓器に関する複数時相の３次元輪郭に基づいて、前記臓器の３次元輪郭上の各頂点の時間的変位量を算出する変位量算出ステップと、
算出した前記各頂点の時間的変位量から前記臓器の輪郭の運動情報を算出する運動算出ステップと、
を有し、
前記変位量算出ステップは、前記複数時相における３次元輪郭の各頂点の法線方向の移動距離を算出し、
前記運動算出ステップは、前記移動距離から前記各頂点の移動距離、または、移動速度を運動情報として算出する
ことを特徴とする医用画像診断方法。 - 生体器官である臓器を撮影した複数時相の画像から前記臓器に関する複数時相の３次元輪郭を算出する輪郭算出機能と、
算出した前記臓器に関する複数時相の３次元輪郭に基づいて、前記臓器の３次元輪郭上の各頂点の時間的変位量を算出する変位量算出機能と、
算出した前記各頂点の時間的変位量から前記臓器の輪郭の運動情報を算出する運動算出機能と、
をコンピュータによって実現し、
前記変位量算出機能は、前記複数時相における３次元輪郭の各頂点の法線方向の移動距離を算出し、
前記運動算出機能は、前記移動距離から前記各頂点の移動距離、または、移動速度を運動情報として算出する
ことを特徴とする医用画像診断方法のプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001120183A JP3725442B2 (ja) | 2001-04-18 | 2001-04-18 | 医用画像診断装置及びその方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001120183A JP3725442B2 (ja) | 2001-04-18 | 2001-04-18 | 医用画像診断装置及びその方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002306483A JP2002306483A (ja) | 2002-10-22 |
| JP3725442B2 true JP3725442B2 (ja) | 2005-12-14 |
Family
ID=18970273
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001120183A Expired - Lifetime JP3725442B2 (ja) | 2001-04-18 | 2001-04-18 | 医用画像診断装置及びその方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3725442B2 (ja) |
Families Citing this family (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3872424B2 (ja) | 2002-12-20 | 2007-01-24 | アロカ株式会社 | 超音波診断装置 |
| US7704208B2 (en) * | 2003-06-03 | 2010-04-27 | Koninklijke Philips Electronics N.V. | Synchronizing a swiveling three-dimensional ultrasound display with an oscillating object |
| US7912528B2 (en) * | 2003-06-25 | 2011-03-22 | Siemens Medical Solutions Usa, Inc. | Systems and methods for automated diagnosis and decision support for heart related diseases and conditions |
| JP4731863B2 (ja) * | 2003-10-15 | 2011-07-27 | 親良 炭 | 変位計測方法及び超音波装置 |
| JP2005161032A (ja) * | 2003-11-10 | 2005-06-23 | Toshiba Corp | 画像処理装置 |
| JP4763989B2 (ja) * | 2004-09-22 | 2011-08-31 | 株式会社東芝 | 磁気共鳴イメージング装置、磁気共鳴データ処理装置及び磁気共鳴データ処理プログラム |
| US20060253024A1 (en) * | 2005-04-26 | 2006-11-09 | Altmann Andres C | Software product for three-dimensional cardiac imaging using ultrasound contour reconstruction |
| WO2006137294A1 (ja) * | 2005-06-21 | 2006-12-28 | National University Corporation Kanazawa University | Ｘ線診断支援装置、プログラム及び記録媒体 |
| JP4996128B2 (ja) * | 2006-04-27 | 2012-08-08 | 株式会社東芝 | 医用画像処理装置及び医用画像処理方法 |
| EP2059907A2 (en) * | 2006-08-29 | 2009-05-20 | Koninklijke Philips Electronics N.V. | Reduction of heart motion artifacts in thoracic ct imaging |
| JP5624258B2 (ja) * | 2007-04-26 | 2014-11-12 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
| JP5414157B2 (ja) * | 2007-06-06 | 2014-02-12 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム |
| JP4966108B2 (ja) * | 2007-06-25 | 2012-07-04 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム |
| JP5523681B2 (ja) * | 2007-07-05 | 2014-06-18 | 株式会社東芝 | 医用画像処理装置 |
| JP5390081B2 (ja) * | 2007-07-25 | 2014-01-15 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
| JP2009112356A (ja) * | 2007-11-02 | 2009-05-28 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置 |
| US8320711B2 (en) * | 2007-12-05 | 2012-11-27 | Biosense Webster, Inc. | Anatomical modeling from a 3-D image and a surface mapping |
| WO2009072397A1 (ja) | 2007-12-07 | 2009-06-11 | Konica Minolta Medical & Graphic, Inc. | 動態撮影システム |
| JP5159299B2 (ja) * | 2007-12-27 | 2013-03-06 | 株式会社東芝 | Ｘ線ｃｔ装置及び医用画像処理装置 |
| US20090177089A1 (en) | 2008-01-04 | 2009-07-09 | Assaf Govari | Three-dimensional image reconstruction using doppler ultrasound |
| JP5558672B2 (ja) * | 2008-03-19 | 2014-07-23 | 株式会社東芝 | 画像処理装置及びｘ線コンピュータ断層撮影装置 |
| US8326007B2 (en) * | 2008-05-12 | 2012-12-04 | General Electric Company | Methods and apparatus for combined 4D presentation of quantitative regional measurements and morphology |
| JP5718548B2 (ja) * | 2008-11-13 | 2015-05-13 | 株式会社東芝 | 超音波診断装置 |
| US8224053B2 (en) * | 2009-03-31 | 2012-07-17 | General Electric Company | Methods and systems for displaying quantitative segmental data in 4D rendering |
| JP2011161220A (ja) * | 2010-01-14 | 2011-08-25 | Toshiba Corp | 画像処理装置、ｘ線コンピュータ断層撮像装置及び画像処理プログラム |
| EP2563212B1 (en) * | 2010-04-28 | 2018-04-04 | Koninklijke Philips N.V. | Visualization of myocardial infarct size in diagnostic ecg |
| DE102010038905A1 (de) * | 2010-08-04 | 2012-02-09 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Bestimmen von Rad- und Karosseriebewegungen eines Fahrzeuges |
| JP5323795B2 (ja) * | 2010-10-12 | 2013-10-23 | 富士フイルム株式会社 | 診断支援装置、診断支援プログラムおよび診断支援方法 |
| JP5651030B2 (ja) * | 2011-02-02 | 2015-01-07 | 日立アロカメディカル株式会社 | 超音波画像処理装置 |
| JP5444408B2 (ja) * | 2012-05-21 | 2014-03-19 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム |
| JP6173886B2 (ja) | 2012-11-21 | 2017-08-02 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理方法 |
| KR101727836B1 (ko) | 2015-12-22 | 2017-04-17 | 한국과학기술원 | 영상 복원 장치 및 방법 및 장치 |
| EP3291175B1 (en) * | 2016-09-05 | 2018-12-12 | RaySearch Laboratories AB | Image processing system and method for interactive contouring of three-dimensional medical data |
| WO2023054442A1 (ja) * | 2021-09-29 | 2023-04-06 | テルモ株式会社 | コンピュータプログラム、情報処理装置及び情報処理方法 |
-
2001
- 2001-04-18 JP JP2001120183A patent/JP3725442B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002306483A (ja) | 2002-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3725442B2 (ja) | 医用画像診断装置及びその方法 | |
| EP1846896B1 (en) | A method, a system and a computer program for integration of medical diagnostic information and a geometric model of a movable body | |
| JP4348383B2 (ja) | 生物学的構造体の、湾曲した解剖学的ｍモード表示を行う装置 | |
| US10660613B2 (en) | Measurement point determination in medical diagnostic imaging | |
| JP4594610B2 (ja) | 超音波画像処理装置及び超音波診断装置 | |
| JP4060615B2 (ja) | 画像処理装置及び超音波診断装置 | |
| JP6441335B2 (ja) | リアルタイム超音波イメージングのトリプレーン画像の自動セグメント化 | |
| JP5271262B2 (ja) | 三次元超音波イメージングのための方法、装置及びコンピュータプログラム | |
| JP2019048152A (ja) | ボリューム測定位相コントラストｍｒｉによる総合的心血管解析 | |
| US6482161B1 (en) | Medical diagnostic ultrasound system and method for vessel structure analysis | |
| US7803113B2 (en) | Ultrasound imaging of rotation | |
| EP3516620B1 (en) | Method and apparatus for correcting dynamic models obtained by tracking methods | |
| US20060045328A1 (en) | Image processing method for displaying information relating to parietal motions of a deformable 3-D object | |
| JP2008253753A (ja) | 心機能表示装置およびそのプログラム | |
| JP2003508139A (ja) | 非剛体運動画像解析 | |
| US9734430B2 (en) | Evaluation system or determination of cardiovascular function parameters | |
| JP2004008304A (ja) | 複数方向の投影映像を用いた３次元形状の生成方法および表示方法 | |
| JP2013192781A (ja) | 医用画像表示制御装置および方法並びにプログラム | |
| WO2006085257A2 (en) | A method, a system and a computer program for segmenting a surface in a multidimensional dataset | |
| JP4276595B2 (ja) | 超音波診断装置 | |
| US7030874B2 (en) | Method of following the three-dimensional deformation of a deformable organ | |
| CN100438831C (zh) | 多视点冠脉造影序列影像的时间配准方法 | |
| JP6135901B2 (ja) | 表示処理プログラム、表示処理方法及び表示処理装置 | |
| Loeckx et al. | Spatiotemporal non-rigid image registration for 3D ultrasound cardiac motion estimation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050322 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050519 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050621 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050816 |
|
| 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: 20050913 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050921 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3725442 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080930 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090930 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090930 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100930 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100930 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110930 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120930 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120930 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130930 Year of fee payment: 8 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |