(57)【要約】【課題】超音波プローブとして従来の１Ｄ超音波プローブを用い、例えば乳腺等の連続構造物と微小石灰化部分等の微小構造物とを正確に区別し、所望の微小構造物抽出画像をリアルタイムで得ることができる超音波診断装置及びプログラムを提供すること。【解決手段】被検体に対して超音波を送信してエコー信号を受信する超音波プローブ12と、前記エコー信号に基づいて前記被検体について複数の断層像を生成する画像生成ユニット25と、生成された前記断層像を記憶する画像メモリ26と、制御プロセッサ28とを具備する超音波診断装置であって、前記制御プロセッサは次の様な制御を行う。すなわち、前画像メモリ26に記憶された断層像と、現時点の前記超音波プローブ12の位置に対応する断層像と、の間の相関を利用して微小構造物を抽出する画像処理を実行させる制御を行う。この微小構造物抽出画像は所定の形態にて表示手段に表示させる。【選択図】      図１
2010000143.tif 000002
【特許請求の範囲】
【請求項１】
被検体を超音波で走査して断層像を取得する超音波診断装置であって、
前記被検体に対して超音波を送信し、当該被検体からのエコー信号を受信する超音波プローブと、
前記超音波プローブによって受信されたエコー信号に基づいて、前記被検体について複数の断層像を生成する画像生成手段と、
前記画像生成手段によって生成された前記断層像を記憶する記憶手段と、
前記記憶手段に記憶された断層像と、現時点の前記超音波プローブの位置に対応する断層像と、の間の相関を利用して微小構造物を抽出する画像処理を実行し、微小構造物抽出画像を生成する画像処理手段と、
前記画像処理手段によって生成された前記微小構造物抽出画像を所定の形態にて表示する表示手段と、
を具備することを特徴とする超音波診断装置。
【請求項２】
前記画像処理手段による前記画像処理とは、前記断層像間の相関を利用して、スペックルパタンを軽減する処理を含む処理であることを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記画像処理手段による前記画像処理とは、前記断層像間の相関を利用して、前記断層像における高輝度部分を抽出する処理を含む処理であることを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記画像処理手段による前記画像処理とは、前記断層像間の空間的に対応する画素を含む参照領域内の最大画素値を差分する処理を含む処理であることを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記超音波プローブの動きを検出するセンサ手段又は前記超音波プローブの動きを推定する推定手段を含み、
前記画像処理手段は、前記センサ手段又は前記推定手段によって前記超音波プローブの動きの停止が検出又は推定されると、当該停止時に取得した前記断層像を前記画像処理に用いないことを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
前記超音波プローブの動きを検出するセンサ手段又は前記超音波プローブの動きを推定する推定手段を含み、
前記画像処理手段は、前記センサ手段又は前記推定手段によって検出又は推定された前記超音波プローブの移動量に基づいて、前記記憶手段に記憶された断層像のうち前記画像処理に用いる断層像を決定することを特徴とする請求項１に記載の超音波診断装置。
【請求項７】
前記超音波プローブの動きを検出するセンサ手段又は前記超音波プローブの動きを推定する推定手段と、
前記画像処理手段は、前記センサ手段又は前記推定手段によって検出又は推定された前記超音波プローブの移動量に基づいて、前記画像処理に用いる断層像のフレーム間隔を決定する断層像フレーム間隔決定手段と、
を含むことを特徴とする請求項１に記載の超音波診断装置。
【請求項８】
被検体を超音波で走査して断層像を取得する超音波診断装置としてコンピュータを機能させるプログラムであって、
超音波プローブによって前記被検体に対して超音波を送信し、当該被検体からのエコー信号を受信する機能と、
前記エコー信号に基づいて、前記被検体について複数の断層像を生成する画像生成機能と、
前記断層像を記憶する記憶機能と、
前記記憶機能によって記憶された断層像と、現時点の前記超音波プローブの位置に対応する断層像と、の間の相関を利用して微小構造物を抽出する画像処理を実行し、微小構造物抽出画像を生成する画像処理機能と、
前記微小構造物抽出画像を所定の形態にて表示する表示機能と、
をコンピュータに実現させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、組織のエコー信号から、生体臓器中の微小な構造物を抽出し表示する超音波診断装置及びプログラムに関する。
【背景技術】
【０００２】
超音波診断は、超音波プローブを体表から当てるだけの簡単な操作で心臓の拍動や胎児の動きの様子をリアルタイム表示で得られる。また、超音波診断は、安全性が高いため繰り返して検査を行うことができる。さらには、システムの規模がＸ線、ＣＴ、ＭＲＩなど他の診断機器に比べて小さく、ベッドサイドへ移動させての検査も容易に行えるなど簡便な診断手法であると言える。この超音波診断において用いられる超音波診断装置は、それが具備する機能の種類によって様々に異なるが、小型なものは片手で持ち運べる程度のものが開発されており、超音波診断はＸ線などのように被曝の影響がなく、産科や在宅医療等においても使用することができる。
【０００３】
このような種々の長所を持つ超音波診断の一つに、乳がんの早期診断がある。乳房組織においては、乳癌の徴候として微小石灰化が発生するケースが多いことが知られている。微小石灰化病変は１個あるいは数個で局所に散在する。石灰は生体組織に比べ硬いため、超音波をよく反射するため画像上高輝度となることが期待される。しかしながら、実際に画像の中から目視する場合には、数百ミクロン程度であっても抽出するのは難しいと言われている。
【０００４】
なお、超音波画像上には、超音波のランダムな干渉に起因するスペックルパタンと呼ばれる干渉縞が発生する場合がある。このスペックルパタンは、肝硬変の診断等に利用される。
【０００５】
一方、このスペックルパタンは、例えば上記乳癌検診の場合には見落とされがちな微小石灰化等の微小構造物と酷似しており、場合によっては紛らわしい画像情報となる。従って、乳癌診断においてはスペックルパタンを除去したいというニーズがある。
【０００６】
このような事情に鑑みて、例えば空間コンパウンド、ＣＦＡＲ（Ｃｏｎｔｒａｓｔ  Ｆａｌｓｅ  Ａｌａｒｍ  Ｒａｔｅ：一定誤警報率）処理、類似度フィルタ等の技術が提案されている。ここで、空間コンパウンドとは、異なる方向からの送受信信号を重畳しスペックルを平滑化するものである。ＣＦＡＲ処理とは、対象画素を周囲の輝度平均で減算し、これを用いて高輝度部分を抽出するものである。類似度フィルタとは、その統計的性質を利用して、スペックルを除去するものである。また、これらのスペックルパタン除去の手法の他、超音波診断の分野ではないが、微小石灰化を自動認識する試みが、主にＸ線診断画像の応用として種々報告されている。
【０００７】
ところで、診断対象の乳腺は乳管などの構造が複雑であり、もとより均質な臓器ではない。従って、従来のフィルタ処理を行うと、微小石灰化が検出されると同時に、乳腺構造も（構造物として）抽出されてしまい、両者を十分に区別することができない。
【０００８】
また、乳管などは微小石灰化に比べて明らかに大きな構造物であるため、フィルタ処理にて残存しても、目視にて弁別が可能となることが期待されるが、それでも弁別が困難となることを、発明者らは研究でしばしば経験している。特に乳腺構造の一部のみが残存する場合は、フィルタ後の画像は点状に見えるため、微小石灰化に類似した画像となる場合がある。
【０００９】
さらに、画像上のスペックルパタンがランダムに変動する場合がある。係る場合には、所定のスペックルリダクション処理を行ったとしても、スペックルが残存してしまい、石灰化部位と区別が困難になる。
【００１０】
このような事情に鑑みて、例えば特許文献１には次のような技術が提案されている。すなわち、特許文献１には、第１の超音波画像及び当該第１の超音波画像の位置に基づいて決定される第２の超音波画像を用いて微小構造物を抽出する画像処理であって、前記第１の超音波画像上の各画素につき、前記第２の画像において空間的に対応する画素を含む参照領域内の最大画素値を差分する微小構造物抽出処理を実行することで、微小構造物抽出画像を生成する画像処理手段と、前記微小構造物抽出画像を所定の形態にて表示する表示手段と、を具備することを特徴とする超音波診断装置が開示されている。
【００１１】
この特許文献１に開示されている超音波診断装置によれば、微細構造物を抽出する信号処理において、断層像一面を利用するのではなく、空間的３次元情報を利用して信号処理を行い、例えば乳腺等の連続構造物と微小石灰化部分等の微小構造物とを正確に区別し、微小構造物を抽出することができる超音波診断装置等が提供される。
【特許文献１】特開２００７−３１３１１４号公報
【発明の開示】
【発明が解決しようとする課題】
【００１２】
ところで、特許文献１に開示されている技術で用いられる超音波プローブは、被検体の３次元領域を超音波走査可能なものである。そのため、特許文献１に開示されている技術で用いられる超音波プローブは、振動子をその配列方向の直交方向に沿って機械的に揺動させ、３次元領域を超音波走査する構成の超音波プローブ（所謂メカ４Ｄ超音波プローブ）や、二次元的に配列された二次元振動素子を用いて電気的制御により３次元領域を超音波走査する構成の超音波プローブ（所謂リアルタイム３Ｄ超音波プローブ）等である。
【００１３】
前記メカ４Ｄ超音波プローブでは、被検体の３次元的走査は前記揺動回路によって行われるため、検査者は超音波プローブ本体を被検体に接触させるだけで、自動的に複数の二次元断層像を取得することができる。また、制御された揺動速度から断面間の正確な距離も検知できる。他方、前記リアルタイム３Ｄ超音波プローブでは、原理的には、従来の二次元断層像を取得するのと同じ時間で、３次元領域を超音波走査することができる。
【００１４】
しかしながら、前記メカ４Ｄ超音波プローブは、その大きさや重さに起因して、微細な構造物を捉える為の走査が困難であるという問題や、リアルタイム性が不十分である等の問題を抱えている。また、前記リアルタイム３Ｄ超音波プローブについては、開発に更なる時間を要することが予想されている。
【００１５】
従って、あくまでも現状においては、超音波プローブとして従来の１Ｄ超音波プローブ（１．５Ｄ超音波プローブも含む）を用いて、臨床的に許容レベルの診断結果を得る（微小構造物の抽出を行うことができる）技術が強く望まれている。換言すれば、超音波プローブとして最も普及している１Ｄ超音波プローブを用い、リアルタイムで所望の微小構造物抽出画像を得ることが可能な技術が望まれている。
【００１６】
本発明は、前記の事情に鑑みて為されたものであり、超音波プローブとして従来の１Ｄ超音波プローブを用い、例えば乳腺等の連続構造物と微小石灰化部分等の微小構造物とを正確に区別し、所望の微小構造物抽出画像をリアルタイムで得ることができる超音波診断装置及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１７】
前記の目的を達成するために、本発明の第１の態様による超音波診断装置は
被検体を超音波で走査して断層像を取得する超音波診断装置であって、
前記被検体に対して超音波を送信し、当該被検体からのエコー信号を受信する超音波プローブと、
前記超音波プローブによって受信されたエコー信号に基づいて、前記被検体について複数の断層像を生成する画像生成手段と、
前記画像生成手段によって生成された前記断層像を記憶する記憶手段と、
前記記憶手段に記憶された断層像と、現時点の前記超音波プローブの位置に対応する断層像と、の間の相関を利用して微小構造物を抽出する画像処理を実行し、微小構造物抽出画像を生成する画像処理手段と、
前記画像処理手段によって生成された前記微小構造物抽出画像を所定の形態にて表示する表示手段と、
を具備することを特徴とする。
【００１８】
前記の目的を達成するために、本発明の第２の態様によるプログラムは、
被検体を超音波で走査して断層像を取得する超音波診断装置としてコンピュータを機能させるプログラムであって、
超音波プローブによって前記被検体に対して超音波を送信し、当該被検体からのエコー信号を受信する機能と、
前記エコー信号に基づいて、前記被検体について複数の断層像を生成する画像生成機能と、
前記断層像を記憶する記憶機能と、
前記記憶機能によって記憶された断層像と、現時点の前記超音波プローブの位置に対応する断層像と、の間の相関を利用して微小構造物を抽出する画像処理を実行し、微小構造物抽出画像を生成する画像処理機能と、
前記微小構造物抽出画像を所定の形態にて表示する表示機能と、
をコンピュータに実現させることを特徴とする。
【発明の効果】
【００１９】
本発明によれば、超音波プローブとして従来の１Ｄ超音波プローブを用い、例えば乳腺等の連続構造物と微小石灰化部分等の微小構造物とを正確に区別し、所望の微小構造物抽出画像をリアルタイムで得ることができる超音波診断装置及びプログラムを提供することができる。
【発明を実施するための最良の形態】
【００２０】
［第１実施形態］
図１は、第１の実施形態に係る超音波診断装置１の構成を示すブロック図である。同図に示すように、本超音波診断装置１は、超音波プローブ１２と、装置本体１１と、入力装置１３と、モニター１４と、を具備する。
【００２１】
超音波プローブ１２は、超音波送信ユニット２１からの駆動信号に基づき超音波を発生し、被検体からの反射波を電気信号に変換する複数の圧電振動子、当該圧電振動子に設けられる整合層、当該圧電振動子から後方への超音波の伝播を防止するバッキング材等を備える。当該超音波プローブ１２から被検体Ｐに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ１２に受信される。このエコー信号の振幅は、反射することになった反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分を依存して、周波数偏移を受ける。
【００２２】
前記装置本体１１は、超音波送信ユニット２１と、超音波受信ユニット２２と、Ｂモード処理ユニット２３と、ドプラ処理ユニット２４と、画像生成ユニット２５と、画像メモリ２６と、画像合成部２７と、制御プロセッサ（ＣＰＵ）２８と、内部記憶部２９と、インターフェース部３０と、を有する。以下、個々の構成要素の機能について説明する。
【００２３】
前記超音波送信ユニット２１は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を備える。パルサ回路では、所定のレート周波数ｆｒ  Ｈｚ（周期；１／ｆｒ秒）で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。トリガ発生回路は、このレートパルスに基づくタイミングで、超音波プローブ１２に駆動パルスを印加する。
【００２４】
なお、前記超音波送信ユニット２１は、制御プロセッサ２８の指示に従って所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に送信駆動電圧の変更については、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は複数の電源ユニットを電気的に切り替える機構によって実現される。
【００２５】
前記超音波受信ユニット２２は、図示していないアンプ回路、Ａ／Ｄ変換器、加算器等を備える。アンプ回路では、超音波プローブ１２を介して取り込まれたエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【００２６】
前記Ｂモード処理ユニット２３は、超音波受信ユニット２２からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。このデータは、画像生成ユニット２５に送信され、反射波の強度を輝度にて表したＢモード画像としてモニター１４に表示される。
【００２７】
前記ドプラ処理ユニット２４は、超音波受信ユニット２２から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。
【００２８】
前記画像生成ユニット２５は、一般的には、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示画像としての超音波診断画像を生成する。
【００２９】
また、画像生成ユニット２５は、スキャンコンバート以外の種々の画像処理を実行する。すなわち、画像生成ユニット２５は、後述する微小構造物抽出処理の他、例えば、スキャンコンバート後の複数の画像フレームを用いた輝度の最大値画像を再生成する方法（平滑化処理）、画像間の差分を用いる方法（差分演算）等を実行する。なお、当該画像生成ユニット２５に入る以前のデータは、「生データ」と呼ばれることがある。
【００３０】
前記画像メモリ（シネメモリ）２６は、例えばフリーズする直前の複数フレームに対応する超音波画像を保存するメモリである。この画像メモリ２６に記憶されている超音波画像を連続表示（シネ表示）することで、超音波動画像を表示することも可能である。
【００３１】
前記画像合成部２７は、画像生成ユニット２５又から受け取った超音波画像を種々のパラメータの文字情報や目盛等と共に合成し、ビデオ信号としてモニター１４に出力する。
【００３２】
前記制御プロセッサ２８は、情報処理装置（計算機）としての機能を持ち、本超音波診断装置本体の動作を制御する。制御プロセッサ２８は、内部記憶部２９から微小構造物抽出機能を実現するための専用プログラム、所定の画像生成・表示等を実行するための制御プログラムを読み出して自身が有するメモリ上に展開し、各種処理に関する演算・制御等を実行する。
【００３３】
前記内部記憶部２９は、所定のスキャンシーケンス、各実施形態に係る微小構造物抽出機能を実現するための専用プログラム、画像生成、表示処理を実行するための制御プログラム、診断情報（患者ＩＤ、医師の所見等）、診断プロトコル、送受信条件、ＣＦＡＲ処理制御プログラム、ボディマーク生成プログラムその他のデータ群が保管されている。また、必要に応じて、画像メモリ２６中の画像の保管などにも使用される。内部記憶部２９のデータは、インターフェース部３０を経由して外部周辺装置へ転送することも可能となっている。
【００３４】
前記インターフェース部３０は、入力装置１３、ネットワーク、新たな外部記憶装置（図示せず）に関するインターフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インターフェース部３０よって、ネットワークを介して他の装置に転送可能である。
【００３５】
前記入力装置１３は、装置本体１１に接続され、オペレータからの各種指示、条件、関心領域（ＲＯＩ）の設定指示、種々の画質条件設定指示等を装置本体１１にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等を有している。例えば、操作者が入力装置１３の終了ボタンやＦＲＥＥＺＥボタンを操作すると、超音波の送受信は終了し、当該超音波診断装置は一時停止状態となる。
【００３６】
前記モニター１４は、画像生成ユニット２５からのビデオ信号に基づいて、生体内の形態学的情報（Ｂモード画像）、血流情報（平均速度画像、分散画像、パワー画像等）、これらの組み合わせを画像として表示する。
【００３７】
なお、本超音波診断装置１には、必要に応じて、位置検出装置１５を設けることが可能である。この装置は、撮影対象（すなわち、診断部位）に対する超音波プローブ１２の位置に関する情報を検出するものである。ここで、超音波プローブ１２の位置に関する情報とは、超音波プローブ１２の絶対位置情報、相対位置情報、超音波プローブ１２の可動前の位置情報と可動速度及び時間、その他の走査時における超音波プローブ１２の位置を特定するための情報である。
【００３８】
図２は、位置検出装置１５の一例を示した図である。同図に示すように、位置検出装置１５は、可動ステージ１５０、駆動部１５１を有している。可動ステージ１５０には、専用アダプタを介して超音波プローブ１２を設置することができる。駆動部１５１は、制御プロセッサ２８からの制御に基づいて、設置された超音波プローブ１２を可動ステージ１５０に沿って移動させる。また、駆動部１５１は、内部にロータリーエンコーダ等を有しており、可動ステージ１５０上における超音波プローブ１２の位置を検出し、その検出結果を制御プロセッサ２８に逐次送信する。
【００３９】
超音波画像の取得に際しては、診断対象である乳房は、水槽１７に水深させた状態で所定の位置に配置され、検査中は動かないよう固定される。また、超音波プローブ１２及び位置検出装置１５は、水槽１７の底面側に配置される。制御プロセッサ２８は、超音波プローブ１２が所定の速度で移動するように駆動部１５１を制御させつつ超音波送受信を実行することで、自走式超音波走査を実行する。超音波プローブ１２からの画像は、実施例１同様に装置本体に送られる。また、駆動部１５１から取得する位置情報は、後述する超音波プローブ位置に関する情報生成にリアルタイムで利用されると共に、各フレームの付帯情報として記入され管理される。
【００４０】
以下、本第１実施形態に係る超音波診断装置及びプログラムの主な特徴の一つである微小構造物の抽出処理について詳細に説明する。
【００４１】
超音波断層像（Ｂモード像）で微小な構造物と思われる輝点（以下、単に輝点）が観察されたとき、実際にはそれが本当に微細石灰化のような構造物であるか、乳腺のような組織構造の一部なのかは、判断しにくい。特に１枚の静止画像では診断することができないと言われている。
【００４２】
しかし両者は、例えば次のような点で異なる。
【００４３】
（相違点１）微細石灰の組成は生体組織よりも硬く、原理的にはより大きな超音波反射信号が得られているはずである。我々の研究結果からも、微細石灰化による輝点の信号レベルは周囲のスペックルパタンの最大値よりもいくらか高いことが判明した。ただし輝度としてモニター表示された場合、この信号レベルの差を目視で判定することは難しいわけである。
【００４４】
（相違点２）微細石灰化は局在する微小構造物であり、一方で乳腺等の生体組織は３次元的に一定の範囲で連続する構造を持つ連続構造物であり、その空間分布の形態において、本質的に大きく異なる。よって奥行き方向の３次元的な連続性も加味すれば、両者の違いが判定できると期待できる。
【００４５】
本第１実施形態では、この点に着目し、例えば乳房、肝臓、膵臓等の診断において、輝度の微小な違い、および空間分布の形態に基づいて両者を区別し、微小構造物が積極的に抽出された画像生成（微小構造物抽出画像）を生成するものである。
【００４６】
なお、本第１実施形態における微小構造物抽出処理は、画像群データをその対象とする。ここで、画像群データとは、複数の二次元画像を有するボリュームデータ、或いは異なる複数の二次元画像から構成されるデータ（完全なボリュームデータを必ずしも構成しない）を意味する。
【００４７】
本第１実施形態においては、説明を具体的にするため、ボリュームデータとしての画像群データを用いた微小構造物抽出処理について説明する。この様な画像群データは、図２に示した自走式走査が可能なデバイスを用いる撮影、或いは、超音波振動素子が１次元的に配列された超音波プローブ（必要に応じて、位置センサが設けられたもの）を用いて、手動によって複数断層を取得する撮影等で取得することができる。
【００４８】
図３は、微小構造物抽出処理を説明するための図である。同図に示すように、画像群データに含まれる複数の二次元画像から、対象画像３１と参照画像３２とを選択する。ここで、前記対象画像３１は、本微小構造物抽出処理の対象とされる画像の一つであり、前記超音波プローブ１２の直下の断層像（ＴｉｍｅＳｍｏｏｔｈ等の処理を施された断層像も含む）である。前記参照画像３２は、前記対象画像３１とは空間的に異なる（例えば、対象画像からｋフレーム離れた）別の断層像であり、本微小構造物抽出処理に利用されるものである。詳細には前記参照画像３２は、前記対象画像３１の前のフレームの画像であり、前記画像メモリ２６に記録されている画像である。これらの画像は、Ｂモード診断同様、超音波プローブ本体から垂直方向の断面であることが望ましい。
【００４９】
図４（ａ）は対象画像の模式図を、図４（ｂ）は参照画像の模式図をそれぞれ示している。本微小構造物抽出処理では、対象画像から参照画像を差し引くことによって、差分画像を生成する。このとき、対象画像上の各画素(xi,yi)（対象画素）に対して、参上画像上に設定される参照領域Ri内に存在する画素の代表値を決定し、これを対象画素(xi,yi)の値から差し引くものとする。
【００５０】
ここで、参照領域とは、参照画像上の座標が対象画素と同一な画素（対応画素）を含むように、任意のサイズで参照画像上に設定されるものである。なお、参照領域Riの代表値は、参照領域Riの特徴を代表するものであればどんなものでもよい。具体例としては、最大値、平均値、中央値等を挙げることができる。本第１実施形態においては、代表値として最大値を採用するものとする。この差分画像の生成は、次式（１）に示すように表すことができる。
【００５１】
Qi(xi,yi) = Pi(xi,yi)  -  MAX [ Pi-k (xi+m,yi+n) ]  ・・・（１）
ここで、Qi(xi,yi)は差分画像の各画素の値、Pi(xi,yi)は対象画像上の対象画素の値、Pi-k (xi,yi)はPi(xi,yi)と空間的に対応する位置にある参照画像上の各画素の値、Pi-k (xi+m,yi+n)は参照領域Ri内の各画素の値、m、nは参照領域の大きさを指定する任意の値、MAX[]は[]内から最大値を選択する演算をそれぞれ示している。従って、図４（ｂ）は、x軸方向及びy軸方向共に±２画素、計２５画素からなる参照領域Riを例示したものとなっている。なお、この減算処理において、結果が負数の場合は全てゼロ（輝度値＝黒）とするのが望ましい。
【００５２】
本第１実施形態における微小構造物抽出処理によって生成される差分画像は、連続構造物及びランダムなスペックルパタンを除去し、微小構造物を好適に映像化するものとなる。これは、以下の理由による。
【００５３】
すなわち、画像群データを構成する二次元超音波画像上に残存する連続構造物は、微小石灰より大きな構造物であり、奥行き方法の連続性を持っていると考えられる。この観点より、超音波断層像のある１点に注目する場合、連続構造物は、例えば図５に示すような信号強度（画像の輝度）（ａ）の変化が緩やかであることが予想される。
【００５４】
他方、微小構造物では、図５に示す信号強度（ｂ）のごとく、特定の画像にのみ含まれることが予想される。従って、連続或いは近接するフレーム画像から差分画像を生成（例えば、画像Fi（対象画像）から画像Fi-ｆ（参照画像）を減算）すると、当該差分画像により、不連続に存在する微小構造物（例えば微小石灰化部分）のみ抽出することができる。
【００５５】
以下、微小構造物抽出処理によって得られた微小構造物抽出画像の表示形態について説明する。なお、以下に述べる各実施例に係る表示形態は、単独で、或いは複数を組み合わせて用いることが可能である。
【００５６】
（表示形態例１）
本表示形態例１においては、微小構造物抽出画像としての差分画像を、当該差分画像の生成に用いた対象画像の取得時の超音波プローブ位置を示す情報と共に表示するものである。超音波プローブ位置を示す情報は、その目的を果たすものであればどの様なものであってもよいが、典型例としては、図６に示すような、ボディマーク上に設定される超音波プローブ１２の模式図を挙げることができる。
【００５７】
この様な超音波プローブ位置を示す情報は、例えば図２に示した位置検出装置１５により検出さる超音波プローブ位置情報に基づいて生成することができる。図６のボディマークの場合、画像合成部２７は、制御プロセッサ２８の制御に従って、超音波プローブ位置が示されたボディマークを生成し、差分画像と合成してモニター１４に送り出す。これにより、図６に示した形態で、超音波プローブ位置を示す情報と共に差分画像を表示することができる。なお、必要に応じて、画像群データを構成する二次元画像の全ての超音波プローブ位置情報に基づいて、超音波プローブ１２の走査範囲或いは既に表示済みの領域を「軌跡」として、ボディマーク上に色別表示するようにしてもよい。
【００５８】
（表示形態例２）
本表示形態例２においては、微小構造物抽出処理によって得られる複数の差分画像（例えば、画像群データに対応する差分画像）を用いて図７に示すようなＭＩＰ処理（Maximan Intensity Projection：最大値投影処理）を実行し、その結果得られるＭＩＰ画像を微小構造物抽出画像として表示するものである。これにより、複数の差分画像に含まれる微小構造物の情報を、１枚のＭＩＰ画像に凝縮することができる。このＭＩＰ画像を電子カルテ等の添付データとして利用することで、管理データサイズの縮減を図ることができる。
【００５９】
（表示形態例３）
本表示形態例３においては、定量的解析を用いて画像群データから一定の信頼性を持つ差分画像を抽出し、これを用いて上記実施例２に係るＭＩＰ処理を実行するものである。すなわち、画像群データに対応する差分画像の各画素について輝度曲線を生成し、これを用いてある期間（例えばフレーム間隔）での時間的変化量及びその標準偏差を算出する。得られた結果のうち著しく異なる値を示す標準偏差（例えば、所定の閾値以上を持つ標準偏差）に対応する画素は、微小構造物である可能性が高いということができる。従って、当該画像を有する差分画像を抽出し、これらを用いたＭＩＰ処理を行うことで、微小構造物の抽出精度を向上させることができる。
【００６０】
（表示形態例４）
本表示形態４においては、微小構造物抽出処理前のＢモード画像、微小構造物抽出画像、差分画像を用いたＭＩＰ処理によって得られるＭＩＰ画像を、重畳表示、Ｄｕａｌ表示、Ｔｒｉｐｌｅｘ表示のいずれかの形態で表示することができる。重畳表示の場合、スペックルパタン除去前のＢモード画像と、除去後の新たな画像は、基本色をそれぞれ変えて重畳することで、両者の識別が可能となる。また、Ｄｕａｌ表示のような並列表示において、異なる種類の画像が同時に表示される各表示形態においては、各画像において同一の位置に対応するようにカーソルが配置される。従って、医師等の観察者は、目的に応じて所望の表示形態、所望のタイミングにて微小構造物抽出画像を表示できると共に、複数種類の画像によって微小構造物を迅速且つ簡単に特定し観察することができる。
【００６１】
以下、本第１実施形態に係る超音波診断装置及びプログラムによる微小構造物抽出処理を、図８に示すフローチャートを参照して説明する。なお、説明を具体的にするため、図８に示す例では、上記表示形態例１に係る表示形態を採用するものとする。同図に示すように、まず、ボリューム走査開始に伴って微小構造物抽出処理へ遷移する指令を受けて、必要なパラメータ群が読み込まれる（ステップＳ１、Ｓ２）。ここで、微小構造物抽出処理に必要なパラメータ群とは、参照画像の数、対象画像との距離、平滑処理（最大値演算）の領域等である。
【００６２】
次に、所定の手法を用いたボリュームスキャンによって乳房に関する画像群データが取得され、メモリに格納される（ステップＳ３）。その後、画像生成ユニット２５は、参照画像毎の参照領域の代表値を算出し（ステップＳ４）、これを用いて既述の微小構造物抽出処理を実行することで、画像群データに対応する複数の差分画像を生成する（ステップＳ５）。得られた差分画像は、モニター１４に例えば超音波プローブ位置情報を有するボディマークと共に表示され、自動的に保存される（ステップＳ６）。
【００６３】
その後、画像フリーズもしくは本発明の映像モードを終了する指令が行われるまで、診断装置は本微小構造物抽出処理を繰り返し実行する。
【００６４】
なお、上述した本第１実施形態に係る超音波診断装置による一連の処理は、プログラム化することで、或いはプログラム化した後当該プログラムを記憶媒体に読み込むことによって、当該医用画像診断装置とは独立したソフトウェア製品単体としての販売、配布も容易になり、また本第１実施形態に係る技術を他のハードウェア上で利用することも可能となる。
【００６５】
以上説明したように、本第１実施形態によれば、超音波プローブとして従来の１Ｄ超音波プローブを用い、例えば乳腺等の連続構造物と微小石灰化部分等の微小構造物とを正確に区別し、所望の微小構造物抽出画像をリアルタイムで得ることができる超音波診断装置及びプログラムを提供することができる。
【００６６】
すなわち、本第１実施形態に係る超音波診断装置及びプログラムによれば、１Ｄ超音波プローブを利用してリアルタイムに、連続構造物と微細構造物との空間的広がり（空間（３Ｄ）情報）の相違を利用して、微細構造物の抽出を行うことができる。換言すれば、リアルタイム性を失わず且つ広く普及している１Ｄ超音波プローブを用いて、断層像内のみでの信号処理で取り除ききれなかった擬陽性を大幅に削減することができる。
【００６７】
具体的には、本第１実施形態に係る超音波診断装置によれば、超音波プローブとして従来の１Ｄ超音波プローブ（１．５Ｄ超音波プローブも含む）を用いて、例えば乳房、肝臓、膵臓等の診断において、不連続に存在する微小構造物を、当該画像と実質的に直交する方向（奥行き方向）関する情報をも用いて、リアルタイムに抽出する。なお、微小構造物抽出処理における最大値平滑化は、対象画像と参照画像との単なる差分では除去しきれない、スペックルパタンの変動や、構造物の断面方向の位置ずれによる残存をも効果的に除去することを可能とする。
【００６８】
図９（ａ）は、通常のＢモード画像を示した図であり、図９（ｂ）は、微小構造物抽出画像を示す図である。
【００６９】
双方を比較した場合、図９（ａ）に示す通常のＢモード画像では、微小石灰化部位に加えて、組織の一部をも映像化されてしまっており、点状の高輝度部位が多く散在している。従って、どの点が微小石灰化部位に対応するのかを、肉眼で識別することはできない。
【００７０】
一方、図９（ｂ）に示す微小構造物抽出画像では、微小石灰化部位のみが抽出され、点状の高輝度部位として映像化されている。従って、３次元的に連続して分布する連続構造物と局在する微小構造物とを弁別し、微小構造物が抽出された微小構造物抽出画像をリアルタイムで生成し表示することができる。医師等は、この微小構造物抽出画像を観察することで、目視ではスペックルパタンとの分別が難しく、且つある特定の断面像にのみ現れるような微小構造物を、短時間で発見することが可能となる。
【００７１】
なお、本第１実施形態に係る超音波診断装置では、微小構造物抽出処理に用いる参照画像とするフレーム及び参照領域のサイズを任意に選択することができる。従って、検査目的や個体差に応じた参照画像とするフレーム及び参照領域のサイズを設定することで、個々の状況に対応させた微小構造物の好適な映像化が可能である。
【００７２】
また、本第１実施形態に係る超音波診断装置では、対象画像取得時の超音波プローブ位置及び走査範囲が設定されたボディマークと共に微小構造物抽出画像を表示する、微小構造物抽出処理で得られた差分画像を用いて生成されるＭＩＰ画像を所定の形態で表示する、微小構造物抽出前後の画像を所定の形態で表示する等、種々の表示形態を採用することができる。
【００７３】
従って、医師等は、所望の表示形態で微小構造物抽出画像を観察、又は種々の表示形態による微小構造物抽出画像を比較することで、目視ではスペックルパタンとの分別が難しく、かつある特定の断面像にのみ現れるような微小構造物を、リアルタイムで発見することが可能となる。
【００７４】
［第２実施形態］
以下、本発明の第２実施形態に係る超音波診断装置及びプログラムについて説明する。なお、前記第１実施形態と重複した説明を避ける為、前記第１実施形態との相違点のみを説明する。
【００７５】
本第２実施形態に係る超音波診断装置では、図１０に示すように、前記超音波プローブ１２に微小位置センサ１０１を装着又は内蔵させる。この微小位置センサ１０１は、前記超音波プローブ１２の位置の移動量を検出し、フレームレート等に基づいて、目的の微細構造物２０１を抽出するのに適当な断層像フレーム間隔を算出する。そして、このように算出された断層像フレーム間隔に基づいて、第１実施形態において説明した微細構造物の抽出処理を行う。
【００７６】
ところで、本第２実施形態においては、前記微小位置センサ１０１によって前記超音波プローブ１２の移動停止が検出された場合、次のような処理を行う。
【００７７】
（超音波プローブ停止時の処理例１）
前記微小位置センサ１０１によって前記超音波プローブ１２の移動停止が検出された場合、当該時点における前記超音波プローブ１２の直下の断層像のみを使用した処理（超音波プローブ１２が同一位置で取得した断層像を複数使用しない処理）に切り換える。
【００７８】
（超音波プローブ停止時の処理例２）
前記微小位置センサ１０１によって前記超音波プローブ１２の移動停止が検出された場合、その停止期間に取得した断層像を削除し（又は利用せずとし）、前記超音波プローブ１２の移動再開時には、当該停止期間の直前の断層像を利用して上述した処理を行う。
【００７９】
以上説明したように、本第２実施形態によれば、前記第１実施形態に係る超音波診断装置及びプログラムと同様の効果を奏する上に、より精度の高い微細構造物抽出の処理を行うことができる超音波診断装置及びプログラムを提供することができる。
【００８０】
［第３実施形態］
以下、本発明の第３実施形態に係る超音波診断装置及びプログラムについて説明する。なお、前記第１実施形態と重複した説明を避ける為、前記第１実施形態との相違点のみを説明する。
【００８１】
本第３実施形態に係る超音波診断装置では、断層像間の相関情報に基づいて、前記超音波プローブ１２の位置の移動量又は当該断層像間の類似性を算出し、この算出結果に基づいて、断層像目的の微細構造物を抽出するのに適当な断層像フレーム間隔を算出する。そして、このように算出された断層像フレーム間隔に基づいて、第１実施形態において説明した微細構造物抽出の処理を行う。
【００８２】
（超音波プローブの移動量を算出する場合）
算出した超音波プローブ１２のフレームレート等に基づいて、目的の微細構造物を抽出するのに適当な断層像フレーム間隔を算出する。なお、前記超音波プローブ１２の移動停止が検出された場合は、前記第２実施形態に係る超音波診断装置における超音波プローブ１２の移動停止時の処理と同様の処理を行うとする。
【００８３】
（断層像間の類似性を算出する場合）
算出した断層像間の類似性に基づいて、現フレームとその類似性の変化が大きくなるフレーム間で、前記第１実施形態において説明した微細構造物抽出の処理を行う。なお、前記画像メモリ２６に保存された断層像に対して、予め指定した所定の枚数以上において類似性が常に高い場合には、超音波プローブ１２の移動が停止していると判断し、前記第２実施形態に係る超音波診断装置における超音波プローブ１２の移動停止時の処理と同様の処理を行うとする。
【００８４】
以上説明したように、本第３実施形態によれば、前記第１実施形態に係る超音波診断装置及びプログラムと同様の効果を奏する上に、より精度の高い微細構造物抽出の処理を行うことができる超音波診断装置及びプログラムを提供することができる。
【００８５】
以上、第１実施形態乃至第３実施形態に基づいて本発明を説明したが、本発明は上述した実施形態に限定されるものではなく、本発明の要旨の範囲内で種々の変形及び応用が可能なことは勿論である。
【００８６】
さらに、上述した実施形態には種々の段階の発明が含まれており、開示される複数の構成要件の適当な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件からいくつかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果の欄で述べられている効果が得られる場合には、この構成要件が削除された構成も発明として抽出され得る。
【図面の簡単な説明】
【００８７】
【図１】本発明の第１実施形態に係る超音波診断装置の構成を示すブロック図。
【図２】位置検出装置の一例を示した図。
【図３】微小構造物抽出処理を説明するための図。
【図４】図４（ａ）は対象画像の模式図を、図４（ｂ）は参照画像の模式図をそれぞれ示す図。
【図５】画素毎の信号強度（画像の輝度）の変化例を示す図。
【図６】微小構造物抽出画像の表示形態の一例を示す図。
【図７】微小構造物抽出画像の表示形態の他の例を示す図。
【図８】本発明の第１実施形態に係る微小構造物抽出処理のフローチャートを示す図。
【図９】（ａ）は通常のＢモード画像を示す図であり、（ｂ）は微小構造物抽出画像を示す図。
【図１０】本発明の第２実施形態に係る超音波診断装置の超音波プローブを示す図。
【符号の説明】
【００８８】
１…超音波診断装置、  １１…装置本体、  １２…超音波プローブ、  １３…入力装置、  １４…モニター、  １５…位置検出装置、  １７…水槽、  ２１…超音波送信ユニット、  ２２…超音波受信ユニット、  ２３…Ｂモード処理ユニット、  ２４…ドプラ処理ユニット、  ２５…画像生成ユニット、  ２６…画像メモリ、  ２７…画像合成部、  ２８…制御プロセッサ、  ２９…内部記憶部、  ３０…インターフェース部、  ３１…対象画像、  ３２…参照画像、  １０１…微小位置センサ、  １５０…可動ステージ、  １５１…駆動部、  ２０１…微細構造物。
【図１】
2010000143.tif 000003
【図２】
2010000143.tif 000004
【図３】
2010000143.tif 000005
【図４】
2010000143.tif 000006
【図５】
2010000143.tif 000007
【図６】
2010000143.tif 000008
【図７】
2010000143.tif 000009
【図８】
2010000143.tif 000010
【図９】
2010000143.tif 000011
【図１０】
2010000143.tif 000012