WO2005013829A1

WO2005013829A1 - 超音波撮像装置および超音波撮像方法

Info

Publication number: WO2005013829A1
Application number: PCT/JP2004/008391
Authority: WO
Inventors: Takashi Azuma; Shin-Ichiro Umemura; Hirotaka Baba
Original assignee: Hitachi Medical Corporation
Priority date: 2003-08-06
Filing date: 2004-06-09
Publication date: 2005-02-17
Also published as: EP1652475B1; JPWO2005013829A1; EP1652475A1; JP5357684B2; JP4426530B2; JP2010042266A; US20060173334A1; US20090118617A1; US8079956B2; CN1805710A; EP1652475A4; US7473226B2

Abstract

画像合成手段を備えた超音波診断装置において、コントラスト分解能の向上と空間分解能の向上を両立させるため、合成する画像間で相関をとり、画像間の変位量を計算。変位量が大きい時は包絡線検波後の信号を合成し、変位量が小さい時はRF信号を合成する。もしくは、変位量に応じてミキシング周波数を可変とし、位置ずれの程度に応じて、空間分解能の向上とコントラスト分可能向上のバランスを取る。

Description

明細書超音波撮像装置およぴ超音波撮像方法技術分野

本発明は、医療診断等に用いられる超音波撮像装置に係り、特に、いろいろな角度から走査して得られる画像信号を合成して 1つの画像として表示する、一般に、コンパウンド走査像と称される高画質の断層像得る技術に関する。背景技術

超音波撮像装置は、被検体の表面に超音波探触子を当て、その探触子から被検体に超音波を送信するとともに、被検体内部からの反射波（ェコー信号）を受信し、そのエコー信号に基づいて被検体の各部の状態を断層像により表示することにより、診断等に供するものである。

図 1 ( A ) は、超音波撮像装置の一般的な構成を示す図である。 1 0 はいわゆるパソコン等の情報処理装置である。 1 1は C P U、 1 2はメモリ、 1 3はハードディスク、 2 1は入力装置、 2 2は表示装置である。これらの装置はバス 1 4で連結され、必要な情報の伝達を行なう。ハードディスク 1 3には、超音波撮像装置の動作に必要な各種のプログラムおよびデータが格納されており、必要に応じて呼び出されてメモリ 1 2 に記憶され、 C P U 1 1の管理の下に、動作に利用される。入力装置 2 1は、キーボードあるいはマウス等であり、超音波撮像装置の使用者が必要な設定情報を入力し、あるいは、操作信号を与えるために使用される。表示装置 1 3には必要な操作画面および撮像結果の画面等が表示される。

2 3および 2 4は超音波の送波ビームフォーマおよび受波ビームフォ一マであり、情報処理装置 1 0のパス 1 4に接続される。さらに、送波ビームフォーマ 2 3および受波ビームフォーマ 2 4は、送受分離器 3 1 を介して、超音波の超音波探触子 3 2に接続される。使用者の操作信号に応じたプログラムに対応した超音波信号が超音波探触子 3 2から送出されるように、送波ビームフォーマ 2 3は情報処理装置 1 0から与えられるディジタル信号に応じた超音波送信信号を作って、送受分離器 3 1 を介して、超音波の超音波探触子 3 2に送る。一方、送信された超音波信号が被検体内部からの反射波として超音波探触子 3 2に受信されると、これに応じた信号が送受分離器 3 1 を介して、受波ビームフォーマ 2 4 に送られる。受波ビームフォーマ 2 4はこれにディジタル化等の必要な信号処理を施して情報処理装置 1 0に送る。情報処理装置 1 0では、受信信号を撮像の目的に対応したプログラムで処理して、結果を表示装置 1 3に表示する。

図 1 ( B ) は、超音波探触子 3 2の一例を示す図であり、 3分割された直線状に配置された振動子 1— 1， 1一 2および 1一 3が n組並列に配列された例である。図では表示を省略したが、これらの振動子には、送信制御回路 2 3からの送信信号を送り、振動子が検出する被検体内部からの反射波に応じた信号を受信制御回路 2 4に送るための配線が設けられる。一方、超音波探触子 3 2の他面には、振動子を被検体に効果的に接触させ、超音波を効率良く被検体との間で送受するための、例えば、金属粉末と髙分子材料の混合物からなる音響整合層 5 とやはり高分子材料からなる音響レンズが設けられている。各振動子に与える信号の大きさおよぴタイミングを制御することにより、超音波探触子 3 2から照射される超音波の焦点深度、照射の向きおよび強度を任意に制御することができ、図に示す F i— Fい F ₂— F ₂のように深さおょぴ向きを制御することができる。また、後述する複数の角度方向から撮像した前記被検体の Bモード画像データを合成するために、例えば、焦点の軸上でビーム方向が A， Bおよび Cの方向になるように制御することができる。なお、超音波探触子 3 2は、使用者が持っための把手も備えているが、図示は省略した。超音波探触子 3 2は、図 1 ( B ) に示す例に限られず、撮像目的に応じて、複数の振動子を等間隔で直線状、曲線状あるいは面状に配列して形成する等、種々の形態で作製される。例えば、リニア走査型の超音波撮像装置は、あらかじめ定められた口径選択機能により、同一時に駆動する振動子群を選択して口径を形成し、その口径を順次移動して被検体や物体の内部を超音波ビームで走査する。また、コンペタス走查型の超音波撮像装置においても同様である。

まず、一方向からだけの撮像の場合、反射率の角度依存性の効果によつて、被検体内の反射部分の輪郭のうち超音波の進行方向に垂直な部分が強調されて画像化され、超音波の進行方向に平行な組織境界部分はほとんど画像化さず、また、高反射体の影の部分が画像化されにくいという欠点があった。これに対して、コンパウンドイメージングのように超音波の進行方向を複数持っている場合、反射体の境界面に垂直に超音波が入射する場合が増えるので、比較的輪郭がつながって見えやすく、陰になる部分も減るというように、この欠点が改善される。

二つの目の長所として、コントラ'スト^解能の劣化の原因となるスぺックルの除去に大きな効果がある。これは、角度が十分に異なる方向から撮像された画像同士では、スペックルパターンが異なるため、本来の被検体内の反射部分を反映した部分は加算の分だけ信号強度が増加するのに比べ、スペックルパターンは加算の平方根に比例してしか信号強度が増加しないため、加算の数を Nとすると Nの平方根に比例してコントラスト分解能が向上する。

ところで、このようなコントラスト分解能の向上は、画像の合成を R F信号を包絡線信号に変換して、インコヒーレント化してから加算するために実現出来ることである。このことは、逆に言えば、インコヒーレント化してから加算しているため、コヒーレント加算のように空間分解能の向上は期待出来ないということになる。

一方、コヒーレント加算という観点、すなわち R F信号の段階で合成する方法として非特許文献 1に開示されているような開口合成法も公知である。開口合成はコヒーレントな段階で複数の異なる方向からの信号を加算するので、実効的に開口が大きくなることと等価であり、空間分解能の向上に寄与する。

これらの、コンパウンド撮像法、開口合成法は、共に、被検体が生体のように、動く対象物である場合には、複数の撮像間の動きが問題になる。これに対して .

、前記の二つの公知例においては、画像の合成を行う前に画像同士の位置ずれを相関法によって推定し、被検体の動きによる画像のぼけを最小にする方法が開示されている。

【特許文献 1】

特表 2 0 0 2— 5 2 6 2 2 5号公報

【非特許文献 1】

ΙΕϋΕ Transaction on Ultrasonics, Ferroelectrics, and Frequency Control, 4 5卷 N o . 4， 1 0 7 7頁力ら 1 0 8 7頁発明の開示

超音波撮像の対象物である人体の臓器は剛体ではないため、変形は局所毎に異なる。したがって、これを撮像範囲の全てに渡って理想的に補正することは困難である。また、開口合成とコンパウンド撮像では、合成する信号間の位置ずれの影響の大きさが異なる。

開口合成は、合成する信号同士間で干渉性があるため、波長の 4分の 1以上ずれると合成によって却って信号が劣化してしまう。一方、コンパウンド撮像の場合は、合成する信号間で干渉性は無いため、包絡線信号の包絡線の幅以内のずれであれば、信号の劣化にはならない。一般に現状の超音波撮像装置においては、包絡線の幅は波長の数倍程度である。

このように、現状の画像合成法は、コンパウンド撮像法と開口合成法の 2種類の方法があり、それぞれ長所が異なるので、二つの長所を兼ね備えた撮像方法が望まれていた。すなわち、コントラスト分解能と空間分解能の向上を同時に実現すること、また被検体の動きによる画像劣化が少ない撮像法を実現することである。

現状の超音波撮像装置においては、超音波探触子 3 2の伝達関数が周波数空間での幅'が制限されていることから、時間軸上での数波程度の幅を持つ。超音波撮像では、超音波探触子 3 2からの超音波の送信おょぴ被検体から反射される超音波の超音波探触子 3 2による受信の送受信で. 超音波探触子 3 2の伝達関数の影響を二回受けるため、包絡線の幅は 2 から 4波程度である。即ち、コヒーレントな場合とインコヒーレントな場合で、 1 0倍程度、動きに対する影響が異なることがわかる。もし合成する画像間のずれの量が、完全に理想的にピクセル毎に捕正できるのであれば、前記公知の二つの方法のいずれかで、動きの影響を補正することが出来るが、現実には、ノイズの影響で捕正値の精度が保障されないことや、断層像と直交する方向にもずれがあることから、局所的なずれを正確に推定することは難しい。

本発明においては、局所的なずれを完全に捕正することは諦めて、大局的なずれ量から、ずれが小さい（例えば波長の 1 / 4より小さい）場合、合成を R Fデータ同士で行ない、ずれが大きい場合は、合成を包絡線検波信号で行なう。ずれの大きさは、異なるフレームの画像間での R Fデータの相互相関に着目して判断することができる。さらに、ヘテロダイン検波の参照周波数を制御することにより、 R Fデータ同士での合成、または、包絡線検波信号での合成と切り分けることができる。図面の簡単な説明

図 1 ( A ) は超音波撮像装置の一般的な構成を示す図、（B ) は、超音波探触子の一例を示す図である。 .

図 2は本発明の超音波撮像装置の実施例 1の信号の処理フローを示す図である。図 3は本発明の超音波撮像装置の実施例 2の信号の処理フローを示す図である。

図 4は本発明の超音波撮像装置の実施例 3の信号の処理フローを示す図である。

図 5 (A) は送受信トランスジユーサ 3 2が受信するビーム進行方向の異なる 3方向 A， Bおよび Cのビームの方向を模式的に示す図、（B) は送受信トランスジユーサ 3 2が受信する R F信号の各ビーム進行方向のビームパターンを示す図、（C)は送受信トランスジユーサ 3 2が受信した R F信号の包絡線検波後のビームパターンを示す図である。

図 6 (A) は、横軸にフレームの画像間でのずれをとり、縦軸にそれらに対応して行なわれる包絡線検波および R Fデータでの加算の選択の例を示す図、（B) は、ヘテロダイン検波の参照周波数 _{ω ι}の設定をするための設定曲線の例を示す図である。

図 7は、図 5 (Β ) に示す、異なる 3方向 A， Bおよび Cのビームを包絡線検波して図 5 (C) に示す包絡線検波後のビームパターンとした後、合成したときの超音波信号強度の分布を示す図である。

図 8 ( A), (B), (C) および（D) は、超音波の中心周波数 ω。を 2 MH z とし、順に、参照波の周波数 _{ω ι}を 2ΜΗ ζ、 1 ΜΗ ζ、 OMH z、 ― 1 MH z としたときのへテロダイン検波後の結果を示すビームパターンの例を示す図である。

図 9 ( A) は図 2のステップ 1 1 5, 1 1 7を経て得られる画像のシミユレーシヨン結果の例、（B) は図 2のステップ 1 1 9， 1 2 0を経て得られる画像のシミユレーション結果の例を示す図である。

図 1 0 (A) — (D) は図 8 (A) 一 (D) に対応する画像例を示す図である。

図 1 1は本発明の超音波撮像装置の実施例 3において使用者が操作する参照波の周波数を設定するための入力装置の例の外観図である。発明を実施するための最良の形態

以下、本発明の実施例について図面を用いて詳細に説明する。本発明の各実施例は図 1 (A), (B) に示すような情報処理装置を持つ超音波撮像装置によって実現される。

〔実施の形態 1〕

図 2は本発明の超音波撮像装置の実施例 1の信号の処理フローを示す図である。

1 0 1は準備段階のステップで、使用する超音波探触子のタイプ、撮像の対象部位を示すコード等を入力装置 2 1 を介して入力する。 1 0 3 は撮像の対象部位の画像合成を行なう撮像フレーム数を指定するステツプで、入力装置 2 1 を介して入力する。画像合成を行なうフレーム数は少なくとも 3以上とする。このフレーム数を 7から 1 5程度まで増やしてもよく、フレームの枚数を増やす程、画像合成の効果および合成による欠点も増えてくるので、撮像対象（撮像対象の動き

) を考慮して決めるのが良い。

1 0 5は撮像の対象部位に対する超音波の送波方向を設定するステツプで、入力装置 2 1を介して入力する。これはステップ 1 0 3で指定した画像合成を行なうフレーム数に対応し、撮像済みのフレーム数に応じて、プログラムにより自動的に設定されるものとしても良い。ここで、超音波の送波方向とは、図 1 (B) で言えば、焦点 F_t— の面上で異なる A， Bおよび Cの方向に超音波を送る向きを変えることである。図 5 (A) では、これを模式的に実線で示す四角形 A、一点鎖線で示す四角形 Bおよび破線で示す四角形 Cのように示した。すなわち四角形 A、 Bとおよぴ Cは、それぞれ、被検体に対し同一焦点深度、異なった角度で超音波ビームを送出、これをリユア走査している様子を模式的に示す。

1 0 7は 1 フレーム分の超音波送受波 · R Fデータ格納のステップである。撮像フレーム数が指定され'（ステップ 1 0 3)、送波方向が設定される（ステップ 1 0 5 ) と C P U 1 1により、撮像の対象部位に応じたプログラムに対応する超音波パルス信号が生成され、これが送波ビームフォーマ 2 3で超音波採触子 3 2の振動子に与える信号に変換されて、送受分離器 3 1を介して超音波探触子 3 2の振動子に与えられる。撮像の対象部位から反射された超音波信号が超音波探触子 3 2の振動子で受信されて電気信号に変換され、送受分離器 3 1を介して受波ビームフォ一マ 2 4に与えられる。ここで、ディジタル信号に変換されて、チヤネル間で整相加算され、 R Fデータとしてメモリ 1 2に記憶される。 1 0 9は撮像済みのフレーム数を評価するステツプで、ステップ 1 0 3で指定したフレーム数だけの超音波送受波 · R Fデータの格納が終わったか否かを検査する。ステップ 1 0 3で指定したフレーム数だけの格納が終わっていない場合は、ステップ 1 0 5に戻り、次の送波方向での 1 フレーム分の超音波送受波 · R Fデータの格納を行なう。ステップ 1 0 8で撮像済みのフレーム数がステップ 1 0 3で指定したフレーム数に達したと評価されたら、データの評価に移る。

1 1 1は異なるフレームの画像間での R Fデータ相互相関を計算するステップである。メモリ 1 2に格納されている各フレームの画像データを順次読み出して、相互相関を計算する。すなわち、各フレームの画像データから、各送波方向の撮像での特定の被検体内の反射部分に着目した形での移動量および変形量が計算される。

1 1 3は異なるフレームの画像間での変形が設定値より大きいか否かを評価するステップである。すなわち、ステップ 1 1 1で計算された、異なるフレームの画像間での R Fデータ相互相関が設定値より大きいときは、.メモリ 1 2に格納されている各フレームの R Fデータを包絡線検波して（ステップ 1 1 5 )、この検波後のデータから、異なるフレームの画像間でのデータ合成を行なう（ステップ 1 1 7 )。ステップ 1 1 1 で計算された、異なるフレームの画像間での R Fデータ相互相関が設定値より小さいときは、メモリ 1 2に格納されている各フレームの R Fデータから異なるフレームの画像間でのデータ合成を行ない（ステップ 1 2 1 ) . この合成後の R Fデータ'を包絡線検波する（ステップ 1 2 1 )。

1 2 3は画像表示ステツプであり、ステップ 1 1 7のデータ、あるいは、ステップ 1 2 1 のデータのいずれかから、所定のプログラムに従つた面像を表示装置 2 2に表示する。 1 2 5は画質の調整の要否の判定ステツプであり、例えば、画像表示ステップ 1 2 3で、画面の一部に「再度データ取得が必要ですか」と言う質問と" Y E S "、" N O " のスイツチの表示をする。そして、使用者が、表示装置 2 2に表示された画像に不満であり、" Y E S " を選択したときに、ステップ 1 0 3に戻って、必要な設定から再度実行する。" N O " を選択したときには、撮像を終了する。再度、撮像を行なうときは、例えば、合成枚数を増加したり、送波方向の変化を細かくしたり等、超音波送受波 · R Fデータ取得の条件を変更して行なう。超音波診断機では、特に、撮像停止ボタンを押さない限り、画像を撮り続け、画面も更新し続けている。そのため、再度、データ取得の必要性を問う画面表示を行わなくても、画像に応じて操作者が設定を変更したり、それを次の撮像に反映する方法もある。

先にも述べたように、 R F信号のようにコヒーレントな信号の場合と、包絡線検波信号のようにィンコヒーレントな信号の場合とで、合成する信号間での、撮像の対象部位の変位の影響が大きく異なる。すなわち、コヒーレントな場合は、波.長の半分のずれがあると、信号は完全に打ち消し合い、 1 / 4以下のずれで漸くその影響を無視しうるようになる。一方インコヒーレントな場合、打ち消しあうことは無いので、包絡線の幅以下のずれであれば、その影響は少ない。

実施例 1 は、この点に着目して、ステップ 1 1 1で異なるフレームの画像間での R Fデータ相互相関を計算し、ステップ 1 1 3で異なるフレームの画像間での変形が設定値より大きいか否かを評.価する。異なるフレームの画像間での撮像の対象部位の変位が設定値より大きい場合には. フレームの画像間でのずれによる影響が小さい、インコヒーレントな信号として処理をすることとし、異なるフレームの画像間での撮像の対象部位の変位が設定値より小さい場合には、コヒーレントな信号として処理をする。

ここで、相互相関による変形の計算方法について詳しく説明を行なう。 1 フレームの深さ方向のデータ数を M、横方向のデータ数を Nとする。また、深さ方向は浅い方から深くなる方に向かって、座標が大きくなり、横方向には走查するラスターの順番に座標が大きくなるように座標系を考える。尚ここでは走査順序は、隣から順に走查することを想定する。フレーム番号 kの中の、深さ方向 i番目、横方向に； j 番目の画素の輝度を、同位置の R F信号の振幅を Α 」として以下説明を行なう。

相関をとる窓の大きさを wと設定すると、画素 i j における、二つのフレーム間相関による変形べクトルは、式（ 1 )、（ 2 ) のように算出される。相互相関 C r (k， k + l )、

C_r(k,k + l) = A ---(i)

2

が最大となる要素番号 i m a X と自己相関 C r ( k， k)

が最大となる要素番号 i m a X 0の差、 V= i m a x O— i m a x lが画素 i j における、 kフレームと k + lフレームの間の深さ方向の変形である。ここで記号 *は相互相関である。ここでは、 wが奇数の場合を記述したが、偶数の場合も、 i 一 (w— 1 ) / 2を i — (w— 2) / 2 に、 i + (w- 1 ) ノ 2を i + 2のように置換すれば同様である。ここで、本来、局所変形は 2次元べクトルであるところを、深さ方向のみの 1次元の問題として扱っているのは、通常、超音波パルスの送受波による点応答は深さ方向に狭く、方位方向に広がっているため、深さ方法は R F波のナイキスト周波数より細かくサンプリングしているのに比べ、方位方向には、粗くデータを取得しているからである。定量的に比べると、音速 1 5 0 0 m / s、深さ方向のサンプリング周波数 5 0 M H zの場合は、深さ方向のサンプリング間隔は 1 5 μ mになる。一方で、方位方向のサンプリング間隔は通常、素子ピッチ程度もしくは、その半分程度であるため、素子幅を中心周波数での波長の半分程度と想定すると、中心周波数が 3 M H zの場合、 1 2 5 /Z m程度であり、深さ方向のサンプリング間隔の 1 0倍程度である。そのため、対象物の動きに対して、フレームレートが十分に早ければ変形を深さ方向に関してのみ測定すれば良い。

相関を取る時に、窓を設定しているのは以下の二つの理由による。第一に、一度に 1ラスタ一分のデータを全て使ってしまうと、生体のように剛体としては扱えない場合には、深さによって変形が異なる場合に不利である。第二に、本発明においては、角度の異なる画像同士で相関を取っているので、 1ラスター全体で比較すると、深さ方向全体で場所を一致させることが出来ない。しかし窓幅を限れば、窓の範囲内では、ラスターをまたがらせないことが可能となる。例えばビーム方向が 5度異なるフレーム間での変形を求める場合、ラスター間隔が 1 2 5 m、深さ方向のサンプリング間隔が 1 5 z mの場合には、 1 2 5 μ / t a n ( 5度） / 1 5 /i m = 9 5なので、窓幅が 9 5サンプリング点以下であれば、方位方向に関しても同じ場所のデータから変位を計算しているといえる。勿論、この場合、窓を移動するに従って、相関をとる相手のラスター番号が変化していくことになる。

これらの方法で画素 i j での変位が求まった後は、ノイズの影響が小さい場合は、そのまま変位量として用いることも可能であるが、ノイズの影響が大きい場合は、空間的なローパスフィルタなどを通して、平均ィ匕した値に置き換えて、処理を行なうことも可能である。

フレームレートに比べ、対象物の動きが早い場合には、方位方向も含めた 2次元相関によって 2次元の変形を計算することや、一度、ラスタ一毎に深さ方向に積分し、別フレームの積分パターン同士の相関によつて方位方向の変形を求め、既に説明した深さ方向の変形と組み合わせて 2次元の変形を計算する方法も可能である。

図 6 ( A ) は、横軸にフレームの画像間でのずれをとり、縦軸にそれらに対応して行なわれる包絡線検波おょぴ R Fデータでの加算の選択の例を示す図である。フレームの画像間でのずれが超音波の中心周波数 ω。の波長の 1 / 4をしきい値とする例である。画像間でのずれが超音波の中心波数 ω。の波長の 1 4を越えているときは、包絡線検波による処理、すなわち、インコヒーレントな信号として処理をする。これに対して、画像間でのずれが超音波の中心周波数 ω。の波長の 1 Ζ 4より小さいときは、 R Fデータによる加算処理をした後で包絡線検波による処理、すなわち、コヒーレントな信号として処理をする。

このように、異なるフレームの画像間での撮像の対象部位のずれに応じて画像生成のための信号を選択することにより、コントラスト分解能と空間分解能の向上を同時に実現し、また被検体の動きによる画像劣化が少ない撮像法を実現する。

(実施例 2 )

図 3は本発明の超音波撮像装置の実施例 2の信号の処理フローを示す図である。図 3において、図 2 と同じ参照符号を付したステップは同じ処理をするステツプである。

実施例 2においても、実施例 1 と同様に、使用する超音波探触子のタイブ、撮像の対象部位を設定し（ステップ 1 0 1 )、画像合成を行なう撮像フレーム数を指定し (ステップ 1 0 3 )、撮像の対象部位に対する超音波の送波方向を設定し（ステップ 1 0 5 )、 1 フレーム分の超音波送受波 · R Fデータを格納し（ステップ 1 0 7 )、撮像済みのフレーム数を評価し（ステップ 1 0 9 )、ステップ 1 0 3で指定したフレーム数に達するまで繰り返す。そして、ステップ 1 1 1で異なるフレームの画像間での R Fデータ相互相関を計算する。

1 1 4はへテロダイン検波の参照周波数 _{ω ι}の設定をするステップであり、 1 1 6は、参照周波数 _{ω ι}の設定に応じたヘテロダイン検波の実行ステップである。 1 1 8は、ステップ 1 1 6の検波の結果を受けて、異なるフレームの画像間での検波後データの合成をするステツプである。ステップ 1 2 3および 1 2 5の画像表示および画質調整に付いては、実施例 1 と同じである。

図 6 ( Β ) は、ステップ 1 1 4でへテロダイン検波の参照周波数 _{ω ι} の設定をするための設定曲線の例を示す図である。横軸に超音波の中心周波数 ω。の波長に対するフレームの画像間でのずれをとり、縦軸にそれに対応して決定される参照周波数に応じたミキシング後の周波数の超音波の中心周波数 ω。に対する比の例を示す図である。ここでも、実施例 1 と同様、フレームの画像間でのずれが超音波の中心周波数 ω。の波長の 1 Ζ 4をしきい値とする。すなわち、画像間でのずれが超音波の中心周波数 ω。の波長の 1 / 4であるときは、ヘテロダイン検波の参照周波数を超音波の中心周波数 ω。と同じ周波数とする。フレームの画像間でのずれが超音波の中心周波数 ω。の波長の 1 Ζ 4 と異なるときは、図に示すミキシング後の周波数になるように参照周波数を決定する。その結果、画像間でのずれが超音波の中心周波数 ω。の波長の 1 / 4であるときのみ、包絡線検波による処理、すなわち、インコヒーレントな信号として処理が行なわれる。一方、画像間でのずれが超音波の中心周波数 ω。の波長の 1 / 4 と異なるときは、その程度に応じた参照周波数が選択されて、 R Fデータによるへテロダイン検波が行なわれた後に包絡線検波による処理、すなわち、コヒーレントな信号として処理が行なわれる。 '

実施例 2では、ステップ 1 1 1で計算された、異なるフレームの画像間での R Fデータ相互相関が所定のしきい値にあるときのみ、参照周波数 _{ω ι}を超音波の中心周波数 ω。と同じ値に設定してホモダイン検波とする。すなわち、超音波受波の R Fデータをコヒーレントな信号として、包絡線検波による画像化をする。異なるフレームの画像間での R Fデータ相互相関がしきい値以外であるときは参照周波数 ω ,を超音波の中心周波数 ω。とは異なる値に設定してヘテロダイン検波とする。

ここで、実施例 2 におけるコヒーレントな信号とインコヒーレントな信号の切り替えを、へテ口ダイン検波部によつて実行する方法についてもう少し具体的に説明する。

図 5 (Α) は送受信トランスジユーサ 3 2が受信するビーム進行方向の異なる 3方向 A， Bおよび Cのビームの方向を模式的に示す図、（B) は送受信トランスジユーサ 3 2が受信する R F信号の各ビーム進行方向のビームパターンを示す図、（C)は送受信トランスジユーサ 3 2が受信した R F信号の包絡線検波後のビームパターンを示す図である。図 5 (B) および（C) を対比して分かるように、包絡線検波後のビームパターンでは、強度信号のみが強く表れ、位相信号が見えにくいものとなつている。

図 7は、図 5 ( B ) に示す、. 異なる 3方向 A， Bおよび Cのビームを包絡線検波して図 5 (C) に示す包絡線検波後のビームパターンとした後、合成したときの超音波信号強度の分布を示す図である。すなわち、横方向に超音波信号の広がりを、縦方向に焦点深度を示す。図 7は、送受信トランスジユーサ 3 2は深さ 5 0 mmの位置に焦点が合わせられており、この焦点位置で ± 2 ° の広がりがあることを示している。縦軸方向に表れている線は、超音波の強度の等高線である。

図 5 (B) に示すような R F信号は時間の関数 F ( t ) として、式（ 3 ) のように記述することが出来る。

F(t)=A(t)sm(Gi_Qt) (3)

ここで、 ω。は超音波の中心周波数、 A ( t ) は包絡線である。

この F ( t ) の包絡線を検波する方法の一つにヘテロダイン検波の特異形としてのホモダイン検波法がある。

F ( t ) に参照波 s i η ( _{ω ι} t ) を乗じた信号を G ( t ) とすると、式（4 ) となる。 .

：、こ: P _v 照の周^ C ω ,超音波 φ，Φ心周波数 ω。と等しくして、 G{t) = -Α{ ιη{{ω₀ - ω₁ )tj+ smt(ro₀ + }J (4)

G(i) = iA(i)2sin(2Q₀i) (5)

ムホモダイン検波とした場合、式（ 5 ) となるので、

低域透過フィルタを通して、二倍高調波を除去すると、包絡線信号が得られる。

すなわち、超音波の中心周波数 ω。と等しい周波数の参照波 s i n (ω 。 t )で実行されるへテロダイン検波の特異形であるホモダイン検波とすれば、実施例 1に記述した包絡線検波によるィンコヒーレント信号により画像合成ができることになる。すなわち、図 6 (B) に例示したように、特定の異なるフレームの画像間のずれ（図では超音波の中心周波数の波長の 1ノ4) でのみ RF信号、すなわち、コヒーレント信号による画像合成とする。そして、その前後のずれに対しては、 R Fデータ相互相関をパラメータとして、図に例示したように、ヘテロダイン検波の参照周波数 _{ω ι}の設定を超音波の中心周波数 ω。と異なったものにする。これにより、コヒーレントからインコヒーレントまで、連続的に変化させることが出来、インコヒーレント信号による画像合成おょぴコヒーレント信号による画像合成を総合化した形の画像合成が実現できる。

図 8は、超音波の中心周波数 ω。を 2 ΜΗ ζ とし、参照波の周波数 _{W l} を 2ΜΗ ζ、 1 ΜΗ ζ、 OMH z、一 1 M H z とした場合を例に、図 7 の超音波信号強度の分布がどのように変化するかを説明する図である。図 8 (A), (B ), (C) および（D) は、順に、参照波の周波数 _{ω ι} を 2 Η ζ、 1 ΜΗ ζ、 0 Μ Η ζ、一 1 ΜΗ ζ としたときのへテロダイン検波後の結果を示すビームパターンの例である。ビームパターンの細い (広がりが小さい）ほど、方位方向の空間分解能が、よく、またサイドローブが強いほど、スペックルパターンが強く生成され、コントラスト分解能が劣化する。

図 8 ( Α) は、参照波の周波数 _{ω ι}が 2 ΜΗ ζ、すなわち、参照波でミキシングし、ローパスフィルタを通した後の周波数（ヘテロダイン検波後の周波数）が 0ΜΗ ζ、であるので、上述した包絡線検波に対応し、実施例 1のステップ 1 1 5による包絡線検波に対応する。この例が、この中では、最もコントラスト分解能がすぐれ、かつ空間分解能が最も悪い。図 8 ( C) は、参照波の周波数 ω iが 0 MH z、すなわち、ヘテロダイン検波後の周波数が 2 MH z、であるので、受信した R F信号を検波することに対応し、実施例 1のステップ 1 1 9による R Fデータの合成後に包絡線検波をすることに対応する。図 8 ( B ) は、参照波の周波数 ω が 1 MH z、すなわち、へテ口ダイン検波後の周波数が 0 MH z、であるので、上述した（A) と（C ) の中間の状態の結果となる。図 8 (D) は、参照波の周波数 ω tがー 1 MH z、すなわち、ヘテロダイン検波後の周波数が 3 MH z、である（参照波が図 7 ( B ) の逆位相）ので、最も空間分解能が優れ、コントラスト分解能が悪い状態となる。

すなわち、ビームパターン形状からも、上記の、コヒーレントからィンコヒーレントに連続的に調節可能であることが分かるが、図 9および図 1 0はこれを合成後の画面の例として示すものである。

図 9 ( A) は図 2のステップ 1 1 5 , 1 1 7を経て得られる画像のシミユレーション結果の例、図 9 ( B ) は図 2のステップ 1 1 9， 1 2 0を経て得られる画像のシミュレーション結果の例である。この例は、人体の一部に超音波を反射しない血管が存在する部分を 5 c m X 5 c mの大きさに画像化したシミュレーション結果の例である。中央下部にある黒い部分が血管に対応し、周辺部が肉体に対応する。また左側と上辺部にある白い点は、点反射体の映像である。図 9 (A)、 (B ) を対比して分かるように、 R F信号の段階で合成する図 9 ( B ) の画像の方が空間分解能が良い。

図 1 0 (A) — (D) は図 8 (A) 一 (D) に対応する画像例であり、図 9 (A)、 ( B ) と同様にして得られたシミュレーション結果の例である。図 1 0 (A) は、参照波の周波数 _{ω ι}を 2 MH z としており、上述した包絡線検波に対応する図 9 (A) と同じである。図 1 0 ( C ) は、参照波の周波数 ω iを 0 MH z としており、受信した R F信号を検波することに対応し、実施例 1のステップ 1 1 9による R Fデータの合成後に包絡線検波をすることに対応する図 9 (B) と同じである。図 1 0 (B) は、図 8 (B) と同様、参照波の周波数 _{ω ι}が 1 MH zであるので、上述した図 1 0 (Α) と（C) の中間の状態の結果となる。図 1 0 (D) は、図 8 (D) と同様、参照波の周波数 _{ω ι}が— 1 ΜΗ ζである（参照波が図 7 (Β) の逆位相）ので、最も空間分解能が優れ、コントラスト分解能が悪い状態となる。

(実施例 3 )

図 4は本発明の超音波撮像装置の実施例 3の信号の処理フローを示す図である。図 4において、図 2または図 3 と同じ参照符号を付したステップは同じ処理をするステップである。

実施例 3は、実施例 2がへテ口ダイン検波の参照周波数を異なるフレームの画像間での R Fデータ相互相関に応じて自動的に選択される（ステツプ 1 1 4) ものとしたのに対して、ステップ 1 0 4で、使用者が自己の考えに基づいてヘテロダイン検波の参照周波数を操作して変更することが出来るようにした点においてのみ異なる。その他の処理は図 3に示す処理フローと同じである.。

図 1 1は、実施例 3においてステツプ 1 0.4で使用者が操作する入力装置の外観図の例である。 5 1は操作ボタンであり、移動ガイド 5 2に沿.つて左右に直線移動が可能とされる。 5 3は表示部である。ここで例示したのは、直線移動可能な、ボリュームタイプの操作ボタンである。例えば、表示 0— 3は、図 1 0 (Α) — (D) に対応する画像が得られる参照周波数を指定する位置である。すなわち、操作ボタン 5 1の操作は、図 6 (Β ) に対応するへテロダイン検波の参照周波数 ω iの設定をすることを意味する。この場合、使用者は、対象物や、対象疾病に応じて、得られる画像がどのように変化するかを経験により積み重ねることができ、ケースに応じて、コントラスト分解能と空間分解能のどちらを優先させるべきかを、使用者の判断に応じて可変とすることが可能になる。上記の説明では、合成する画像の枚数を 3枚で説明したが、この枚数を大きくして、その一部により、コヒーレント画像の合成を行ない、残りでインコヒーレント画像の合成を行ない、これらの両方を混ぜて画像化することも可能となる。

本発明は、上記にある特定の実施の形態に限定されるものでなく、その技術思想の範囲を逸脱しない範囲で様々な変形が可能である。

産業上の利用可能性

本発明は、動きが少ない場合は合成によって、空間分解能が向上し、動きが大きい場合でもコントラスト分解能が向上する。別の言い方をすれば、合成する画像間の位置ずれに応じて、コントラスト分解能と空間分解能のパランスのとれた画像を得ることが出来る。更には、撮像対象に応じた、もしくは使用者の選択によって、コントラスト分解能と空間分解能の調節がとれた画像を得ることが出来る。

Claims

請求の範囲

1 . 被検体に与えるべき送波超音波ビームを形成するための超音波パルス信号を送出し、前記被検体から反射された超音波パルス信号から形成された受波超音波ビームを受信する情報処理装置と、前記超音波パルス信号に応じて送波超音波ビームを形成するための送波超音波ビームフォーマと、該送波超音波ビームフォーマの出力を送受分離器を介して与えられて前記被検体に送波超音波ビームを送出する超音波探触子と , 該超音波探触子に受信される前記被検体から反射された超音波パルス信号を送受分離器を介して与えられて前記受波超音波ビームを形成する受波ビームフォーマとよりなり、前記情報処理装置は被検体に対して複数の角度方向から撮像した前記被検体の Bモード画像データを合成し、合成画像を表示するためのプログラムを備えるとともに、前記異なる角度 • での各フレームの前記受波超音波ビームの間で相関を計算して、前記相関の結果に応じて前記受波超音波ビームを R F信号を用いて前記合成後に包絡線検波を行なうか、包絡線検波後の絶対値信号を用いて前記合成を行なうか、を決定することを特徴とする超音波撮像装置。

2 . 請求項 1に記載の超音波撮像装置において、前記 R F信号を用いた合成後の包絡線検波か、包絡線検波後の絶対値信号を用いた合成かは前記相関の結果として得られる前記被検体の変位量が、前記超音波パルスの送信及び受信の条件により定まる空間分解能よりも小さい場合、 R F信号を用いた合成を選択し、前記変位量が、前記空間分解能よりも大きい場合、包絡線検波後の絶対値信号を用いた合成を選択することを特徴とする超音波撮像装置。

3 . 被検体に与えるべき送波超音波ビームを形成するための超音波パルス信号を送出し、前記被検体から反射された超音波パルス信号から形成された受波超音波ビームを受信する情報処理装置と、前記超音波パルス信号に応じて送波超音波ビームを形成するための送波超音波ビームフォーマと、該送波超音波ビームフォーマの出力を送受分離器を介して与えられて前記被検体に送波超音波ビームを送出する超音波探触子と . 該超音波探触子に受信される前記被検体から反射された超音波パルス信号を送受分離器を介して与えられて前記受波超音波ビームを形成する受波ビームフォーマとよりなり、前記情報処理装置は被検体に対して複数の角度方向から撮像した前記被検体の Bモード画像データを合成し、合成画像を表示するためのプログラムを備えるとともに、前記異なる角度での各フレームの前記受波超音波ビームの間で相関を計算して、前記相関の結果に応じて前記受波超音波ビームのヘテロダイン検波の参照周波数を選択する超音波撮像装置であって、前記相関が所定の値であるとき、前記参照周波数を前記超音波パルスの中心周波数と等しくすることを特徴とする超音波撮像装置。

4 . 被検体に与えるべき送波超音波ビームを形成するための超音波パルス信号を送出し、前記被検体から反射された超音波ルス信号から形成された受波超音波ビームを受信する情報処理装置と、前記超音波パルス信号に応じて送波超音波ビームを形成するための送波超音波ビームフォーマと、該送波超音波ビームフォーマの出力を送受分離器を介して与えられて前記被検体に送波超音波ビームを送出する超音波探触子と、該超音波探触子に受信される前記被検体から反射された超音波パルス信号を送受分離器を介して与えられて前記受波超音波ビームを形成する受波ビームフォーマとよりなり、前記情報処理装置は被検体に対して複数の角度方向から撮像した前記被検体の Bモード画像データを合成し、合成画像を表示するためのプログラムを備えるとともに、前記異なる角度での各フレームの前記受波超音波ビームを指定された参照周波数でへテ口ダイン検波を行なう超音波撮像装置であって、前記参照周波数の指定が予め操作者が行なう操作により設定されることを特徴とする超音波撮

5 . 超音波探触子から超音波パルスを被検体に送信すること、前記被検体から反射される超音波パルスを受信すること、受信した超音波パルスから複数の角度方向から撮像した前記被検体の Bモード画像データを R F信号を用いて包絡線検波を行なうか、包絡線検波後の絶対値信号を用いて合成を行なうかのいずれかを選択して合成すること、合成画像を表示すること、よりなる超音波撮像方法であって、前記複数の角度方向から撮像した前記被検体の Bモード画像データの各フレームの相互相間を計算して前記被検体内の反射点の相対的な変位量を計算するとともに、該相対的な変位量に応じて前記選択を行なうことを特徴とする超音波撮像方法。

6 . 請求項 5に記載の超音波撮像方法において、前記相互相間から得られる変位量が、前記超音波パルスの送受及び受信の条件により定まる空間分解能よりも小さい場合、前記 R F信号を用いて前記合成後に包絡線検波を行ない、前記相互相間から得られる変位量が、前記空間分解能よりも大きい場合、前記包絡線検波後の絶対値信号を用いて前記合成を行なうことを特徴とする超音波撮像方法。

7 . 請求項 5に記載の超音波撮像方法において、前記相互相関の計算結果に応じて前記受波超音波ビームのへテ口ダイン検波の参照周波数を選択する超音波撮像方法であって、前記相互相関が所定の値であるとき、前記参照周波数を前記超音波パルスの中心周波数と等しくすることを特徴とする超音波撮像方法。

8 . 超音波探触子から超音波パルスを被検体に送信すること、前記被検体から反射される超音波パルスを受信すること、受信した受波超音波ビームを指定された参照周波数でヘテロダイン検波を行なうこと、超音波パルスから複数の角度方向から撮像した前記被検体の Bモード画像データを合成すること、合成画像を表示することよりなる超音波撮像方法であって、前記参照周波数が撮像の最初の段階で操作者が行なう操作により設定されることを特徴とする超音波撮像方法。