(57)【特許請求の範囲】
【請求項１】
異なる撮像装置で撮像された被検体の被検部位の複数の画像データを整合処理する画像診断支援システムであって、
静磁場空間に被検体を収容し、磁気共鳴信号を得るパルスシーケンスでデータを収集し、収集したデータに基づいて断層画像を生成する磁気共鳴撮像装置と、
被検体に超音波を送波し超音波断層画像を生成する超音波撮像装置と、
前記超音波撮像装置による被検部位の超音波断層画像データを用いて、前記超音波断層画像における臓器の境界などの特徴点又は特徴線を抽出し、断層画像上で該特徴情報と整合させるように、前記磁気共鳴撮像装置で事前に取得された磁気共鳴断層画像における血管、腫瘍、体表及び生検針体のうちの少なくとも１つを参照しながら該磁気共鳴断層画像を変形させて、磁気共鳴断層画像を補正する画像整合装置とを有する画像診断支援システム。
【請求項２】
前記画像整合装置により補正された磁気共鳴断層画像を表示する表示装置を有する請求項１に記載の画像診断支援システム。
【請求項３】
前記画像整合装置は、前記磁気共鳴撮像装置で事前に取得された磁気共鳴断層画像データを記憶する記憶装置を有し、該記憶装置から読み出した磁気共鳴断層画像に対して前記補正処理を行う請求項１又は請求項２に記載の画像診断支援システム。
【請求項４】
前記画像整合装置は、前記磁気共鳴撮像装置で事前に取得された磁気共鳴断層画像データを記憶する記憶装置と、前記超音波撮像装置による被検部位の超音波断層画像データから超音波断層画像における臓器の境界などの特徴点又は特徴線を抽出する特徴抽出部と、事前に取得され前記記憶装置に保持された磁気共鳴断層画像を読み出し、断層画像上で前記特徴抽出部で抽出された超音波断層画像の特徴情報と整合させるように、読み出した磁気共鳴断層画像における血管、腫瘍、体表及び生検針体のうちの少なくとも１つを参照しながら該磁気共鳴断層画像を変形させて、超音波断層画像との矛盾をなくした補正磁気共鳴断層画像を生成する磁気共鳴画像補正部とを有する請求項１又は請求項２に記載の画像診断支援システム。
【請求項５】
前記磁気共鳴撮像装置は、励起パルスで被検体内のスピンを励起することによって生じる磁気共鳴信号を、エコーデータとして２次元フーリエ空間に収集し、収集したエコーデータを２次元逆フーリエ変換することにより、断層画像を再構成して前記画像整合装置に出力する請求項１から請求項４のいずれかに記載の画像診断支援システム。
【請求項６】
前記超音波撮像装置は、超音波プローブを通して得られたエコー受信信号に基づいてＢモード断層画像データをリアルタイムに生成して前記画像整合装置に出力する請求項１から請求項５のいずれかに記載の画像診断支援システム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、磁気共鳴撮像（ＭＲＩ：Magnetic Resonance Imaging）装置による画像（以下、ＭＲ画像という）、および超音波撮像装置による画像（以下、ＵＳ画像という）を合成して、医師の診断を支援する画像診断支援システムおよびその合成のための画像処理方法に関するものである。
【０００２】
【従来の技術】
ＭＲＩ装置では、静磁場空間に被検体を収容し、磁気共鳴を利用して被検体の被検部位を撮像する。
具体的には、励起パルスで被検体内のスピン（ｓｐｉｎ）を励起し、それによって生じる磁気共鳴信号を、たとえばスピンエコー（ｓｐｉｎｅｃｈｏ）またはグラディエントエコー（ｇｒａｄｉｅｎｔｅｃｈｏ）として２次元フーリエ空間に収集する。
磁気共鳴信号には、いわゆるビュー（ｖｉｅｗ）毎に異なる位相エンコードを付与し、２次元フーリエ空間において位相軸上の位置が異なる複数のビューのエコーデータをそれぞれ収集する。
そして、収集した全ビューのエコーデータを２次元逆フーリエ変換することにより、画像を再構成し、再構成画像を表示装置に表示する。
【０００３】
また、超音波撮像装置では、超音波プローブを被検体の被検部位に当接させて被検体に超音波を送波し、たとえば被検体内での非線形効果による高調波エコーに基づいて画像を生成し、表示装置に表示する。
【０００４】
これらのＭＲＩ装置および超音波撮像装置は、同じ電磁波シールド内で使用可能であるが、それぞれに装置が独立に画像を収集して表示するような形態で利用され、それぞれの画像は個別に医師等の検査時や手術中の診断の判定に用いられていた。
【０００５】
【発明が解決しようとする課題】
ところで、このような利用形態では、検査中や手術中に個別に装置を使い分ける必要があり、操作が煩雑であることから、ＭＲＩ装置の撮影断層像に応じた超音波撮影断面を容易に収集し、両画像を合成表示して診断に有用な画像を得ることが可能な超音波およびＭＲＩ複合診断装置が提案されている（たとえば、特開平９−２４０３４号公報、特開平９−２４０３５号公報  参照）。
【０００６】
この超音波およびＭＲＩ複合診断装置は、ＭＲ画像とＵＳ画像を合成表示するのであるが、ＭＲ画像とＵＳ画像には、以下に示すような長所と短所を有することから、両画像を単に合成するだけでは、両者の画像情報の矛盾等から真に有用な画像をリアルタイムに得ることは困難である。
【０００７】
すなわち、ＭＲＩ装置は、一般に軟部組織のコントラスト分解能に優れており、高空間分解能画像を得ることができるが、撮像時間が超音波診断装置に比べて長いという欠点がある。
超音波診断装置は、リアルタイムイメージングが可能であるが、ＭＲ画像に比べて画像が粗いとい欠点がある。
【０００８】
また、上記提案された超音波およびＭＲＩ複合診断装置において、たとえばＭＲＩ装置で形態的な情報を取得し、超音波撮像装置のドップラモードの画像を重ね合わせる場合を想定すると、生体を撮像対象とすることから、呼吸や突発的な運動により相互の位置ずれが起きる可能性が容易に推察できる。
【０００９】
本発明は、かかる事情に鑑みてなされたものであり、その目的は、ＭＲＩ装置と超音波撮像装置の画像情報から真に有用な合成画像をリアルタイムに得ることは可能な画像診断支援システムおよび画像処理方法を提供することにある。
【００１０】
【課題を解決するための手段】
上記目的を達成するため、本発明の第１の観点は、異なる撮像装置で撮像された被検体の被検部位の複数の画像データを整合処理する画像診断支援システムであって、静磁場空間に被検体を収容し、磁気共鳴信号を得るパルスシーケンスでデータを収集し、収集したデータに基づいて画像を生成する磁気共鳴撮像装置と、被検体に超音波を送波し画像を生成する超音波撮像装置と、上記超音波撮像装置による被検部位の超音波画像データを用いて、上記超音波画像の特徴を抽出し、当該特徴情報と整合させるように、上記磁気共鳴撮像装置で事前に取得された磁気共鳴画像を変形させて、磁気共鳴画像を補正する画像整合装置とを有する。
【００１１】
第１の観点では、上記画像整合装置により補正された磁気共鳴画像を表示する表示装置を有する。
【００１２】
また、第１の観点では、上記画像整合装置は、上記磁気共鳴撮像装置で事前に取得された磁気共鳴画像データを記憶する記憶装置を有し、上記画像整合装置は、上記記憶装置から読み出した磁気共鳴画像に対して上記補正処理を行う。
【００１３】
また、第１の観点では、上記画像整合装置は、上記磁気共鳴撮像装置で事前に取得された磁気共鳴画像データを記憶する記憶装置と、上記超音波撮像装置による被検部位の超音波画像データから超音波画像の特徴を抽出する特徴抽出部と、事前に取得され上記記憶装置に保持された磁気共鳴画像を読み出し、上記特徴抽出部で抽出された超音波画像の特徴情報と整合させるように、読み出した磁気共鳴画像を変形させて、超音波画像との矛盾をなくした補正磁気共鳴画像を生成する磁気共鳴画像補正部とを有する。
【００１４】
また、第１の観点では、上記磁気共鳴撮像装置は、励起パルスで被検体内のスピンを励起することよって生じる磁気共鳴信号を、エコーデータとして２次元フーリエ空間に収集し、収集したエコーデータを２次元逆フーリエ変換することにより、画像を再構成して上記画像整合装置に出力する。
【００１５】
また、第１の観点では、上記超音波撮像装置は、超音波プローブを通して得られたエコー受信信号に基づいてＢモード画像データをリアルタイムに生成して上記画像整合装置に出力する。
【００１６】
本発明の第２の観点は、異なる撮像装置で撮像された被検体の被検部位の複数の画像データを整合処理する画像処理方法であって、磁気共鳴撮像装置における静磁場空間に被検体を収容し、磁気共鳴信号を得るパルスシーケンスでデータを収集し、収集したデータに基づいて画像を生成することにより磁気共鳴画像データを事前に取得しておき、超音波撮像装置により、被検体に超音波を送波して超音波画像データをリアルタイムに生成し、生成された超音波画像の特徴を抽出し、上記抽出した特徴情報と整合させるように、上記事前に取得された磁気共鳴画像を変形させて、磁気共鳴画像を補正する。
【００１７】
第２の観点では、上記補正された磁気共鳴画像を表示する。
【００１８】
また、第２の観点では、上記磁気共鳴画像データは、励起パルスで被検体内のスピンを励起することよって生じる磁気共鳴信号を、エコーデータとして２次元フーリエ空間に収集し、収集したエコーデータを２次元逆フーリエ変換することにより、画像を再構成して生成する。
【００１９】
また、第２の観点では、上記超音波画像データは、超音波プローブを通して得られたエコー受信信号に基づいてＢモード画像データをリアルタイムに形成して生成する。
【００２０】
本発明によれば、まず、磁気共鳴撮像装置において、静磁場空間に被検体が収容され、励起パルスで被検体内のスピンを励起することよって生じる磁気共鳴信号が、エコーデータとして２次元フーリエ空間に収集される。そして、収集したエコーデータが２次元逆フーリエ変換され、これにより、画像が再構成されて画像整合装置に出力される。
画像整合装置では、磁気共鳴撮像装置による高分解能、高コントラスト分解能の磁気共鳴画像データがたとえば記憶装置に一端記憶される。
次に、超音波撮像装置の超音波プローブを通して得られたエコー受信信号に基づいてＢモード画像データがリアルタイムに形成され、生成された超音波画像データが画像整合装置に出力される。
画像整合装置では、超音波撮像装置によりリアルタイムに得られる超音波画像データを受けて、特徴抽出部により超音波画像から、たとえば被検部位の臓器の境界などの特徴点（線）が抽出され、磁気共鳴画像補正部に出力される。
磁気共鳴画像補正部では、超音波画像の特徴から、その情報と整合させるように記憶装置に格納されている事前に取得された磁気共鳴画像が変形されて、超音波画像と矛盾をなくした補正磁気共鳴画像が生成される。
この補正され、リアルタイムに生成された超音波画像と整合され、両画像間に矛盾がなく、しかも高分解能、高コントラスト分解能の磁気共鳴画像が表示装置に表示される。
【００２１】
【発明の実施の形態】
以下、本発明の実施形態に係る画像診断支援システムについて図面に関連付けて説明する。
【００２２】
図１は、本発明に係る画像診断支援システムの一実施形態を示す構成図である。
本画像診断支援システム１は、図１に示すように、ＭＲＩ装置２、超音波撮像装置３、画像整合装置４、および表示装置５を有している。
なお、本実施形態では、超音波撮像装置３は、被検体内での非線形効果による高調波エコーに基づいて画像を生成する装置を例に説明するが、本発明は他の超音波撮像装置にも適用することができる。
【００２３】
ＭＲＩ装置２は、静磁場空間に被検体を収容し、磁気共鳴を利用して被検体の被検部位に、励起パルスで被検体内のスピン（ｓｐｉｎ）を励起し、それによって生じる磁気共鳴信号を、たとえばスピンエコー（ｓｐｉｎｅｃｈｏ）またはグラディエントエコー（ｇｒａｄｉｅｎｔｅｃｈｏ）として２次元フーリエ空間に収集し、収集した全ビューのエコーデータを２次元逆フーリエ変換することにより、画像を再構成し、再構成画像データを信号Ｓ２として画像整合装置４に供給するとともに、表示部に表示する。
【００２４】
図２は、本実施形態に係るＭＲＩ装置２の構成例を示す図である。
【００２５】
本実施形態に係るＭＲＩ装置２は、図２に示すように、マグネットからの放射電磁波の洩漏や外乱電磁波の進入を防止する閉空間を形成した図示しないスキャンルームに配設される本体装置２Ａ、およびたとえばスキャンルームに隣接して設けられた操作ルーム内のオペレータＯＰが操作等するオペレータコンソール２Ｂを主構成要素として有している。
【００２６】
本体装置２Ａは、図２に示すように、マグネットシステム２１、ＲＦ駆動部２２、勾配駆動部２３、データ収集部２４、制御部２５、およびクレードル２６を有している。
【００２７】
マグネットシステム２１は、図２に示すように、上下の磁石に挟まれた内部空間（ギャップ：ｇａｐ）２１１を有し、ギャップ２１１内には、クッションを介して被検体６を載せたクレードル２６が図示しない搬送部によって搬入される。
【００２８】
マグネットシステム２１内には、図２に示すように、ギャップ２１１内のマグネットセンタ（走査する中心位置）の周囲に、主磁場マグネット部２１２ａ，２１２ｂ、勾配コイル部２１３ａ，２１３ｂ、およびＲＦコイル部２１４ａ，２１４ｂが配置されている。
【００２９】
主磁場マグネット部２１２ａ，２１２ｂ、勾配コイル部２１３ａ，２１３ｂ、およびＲＦコイル部２１４ａ，２１４ｂのそれぞれは、検査時に被検体６が位置するギャップ２１１内の空間を挟んで対向する１対のコイルからなる。
【００３０】
主磁場マグネット部２１２ａ，２１２ｂは、ギャップ２１１内に静磁場を形成する。静磁場の方向は、たとえば概ね被検体６の体軸方向と直交する方向である。すなわち、垂直磁場を形成する。主磁場マグネット部２１２ａ，２１２ｂを構成する一対の主磁場マグネットは、たとえば永久磁石などを用いて構成される。
【００３１】
勾配コイル部２１３ａ，２１３ｂは、ＲＦコイル部２１４ａ，２１４ｂが受信する磁気共鳴信号に３次元の位置情報を持たせるために、主磁場マグネット部２１２ａ，２１２ｂが形成した静磁場の強度に勾配を付加する勾配磁場を発生する。勾配コイル部２１３ａ，２１３ｂが発生する勾配磁場は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト（ｒｅａｄ  ｏｕｔ）勾配磁場およびフェーズエンコード（ｐｈａｓｅ  ｅｎｃｏｄｅ）勾配磁場の３種類であり、これら３種類の勾配磁場に対応して勾配コイル部２１３は３系統の勾配コイルを有する。
【００３２】
ＲＦコイル部２１４ａ，２１４ｂは、主磁場マグネット部２１２が形成した静磁場空間内で被検体６の体内にスピンを励起するための高周波磁場を形成する。ここで、高周波磁場を形成することをＲＦ励起信号の送信という。ＲＦコイル部２１４は、被検体６の体内に励起されたスピンが生じる電磁波を磁気共鳴信号として受信する。
ＲＦコイル部２１４は、図示しない送信用コイルおよび受信用コイルを有する。送信用コイルおよび受信用コイルは、同じコイルを兼用するかあるいはそれぞれ専用のコイルを用いる。
【００３３】
なお、本実施形態に場合、ＲＦコイル部２１４ａ，２１４ｂは、ＲＦ駆動部２２によるプロトコル対応の駆動信号ＤＲ１を受けて高周波磁場を形成する。
磁気共鳴撮影処理においては、１繰り返し時間（ＴＲ；ｒｅｐｅｔｉｔｉｏｎｔｉｍｅ）毎に用いるパルスシーケンス（スキャンシーケンス）の数は、被検部位毎に対応して設定されたプロトコルによって異なる。
たとえば頭部等の被検部位に応じたプロトコル毎に、それぞれ異なる回数、たとえば６４回〜５１２回繰り返されて、６４ビューから５１２ビューのビューデータが得られる。
【００３４】
ＲＦ駆動部２２は、制御部２５の指示に基づいたプロトコル対応の駆動信号ＤＲ１をＲＦコイル部２１４ａ，２１４ｂに与えてＲＦ励起信号を発生させて、被検体６の体内のスピンを励起する。
【００３５】
勾配駆動部２３は、制御部２５の指示に基づいたプロトコル対応の駆動信号ＤＲ２を勾配コイル部２１３ａ，２１３ｂに与えて勾配磁場を発生させる。
勾配駆動部２３は、勾配コイル部２１３の３系統の勾配コイルに対応して、図示しない３系統の駆動回路を有する。
【００３６】
データ収集部２４は、ＲＦコイル部２１４ａ，２１４ｂが受信した受信信号を取り込み、それをビューデータ（ｖｉｅｗ  ｄａｔａ）として収集して、オペレータコンソール２Ｂのデータ処理部２７に出力する。
【００３７】
制御部２５は、オペレータコンソール２Ｂのデータ処理部２７から送られてくる被検体６の被検部位に対応した実行すべきプロトコルに即して、あらかじめ決められた繰り返し時間ＴＲ内において所定のパルスシーケンスが所定回数繰り返される駆動信号ＤＲ１をＲＦコイル部２１４に印加するようにＲＦ駆動部２２を制御する。
同様に、制御部２５は、実行すべきプロトコルに即して、１ＴＲ内に、所定のパターンのパルス信号を勾配コイル２１３ａ，２１３ｂに印加するように勾配駆動部２３を制御する。
また、制御部２５は、ＲＦコイル部２１４ａ，２１４ｂが受信した受信信号を取り込み、それをビューデータ（ｖｉｅｗ  ｄａｔａ）として収集して、オペレータコンソール２Ｂのデータ処理部２７に出力するように、データ収集部２４を制御する。
【００３８】
この磁気共鳴撮像用パルスシーケンスは、いわゆるスピンエコー（ＳＥ：Ｓｐｉｎ  Ｅｃｈｏ）法、グラディエントエコー（ＧＲＥ：ＧＲａｄｉｅｎｔ  Ｅｃｈｏ）法、ファーストスピンエコー（ＦＳＥ：Ｆａｓｔ  Ｓｐｉｎ  Ｅｃｈｏ）法、ファーストリカバリＳＥ（Ｆａｓｔ  Ｒｅｃｏｖｅｒｙ  Ｓｐｉｎ  Ｅｃｈｏ）法、エコープラナー・イメージング（ＥＰＩ：Ｅｃｈｏ  Ｐｌａｎａｒ  Ｉｍａｇｉｎｇ）法等、各撮像方法によって異なる。
【００３９】
ここで、各撮像方法のパルスシーケンスのうち、ＳＥ法のパルスシーケンスについて、図３に関連付けて説明する。
図３（ａ）はＳＥ法におけるＲＦ励起用の９０°パルスおよび１８０°パルスのシーケンスであり、ＲＦ駆動部２２がＲＦコイル部２１４に印加する駆動信号ＤＲ１に相当する。
図３（ｂ）、（ｃ）、（ｄ）、および（ｅ）は、それぞれスライス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコード勾配Ｇｐ、およびスピンエコーＭＲのシーケンスであり、スライス勾配Ｇｓ、リードアウト勾配Ｇｒ、およびフェーズエンコード勾配Ｇｐのパルスは、勾配駆動部２３が勾配コイル部２１３に印加する駆動信号ＤＲ２に相当する。
【００４０】
図３（ａ）に示すように、ＲＦ駆動部２２によりＲＦコイル部２１４ａ，２１４ｂに対して９０°パルスが印加され、スピンの９０°励起が行われる。このとき、図３（ｂ）に示すように、勾配駆動部２３により勾配コイル部２１３ａ，２１３ｂに対してスライス勾配パルスＧｓが印加され、所定のスライスについて選択励起が行われる。
図３（ａ）に示すように、９０°励起から所定の時間後に、ＲＦ駆動部２２によりＲＦコイル部２１４ａ，２１４ｂに対して１８０°パルスが印加され、１８０°励起、すなわちスピン反転が行われる。このときも、図３（ｂ）に示すように、勾配駆動部２３により勾配コイル部２１３ａ，２１３ｂに対してスライス勾配パルスＧｓが印加され、同じスライスについて選択的な反転が行われる。
【００４１】
図３（ｃ）および（ｄ）に示すように、９０°励起とスピン反転の間の期間に、勾配駆動部２３により勾配コイル部２１３ａ，２１３ｂに対してリードアウト勾配パルスＧｒ、およびフェーズエンコード勾配パルスＧｐが印加される。
そして、リードアウト勾配パルスＧｒによりスピンのディフェーズが行われ、フェーズエンコード勾配パルスＧｐによりスピンのフェーズエンコードが行われる。
【００４２】
スピン反転後、図３（ｃ）に示すように、勾配駆動部２３により勾配コイル部２１３ａ，２１３ｂに対してリードアウト勾配パルスＧｒが印加されて、リフェーズされて、図３（ｅ）に示すように、スピンエコーＭＲが発生される。
このスピンエコーＭＲは、データ収集部２４によりビューデータとして収集される。
【００４３】
制御部２５は、このようなパルスシーケンスで、実行プロトコルに応じて、周期ＴＲでたとえば６４〜５１２回繰り返すように、ＲＦ駆動部２２、勾配駆動部２３、およびデータ収集部２４を制御する。
また、制御部２５は、繰り返しのたびに、フェーズエンコード勾配パルスＧｐを変更し、毎回異なるフェーズエンコードを行うように、制御を行う。
【００４４】
オペレータコンソール２Ｂは、図２に示すように、データ処理部２７、操作部２８、および表示部２９を有している。
【００４５】
データ処理部２７は、データ収集部２４から取り込んだデータをメモリに記憶する。メモリ内にはデータ空間が形成される。メモリに形成されるデータ空間は、２次元フーリエ空間を構成する。データ処理部２７は、これら２次元フーリエ空間のデータを２次元逆フーリエ変換、すなわちフーリエ周波数空間から実空間への変換を行って、被検体６の画像を生成（再構成）する。そして、データ処理部２７は、再構成画像を信号Ｓ２として画像整合装置４に供給する。なお、２次元フーリエ空間をｋスペースともいう。
【００４６】
データ処理部２７には、制御部２５が接続されており、制御部２５の上位にあってそれを統括する。
データ処理２７には、また、操作部２８、および表示部２９が接続されている。
【００４７】
操作部２８は、ポインティングデバイスを備えたキーボードやマウス等により構成され、オペレータＯＰの操作に応じた操作信号をデータ処理部２７に出力する。また、操作部２８からは、たとえば上述した実行すべきプロトコルの入力が行われる。データ処理部２７は、操作部２８から入力されたプロトコルに関する情報（プロトコル番号等）を制御部２５に供給する。
【００４８】
表示部２９は、グラフィックディスプレイ等により構成され、操作部２８からの操作信号に応じて、本体装置２Ａの動作状態に応じた所定の情報を表示する。
【００４９】
超音波撮像装置３は、超音波プローブを被検体の被検部位に当接させて被検体に超音波を送波し、被検体内での非線形効果による高調波エコーに基づいて画像を生成し、生成した画像データを信号Ｓ３として画像整合装置４に供給するとともに、表示部に表示する。
【００５０】
図４は、本実施形態に係る超音波撮像装置３の構成例を示す図である。
【００５１】
超音波撮像装置３は、図４に示すように、超音波プローブ３１、送受信部３２、Ｂモード処理部３３、画像処理部３４、表示部３５、制御部３６、および操作部３７を有している。
【００５２】
超音波プローブ３１は、被検体６に当接されて超音波の送受波に使用される。超音波プローブ３１は、図示しない超音波トランスデューサアレイ(transducerarray)を有する。超音波トランスデューサアレイは、複数の超音波トランスデューサをアレイ状に配列して構成される。個々の超音波トランスデューサは、たとえばＰＺＴ（チタン(Ti)酸ジルコン(Zr)酸鉛）セラミックス(ceramics)等の圧電材料で構成される。
超音波プローブ３１は、ケーブル３８により送受信部３２に接続されている。
【００５３】
送受信部３２は、超音波ブローブ３１に駆動信号を与えて、超音波を送波させ、また、超音波プローブ３１が受波したエコーを受信する。
【００５４】
図５は、図４の送受信部３２の構成例を示す図である。
送受信部３２は、図５に示すように、送波タイミング発生回路３２１、送波ビームフォーマ３２２、送受切換回路３２３、および受信ビームフォーマ３２４を有している。
【００５５】
送波タイミング発生回路３２１は、送波タイミング信号Ｓ３２１を周期的に発生して送波ビームフォーマ３２２に供給する。
送波ビームフォーマ３２２は、送波タイミング信号Ｓ３２１に基づいて、送波ビームフォーミング信号、すなわち、超音波トランスデューサアレイ中の送波アパーチャ(aperture)を構成する複数の超音波トランスデューサを時間差をもって駆動する複数の駆動信号Ｓ３２２を発生し、送受切換回路３２３に出力する。
【００５６】
送受切換回路３２３は、複数の駆動信号Ｓ３２２を超音波プローブ３１の超音波トランスデューサアレイにケーブル３８を介して送信する。アレイ中の送波アパーチャを構成する複数の超音波トランスデューサは、複数の駆動信号の時間差に対応した位相差を持つ複数の超音波をぞれぞれ発生する。それら超音波の波面合成により超音波ビームが形成される。
超音波ビームの送波は、送波タイミング発生回路３２１が発生する送波タイミング信号Ｓ３２１により、所定の時間間隔で繰り返し行われる。
そして、超音波ビームの方位は送波ビームフォーマ３２２によって順次変更される。
それにより、被検体６の内部が、超音波ビームが形成する音線によって走査される。すなわち被検体６の内部が音線順次で走査される。
【００５７】
また、送受切換回路３２３は、超音波トランスデューサアレイ中の受波アパーチャが受波した複数のエコー信号を受波ビームフォーマ３２４に入力する。
受波ビームフォーマ３２４は、複数の受波エコーに時間差を付与して位相を調整し次いでそれら加算して、音線に沿ったエコー受信信号の形成、すなわち、受波のビームフォーミングを行う。受波ビームフォーマ３２４により、受波の音線も送波に合わせて走査される。
受波ビームフォーマ３２４は、各音線毎のエコー受信信号をＢモード処理部３３に出力する。
【００５８】
以上の、送波タイミング発生回路３２１、送波ビームフォーマ３２２、送受切換回路３２３、および受信ビームフォーマ３２４は、制御部３６の制御信号ＣＴＬによって制御される。
【００５９】
Ｂモード処理部３３は、送受信部３２の受波ビームフォーマ３２４による各音線毎のエコー受信信号を受けて、Ｂモード画像データを形成する。
【００６０】
図６は、図４のＢモード処理部３３の構成例を示す図である。
Ｂモード処理部３３は、図６に示すように、基本波処理部３３１および高調波処理部３３２を有している。基本波処理部３３１および高調波処理部３３２には、受波ビームフォーマ３２４の出力信号が共通に入力される。
【００６１】
基本波処理部３３１は、基本波エコー、すなわち送波超音波の中心周波数と同じ周波数を持つエコー受信信号を通過させる図示しないフィルタを有する。
高調波処理部３３２は、高調波エコー、すなわち送波超音波の中心周波数のたとえば２次の高調波（第２高調波）と同じ周波数を持つエコー受信信号を通過させる図示しないフィルタを有する。なお、このフィルタは、必要に応じて３次またはそれ以上の高次の高調波に対応するものとしても良いのは勿論である。
【００６２】
基本波処理部３３１は、入力信号につき、基本波エコーを対数増幅および包路線検波することにより、音線上の個々の反射点でのエコーの強度を表す信号、すなわちＡスコープ信号を得て、このＡスコープ信号の各瞬時の振幅をそれぞれ輝度値として、Ｂモード画像データを形成する。すなわち基本波処理部３３１は基本波エコーに基づくＢモード画像データを生成し、画像処理部３４に出力する。
【００６３】
高調波処理部３３２は、入力信号につき、第２高調波エコーを対数増幅および包絡線検波することにより、音線上の個々の反射点でのエコーの強度を表す信号すなわちＡスコープ信号を得て、このＡスコープ信の各瞬時の振幅をそれぞれ輝度値として、Ｂモード画像データを形成する。すなわち高調波処理部３３２は、第２高調波エコーに基づくＢモード画像データをそれぞれ生成し、画像処理部３４に出力する。
【００６４】
画像処理部３４は、Ｂモード処理部３３から入力される複数系統のＢモード画像データに基づいて複数のＢモード画像をそれぞれ生成し、生成したＢモード画像データを信号Ｓ３として画像処理装置４に出力する。
【００６５】
図７は、図４の画像処理部３４の構成例を示す図である。
画像処理部３４は、図７に示すように、バス３４１によって接続された音線データメモリ３４２、ディジタル・スキャンコンバータ３４３、画像メモリ３４４、および画像処理プロセッサ３４５を有している。
【００６６】
Ｂモード処理部３３から音線毎に入力された基本波エコーおよび第２高調波エコーによるＢモード画像データは、音線データメモリ３４２にそれぞれ記憶される。音線データメモリ３４２内にはそれぞれの音線デー夕空間が形成される。
【００６７】
ディジタル・スキャンコンバータ３４３は、走査変換により音線データ空間のデータを物理空間のデータに変換する。
ディジタル・スキャンコンバータ３４３によって変換された画像データは、画像メモリ３４４に記憶される。すなわち、画像メモリ３４４は、物理空間の画像データを記憶する。
画像処理プロセッサ３４５は、音線データメモリ３４２および画像メモリ３４４のデータについてそれぞれ所定のデータ処理を施す。
【００６８】
また、画像処理部３４には、表示部３５が接続されている。
表示部３５は、画像処理部３４から画像信号が与えられ、それに基づいて画像を表示する。表示部３５は、たとえばカラー画像が表示可能なグラフィックディスプレイ等によって構成される。
【００６９】
制御部３６は、送受信部３２、Ｂモード処理部３３、画像処理部３４、および表示部３５の各部に制御信号ＣＴＬを与えてその動作を制御する。
また、制御部３６には、被制御の各部から各種の報知信号が入力される。
そして、制御部３６による制御の下で、超音波撮像が遂行される。
さらに、制御部３６には操作部３７が接続されている。
操作部３７は操作者によって操作され、制御部３６に所望の指令や情報を入力する。操作部３７は、たとえばキーボードやその他の操作具を備えた操作パネルで構成される。
【００７０】
画像整合装置４は、たとえばワークステーション等により構成され、ＭＲＩ装置２により事前に取得された、たとえば図８（ａ）に示すような、血管▲１▼、腫瘍▲２▼、体表▲３▼、並びに検査のために被検体内に挿入された生検針体▲４▼を含むの臓器に関する高分解能、高コントラスト分解能のＭＲ画像信号Ｓ２を受けて一旦記憶装置等に蓄積しておき、超音波撮像装置３によりリアルタイムに得られるＵＳ画像信号Ｓ３を受けて、図８（ｂ）に示し詳細は後述するように、ＵＳ画像から臓器の境界などの特徴点（線）を抽出し、ＵＳ画像の特徴から、その情報と整合させるように（つじつまがあうように）、事前に取得されメモリに保持したＭＲ画像を、図８（ｃ）に示すように変形させて、ＭＲ画像とＵＳ画像の画像情報間に矛盾をなくした補正ＭＲ画像を生成し、表示装置５に表示させる。
すなわち、本実施形態のように、同一対象を時間をおいて撮像した場合、それらの複数枚の画像を相互に比較するためには、正確に重ね合せをする必要がある。すなわち、画像整合装置４は、位置合わせという操作を行う。
位置合わせに際しては、ＭＲ画像とＵＳ画像のでの対応関係を明らかにしておく必要がある。
たとえばいくつかの基準点が共通に求められる場合には、これらの座標を用いて後述する座標変換方法に従って座標変換を行い、一方の画像を他方の画像に重ね合わせることができる。
【００７１】
図９は、画像整合装置４の構成例を示す図である。
画像整合装置４は、図９に示すように、事前に取得した高分解能、高コントラスト分解能のＭＲ画像信号Ｓ２を記憶するハードディスク装置等からなる記憶装置４１、超音波撮像装置３によるＵＳ画像信号Ｓ３からＵＳ画像の特徴を抽出する特徴抽出部４２、および特徴抽出部４２で抽出されたＵＳ画像の特徴から、その情報と整合させるように、記憶装置４１に保持したＭＲ画像を変形させて、ＭＲ画像とＵＳ画像の画像情報間に矛盾をなくした補正ＭＲ画像を生成するＭＲ画像補正部４３を有する。
【００７２】
なお、臓器の特徴点（線）は、図８（ａ）に示すような、血管▲１▼、腫瘍▲２▼、体表▲３▼等に基づいて検出する。
特徴抽出部４２は、たとえば位置情報および輝度データに基づいて、抽出体対象領域の形状や位置等の特徴を検出する。
なお、腫瘍▲２▼の特徴パラメータは、形状、形態、位置、エコー、およびテクスチャの５つに大別できる。
たとえば腫瘤の形状や形態の特徴パラメータをもとめるためには、輪郭抽出処理を行う。
エコーに関する特徴は、たとえば腫瘤内部、外側、後部、後部外側陰影、および同じ深さの濃度平均値と分散を求め、それらの差や比を一つの特徴パラメータとする。
また、テクスチャは、腫瘤の種類により内部組成が異なることから違いがある。テクスチャの特徴パラメータとしては、一様性、コントラスト、エントロピー、濃度相関が用いられる。テクスチャを定量化するには、たとえばフーリエ変換法、同時共起行列法、フラクタル法などが用いられる。
また、特徴パラメータである位置パラメータは、たとえば皮膚から臓器の体表▲３▼までの距離等が用いられる。
【００７３】
ＭＲ画像補正部４３がＭＲ画像を矛盾がないようにＵＳ画像に整合させる変形処理においては、以下に示すような位置合わせが行われる。
【００７４】
すなわち、上述したように本実施形態のように、同一対象を時間をおいて撮像したり、あるいは異なる撮像系や波長を用いて撮像した場合、それらの複数枚の画像を相互に比較するためには、正確に重ね合せをする必要がある。
このような操作を位置合わせという。
位置合わせに際しては、画像間での対応関係を明らかにしておく必要がある。いくつかの基準点が共通に求められる場合には、これらの座標を用いて後述する座標変換方法に従って座標変換を行い、一方の画像を他方の画像に重ね合わせることができる。
基準点が明確でない場合には、一方の画像の小領域が他方の画像内のどの領域に最も似ているかを検出することによって２枚の画像間での対応付けを行う。
基本的には、２枚の画像が平行移動で重なる関係にある場合に有効である。若干のひずみがあっても適用可能であるが、たとえば大きさが異なるなどの場合には上述したように、画像の特徴を抽出して特徴記述間のマッチングを行う。
領域域間での類似度を判定して対応点（領域）の検出を行う方法は、テンプレートマッチングと呼ばれる。
テンプレートマッチング法では  基準とする画像内の小領域をテンプレートとし、他方の画像内でこのテンプレートと最も類似した領域を探し出す（マッチング）ことが行われる。
【００７５】
座標変換
幾何学的変換（変形）は、ｘ−ｙ座標系で表現された画像ｆ（ｘ，ｙ）を、別の座標系であるｕ−ｖ座標系を用いた画像ｇ（ｕ，ｖ）に変換する座標変換の問題として扱うことができる。
すなわち、図１０（ａ），（ｂ）に示すように、画像２中の座標（ｕ，ｖ）にある画素P ’が、画像１中の座標（ｘ，ｙ）にある画素P に対応するとする。このとき、座標変換を表現する関数ｐ（ｘ，ｙ）、ｑ（ｘ，ｙ）を用いて、次のように表現できる。
【００７６】
【数１】
ｕ＝ｐ（ｘ，ｙ）
ｖ＝ｑ（ｘ，ｙ）                                    …（１）
【００７７】
関数ｐ，ｑとしては、次の与え方がある。
(1) 平行移動、回転拡大・縮小などあらかじめ変換式が与えられている場合、
(2) 基準となる画像( たとえば正方格子状のパターン) を撮像し、撮像系におけるひずみ特性を解析的に求める場合、
(3) 互いに位置を合わせようとする画像相互で対応する点を指定し、これらの対応関係から変換式を推定する場合、
である。
【００７８】
(3) の場合、変換式としてはたとえば次式で表現される多項式表現が用いられる。
【００７９】
【数２】
0003878456.tif 000002
【００８０】
画像１と２とで対応する点の組（ｘk ，ｙk ），（ｕk ，ｖk ）を複数求め、式（２）に代入してａij、ｂijに関する連立方程式をたてる。これを最小２乗法を用いて解くことにより係数ａij、ｂijの値が求まり変換式が確定する。
【００８１】
次に、上記構成による動作を説明する。
【００８２】
先ず、クッションを介してクレードル２６上に載せられた被検体６が、図示しない搬送部によって、本体装置２０のマグネットシステム２１のギャップ２１１内に搬入される。
【００８３】
次に、被検体６の被検部位、たとえば腹部をギャップ２１１内のマグネットセンタに位置させる。このとき、マグネットセンタを含むギャップ２１１内の所定の領域には、主磁場マグネット部２１２による静磁場が形成されている。
【００８４】
そして、オペレータＯＰにより、被検部位に対応したプロトコル情報が操作部２８から入力される。
操作部２８から入力されたプロトコルに関する情報（プロトコル番号等）がデータ処理部３１により制御部２５に供給される。
【００８５】
制御部２５では、オペレータコンソール２Ｂのデータ処理部２７により実行すべきプロトコルの指定があると、オペレータコンソール２Ｂのデータ処理部２７から送られてくる被検体６の被検部位に対応した実行すべきプロトコルに即して、あらかじめ決められた繰り返し時間ＴＲ内において所定のパルスシーケンスが所定回数繰り返される駆動信号ＤＲ１をＲＦコイル部２１４に印加するようにＲＦ駆動部２２が制御され、実行すべきプロトコルに即して、１ＴＲ内に、所定のパターンのパルス信号を勾配コイル２１３に印加するように勾配駆動部２３が制御される。
【００８６】
ＲＦ駆動部２２では、制御部２５の指示に基づいたプロトコル対応の駆動信号ＤＲ１がＲＦコイル部２１４に印加され、勾配駆動部２３では、制御部２５の指示に基づいたプロトコル対応の駆動信号ＤＲ２が勾配コイル部２１３に印加される。
【００８７】
そして、ＲＦコイル部２１４より発生された高周波磁場に基づき、被検部位のスピンが励起され、励起信号の送信を打ち切った後に外部に放射される電磁波が受信コイルで受信される。
【００８８】
これにより、被検体６の被検部位で励起されたスピンが生じる電磁波が磁気共鳴信号として取り出され、これがデータ収集部２４で収集され、検査結果のデータとしてオペレータコンソール２Ｂのデータ処理部２７に出力される。
すなわち、被検部位の撮像が行われる。
【００８９】
データ処理部２７では、データ収集部２４から入力したデータがメモリに記憶され、メモリ内にデータ空間が形成される。データ処理部２７では、これら２次元フーリエ空間のデータを２次元逆フーリエ変換して被検体６の被検部位の画像が生成（再構成）される。
そして、データ処理部２７から再構成画像が信号Ｓ２として画像整合装置４に供給される。
【００９０】
画像整合装置４では、信号Ｓ２として入力したＭＲ画像が記憶装置４１に一旦格納される。
【００９１】
次に、超音波撮像装置３により、ＭＲＩ装置２で撮像した被検体の同一被検部位が撮像される。
【００９２】
すなわち、超音波プローブ３１が被検体６の所望の個所に当接され、操作部３７を操作して撮像が行われる。
撮像は、制御部３６による制御の下で遂行される。
具体的には、たとえばセクタスキャンにより、各音線ごとに超音波ビームが送波され、そのエコーが送受信部３２で受信され、各音線のエコー受信信号に基づき、Ｂモード処理部３３でＢモード画像データが形成される。Ｂモード画像データは、基本波エコーに基づくものと第２高調波エコーに基づくものとがそれぞれ形成され、画像処理部３４の音線データメモリ３４２に記憶される。
画像処理部３４では、画像処理プロセッサ３４５により音線データメモリ３４２の複数系統のＢモード画像データが、デイジタル・スキャンコンバータ３４３で走査変換されて、それぞれ画像メモリ３４４に書き込まれる。
ここで、操作部３７を操作して、これらのＢモード画像を表示部３５に表示させる。そして、表示された基本波エコー像と第２高調波エコー像とを観察し、両画像の比較対照等により診断（検査）が行われる。第２高調波エコー像は、被検体の体表から始まる画像を含むので、基本波エコー像との比較対照を行うのに都合が良い。
また、画像処理部３４からは、生成したＢモード画像データが信号Ｓ３として画像処理装置４に出力される。
【００９３】
画像整合装置４では、上述したように、ＭＲＩ装置２により事前に取得された所定の臓器に関する高分解能、高コントラスト分解能のＭＲ画像データが記憶装置４１に格納されている。
そして、画像整合装置４では、超音波撮像装置３によりリアルタイムに得られるＵＳ画像信号Ｓ３を受けて、特徴抽出部４２によりＵＳ画像から臓器の境界などの特徴点（線）が抽出され、ＭＲ画像補正部４３に出力される。
ＭＲ画像補正部４３では、ＵＳ画像の特徴から、その情報と整合させるように記憶装置４１に格納されている事前に取得されたＭＲ画像が変形されて、ＭＲ画像とＵＳ画像の画像情報間に矛盾をなくした補正ＭＲ画像が生成される。
この補正され、リアルタイムに生成された超音波画像と整合され、両画像間に矛盾がなく、しかも高分解能、高コントラスト分解能のＭＲ画像が表示装置５に表示される。
【００９４】
そして、被検体６の被検部位のデータ収集が完了すると、図示しない搬送部によって、クレードル２６と共に被検体６がギャップ２１１の外に搬出される。
【００９５】
以上説明したように、本実施形態によれば、超音波撮像装置３によりリアルタイムに得られるＵＳ画像信号Ｓ３を受けて、特徴抽出部４２によりＵＳ画像から臓器の境界などの特徴点（線）を抽出し、ＵＳ画像の特徴から、その情報と整合させるように記憶装置４１に格納されている事前に取得されたＭＲ画像を変形させて、ＭＲ画像とＵＳ画像の画像情報間に矛盾をなくした補正ＭＲ画像を生成し、補正ＭＲ画像が表示装置５に表示する画像整合装置４を設けたので、ＭＲＩ装置と超音波撮像装置の画像情報から真に有用な高分解能、高コントラスト分解能のＭＲ画像をリアルタイムに得ることができる利点がある。
【００９６】
【発明の効果】
以上説明したように、本発明によれば、ＭＲＩ装置と超音波撮像装置の画像情報から真に有用な合成画像をリアルタイムに得ることができる利点がある。
【図面の簡単な説明】
【図１】本発明に係る画像診断支援システムの一実施形態を示す構成図である。
【図２】本実施形態に係るＭＲＩ装置の構成例を示す図である。
【図３】スピンエコー法のパルスシーケンスについて説明するためのタイミングチャートである。
【図４】本実施形態に係る超音波撮像装置の構成例を示す図である。
【図５】図４の送受信部の構成例を示す図である。
【図６】図４のＢモード処理部の構成例を示す図である。
【図７】図４の画像処理部の構成例を示す図である。
【図８】本発明に係るＭＲＩ装置によるＭＲ画像、超音波撮像装置によるＵＳ画像、および補正後のＭＲ画像を示す図である。
【図９】本発明に係る画像整合装置の構成例を示す図である。
【図１０】座標変換を説明するための図である。
【符号の説明】
１…画像診断支援システム、２…ＭＲＩ装置、２Ａ…本体装置、２１…マグネットシステム、２１１…ギャップ、２１２…主磁場マグネット部、２１３…勾配コイル部、２１４…ＲＦコイル部、２２…ＲＦ駆動部、２３…勾配駆動部、２４…データ収集部、２５…制御部、２６…クレードル、２Ｂ…オペレータコンソール、２７…データ処理部、２８…操作部、２９…表示部、３…超音波撮像装置、３１…超音波プローブ、３２…送受信部、３２１…送波タイミング発生回路、３２２…送波ビームフォーマ、３２３…送受切換回路、３２４…受波ビームフォーマ、３３…Ｂモード処理部、３３１…基本波処理部、３３２…高調波処理部、３４…画像処理部、３４１…バス、３４２…音線データメモリ、３４３…ディジタル・スキャンコンバータ、３４４…画像メモリ、３４５…画像プロセッサ、３５…表示部、３６…制御部、３７…操作部、４…画像整合装置、４１…記憶装置、４２…特徴抽出部、４３…ＭＲ画像補正部、５…表示装置、６…被検体。
【図１】
0003878456.tif 000003
【図２】
0003878456.tif 000004
【図３】
0003878456.tif 000005
【図４】
0003878456.tif 000006
【図５】
0003878456.tif 000007
【図６】
0003878456.tif 000008
【図７】
0003878456.tif 000009
【図８】
0003878456.tif 000010
【図９】
0003878456.tif 000011
【図１０】
0003878456.tif 000012