(57)【要約】      （修正有）【課題】画像合成手段を備えた超音波診断装置において、コントラスト分解能の向上と空間分解能の向上を両立させる。【解決手段】画像合成手段を備えた超音波診断装置において、コントラスト分解能の向上と空間分解能の向上を両立させるため、合成する画像間で相関をとり、画像間の変位量を計算。変位量が大きい時は包絡線検波後の信号を合成し、変位量が小さい時はＲＦ信号を合成する。もしくは、変位量に応じてミキシング周波数を可変とし、位置ずれの程度に応じて、空間分解能の向上とコントラスト分可能向上のバランスを取る。【選択図】図２
2010042266.tif 000002
【特許請求の範囲】
【請求項１】
被検体に与えるべき送波超音波ビームを形成するための超音波パルス信号を送出し、前記被検体から反射された超音波パルス信号から形成された受波超音波ビームを受信する情報処理装置と、前記超音波パルス信号に応じて送波超音波ビームを形成するための送波超音波ビームフォーマと、該送波超音波ビームフォーマの出力を送受分離器を介して与えられて前記被検体に送波超音波ビームを送出する超音波探触子と、該超音波探触子に受信される前記被検体から反射された超音波パルス信号を送受分離器を介して与えられて前記受波超音波ビームを形成する受波ビームフォーマとよりなり、
前記情報処理装置は被検体に対して複数の角度方向から撮像した前記被検体のＢモード画像データを合成し、合成画像を表示するためのプログラムを処理する処理部を備えるとともに、前記異なる角度での各フレームの前記受波超音波ビームを指定された参照周波数でヘテロダイン検波を行なう超音波撮像装置であって、前記参照周波数の指定が予め操作者が行なう操作により設定されることを特徴とする超音波撮像装置。
【請求項２】
超音波探触子から超音波パルスを被検体に送信した結果、前記被検体から反射される超音波パルスから形成された受信超音波ビームに基づき、複数の角度方向から撮像した前記被検体のＢモード画像データを合成し、合成画像を表示する情報処理装置による超音波撮像方法であって、
前記受波超音波ビームをヘテロダイン検波するための参照周波数の指定が可能であることを特徴とする超音波撮像方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療診断等に用いられる超音波撮像装置に係り、特に、いろいろな角度から走査して得られる画像信号を合成して１つの画像として表示する、一般に、コンパウンド走査像と称される高画質の断層像を得る技術に関する。
【背景技術】
【０００２】
超音波撮像装置は、被検体の表面に超音波探触子を当て、その探触子から被検体に超音波を送信するとともに、被検体内部からの反射波（エコー信号）を受信し、そのエコー信号に基づいて被検体の各部の状態を断層像により表示することにより、診断等に供するものである。
【０００３】
図１（Ａ）は、超音波撮像装置の一般的な構成を示す図である。１０はいわゆるパソコン等の情報処理装置である。１１はＣＰＵ、１２はメモリ、１３はハードディスク、２１は入力装置、２２は表示装置である。これらの装置はバス１４で連結され、必要な情報の伝達を行なう。ハードディスク１３には、超音波撮像装置の動作に必要な各種のプログラムおよびデータが格納されており、必要に応じて呼び出されてメモリ１２に記憶され、ＣＰＵ１１の管理の下に、動作に利用される。入力装置２１は、キーボードあるいはマウス等であり、超音波撮像装置の使用者が必要な設定情報を入力し、あるいは、操作信号を与えるために使用される。表示装置１３には必要な操作画面および撮像結果の画面等が表示される。
【０００４】
２３および２４は超音波の送波ビームフォーマおよび受波ビームフォーマであり、情報処理装置１０のバス１４に接続される。さらに、送波ビームフォーマ２３および受波ビームフォーマ２４は、送受分離器３１を介して、超音波の超音波探触子３２に接続される。使用者の操作信号に応じたプログラムに対応した超音波信号が超音波探触子３２から送出されるように、送波ビームフォーマ２３は情報処理装置１０から与えられるディジタル信号に応じた超音波送信信号を作って、送受分離器３１を介して、超音波の超音波探触子３２に送る。一方、送信された超音波信号が被検体内部からの反射波として超音波探触子３２に受信されると、これに応じた信号が送受分離器３１を介して、受波ビームフォーマ２４に送られる。受波ビームフォーマ２４はこれにディジタル化等の必要な信号処理を施して情報処理装置１０に送る。情報処理装置１０では、受信信号を撮像の目的に対応したプログラムで処理して、結果を表示装置１３に表示する。
【０００５】
図１（Ｂ）は、超音波探触子３２の一例を示す図であり、３分割された直線状に配置された振動子１−１，１−２および１−３がｎ組並列に配列された例である。図では表示を省略したが、これらの振動子には、送信制御回路２３からの送信信号を送り、振動子が検出する被検体内部からの反射波に応じた信号を受信制御回路２４に送るための配線が設けられる。一方、超音波探触子３２の他面には、振動子を被検体に効果的に接触させ、超音波を効率良く被検体との間で送受するための、例えば、金属粉末と高分子材料の混合物からなる音響整合層５とやはり高分子材料からなる音響レンズが設けられている。各振動子に与える信号の大きさおよびタイミングを制御することにより、超音波探触子３２から照射される超音波の焦点深度、照射の向きおよび強度を任意に制御することができ、図に示すＦ1−Ｆ1、Ｆ2−Ｆ2のように深さおよび向きを制御することができる。また、後述する複数の角度方向から撮像した前記被検体のＢモード画像データを合成するために、例えば、焦点Ｆ1−Ｆ1の軸上でビーム方向がＡ，ＢおよびＣの方向になるように制御することができる。なお、超音波探触子３２は、使用者が持つための把手も備えているが、図示は省略した。
【０００６】
超音波探触子３２は、図１（Ｂ）に示す例に限られず、撮像目的に応じて、複数の振動子を等間隔で直線状、曲線状あるいは面状に配列して形成する等、種々の形態で作製される。例えば、リニア走査型の超音波撮像装置は、あらかじめ定められた口径選択機能により、同一時に駆動する振動子群を選択して口径を形成し、その口径を順次移動して被検体や物体の内部を超音波ビームで走査する。また、コンベクス走査型の超音波撮像装置においても同様である。
【０００７】
まず、一方向からだけの撮像の場合、反射率の角度依存性の効果によって、被検体内の反射部分の輪郭のうち超音波の進行方向に垂直な部分が強調されて画像化され、超音波の進行方向に平行な組織境界部分はほとんど画像化さず、また、高反射体の影の部分が画像化されにくいという欠点があった。これに対して、コンパウンドイメージングのように超音波の進行方向を複数持っている場合、反射体の境界面に垂直に超音波が入射する場合が増えるので、比較的輪郭がつながって見えやすく、陰になる部分も減るというように、この欠点が改善される。
【０００８】
二つの目の長所として、コントラスト分解能の劣化の原因となるスペックルの除去に大きな効果がある。これは、角度が十分に異なる方向から撮像された画像同士では、スペックルパターンが異なるため、本来の被検体内の反射部分を反映した部分は加算の分だけ信号強度が増加するのに比べ、スペックルパターンは加算の平方根に比例してしか信号強度が増加しないため、加算の数をＮとするとＮの平方根に比例してコントラスト分解能が向上する。
【０００９】
ところで、このようなコントラスト分解能の向上は、画像の合成をＲＦ信号を包絡線信号に変換して、インコヒーレント化してから加算するために実現出来ることである。このことは、逆に言えば、インコヒーレント化してから加算しているため、コヒーレント加算のように空間分解能の向上は期待出来ないということになる。
【００１０】
一方、コヒーレント加算という観点、すなわちＲＦ信号の段階で合成する方法として非特許文献１に開示されているような開口合成法も公知である。開口合成はコヒーレントな段階で複数の異なる方向からの信号を加算するので、実効的に開口が大きくなることと等価であり、空間分解能の向上に寄与する。
【００１１】
これらの、コンパウンド撮像法、開口合成法は、共に、被検体が生体のように、動く対象物である場合には、複数の撮像間の動きが問題になる。これに対して、前記の二つの公知例においては、画像の合成を行う前に画像同士の位置ずれを相関法によって推定し、被検体の動きによる画像のぼけを最小にする方法が開示されている。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特表２００２−５２６２２５号公報
【非特許文献】
【００１３】
【非特許文献１】IEEE Transaction on Ultrasonics, Ferroelectrics, and Frequency Control, ４５巻Ｎｏ．４，１０７７頁から１０８７頁
【発明の概要】
【００１４】
超音波撮像の対象物である人体の臓器は剛体ではないため、変形は局所毎に異なる。したがって、これを撮像範囲の全てに渡って理想的に補正することは困難である。また、開口合成とコンパウンド撮像では、合成する信号間の位置ずれの影響の大きさが異なる。
【００１５】
開口合成は、合成する信号同士間で干渉性があるため、波長の４分の１以上ずれると合成によって却って信号が劣化してしまう。一方、コンパウンド撮像の場合は、合成する信号間で干渉性は無いため、包絡線信号の包絡線の幅以内のずれであれば、信号の劣化にはならない。一般に現状の超音波撮像装置においては、包絡線の幅は波長の数倍程度である。
【００１６】
このように、現状の画像合成法は、コンパウンド撮像法と開口合成法の２種類の方法があり、それぞれ長所が異なるので、二つの長所を兼ね備えた撮像方法が望まれていた。すなわち、コントラスト分解能と空間分解能の向上を同時に実現すること、また被検体の動きによる画像劣化が少ない撮像法を実現することである。
【００１７】
現状の超音波撮像装置においては、超音波探触子３２の伝達関数が周波数空間での幅が制限されていることから、時間軸上での数波程度の幅を持つ。超音波撮像では、超音波探触子３２からの超音波の送信および被検体から反射される超音波の超音波探触子３２による受信の送受信で、超音波探触子３２の伝達関数の影響を二回受けるため、包絡線の幅は２から４波程度である。即ち、コヒーレントな場合とインコヒーレントな場合で、１０倍程度、動きに対する影響が異なることがわかる。もし合成する画像間のずれの量が、完全に理想的にピクセル毎に補正できるのであれば、前記公知の二つの方法のいずれかで、動きの影響を補正することが出来るが、現実には、ノイズの影響で補正値の精度が保障されないことや、断層像と直交する方向にもずれがあることから、局所的なずれを正確に推定することは難しい。
【００１８】
本発明においては、局所的なずれを完全に補正することは諦めて、大局的なずれ量から、ずれが小さい（例えば波長の１／４より小さい）場合、合成をＲＦデータ同士で行ない、ずれが大きい場合は、合成を包絡線検波信号で行なう。ずれの大きさは、異なるフレームの画像間でのＲＦデータの相互相関に着目して判断することができる。さらに、ヘテロダイン検波の参照周波数を制御することにより、ＲＦデータ同士での合成、または、包絡線検波信号での合成と切り分けることができる。
【図面の簡単な説明】
【００１９】
【図１】（Ａ）は超音波撮像装置の一般的な構成を示す図、（Ｂ）は、超音波探触子の一例を示す図である。
【図２】本発明の超音波撮像装置の実施例１の信号の処理フローを示す図である。
【図３】本発明の超音波撮像装置の実施例２の信号の処理フローを示す図である。
【図４】本発明の超音波撮像装置の実施例３の信号の処理フローを示す図である。
【図５】（Ａ）は送受信トランスジューサ３２が受信するビーム進行方向の異なる３方向Ａ，ＢおよびＣのビームの方向を模式的に示す図、（Ｂ）は送受信トランスジューサ３２が受信するＲＦ信号の各ビーム進行方向のビームパターンを示す図、（Ｃ）は送受信トランスジューサ３２が受信したＲＦ信号の包絡線検波後のビームパターンを示す図である。
【図６】（Ａ）は、横軸にフレームの画像間でのずれをとり、縦軸にそれらに対応して行なわれる包絡線検波およびＲＦデータでの加算の選択の例を示す図、（Ｂ）は、ヘテロダイン検波の参照周波数ω1の設定をするための設定曲線の例を示す図である。
【図７】図５（Ｂ）に示す、異なる３方向Ａ，ＢおよびＣのビームを包絡線検波して図５（Ｃ）に示す包絡線検波後のビームパターンとした後、合成したときの超音波信号強度の分布を示す図である。
【図８】（Ａ），（Ｂ），（Ｃ）および（Ｄ）は、超音波の中心周波数ω0を２ＭＨｚとし、順に、参照波の周波数ω1を２ＭＨｚ、１ＭＨｚ、０ＭＨｚ、−１ＭＨｚとしたときのヘテロダイン検波後の結果を示すビームパターンの例を示す図である。
【図９】（Ａ）は図２のステップ１１５,１１７を経て得られる画像のシミュレーション結果の例、（Ｂ）は図２のステップ１１９,１２０を経て得られる画像のシミュレーション結果の例を示す図である。
【図１０】（Ａ）−（Ｄ）は図８（Ａ）−（Ｄ）に対応する画像例を示す図である。
【図１１】本発明の超音波撮像装置の実施例３において使用者が操作する参照波の周波数を設定するための入力装置の例の外観図である。
【発明を実施するための形態】
【００２０】
以下、本発明の実施例について図面を用いて詳細に説明する。本発明の各実施例は図１（Ａ），（Ｂ）に示すような情報処理装置を持つ超音波撮像装置によって実現される。
【実施例１】
【００２１】
図２は本発明の超音波撮像装置の実施例１の信号の処理フローを示す図である。
【００２２】
１０１は準備段階のステップで、使用する超音波探触子のタイプ、撮像の対象部位を示すコード等を入力装置２１を介して入力する。１０３は撮像の対象部位の画像合成を行なう撮像フレ―ム数を指定するステップで、入力装置２１を介して入力する。画像合成を行なうフレーム数は少なくとも３以上とする。このフレーム数を７から１５程度まで増やしてもよく、フレームの枚数を増やす程、画像合成の効果および合成による欠点も増えてくるので、撮像対象（撮像対象の動き）を考慮して決めるのが良い。
【００２３】
１０５は撮像の対象部位に対する超音波の送波方向を設定するステップで、入力装置２１を介して入力する。これはステップ１０３で指定した画像合成を行なうフレーム数に対応し、撮像済みのフレーム数に応じて、プログラムにより自動的に設定されるものとしても良い。ここで、超音波の送波方向とは、図１（Ｂ）で言えば、焦点Ｆ1−Ｆ1の面上で異なるＡ，ＢおよびＣの方向に超音波を送る向きを変えることである。図５（Ａ）では、これを模式的に実線で示す四角形Ａ、一点鎖線で示す四角形Ｂおよび破線で示す四角形Ｃのように示した。すなわち四角形Ａ、ＢとおよびＣは、それぞれ、被検体に対し同一焦点深度、異なった角度で超音波ビームを送出、これをリニア走査している様子を模式的に示す。
【００２４】
１０７は１フレーム分の超音波送受波・ＲＦデータ格納のステップである。撮像フレ―ム数が指定され（ステップ１０３）、送波方向が設定される（ステップ１０５）とＣＰＵ１１により、撮像の対象部位に応じたプログラムに対応する超音波パルス信号が生成され、これが送波ビームフォーマ２３で超音波探触子３２の振動子に与える信号に変換されて、送受分離器３１を介して超音波探触子３２の振動子に与えられる。撮像の対象部位から反射された超音波信号が超音波探触子３２の振動子で受信されて電気信号に変換され、送受分離器３１を介して受波ビームフォーマ２４に与えられる。ここで、ディジタル信号に変換されて、チャネル間で整相加算され、ＲＦデータとしてメモリ１２に記憶される。１０９は撮像済みのフレーム数を評価するステップで、ステップ１０３で指定したフレーム数だけの超音波送受波・ＲＦデータの格納が終わったか否かを検査する。ステップ１０３で指定したフレーム数だけの格納が終わっていない場合は、ステップ１０５に戻り、次の送波方向での１フレーム分の超音波送受波・ＲＦデータの格納を行なう。ステップ１０８で撮像済みのフレーム数がステップ１０３で指定したフレーム数に達したと評価されたら、データの評価に移る。
【００２５】
１１１は異なるフレームの画像間でのＲＦデータ相互相関を計算するステップである。メモリ１２に格納されている各フレームの画像データを順次読み出して、相互相関を計算する。すなわち、各フレームの画像データから、各送波方向の撮像での特定の被検体内の反射部分に着目した形での移動量および変形量が計算される。
【００２６】
１１３は異なるフレームの画像間での変形が設定値より大きいか否かを評価するステップである。すなわち、ステップ１１１で計算された、異なるフレームの画像間でのＲＦデータ相互相関が設定値より大きいときは、メモリ１２に格納されている各フレームのＲＦデータを包絡線検波して（ステップ１１５）、この検波後のデータから、異なるフレームの画像間でのデータ合成を行なう（ステップ１１７）。ステップ１１１で計算された、異なるフレームの画像間でのＲＦデータ相互相関が設定値より小さいときは、メモリ１２に格納されている各フレームのＲＦデータから異なるフレームの画像間でのデータ合成を行ない（ステップ１２１）、この合成後のＲＦデータを包絡線検波する（ステップ１２１）。
【００２７】
１２３は画像表示ステップであり、ステップ１１７のデータ、あるいは、ステップ１２１のデータのいずれかから、所定のプログラムに従った画像を表示装置２２に表示する。１２５は画質の調整の要否の判定ステップであり、例えば、画像表示ステップ１２３で、画面の一部に「再度データ取得が必要ですか」と言う質問と”ＹＥＳ”、”ＮＯ”のスイッチの表示をする。そして、使用者が、表示装置２２に表示された画像に不満であり、”ＹＥＳ”を選択したときに、ステップ１０３に戻って、必要な設定から再度実行する。”ＮＯ”を選択したときには、撮像を終了する。再度、撮像を行なうときは、例えば、合成枚数を増加したり、送波方向の変化を細かくしたり等、超音波送受波・ＲＦデータ取得の条件を変更して行なう。超音波診断機では、特に、撮像停止ボタンを押さない限り、画像を撮り続け、画面も更新し続けている。そのため、再度、データ取得の必要性を問う画面表示を行わなくても、画像に応じて操作者が設定を変更したり、それを次の撮像に反映する方法もある。
【００２８】
先にも述べたように、ＲＦ信号のようにコヒーレントな信号の場合と、包絡線検波信号のようにインコヒーレントな信号の場合とで、合成する信号間での、撮像の対象部位の変位の影響が大きく異なる。すなわち、コヒーレントな場合は、波長の半分のずれがあると、信号は完全に打ち消し合い、１／４以下のずれで漸くその影響を無視しうるようになる。一方インコヒーレントな場合、打ち消しあうことは無いので、包絡線の幅以下のずれであれば、その影響は少ない。
【００２９】
実施例１は、この点に着目して、ステップ１１１で異なるフレームの画像間でのＲＦデータ相互相関を計算し、ステップ１１３で異なるフレームの画像間での変形が設定値より大きいか否かを評価する。異なるフレームの画像間での撮像の対象部位の変位が設定値より大きい場合には、フレームの画像間でのずれによる影響が小さい、インコヒーレントな信号として処理をすることとし、異なるフレームの画像間での撮像の対象部位の変位が設定値より小さい場合には、コヒーレントな信号として処理をする。
【００３０】
ここで、相互相関による変形の計算方法について詳しく説明を行なう。１フレームの深さ方向のデータ数をＭ、横方向のデータ数をＮとする。また、深さ方向は浅い方から深くなる方に向かって、座標が大きくなり、横方向には走査するラスターの順番に座標が大きくなるように座標系を考える。尚ここでは走査順序は、隣から順に走査することを想定する。フレーム番号ｋの中の、深さ方向ｉ番目、横方向にｊ番目の画素の輝度をIki,j、同位置のＲＦ信号の振幅をＡki,jとして以下説明を行なう。
【００３１】
相関をとる窓の大きさをｗと設定すると、画素ｉｊにおける、二つのフレーム間相関による変形ベクトルＶki,jは、式（１）、（２）のように算出される。相互相関Ｃｒ（ｋ，ｋ＋１）が最大となる要素番号ｉｍａｘ１と自己相関Ｃｒ（ｋ，ｋ）が最大となる要素番号ｉｍａｘ０の差、Ｖ＝ｉｍａｘ０−ｉｍａｘ１が画素ｉｊにおける、ｋフレームとｋ＋１フレームの間の深さ方向の変形である。
【００３２】
【数１】
2010042266.tif 000003
【００３３】
【数２】
2010042266.tif 000004
【００３４】
ここで記号＊は相互相関である。ここでは、ｗが奇数の場合を記述したが、偶数の場合も、ｉ−（ｗ一１）／２をｉ−（ｗ−２）／２に、ｉ＋（ｗ−１）／２をｉ＋ｗ／２のように置換すれば同様である。
【００３５】
ここで、本来、局所変形は２次元ベクトルであるところを、深さ方向のみの１次元の問題として扱っているのは、通常、超音波パルスの送受波による点応答は深さ方向に狭く、方位方向に広がっているため、深さ方法はＲＦ波のナイキスト周波数より細かくサンプリングしているのに比べ、方位方向には、粗くデータを取得しているからである。定量的に比べると、音速１５００ｍ／ｓ、深さ方向のサンプリング周波数５０ＭＨｚの場合は、深さ方向のサンプリング間隔は１５μｍになる。一方で、方位方向のサンプリング間隔は通常、素子ピッチ程度もしくは、その半分程度であるため、素子幅を中心周波数での波長の半分程度と想定すると、中心周波数が３ＭＨｚの場合、１２５μｍ程度であり、深さ方向のサンプリング間隔の１０倍程度である。そのため、対象物の動きに対して、フレームレートが十分に早ければ変形を深さ方向に関してのみ測定すれば良い。
【００３６】
相関を取る時に、窓を設定しているのは以下の二つの理由による。第一に、一度に１ラスタ−分のデータを全て使ってしまうと、生体のように剛体としては扱えない場合には、深さによって変形が異なる場合に不利である。第二に、本発明においては、角度の異なる画像同士で相関を取っているので、１ラスター全体で比較すると、深さ方向全体で場所を一致させることが出来ない。しかし窓幅を限れば、窓の範囲内では、ラスターをまたがらせないことが可能となる。例えばビーム方向が５度異なるフレーム間での変形を求める場合、ラスタ一間隔が１２５μｍ、深さ方向のサンプリング間隔が１５μｍの場合には、１２５μｍ／ｔａｎ（５度）／１５μｍ＝９５なので、窓幅が９５サンプリング点以下であれば、方位方向に関しても同じ場所のデータから変位を計算しているといえる。勿論、この場合、窓を移動するに従って、相関をとる相手のラスター番号が変化していくことになる。
【００３７】
これらの方法で画素ｉｊでの変位が求まった後は、ノイズの影響が小さい場合は、そのまま変位量として用いることも可能であるが、ノイズの影響が大きい場合は、空間的なローパスフイルタなどを通して、平均化した値に置き換えて、処理を行なうことも可能である。
【００３８】
フレームレートに比べ、対象物の動きが早い場合には、方位方向も含めた２次元相関によって２次元の変形を計算することや、一度、ラスター毎に深さ方向に積分し、別フレームの積分パターン同士の相関によって方位方向の変形を求め、既に説明した深さ方向の変形と組み合わせて２次元の変形を計算する方法も可能である。
【００３９】
図６（Ａ）は、横軸にフレームの画像間でのずれをとり、縦軸にそれらに対応して行なわれる包絡線検波およびＲＦデータでの加算の選択の例を示す図である。フレームの画像間でのずれが超音波の中心周波数ω0の波長の１／４をしきい値とする例である。画像間でのずれが超音波の中心周波数ω0の波長の１／４を越えているときは、包絡線検波による処理、すなわち、インコヒーレントな信号として処理をする。これに対して、画像間でのずれが超音波の中心周波数ω0の波長の１／４より小さいときは、ＲＦデータによる加算処理をした後で包絡線検波による処理、すなわち、コヒーレントな信号として処理をする。
【００４０】
このように、異なるフレームの画像間での撮像の対象部位のずれに応じて画像生成のための信号を選択することにより、コントラスト分解能と空間分解能の向上を同時に実現し、また被検体の動きによる画像劣化が少ない撮像法を実現する。
【実施例２】
【００４１】
図３は本発明の超音波撮像装置の実施例２の信号の処理フローを示す図である。図３において、図２と同じ参照符号を付したステップは同じ処理をするステップである。
【００４２】
実施例２においても、実施例１と同様に、使用する超音波探触子のタイプ、撮像の対象部位を設定し（ステップ１０１）、画像合成を行なう撮像フレ―ム数を指定し（ステップ１０３）、撮像の対象部位に対する超音波の送波方向を設定し（ステップ１０５）、１フレーム分の超音波送受波・ＲＦデータを格納し（ステップ１０７）、撮像済みのフレーム数を評価し（ステップ１０９）、ステップ１０３で指定したフレーム数に達するまで繰り返す。そして、ステップ１１１で異なるフレームの画像間でのＲＦデータ相互相関を計算する。
【００４３】
１１４はヘテロダイン検波の参照周波数ω1の設定をするステップであり、１１６は、参照周波数ω1の設定に応じたヘテロダイン検波の実行ステップである。１１８は、ステップ１１６の検波の結果を受けて、異なるフレームの画像間での検波後データの合成をするステップである。ステップ１２３および１２５の画像表示および画質調整に付いては、実施例１と同じである。
【００４４】
図６（Ｂ）は、ステップ１１４でヘテロダイン検波の参照周波数ω1の設定をするための設定曲線の例を示す図である。横軸に超音波の中心周波数ω0の波長に対するフレームの画像間でのずれをとり、縦軸にそれに対応して決定される参照周波数に応じたミキシング後の周波数の超音波の中心周波数ω0に対する比の例を示す図である。ここでも、実施例１と同様、フレームの画像間でのずれが超音波の中心周波数ω0の波長の１／４をしきい値とする。すなわち、画像間でのずれが超音波の中心周波数ω0の波長の１／４であるときは、ヘテロダイン検波の参照周波数を超音波の中心周波数ω0と同じ周波数とする。フレームの画像間でのずれが超音波の中心周波数ω0の波長の１／４と異なるときは、図に示すミキシング後の周波数になるように参照周波数を決定する。その結果、画像間でのずれが超音波の中心周波数ω0の波長の１／４であるときのみ、包絡線検波による処理、すなわち、インコヒーレントな信号として処理が行なわれる。一方、画像間でのずれが超音波の中心周波数ω0の波長の１／４と異なるときは、その程度に応じた参照周波数が選択されて、ＲＦデータによるヘテロダイン検波が行なわれた後に包絡線検波による処理、すなわち、コヒーレントな信号として処理が行なわれる。
【００４５】
実施例２では、ステップ１１１で計算された、異なるフレームの画像間でのＲＦデータ相互相関が所定のしきい値にあるときのみ、参照周波数ω1を超音波の中心周波数ω0と同じ値に設定してホモダイン検波とする。すなわち、超音波受波のＲＦデータをコヒーレントな信号として、包絡線検波による画像化をする。異なるフレームの画像間でのＲＦデータ相互相関がしきい値以外であるときは参照周波数ω1を超音波の中心周波数ω0とは異なる値に設定してヘテロダイン検波とする。
【００４６】
ここで、実施例２におけるコヒーレントな信号とインコヒーレントな信号の切り替えを、ヘテロダイン検波部によって実行する方法についてもう少し具体的に説明する。
【００４７】
図５（Ａ）は送受信トランスジューサ３２が受信するビーム進行方向の異なる３方向Ａ，ＢおよびＣのビームの方向を模式的に示す図、（Ｂ）は送受信トランスジューサ３２が受信するＲＦ信号の各ビーム進行方向のビームパターンを示す図、（Ｃ）は送受信トランスジューサ３２が受信したＲＦ信号の包絡線検波後のビームパターンを示す図である。図５（Ｂ）および（Ｃ）を対比して分かるように、包絡線検波後のビームパターンでは、強度信号のみが強く表れ、位相信号が見えにくいものとなっている。
【００４８】
図７は、図５（Ｂ）に示す、異なる３方向Ａ，ＢおよびＣのビームを包絡線検波して図５（Ｃ）に示す包絡線検波後のビームパターンとした後、合成したときの超音波信号強度の分布を示す図である。すなわち、横方向に超音波信号の広がりを、縦方向に焦点深度を示す。図７は、送受信トランスジューサ３２は深さ５０ｍｍの位置に焦点が合わせられており、この焦点位置で±２°の広がりがあることを示している。縦軸方向に表れている線は、超音波の強度の等高線である。
【００４９】
図５（Ｂ）に示すようなＲＦ信号は時間の関数Ｆ（ｔ）として、式（３）のように記述することが出来る。
【００５０】
【数３】
2010042266.tif 000005
【００５１】
ここで、ω0は超音波の中心周波数、Ａ（ｔ）は包絡線である。
【００５２】
このＦ（ｔ）の包絡線を検波する方法の一つにヘテロダイン検波の特異形としてのホモダイン検波法がある。
【００５３】
Ｆ（ｔ）に参照波ｓｉｎ（ω1ｔ）を乗じた信号をＧ（ｔ）とすると、式（４）となる。
【００５４】
【数４】
2010042266.tif 000006
【００５５】
ここで、参照波の周波数ω1を超音波の中心周波数ω0と等しくして、ホモダイン検波とした場合、式（５）となるので、
【００５６】
【数５】
2010042266.tif 000007
【００５７】
低域透過フィルタを通して、二倍高調波を除去すると、包絡線信号が得られる。
【００５８】
すなわち、超音波の中心周波数ω0と等しい周波数の参照波ｓｉｎ（ω0ｔ）で実行されるヘテロダイン検波の特異形であるホモダイン検波とすれば、実施例１に記述した包絡線検波によるインコヒーレント信号により画像合成ができることになる。すなわち、図６（Ｂ）に例示したように、特定の異なるフレームの画像間のずれ（図では超音波の中心周波数の波長の１／４）でのみＲＦ信号、すなわち、コヒーレント信号による画像合成とする。そして、その前後のずれに対しては、ＲＦデータ相互相関をパラメータとして、図に例示したように、ヘテロダイン検波の参照周波数ω1の設定を超音波の中心周波数ω0と異なったものにする。これにより、コヒーレントからインコヒーレントまで、連続的に変化させることが出来、インコヒーレント信号による画像合成およびコヒーレント信号による画像合成を総合化した形の画像合成が実現できる。
【００５９】
図８は、超音波の中心周波数ω0を２ＭＨｚとし、参照波の周波数ω1を２ＭＨｚ、１ＭＨｚ、０ＭＨｚ、−１ＭＨｚとした場合を例に、図７の超音波信号強度の分布がどのように変化するかを説明する図である。
【００６０】
図８（Ａ），（Ｂ），（Ｃ）および（Ｄ）は、順に、参照波の周波数ω1を２Ｈｚ、１ＭＨｚ、０ＭＨｚ、−１ＭＨｚとしたときのヘテロダイン検波後の結果を示すビームパターンの例である。ビームパターンの細い（広がりが小さい）ほど、方位方向の空間分解能が、よく、またサイドローブが強いほど、スペックルパターンが強く生成され、コントラスト分解能が劣化する。
【００６１】
図８（Ａ）は、参照波の周波数ω1が２ＭＨｚ、すなわち、参照波でミキシングし、ローパスフイルタを通した後の周波数（ヘテロダイン検波後の周波数）が０ＭＨｚ、であるので、上述した包絡線検波に対応し、実施例１のステップ１１５による包絡線検波に対応する。この例が、この中では、最もコントラスト分解能がすぐれ、かつ空間分解能が最も悪い。図８（Ｃ）は、参照波の周波数ω1が０ＭＨｚ、すなわち、ヘテロダイン検波後の周波数が２ＭＨｚ、であるので、受信したＲＦ信号を検波することに対応し、実施例１のステップ１１９によるＲＦデータの合成後に包絡線検波をすることに対応する。図８（Ｂ）は、参照波の周波数ω1が１ＭＨｚ、すなわち、ヘテロダイン検波後の周波数が０ＭＨｚ、であるので、上述した（Ａ）と（Ｃ）の中間の状態の結果となる。図８（Ｄ）は、参照波の周波数ω1が−１ＭＨｚ、すなわち、ヘテロダイン検波後の周波数が３ＭＨｚ、である（参照波が図７（Ｂ）の逆位相）ので、最も空間分解能が優れ、コントラスト分解能が悪い状態となる。
【００６２】
すなわち、ビームパターン形状からも、上記の、コヒーレントからインコヒーレントに連続的に調節可能であることが分かるが、図９および図１０はこれを合成後の画面の例として示すものである。
【００６３】
図９（Ａ）は図２のステップ１１５,１１７を経て得られる画像のシミュレーション結果の例、図９（Ｂ）は図２のステップ１１９,１２０を経て得られる画像のシミュレーション結果の例である。この例は、人体の一部に超音波を反射しない血管が存在する部分を５ｃｍ×５ｃｍの大きさに画像化したシミュレーション結果の例である。中央下部にある黒い部分が血管に対応し、周辺部が肉体に対応する。また左側と上辺部にある白い点は、点反射体の映像である。図９（Ａ）、（Ｂ）を対比して分かるように、ＲＦ信号の段階で合成する図９（Ｂ）の画像の方が空間分解能が良い。
【００６４】
図１０（Ａ）−（Ｄ）は図８（Ａ）−（Ｄ）に対応する画像例であり、図９（Ａ）、（Ｂ）と同様にして得られたシミュレーション結果の例である。図１０（Ａ）は、参照波の周波数ω1を２ＭＨｚとしており、上述した包絡線検波に対応する図９（Ａ）と同じである。図１０（Ｃ）は、参照波の周波数ω1を０ＭＨｚとしており、受信したＲＦ信号を検波することに対応し、実施例１のステップ１１９によるＲＦデータの合成後に包絡線検波をすることに対応する図９（Ｂ）と同じである。図１０（Ｂ）は、図８（Ｂ）と同様、参照波の周波数ω1が１ＭＨｚであるので、上述した図１０（Ａ）と（Ｃ）の中間の状態の結果となる。図１０（Ｄ）は、図８（Ｄ）と同様、参照波の周波数ω1が−１ＭＨｚである（参照波が図７（Ｂ）の逆位相）ので、最も空間分解能が優れ、コントラスト分解能が悪い状態となる。
【実施例３】
【００６５】
図４は本発明の超音波撮像装置の実施例３の信号の処理フローを示す図である。図４において、図２または図３と同じ参照符号を付したステップは同じ処理をするステップである。
【００６６】
実施例３は、実施例２がヘテロダイン検波の参照周波数を異なるフレームの画像間でのＲＦデータ相互相関に応じて自動的に選択される（ステップ１１４）ものとしたのに対して、ステップ１０４で、使用者が自己の考えに基づいてヘテロダイン検波の参照周波数を操作して変更することが出来るようにした点においてのみ異なる。その他の処理は図３に示す処理フローと同じである。
【００６７】
図１１は、実施例３においてステップ１０４で使用者が操作する入力装置の外観図の例である。５１は操作ボタンであり、移動ガイド５２に沿って左右に直線移動が可能とされる。５３は表示部である。ここで例示したのは、直線移動可能な、ボリュームタイプの操作ボタンである。例えば、表示０−３は、図１０（Ａ）−（Ｄ）に対応する画像が得られる参照周波数を指定する位置である。すなわち、操作ボタン５１の操作は、図６（Ｂ）に対応するヘテロダイン検波の参照周波数ω1の設定をすることを意味する。この場合、使用者は、対象物や、対象疾病に応じて、得られる画像がどのように変化するかを経験により積み重ねることができ、ケースに応じて、コントラスト分解能と空間分解能のどちらを優先させるべきかを、使用者の判断に応じて可変とすることが可能になる。
【００６８】
上記の説明では、合成する画像の枚数を３枚で説明したが、この枚数を大きくして、その一部により、コヒーレント画像の合成を行ない、残りでインコヒーレント画像の合成を行ない、これらの両方を混ぜて画像化することも可能となる。
【００６９】
本発明は、上記にある特定の実施の形態に限定されるものでなく、その技術思想の範囲を逸脱しない範囲で様々な変形が可能である。
【産業上の利用可能性】
【００７０】
本発明は、動きが少ない場合は合成によって、空間分解能が向上し、動きが大きい場合でもコントラスト分解能が向上する。別の言い方をすれば、合成する画像間の位置ずれに応じて、コントラスト分解能と空間分解能のバランスのとれた画像を得ることが出来る。更には、撮像対象に応じた、もしくは使用者の選択によって、コントラスト分解能と空間分解能の調節がとれた画像を得ることが出来る。
【符号の説明】
【００７１】
１０…情報処理装置
１１…ＣＰＵ
１２…メモリ
１３…ハードディスク
１４…バス
２１…入力装置
２２…表示装置
２３…送波ビームフォーマ
２４…受波ビームフォーマ
３１…送受分離器
３２…超音波探触子。
【図１】
2010042266.tif 000008
【図２】
2010042266.tif 000009
【図３】
2010042266.tif 000010
【図４】
2010042266.tif 000011
【図５】
2010042266.tif 000012
【図６】
2010042266.tif 000013
【図７】
2010042266.tif 000014
【図８】
2010042266.tif 000015
【図９】
2010042266.tif 000016
【図１０】
2010042266.tif 000017
【図１１】
2010042266.tif 000018
【手続補正書】
【提出日】平成21年10月7日(2009.10.7)
【手続補正１】
【補正対象書類名】図面
【補正対象項目名】図３
【補正方法】変更
【補正の内容】
【図３】
2010042266.tif 000019
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図４
【補正方法】変更
【補正の内容】
【図４】
2010042266.tif 000020