JP2009011628A - 画像診断支援装置 - Google Patents
画像診断支援装置 Download PDFInfo
- Publication number
- JP2009011628A JP2009011628A JP2007178228A JP2007178228A JP2009011628A JP 2009011628 A JP2009011628 A JP 2009011628A JP 2007178228 A JP2007178228 A JP 2007178228A JP 2007178228 A JP2007178228 A JP 2007178228A JP 2009011628 A JP2009011628 A JP 2009011628A
- Authority
- JP
- Japan
- Prior art keywords
- image
- deformation
- ultrasonic
- displacement
- organ
- 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.)
- Pending
Links
- 238000003745 diagnosis Methods 0.000 title claims description 11
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 210000000056 organ Anatomy 0.000 claims abstract description 31
- 239000000523 sample Substances 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000002591 computed tomography Methods 0.000 claims description 8
- 230000033001 locomotion Effects 0.000 claims description 6
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 6
- 238000003825 pressing Methods 0.000 claims description 6
- 238000002059 diagnostic imaging Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 10
- 239000000126 substance Substances 0.000 abstract 1
- 238000001356 surgical procedure Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 11
- 206010028980 Neoplasm Diseases 0.000 description 2
- 206010052428 Wound Diseases 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000002357 laparoscopic surgery Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 208000004550 Postoperative Pain Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002350 laparotomy Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】人体モデル画像や臓器モデル画像を使った手術シミュレーションにおいて、臓器の変形を、極力、被検体の実体に近づける。
【解決手段】画像支援装置は、先ず被検体の臓器等の変形のしやすさ、しにくさを、実測結果から画像処理によって求める。この実測では、先ず超音波深触子の押圧の有無の2つの状態で、同一部位から超音波計測を行い2つの超音波断層像を得、この2つの超音波断層像相互の画素単位又は所定の画素数(M×M)ブロック単位の押圧による変位を求める。次に、この変位から画素単位又はブロック単位の歪又は弾性率を求める。かくして得た弾性画像を、ボリューム画像の画素対応に埋め込む。この結果を用いて変形の様子を計算により求め、変形再現を行う。
【選択図】図2
【解決手段】画像支援装置は、先ず被検体の臓器等の変形のしやすさ、しにくさを、実測結果から画像処理によって求める。この実測では、先ず超音波深触子の押圧の有無の2つの状態で、同一部位から超音波計測を行い2つの超音波断層像を得、この2つの超音波断層像相互の画素単位又は所定の画素数(M×M)ブロック単位の押圧による変位を求める。次に、この変位から画素単位又はブロック単位の歪又は弾性率を求める。かくして得た弾性画像を、ボリューム画像の画素対応に埋め込む。この結果を用いて変形の様子を計算により求め、変形再現を行う。
【選択図】図2
Description
本発明は、MRI装置やCT装置などの画像診断装置で撮像された被検体のボリューム画像データを表示し、ボリューム画像データ内にある被検体の臓器の形状を変形させることが可能な画像診断支援装置に関する。
近年、計算機性能の向上によりコンピュータグラフィックスを用いたシミュレーションシステムが開発されている。医療分野においても、MRI装置やCT装置などで撮像した被検体のボリューム画像データを用いて人体及び臓器モデルを作成し、作成したモデルの変形及び切除などの手術シミュレーションを行うシステムが研究されている。応用例として腹腔鏡手術支援シミュレーションがある。腹腔鏡手術は腹部に5〜10mmの穴を開け、その穴から腹腔鏡や術具などを挿入して腹腔鏡を見ながら腹部内部の患部の切開や摘出、縫合を行うものである。開腹せずに小さな傷で済むため、術後の痛みが少ないことや傷口が目立たないこと、入院期間が短期間であることなどの利点があり、今後ますます盛んになる手術法であると考えられる。このような背景の中で、医用画像を用いた手術シミュレーションシステムは手術計画や手術練習用などで有用になると期待されている。
手術シミュレーションシステムを医療現場に適用するにあたって、医用画像から作成した人体または臓器モデルの変形をいかに実物に近づけるかが課題となる。
特許文献1では、臓器モデルをポリゴンデータで表現し、ポリゴンデータを形成する格子点間の接続線にバネモデルで定義した弾性情報を付与することで臓器の変形を実現している。
特開平10−111958
物体の運動はほぼ等価的にバネモデルで表すことが可能であるため、バネモデルは臓器変形の際のモデルとして適当であると言える。しかし組織毎に弾性率を示すバネ定数が存在するため、異なる組織で形成される臓器の場合には個々にバネ定数を設定する必要がある。人体組織のバネ定数を実際に測定するのは困難であり、実情は文献値や経験的に決定したバネ定数を用いて臓器変形を行っている。このため臓器内に例えば硬さの異なる腫瘍等が存在する場合、腫瘍をも周囲の組織と同様のバネ定数を用いることになってしまい、精度のよいシミュレーションを行うことができない。
本発明はこのような事情を鑑みてなされたもので、超音波診断装置を用いて組織の弾性率を求め、求めた弾性率をボリューム画像データに対応させることで、より精度の高い臓器変形シミュレーションを実現することを目的とする。
前記目的を達成するために本発明に係る画像診断支援装置は、
磁気共鳴撮像(MRI)装置またはX線コンピュータ断層(CT)装置などの画像診断装置で撮像された被検体のボリューム画像データを読み込む画像読み込み手段と、読み込んだ画像をモニタなどの表示装置に表示する画像表示手段とを有する画像診断支援装置において、
超音波診断装置を用いて得られた時系列で且つ外部からの押圧の有無による2つの超音波断層画像間で演算を行って超音波断層画像上の押圧の有無による各点の移動量または変位を計測する変位計測手段と、前記変位計測手段で求めた変位から超音波断層画像上の各点の組織の変形のしやすさ、しにくさを表わすデータを演算する弾性演算手段と、前記変位計測手段を求めるのに使用した超音波断層画像のスキャン断面の座標を求めるスキャン断面算出手段と、前記画像読込手段によって読み込んだボリューム画像データから、前記スキャン断面算出手段によって求めた超音波断層画像のスキャン断面の座標に対応する画素の座標を算出し、その画素に前記弾性演算手段によって求めた組織の変形のしやすさ、しにくさを表わすデータを対応させる弾性画像作成手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータによって色を変え、ボリューム画像データ上に付加してモニタなどの表示装置に表示する弾性画像表示手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータを用いて画像内に映る臓器などの形状を変形させる臓器変形手段と、
を備えたことを特徴とする。
磁気共鳴撮像(MRI)装置またはX線コンピュータ断層(CT)装置などの画像診断装置で撮像された被検体のボリューム画像データを読み込む画像読み込み手段と、読み込んだ画像をモニタなどの表示装置に表示する画像表示手段とを有する画像診断支援装置において、
超音波診断装置を用いて得られた時系列で且つ外部からの押圧の有無による2つの超音波断層画像間で演算を行って超音波断層画像上の押圧の有無による各点の移動量または変位を計測する変位計測手段と、前記変位計測手段で求めた変位から超音波断層画像上の各点の組織の変形のしやすさ、しにくさを表わすデータを演算する弾性演算手段と、前記変位計測手段を求めるのに使用した超音波断層画像のスキャン断面の座標を求めるスキャン断面算出手段と、前記画像読込手段によって読み込んだボリューム画像データから、前記スキャン断面算出手段によって求めた超音波断層画像のスキャン断面の座標に対応する画素の座標を算出し、その画素に前記弾性演算手段によって求めた組織の変形のしやすさ、しにくさを表わすデータを対応させる弾性画像作成手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータによって色を変え、ボリューム画像データ上に付加してモニタなどの表示装置に表示する弾性画像表示手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータを用いて画像内に映る臓器などの形状を変形させる臓器変形手段と、
を備えたことを特徴とする。
また前記スキャン断面算出手段は、超音波診断装置の超音波探触子の位置及び傾きを検出する位置センサを用い、位置センサの出力に基づいて超音波断層画像のスキャン面座標を算出することを特徴とする。
本発明によれば、超音波診断装置を用いて被検体内の組織の変形のしやすさ、しにくさを表わすデータを求め、求めたデータをボリューム画像データに対応させることで、より精度の高い臓器変形シミュレーションを実現することができる。
以下、添付図面に従って本発明に係る画像診断支援装置の好ましい実施形態について説明する。
図1は本発明の一実施形態による画像診断支援装置の概要構成図である。本画像診断支援装置は、機械的又は電子的にビーム操作を行って被検体に超音波を送信及び受信する超音波探触子10と、探触子に付随して設けられた探触子の空間的な位置や傾きを検出する位置センサ10Aと、CT装置やMRI装置で撮影されたボリューム画像データの読み込みや探触子で受信した受波から超音波断層画像を作成する画像処理装置11及び画像表示装置12、外部入力装置13、マウス14、画像表示装置の画面に表示されたソフトウェアスイッチを操作するためのコントローラ24とより成る。
図1は本発明の一実施形態による画像診断支援装置の概要構成図である。本画像診断支援装置は、機械的又は電子的にビーム操作を行って被検体に超音波を送信及び受信する超音波探触子10と、探触子に付随して設けられた探触子の空間的な位置や傾きを検出する位置センサ10Aと、CT装置やMRI装置で撮影されたボリューム画像データの読み込みや探触子で受信した受波から超音波断層画像を作成する画像処理装置11及び画像表示装置12、外部入力装置13、マウス14、画像表示装置の画面に表示されたソフトウェアスイッチを操作するためのコントローラ24とより成る。
画像処理装置11は主として各構成要素の動作を制御する中央処理装置(CPU)20と、画像処理装置の制御プログラムが格納された主メモリ21と、画像データを格納するデータ記録装置22と、被検体の画像データを一時的に記憶する表示メモリ23と画像表示装置の画面に表示されたソフトスイッチを操作するためのコントローラ24より成る。
画像表示装置12は、表示メモリからの画像データに基づいて画像を表示する。ポインティングデバイスとしてのマウス14は、画像表示装置上のソフトスイッチを操作する。外部入力装置13はキーボード等であり、各種の入力やパラメータ設定等を行う。データ記録装置22は、磁気ディスク等の記憶装置や、取り出し可能な外部メディアに対してデータの書き込みや読み出しを行う装置でもよい。
画像表示装置12は、表示メモリからの画像データに基づいて画像を表示する。ポインティングデバイスとしてのマウス14は、画像表示装置上のソフトスイッチを操作する。外部入力装置13はキーボード等であり、各種の入力やパラメータ設定等を行う。データ記録装置22は、磁気ディスク等の記憶装置や、取り出し可能な外部メディアに対してデータの書き込みや読み出しを行う装置でもよい。
図2は以上のように構成された画像診断支援装置によって、超音波診断装置を用いて組織の変形のしやすさ・しにくさを表わすデータを求め、求めたデータをボリューム画像データに対応させることで、より精度の高い臓器変形シミュレーションを実現する処理の一例を示すフローチャートである。以下、図2に示すフローチャートの各記号(1)〜(7)に従って本発明の処理の詳細を説明する。
(1)画像読込
CT装置やMRI装置などの画像診断装置によって撮影された被検体のボリューム画像データを画像記録装置又は画像データベースからオンラインまたはオフラインで読み出し、主メモリに展開する。ここでボリューム画像データは被検体を撮影した位置的に連続した数枚から数百枚の断層画像のことを指す。
また読み込んだボリューム画像データを表示メモリに展開し、モニタなどの画像表示装置に表示する。表示画像はボリューム画像データ内の1断面分のデータでも良いし、ボリューム画像データからボリュームレンダリング法又はサーフェスレンダリング法などの擬似3次元画像処理を施した画像でも良い。
CT装置やMRI装置などの画像診断装置によって撮影された被検体のボリューム画像データを画像記録装置又は画像データベースからオンラインまたはオフラインで読み出し、主メモリに展開する。ここでボリューム画像データは被検体を撮影した位置的に連続した数枚から数百枚の断層画像のことを指す。
また読み込んだボリューム画像データを表示メモリに展開し、モニタなどの画像表示装置に表示する。表示画像はボリューム画像データ内の1断面分のデータでも良いし、ボリューム画像データからボリュームレンダリング法又はサーフェスレンダリング法などの擬似3次元画像処理を施した画像でも良い。
(2)超音波スキャン
被検体へ超音波探触子10を当て、2回にわたり超音波スキャンを行う。これらのスキャンは、いずれか一方が、探触子を、通常の計測時と同じく被検体に単に当接させた状態で、他方が探触子を被検体に外力を加えた(押圧した)状態で行われる。外力を加える理由は、それによって臓器に変形を与えるためである。各超音波スキャンに際し、被検体に接触した探触子10から超音波を送信し、さらに診断部位からの反射エコー信号を探触子10で受信する。この超音波送受信を繰り返して得られた各反射エコー信号を整相、A/D変換し、時系列的に連続する複数の断層画像データとして画像処理装置の主メモリに記憶する。
尚、外力を人間によって加えるとしたが、人間の力はおおまかであることから、押圧機具によって、押圧力を表示器に表示させると共に押圧を定量化させて実行するやり方もある。
被検体へ超音波探触子10を当て、2回にわたり超音波スキャンを行う。これらのスキャンは、いずれか一方が、探触子を、通常の計測時と同じく被検体に単に当接させた状態で、他方が探触子を被検体に外力を加えた(押圧した)状態で行われる。外力を加える理由は、それによって臓器に変形を与えるためである。各超音波スキャンに際し、被検体に接触した探触子10から超音波を送信し、さらに診断部位からの反射エコー信号を探触子10で受信する。この超音波送受信を繰り返して得られた各反射エコー信号を整相、A/D変換し、時系列的に連続する複数の断層画像データとして画像処理装置の主メモリに記憶する。
尚、外力を人間によって加えるとしたが、人間の力はおおまかであることから、押圧機具によって、押圧力を表示器に表示させると共に押圧を定量化させて実行するやり方もある。
(3)スキャン断面算出
位置センサ10Aを用いてスキャン中の探触子10の断面の座標を求める。位置センサ10Aは図3に示すように探触子に取り付けられ、探触子10の位置及び傾き、言い換えれば超音波スキャン断面の位置及び傾きを検出できるようになっている。図3において、探触子10の平面状の面に振動子が配列されており、その振動子配列の中心点と位置センサ10Aの中心点との距離は正確に求められている。なお、探触子と位置センサの関係は同図に限らず、図4のように構成することができる。すなわち、位置センサに関連付けて棒状のポインタを取り外し可能に設け、その先端の点を中心点に対応した基準点とすることもできる。
位置センサ10Aを用いてスキャン中の探触子10の断面の座標を求める。位置センサ10Aは図3に示すように探触子に取り付けられ、探触子10の位置及び傾き、言い換えれば超音波スキャン断面の位置及び傾きを検出できるようになっている。図3において、探触子10の平面状の面に振動子が配列されており、その振動子配列の中心点と位置センサ10Aの中心点との距離は正確に求められている。なお、探触子と位置センサの関係は同図に限らず、図4のように構成することができる。すなわち、位置センサに関連付けて棒状のポインタを取り外し可能に設け、その先端の点を中心点に対応した基準点とすることもできる。
このような位置センサを用いるには、先ず位置センサの初期化を実施する。位置センサの位置が空間座標上のどの位置にあるかを数値化するためには、位置センサの基準点(原点)を設定する必要がある。位置センサの基準点設定では図5に示すように被検体の体表位置に探触子を当てて、その位置を基準点とする。ここで設定する基準点は、後に説明する弾性画像作成処理にてボリューム画像データとの整合をとるため、ボリューム画像データに映っており且つ、剣状突起などの外観から分かりやすい体表上の特徴点を選択することが望ましい。
基準点の設定をした後、超音波スキャンをするために移動させた探触子の空間的な位置座標及び傾きを求める。例えば図6に示すように探触子の位置を基準点から現在点に移動させてスキャンした場合、X方向に10mm、Y方向に50mm、Z方向に20mm移動ならば現在点の座標は(X,Y,Z)=(10,50,20)となる。さらに基準点を設定した際のスキャン断面がZ軸に平行でX,Y軸に垂直な面である場合、図7に示すようになるが、これに対しスキャン面がX,Y,Z軸に任意の傾きを待つ場合はスキャン面の傾きは図8に示すように(α,β,γ)で表す。
(4)変位計測
超音波スキャンによって得た時系列的に連続した2枚の超音波断層画像から画素毎に画像間の位置的変位を求める。この位置的変位は例えばブロック・マッチング法によって行う。以下図9を参照して説明する。
超音波スキャンによって得た時系列的に連続した2枚の超音波断層画像から画素毎に画像間の位置的変位を求める。この位置的変位は例えばブロック・マッチング法によって行う。以下図9を参照して説明する。
図9において、現画像をN枚目、現画像より以前の画像(前画像)をN−1枚目とする。そして各画像をM×M画素からなる複数のブロックに分割する。ここで画像上の臓器は前画像の臓器Aの位置から現画像の臓器Bの位置に移動又は変形したとする。このような状態で、N枚目画像のM×Mからなる各ブロックの左上の画素座標を(Pu,Pv)とし、その座標の画素値をXn(Pu,Pv)とする。そしてN−1枚目画像でブロックを(i,j)画素ずらした時のブロックの左上の画素座標を(Pu+i,Pv+j)とし、その座標の画素値をXn−1(Pu+i,Pv+j)とする。
次に対象ブロックについてN−1枚目とN枚目の画素値の差分をとり、その絶対値和Sを次式で算出する。
次に対象ブロックについてN−1枚目とN枚目の画素値の差分をとり、その絶対値和Sを次式で算出する。
上記の演算を分割した各ブロックに対して行い、ブロック毎に変位ΔLを求める。分割するブロックが大きいと求められる変位は大まかなものになるが、分割するブロック数を多くすればより細かい領域で変位が求められる。しかし、求まる変位はブロック単位である。ブロック単位が細かくなる程、画素単位の変位に近づき、その近似値となる。理想的には、画素単位の変位が求まることが好ましいが、ブロック単位を2×2とか3×3とかにすることで、画素単位の近似変位値を求めることができる。またブロック単位の変位を、その構成画素対応に補間させ、画素対応の変位を求めることができる。
ブロック・マッチングは上記で説明した以外にも、探触子で受信する受波波形に対しても適用できる。
ブロック・マッチングは上記で説明した以外にも、探触子で受信する受波波形に対しても適用できる。
(5)歪または弾性率の算出
変位計測で求めた超音波断層画像の各画素単位(又はブロック単位。以下同じ)の変位から歪を求める。歪は変位計測で求めた各画素の変位を空間的に微分することで求められる。画像における空間微分は隣接する画素との差分と等価であることから、隣接する画素との変位差分で求める。
変位計測で求めた超音波断層画像の各画素単位(又はブロック単位。以下同じ)の変位から歪を求める。歪は変位計測で求めた各画素の変位を空間的に微分することで求められる。画像における空間微分は隣接する画素との差分と等価であることから、隣接する画素との変位差分で求める。
例えば、図10(a)に示すように周辺組織とは硬さの異なる組織A(ここで組織Aは周辺組織よりも硬いとする。)が存在し、図10(a)のPの列に対して、探触子によって圧迫する前と後の2枚の時系列的な超音波画像から変位を求めると図10(b)のように、圧迫側の変位は大きく、深くなるに従って小さい変位となる。それ故に、Pの列に対して図10(b)の如く、算出変位は折れ線状になる。この場合、変位は周辺組織より硬い組織のほうが小さい。さらに図10(b)を深さ方向に変位に対して空間微分すると、図10(c)のような歪が求まる。
列Pは、一例であって、実際は画像の左端から右端までのすべての列に対して同様な処理を行って、全画素について歪を求める。また、列だけではなく、行方向に対しても同様な処理を行って、行と列とのそれぞれの歪を画素単位に求める例もある。この場合、行と列とを、変化させる方向に応じて選択的に使い分ける例もあれば、両者の平均値をとって1つの弾性率にする例もある。
以上、画像全体の歪を求める例を説明したが、特定の関心領域のみに絞って、歪を求める例もある。
なお、上記実施形態では、組織の変形のしやすさ、しにくさを表わす指標として歪を用いたが、弾性率を用いることもできる。組織の弾性率を求めるためには、探触子10の被検体との接触面に圧力センサを設け、2つの画像取得時の圧力差で前記歪値を割算するようにすれば良い。
以上、画像全体の歪を求める例を説明したが、特定の関心領域のみに絞って、歪を求める例もある。
なお、上記実施形態では、組織の変形のしやすさ、しにくさを表わす指標として歪を用いたが、弾性率を用いることもできる。組織の弾性率を求めるためには、探触子10の被検体との接触面に圧力センサを設け、2つの画像取得時の圧力差で前記歪値を割算するようにすれば良い。
(6)歪画像又は弾性率画像の作成
読み込んだボリューム画像データと、スキャン断面算出処理にて求めた超音波探触子のスキャン断面座標、歪又は弾性率算出処理で求めた超音波断層画像の各画素の歪又は弾性率を用いて、超音波断層画像の各画素がボリューム画像データのどの画素に対応するかを求め、歪又は弾性率とボリューム画像データを対応づける。
読み込んだボリューム画像データと、スキャン断面算出処理にて求めた超音波探触子のスキャン断面座標、歪又は弾性率算出処理で求めた超音波断層画像の各画素の歪又は弾性率を用いて、超音波断層画像の各画素がボリューム画像データのどの画素に対応するかを求め、歪又は弾性率とボリューム画像データを対応づける。
超音波断層画像の各画素が空間座標上でどの位置になるかは、スキャン断面算出処理で得た探触子の位置と傾きから求められるので、ボリューム画像データと超音波断層画像の位置を対応づけるには、スキャン断面がボリューム画像データのどの画素に対応しているかを求めればよい。
スキャン断面の空間座標に対応するボリューム画像データの画素を求めるには、先ずボリューム画像データ上でスキャン断面算出処理にて設定した初期化位置を画像上で指定する。初期化位置の指定では、例えば図11に示すようにボリューム画像データのアキシャル面、サジタル面、コロナル面を表示し、これら3つの断面画像からスキャン断面算出処理で設定した基準点を指定する。
ボリューム画像データ上で基準点を設定したら、基準点を原点としてスキャン断面算出処理と同様の座標系をとり(図6参照)、スキャン断面座標に対応するボリューム画像データの画素座標を求める。
求めた画素座標に歪又は弾性率算出処理で求めた歪又は弾性率を対応させる。対応の方法としては、例えばボリューム画像データと同じサイズで画素値が0で初期化されたマスクデータを用意し、そのマスクデータの画素に歪又は弾性率を代入する。
このようにして歪又は弾性率を求めた超音波断層画像上の位置に対応するボリューム画像データの画素に弾性率を対応させていく。対応していない位置があれば、補間でこの対応していないボリューム画素位置に、周囲の歪又は弾性率から補間した歪又は弾性率を埋め込む。
このようにして歪又は弾性率を求めた超音波断層画像上の位置に対応するボリューム画像データの画素に弾性率を対応させていく。対応していない位置があれば、補間でこの対応していないボリューム画素位置に、周囲の歪又は弾性率から補間した歪又は弾性率を埋め込む。
(7)臓器変形
弾性画像を用いてボリューム画像データ内にある臓器モデルの変形を行う。
臓器変形では、図12に示すようにボリューム画像データ内の臓器モデルが存在する各画素を制御点とし、さらに隣接する制御点とはばねモデルで定義される接続線によって関連付けられているものとする。
臓器モデルの変形は、この制御点の移動によって実現する。制御点の空間上の位置が元の位置から変化することで臓器モデルの変形が実現される。さらに制御点は隣接する接続線とばねモデルで定義した接続線によって接続されているので、ある制御点を移動させると、接続されている制御点もばねモデルのばね定数に応じた移動を行う。このばね定数を求めた歪又は弾性率にすることで実際の臓器の特徴をあわせることが可能となる。こうした変形の技術は、例えば前述した特許文献1に詳しい。
なお図12に示す制御点のモデルは説明を容易にするため2次元としたが、3次元モデルでも適用可能である。
弾性画像を用いてボリューム画像データ内にある臓器モデルの変形を行う。
臓器変形では、図12に示すようにボリューム画像データ内の臓器モデルが存在する各画素を制御点とし、さらに隣接する制御点とはばねモデルで定義される接続線によって関連付けられているものとする。
臓器モデルの変形は、この制御点の移動によって実現する。制御点の空間上の位置が元の位置から変化することで臓器モデルの変形が実現される。さらに制御点は隣接する接続線とばねモデルで定義した接続線によって接続されているので、ある制御点を移動させると、接続されている制御点もばねモデルのばね定数に応じた移動を行う。このばね定数を求めた歪又は弾性率にすることで実際の臓器の特徴をあわせることが可能となる。こうした変形の技術は、例えば前述した特許文献1に詳しい。
なお図12に示す制御点のモデルは説明を容易にするため2次元としたが、3次元モデルでも適用可能である。
10 深触子
10A 位置センサ
11 画像処理装置
12 画像表示装置
10A 位置センサ
11 画像処理装置
12 画像表示装置
Claims (2)
- 磁気共鳴撮像(MRI)装置またはX線コンピュータ断層(CT)装置などの画像診断装置で撮像された被検体のボリューム画像データを読み込む画像読み込み手段と、読み込んだ画像をモニタなどの表示装置に表示する画像表示手段とを有する画像診断支援装置において、
超音波診断装置を用いて得られた時系列で且つ外部からの押圧の有無による2つの超音波断層画像間で演算を行って超音波断層画像上の押圧の有無による各点の移動量または変位を計測する変位計測手段と、前記変位計測手段で求めた変位から超音波断層画像上の各点の組織の変形のしやすさ、しにくさを表わすデータを演算する弾性演算手段と、前記変位計測手段を求めるのに使用した超音波断層画像のスキャン断面の座標を求めるスキャン断面算出手段と、前記画像読込手段によって読み込んだボリューム画像データから、前記スキャン断面算出手段によって求めた超音波断層画像のスキャン断面の座標に対応する画素の座標を算出し、その画素に前記弾性演算手段によって求めた組織の変形のしやすさ、しにくさを表わすデータを対応させる弾性画像作成手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータによって色を変え、ボリューム画像データ上に付加してモニタなどの表示装置に表示する弾性画像表示手段と、前記弾性画像作成手段によって作成した弾性画像の組織の変形のしやすさ、しにくさを表わすデータを用いて画像内に映る臓器などの形状を変形させる臓器変形手段と、
を備えたことを特徴とする画像診断支援装置 - 前記スキャン断面算出手段は、超音波診断装置の超音波探触子の位置及び傾きを検出する位置センサを用い、位置センサの出力に基づいて超音波断層画像のスキャン面座標を算出することを特徴とした請求項1に記載の画像診断支援装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007178228A JP2009011628A (ja) | 2007-07-06 | 2007-07-06 | 画像診断支援装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007178228A JP2009011628A (ja) | 2007-07-06 | 2007-07-06 | 画像診断支援装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009011628A true JP2009011628A (ja) | 2009-01-22 |
Family
ID=40353254
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007178228A Pending JP2009011628A (ja) | 2007-07-06 | 2007-07-06 | 画像診断支援装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2009011628A (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010194146A (ja) * | 2009-02-26 | 2010-09-09 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置及びその制御プログラム |
KR20130079987A (ko) * | 2012-01-03 | 2013-07-11 | 삼성전자주식회사 | 장기 변형모델을 추정하는 방법, 장치 및 의료영상시스템 |
JP2015128546A (ja) * | 2014-01-09 | 2015-07-16 | 富士通株式会社 | 臓器位置推定装置、臓器位置推定装置の制御方法および臓器位置推定装置の制御プログラム |
US9558587B2 (en) | 2012-07-10 | 2017-01-31 | Samsung Electronics Co., Ltd. | Method and apparatus for creating model of patient specified target organ |
-
2007
- 2007-07-06 JP JP2007178228A patent/JP2009011628A/ja active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010194146A (ja) * | 2009-02-26 | 2010-09-09 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置及びその制御プログラム |
KR20130079987A (ko) * | 2012-01-03 | 2013-07-11 | 삼성전자주식회사 | 장기 변형모델을 추정하는 방법, 장치 및 의료영상시스템 |
KR101993381B1 (ko) | 2012-01-03 | 2019-06-26 | 삼성전자주식회사 | 장기 변형모델을 추정하는 방법, 장치 및 의료영상시스템 |
US9558587B2 (en) | 2012-07-10 | 2017-01-31 | Samsung Electronics Co., Ltd. | Method and apparatus for creating model of patient specified target organ |
JP2015128546A (ja) * | 2014-01-09 | 2015-07-16 | 富士通株式会社 | 臓器位置推定装置、臓器位置推定装置の制御方法および臓器位置推定装置の制御プログラム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105407811B (zh) | 用于超声图像的3d获取的方法和系统 | |
JP6367425B2 (ja) | 超音波診断装置 | |
JP5203605B2 (ja) | 超音波診断装置 | |
JP5400343B2 (ja) | 超音波による分娩の診察のための方法及び装置 | |
US20110313291A1 (en) | Medical image processing device, medical image processing method, medical image diagnostic apparatus, operation method of medical image diagnostic apparatus, and medical image display method | |
KR101100464B1 (ko) | 부 관심영역에 기초하여 3차원 초음파 영상을 제공하는 초음파 시스템 및 방법 | |
EP1795129B1 (en) | Apparatus and method for displaying an ultrasound image | |
JP6835587B2 (ja) | 医用4dイメージングにおける動き適応型可視化 | |
JP2009011628A (ja) | 画像診断支援装置 | |
JP4515799B2 (ja) | 超音波診断装置 | |
JP5075830B2 (ja) | 超音波診断装置 | |
JP4247265B2 (ja) | 超音波診断装置 | |
JP4608458B2 (ja) | 超音波診断装置 | |
JP7033692B6 (ja) | スマート剪断波エラストグラフィのための超音波システムと方法 | |
JP7309850B2 (ja) | 異方性組織の誘導せん断波エラストグラフィのための超音波システム及び方法 | |
JP2010274114A (ja) | 歪み画像データをスケーリングするためのシステムおよび方法 | |
CN115243621A (zh) | 三维超声成像数据的背景多平面重建以及相关联的设备、系统和方法 | |
JP4908915B2 (ja) | 超音波診断装置 | |
JP2010246909A (ja) | 圧力情報を提供する超音波システムおよび方法 | |
CN113301855A (zh) | 用于各向异性组织的剪切波表征的超声系统和方法 | |
JP2007275395A (ja) | 超音波診断装置 | |
WO2019082574A1 (ja) | 超音波診断装置及びその動作方法 | |
JP2017080119A (ja) | 超音波診断装置 | |
CN111970974A (zh) | 提供三维超声图像 |