JPWO2011089898A1 - 変位推定方法、変位推定装置 - Google Patents

Abstract

超音波信号を用いて、変位を反復的に推定する変位推定方法であって、少なくとも１つの超音波信号を走査して、媒体に送信するステップと、走査された前記媒体から反射した前記超音波信号を受信するステップと、ウインドウサイズを計算するステップと、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、計算された前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む変位推定方法が開示される。

Description

本発明は、超音波信号を用いて変位を推定する方法を導入するものである。これは、最終的な結果、または、更なる処理のための中間ステップの何れであっても、変位推定が必要な適用例において、用いることができ、医療用および工業用超音波機器に用いることができる。

超音波装置は、高周波の機械的音波を、走査対象の媒体に対して送受信することにより、媒体を非侵襲で診断するものである。

そのような装置のトランスデューサは、超音波を、走査対象の媒体に送信する。

超音波は、散乱および反射することにより、その媒体の下層構造と相互作用する。

なお、下層構造は、例えば、媒体の表面よりも、媒体の内部側における構造などをいう。下層構造は、例えば、人間の身体における血管などである。

散乱および反射された超音波は、下層構造の有益な情報を含み、トランスデューサによって受信される。その後、超音波装置によって処理され、例えば画像情報としてユーザに提示される。

図８は、超音波ＲＦ信号を示す図である。

超音波装置が反射波から得る最も基礎的なデータの種類の１つとして、ＲＦ（Radiofrequency Signal）信号が挙げられる（図８）。

これは、受信波を、アナログからデジタルに直接変換したものである。

異なる適用例では、超音波ＲＦ信号から、輝度モード（Ｂモード）画像、ドップラー画像等、他の種類のデータを導出できる。

多くの適用例のうちの１つは、走査された媒体の下層構造の動きを分析することである。

ドップラー効果は、構造の動きの方向および強度を推測する簡易な方法として利用される。

しかし、ドップラー効果の精度は、非常に限られたものである。

医療用超音波における、血流のモニタリングのような、精度が厳密に要求されない適用例には、ドップラー超音波は適している。

しかし、構造が小さく、かつ動きが極めて小さい血管などが対象物であり、高水準の精度が必要とされる場合には、より感度の高い技術が必要である。

最近の新しい適用例として、構造の弾性を推測するために、構造の変位を用いる超音波弾性率計測法がある。

弾性を正確に推定するためには、構造の変位を、受信した超音波信号から、正確に推定することが必要である。

超音波を用いた変位推定には、解像度（つまり、極めて小さな変位を推定する能力）とは別に、高水準の精度も必要である。

多くの先行技術において、構造の変位は、超音波Ｂモード画像から推定される。

しかし、推定された変位の質は、Ｂモード画像の画質と解像度とに大きく依存する。

殆どの超音波装置において、Ｂモード画像の解像度では、マイクロメートル単位の変位を推定できない。

また、受信したＲＦ信号から、直接的に変位を推定することに焦点を当てた技術もある。

最も一般的な技術の１つは、［１］のような相互相関の技術である。

しかしながら、相互相関の技術では、計算量が大きくなり、複数のサンプリング点に対応する変位しか推定できない。

マイクロメートル単位の変位は、通常、サンプリング間隔のひと間隔のごく一部に相当するに過ぎない。

このため、相互相関を用いては推定できない。

そのような状況において、信号補間に基づいて変位を推定する技術がある。

しかし、この方法は、処理時間を増大させ、また、推定品質は、補間方法に依存する。

自己相関は、受信したＲＦ信号の直交復調信号（ベースバンド信号としても知られる）の位相情報に依存する。

この方法は、［２］に記載の、ＲＦ信号のサブサンプルに対応する変位を推定できるという利点を有する。

しかし、この方法では、ノイズが発生しやすく、また、振幅変調効果に影響される。

つまり、推定された変位は、信号電力の高い領域よりに、バイアスがかかっている。

ノイズの影響を克服するために、より大きなセットのサンプルを選択して、自己相関を実行することもできるが、この方法では、より詳細な変位を推定する能力が落ちる。

上記技術から、推定された変位の不正確性を克服する方法がいくつか開発されている。

［３］のように、「ｃｏａｒｓｅ ｔｏ ｆｉｎｅ」アプローチは、異なる推定段階に、異なるウインドウ領域を利用することで、第１段階では、変位を粗く推定し、第２段階では、より細かく推定して精度を向上させる。

しかし、推定品質を評価する方法は提供されていない。

［１］のように、第１段階の結果の精度を向上させるために、第２推定段階を実行し、第１段階の推定結果と組み合わせることで、第１段階の推定結果に基づく信号ワーピングを用いる技術もある。

［２］のように、他の変位修正方法も利用できる。

これらの方法は、発散を発生させるリスクが高いため、段階数を制限しなければならない。

また、収束に導き、かつＲＦ信号電力の不均一な分布によって起こるノイズおよびバイアスの影響を克服する方法はない。

変位推定品質を示すことができ、推定段階数を制限する必要なく、反復推定の収束を導き、ごく小さな変位に対しても精度の高い、より包括的な反復推定方法を提供するため、これらの方法を改善する必要がある。

このように、従来例として、特許文献１〜６の技術が知られる。

米国特許第６２７７０７４号明細書 米国特許第６５０６１５８号明細書（米国特許出願公開第２００１／００３４４８５号明細書） 米国特許出願公開第２００８／００１９６０９号明細書 米国特許出願公開第２００９／０２２１９１６号明細書 米国特許第６２７０４５９号明細書 国際公開第２００８／０３８６１５号

単一段階の変位推定方法に関し、現在の方法は、精度および変位解像度の点で課題を有する。

相互相関に基づく方法は、信号補間なしでは、サブサンプリング間隔に対応する小さな変位を推定することができない。

しかしながら、この方法では、ノイズの影響が少ない。

よって、相互相関ベースの方法では、精度は高いが、変位解像度は低い。

相互相関方法の変位解像度を向上させる目的で、信号補間を適用することができる。

しかし、これは、推定対象の変位が、マイクロメートル単位の場合に、処理電力を大幅に増大させる。

また、この方法の精度は、補間アルゴリズムの質に依存する。

自己相関に基づく方法では、小さな変位を推定することができるが、非常に、ノイズを発生しやすい。

よって、この方法では、変位解像度は高いが、精度は低い。

よって、精度を向上させるため、２段階の推定方法が利用できる。

しかし、これらの方法は、包括的な反復手法の欠如、変位の質を評価する方法の欠如、変位の精度を制限する要素、つまり、ノイズおよび振幅変調効果を解消する方法の欠如という３つの観点で、課題を有する。

包括的な反復手法には、反復推定の収束を導く技術を組み込む必要がある。

本発明の目的は、組織変位を精度よく推定できる変位推定装置（変位推定方法）を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。

なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。

上記の課題を解決するために、本発明の変位推定方法は、超音波信号を用いて、変位を反復的に推定する変位推定方法であって、少なくとも１つの超音波信号を走査して、媒体に送信するステップと、走査された前記媒体から反射した前記超音波信号を受信するステップと、ウインドウサイズを計算するステップと、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む変位推定方法である。

本発明において、発明者は、反復推定手法を用いる超音波を用いて、変位を推定する方法を教示する。

発明者の実験によると、推定ウインドウサイズが同じままか、または、全ての反復推定ラウンドに対して、適切に構成されていない場合には、推定された変位は、収束しない。

よって、この収束に役立つウインドウ計算方法を導入する。

ウインドウ計算方法は、反復推定の収束を導くために、ウインドウサイズを計算することの部分と、信号電力の不均一な分布によって起こるノイズおよびバイアスの影響を克服するために、ウインドウ境界を計算することの部分との２つの部分を含む。

ウインドウサイズの計算により、反復推定ラウンド毎に、異なるウインドウサイズを決定する。

これは、後続のラウンドにおいて、より詳細な推定結果を得るためと、推定を確実に収束させるためという２つの目的を有する。

つまり、不適切なウインドウサイズでの処理がされて、収束が得られず、正しい変位量が得られないことが回避されて、適切なウインドウサイズでの処理がされて、確実に収束が得られ、より正確な変位量が得られる。

ウインドウ境界の計算では、信号電力を利用し、各推定位置での、推定ウインドウの境界を決定する。

従来の推定方法において、ウインドウは通常、その位置を中心に対称である。

これは、不均一な信号電力の分布によって、バイアスが起こる一因である。

信号電力を用いて、ウインドウ境界を決定することにより、そのような影響を克服できる。

さらなる収束に導くため、推定された変位の質の評価を導入し、反復推定方法が、常に最良の結果に収束するよう保証する。

本発明は、収束を導く方法、精度を制限する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて、変位推定の精度を保証するものである。

図１０は、シミュレーションに、バックトラックアルゴリズムを適用した後の、変位推定の改善結果を示す図である。

図１１は、ファントム実験における、同じ結果を示す図である。

図１０は、推定の第１ラウンドにおける、推定された変位（変位１００１）が、シミュレーションされた変位プロファイル（プロファイル１０００）から逸脱していることを示す図である。

しかしながら、収束を導くウィンドウイング方式を用いるバックトラックアルゴリズムの最終ラウンドでは、推定された変位（変位１００２）は、シミュレーションされたプロファイルと、非常によく適合する。

図１１は、グローバル変位が生成されるファントム実験による推定結果を示す図である。

例示のため、単線が選択される（プロファイル１１００を参照）。

グローバル変位が生成されるので、一定の変位プロファイル（プロファイル１１０１）が期待される。

推定の第１ラウンドでは、結果（データ１１０２を参照）は、期待される変位に、あまり適合しない。

しかしながら、最終ラウンドの推定結果（データ１１０３を参照）は、第１ラウンドの推定結果よりも、グローバルな動きを、ずっとよく表している。

また、適切なウインドウサイズでの処理ができる変位推定装置を提供できる。

本発明の目的は、組織変位を精度よく推定できる変位推定装置（変位推定方法）を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。

なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。

図１は、本発明の変位推定方法を示す図である。 図２は、本発明のウインドウ計算方法を示す図である。 図３は、ウインドウサイズ計算の例を示す図である。 図４は、対称ウインドウおよびエネルギー均等ウインドウの例を示す図である。 図５は、本発明のＲＦ信号ワーピング方法を示す図である。 図６は、ＲＦ信号ワーピングおよびＲＦサンプル遅延を示す図である。 図７は、本発明の収束に導く方法を示す図である。 図８は、超音波ＲＦ信号を示す図である。 図９は、自己相関を用いた変位推定を示す図である。 図１０は、シミュレーションに本発明を用いた変位推定結果の例を示す図である。 図１１は、ファントム実験に本発明を用いた変位推定結果の例を示す図である。 図１２は、下層構造の変位、およびＲＦ信号への影響を示す図である。 図１３は、本装置のブロック図である。

以下、図面を参照して詳しく説明される。

実施形態の変位推定方法は、超音波信号（図１３の超音波信号１２０１ｓ）を用いて、変位（図１２の変位１２０７を参照）を反復的に推定する変位推定方法（図１３の変位推定装置１Ｘ参照）である。

つまり、例えば、変位の大きさの特定が複数回行われて、特定の処理が複数回行われることにより、比較的精度のよい、変位の大きさが特定される。

そして、少なくとも１つの超音波信号を走査して、媒体（媒体１２００Ｍ（図１２）、および、媒体１２００Ｍにある、下層構造（変位の測定の対象）１２０３ｘ（図１２、１３）を参照）に送信する（図１３の超音波処理部１Ｘａにある送信部１Ｘ１を参照）。

つまり、例えば、ライン方向（図８の方向８１Ｌ）の複数の位置へと超音波信号が送信されて、それら複数の位置の走査が行われる。

そして、走査された前記媒体から反射した前記超音波信号（超音波信号１２０１ｓ）を受信する（受信部１Ｘ２を参照）。

なお、超音波処理部１Ｘａは、例えば、プローブ（探触子）などである。

そして、ウインドウサイズ（図１３のウインドウサイズ１Ｘ３Ｍ：例えば、図３の（ａ）欄、（ｂ）欄のウインドウサイズ（winSize）など）を計算する（図１３の情報処理部１Ｘｂに含まれるサイズ計算部１Ｘ３を参照）。

つまり、それぞれのラウンドにおけるウインドウサイズが決定される。

なお、サイズ計算部１Ｘ３は、例えば、ウインドウ計算部１０２（図１）を含んでもよいし、ウインドウ計算部１０２の少なくとも一部などでもよい。

すなわち、例えば、第ｎ−１のラウンドでの相関値（図１３の相関値１Ｘ７Ｍ）に基づいて、第ｎのラウンドでのウインドウサイズが決定される（ｎ≧２）。

そして、計算された前記ウインドウサイズに基づいて、ウインドウの境界（境界の情報１Ｘ４Ｍ：図４の（ａ）欄、（ｂ）欄に示される、ウインドウの始端（ｗｉｎＳｔａｒｔ）、および、終端（ｗｉｎＥｎｄ）などを参照）を計算する（境界計算部１Ｘ４参照）。

なお、境界計算部１Ｘ４は、例えば、ウインドウ計算部１０２（図１）の少なくとも一部等でもよい。

つまり、それら始端および終端の間の幅が、計算されたウインドウサイズである始端および終端が特定される。

そして、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位（変位１ＸｂＭ）を推定する（推定部１Ｘ５参照）。

なお、推定部１Ｘ５は、例えば、例えば、変位推定部１０１（図１）の少なくとも一部等でもよい。

つまり、例えば、本装置からの距離が互いに異なる複数の位置における変位が特定される。

そして、推定された前記変位に基づいて、前記超音波信号をワーピングする（ワーピング部１Ｘ６参照）。

なお、ワーピング部１Ｘ６は、例えば、ＲＦ信号ワーピング部１００（図１）の少なくとも一部等でもよい。

そして、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値（相関値１Ｘ７Ｍ）が大きくなるように、この変位推定方法の収束を導く（収束制御部１Ｘ７参照）。

なお、収束制御部１Ｘ７は、例えば、収束ガイド部１０３の少なくとも一部等でもよい。

つまり、ワーピングされた超音波信号から算出される相関値が最大であるラウンドで推定された、変位の大きさが、最も精度の高い大きさとして特定される。

そして、計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される。

つまり、例えば、それぞれのラウンドで計算されるウインドウサイズは、そのラウンドよりも前のラウンドで計算されるウインドウサイズよりも小さいサイズであり、計算されるウインドウサイズが、小さいサイズへと変更される。

そして、前記ウインドウサイズは、全ての（深度の）ウインドウにおける信号エネルギーが等しくなるように計算される。

そして、超音波信号の相関値（第ｎ−１ラウンドでの相関値）が大きくなるほど、より小さい前記ウインドウサイズ（第ｎラウンドでのウインドウサイズ）が計算される。ここで、相関値など各ラウンドにおける収束度合いの評価値とウインドウサイズとの関係は、生体での実験値に基づく所定の関係式、あるいは、両者の対応関係を示すテーブルなどを参照して決定できる。

そして、第ｎ−１ラウンドでの相関値に対して、決定された対応関係が対応させるウインドウサイズが、第ｎラウンドでのウインドウサイズと特定される。

つまり、例えば、本装置において、次の動作が行われる。

すなわち、測定対象物（例えば新生血管を伴う悪性腫瘍など）である下層構造１２０３ｘ（図１２）へと、超音波信号１２０１ｓが送信され、送信された超音波信号１２０１ｓが受信されてもよい。

特定される遅延時間１２０８は、第１のパルスの位置（図１２の（ｂ）欄のＲＦ信号１２０４を参照）と、第２のパルスの位置（ＲＦ信号１２０５を参照）との間における遅延時間である。

ここで、第１のパルスは、受信された超音波信号１２０１ｓにおける、下層構造１２０３ｘが変位１２０７をする前のパルスである。

そして、第２のパルスは、測定対象物が変位１２０７をした後のパルスである。

そして、特定された遅延時間１２０８に基づいて変位１２０７を決定する。

そして、特定された、変位１２０７の大きさでの、下層構造１２０３ｘの動きが、予め定められた動きであるか否かが判定されてもよい。

そして、予め定められた動きと判定される場合に、当該下層構造１２０３ｘが、悪性腫瘍（癌）と判定され、当該動きではないと判定される場合に、良性腫瘍、あるいは、正常組織と判定されてもよい。

なお、このような、癌か否かの判定は、例えば、収束制御部１ｘ７（図１３）により行われてもよい。

このような本装置において、より具体的には、例えば、後で詳しく説明されるように、第ｎ−１のラウンドで特定された相関値が大きいほど、第ｎのラウンドにおいて、より大きいウインドウサイズでの処理がされてもよい。

なお、ここで、ウインドウサイズによる処理は、例えば、第ｎのラウンドにおける相関値が、受信される超音波信号１２０１ｓにおける、当該ウインドウサイズのウインドウの部分のデータから算出される処理でもよい。

本処理により、各ラウンドにおいて超音波の受信信号の信号レベルに応じた適切なウインドウが設定できるため、収束の精度が高まると共に、少ない反復回数で収束が得られ、高精度な変位推定が少ない処理量で実現できる。

以下の実施形態は、単に、様々な発明ステップの原理を説明するものである。

ここに説明する具体例の様々な変形は、当業者には明らかであろう。

したがって、本発明の技術範囲は、本明細書に記載の、具体的かつ例示的な内容によってではなく、特許請求項の範囲によってのみ限定されるものとする。

図１は、変位推定方法を示す図である。

本発明の主な実施形態は、図１に示される。

ＲＦＳｉｇ（ｄ，ｌ，ｆ）は、超音波部から得たＲＦ信号であり、ｄは深度方向（図８の方向８１Ｄ）、ｌはライン方向（方向８１Ｌ）、ｆはフレーム方向（方向８１Ｆ）を表す。

これらに関することは、さらに、図８に示される。

そのようなＲＦ信号の１フレーム内の１ライン（図示される８０１を参照）は、超音波トランスデューサ８００が、当該ラインの位置で走査方向（方向８０２）に向かってパルスを送信するメカニズムによって生成される。

ここで、当該パルスは、反射および散乱しながら、進路に沿って、下層構造と相互作用し、当該超音波トランスデューサ８００は、反射および散乱した信号を受信する。

そして、その信号は、超音波部によって、対応するラインに変換される。

同じメカニズムにより、複数のラインが様々な位置で生成され、複数のフレームが様々なタイムインスタンスに生成される。

図１２の（ａ）欄は、下層構造の変位を示し、（ｂ）欄は、ＲＦ信号への影響を示す。

走査された媒体の下層構造の変位が、ＲＦ信号の遅延として反射されるメカニズムが、図１２により示される。

例示を目的として、１ラインの１パルスが送信されるとし、下層構造には、反射境界が１つだけあるとする。

第１のタイムインスタンス（（ａ）欄の左側部分を参照）において、超音波トランスデューサ１２００は、走査された媒体に向けてパルスを送信する（パルス１２０１を参照）。

このパルスは、下層構造１２０３（１２０３ｘ）の反射境界に到達するまで、媒体を通して伝搬する。

この境界により、パルスはトランスデューサに向かって反射する（反射したパルス１２０２を参照）。

トランスデューサは、このパルスを受信し、ＲＦ信号１２０４（図１２の（ｂ）欄の表の第１行を参照）に変換する。

ＲＦ信号１２０４の全体における、パルスの位置は、送信されたパルス１２０１が、下層構造１２０３に向かって進み、反射して戻ってくるのに必要な時間を示す。

第２のタイムインスタンス（図１２の（ａ）欄の右側を参照）において、ＲＦ信号１２０５（図１２の（ｂ）欄の表の第３行）は、同一メカニズムにより作成されるが、下層構造は、新たな位置に変位している（変位した後の下層構造１２０６を参照）。

変位ｘ（変位１２０７）により、送信されたパルス（パルス１２０１）が、伝搬され、下層構造の境界で反射するのに要する時間が長くなり、その結果、ＲＦ信号の遅延が起こる（遅延時間１２０８を参照）。

なお、この例において、変位ｘの向きは、図１２に示される方向１２０３ｄの向きである。

この遅延時間１２０８の遅延値が分かると、変位（変位１２０７）を推測することができる。

そして、複数の下層構造を有する媒体では、複数の反射パルスが生成される。

それら複数の下層構造が変位する量が互いに異なることが原因で、ＲＦ信号における、反射パルスの遅延は、異なるものになる。

これらのパルスが重複し、かつ相互作用することにより、各下層構造の正確な変位を推定することが難しくなる。

本発明は、反復して推定を行うことにより、下層構造の変位を、正確に推定する方法を提示するものである。

図１に示される主な実施形態によると、本発明は、以下の主なブロックを備える。

つまり、本装置では、ＲＦ信号ワーピング部１００（図５参照）と、変位推定部１０１（図９参照）と、ウインドウ計算部１０２（図２参照）と、収束ガイド部１０３（図７参照）と、出力スイッチ１０４とを備える。

ＲＦ信号ワーピング部１００は、前回の推定ラウンド（第ｎ−１ラウンド）の変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）に基づいて、対象推定ラウンド（第ｎラウンド）においてのワーピングをする。このワーピングでは、対象推定ラウンドについて、ＲＦ信号の、選択されたフレーム内の、選択されたラインをワーピングする。

なお、これは、第１ラウンドでは、ゼロに設定される。

このブロックの出力は、推定ラウンドに用いられる、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）である。

なお、第１の推定ラウンド（つまり、ｒｏｕｎｄ ＝ １）では、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，０）をゼロに初期化したことにより、ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，１）と、超音波部から得たＲＦ信号であるＲＦＳｉｇ（ｄ，ｌ，ｆ） （先述）は、実質的に、互いに同一である。

変位を推定するため、予め定義されたフレームのセットを、このブロックに対する入力として選択することができる。

例示を目的として、２つのフレームが選択され、これらはｆ１およびｆ２で表される。

ＲＦ信号ワーピングの目的は、先行するラウンドで推定された変位に基づいて、１つのＲＦ信号を、他のＲＦ信号と適合するように、修正することである。

そして、その後、修正されたＲＦ信号から、残留変位を推定し、先行するラウンドの変位に追加することで、新たな変位を作成する。

いくつかの数（予め定められた数）の推定ラウンドの後、推定された変位は、修正されたＲＦ信号にほぼ一致するものとなり、残留変位はゼロに収束する。

変位推定部１０１は、ＲＦ信号の時間遅延から導出することで、ＲＦ信号ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）から、変位を推定する。

このタスクの実行に好ましい方法は、各深度における推定ウインドウが、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）によって記述される、自己相関の方法であるが、これに限定されるものではない。

ウインドウ計算部１０２は、各深度での、選択されたラインに用いられる推定ウィンドウパラメータｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）を計算する。

ウインドウ計算部１０２は、ＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を取得し、かつ、対象ラウンド数ｃｕｒＲｏｕｎｄと、ＲＦ信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とを、収束ガイド部１０３からの入力として取得する。

収束ガイド部１０３は、本技術の、収束を導く処理を行う。

つまり、推定された変位の質を決定するために、ワープされたＲＦ信号間の差分が計算される。

この差分が最小である（対象ラウンドを含む）ラウンドが、最も正確な変位推定結果を表していると判定される。

このブロックの出力は、対象推定ラウンド数ｃｕｒＲｏｕｎｄと、ＲＦ信号差分が最小であるラウンド数ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）と、推定された変位の質の測定値としてのＲＦ信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とである。

ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が小さいほど、質がよいことを示し、推定ラウンド全体にわたって、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が減少していると、収束を示す。

なお、ワープされたＲＦ信号間の差分を計算する方法として、「二乗差分和」および「絶対差分」の方法が考えられるが、これらに限定されるものではない。

出力スイッチ１０４は、入力として、ｃｕｒＲｏｕｎｄおよびｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）を取得する。

対象ラウンド数が最大値に到達した場合、最も正確な推定結果が、最終出力として選択される。

そうでなければ、その最も正確な推定結果は、後続の推定ラウンドに用いるために、ＲＦ信号ワーピング部１００の入力として選択され、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）およびｃｕｒＲｏｕｎｄが、入力として変位推定部１０１に送られる。

次の段落に、図１に示される主な実施形態に基づく具体的な実施形態をいくつか示す。

図５は、ＲＦ信号ワーピング方法を示す図である。

本発明の具体的な実施形態として、主な実施形態のＲＦ信号ワーピング部１００は、図５に示される方法によって実現される。

変位の結果、ＲＦ信号は移動され、かつ伸長される。

よって、先行するラウンドの、最も正確な変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））から、ＲＦＳｉｇ（ｄ，ｌ，ｆ１）を、ＲＦＳｉｇ（ｄ，ｌ，ｆ２）に適合させるサンプルにおける遅延値（図１２の遅延時間１２０８を参照）が、遅延値計算部５００において、以下のように算出される。

ここで、ｆｓは、ＲＦ信号のサンプリング周波数であり、ｃは、走査された媒体における音速である。

走査された媒体における変位は、結果として、得られたＲＦ信号における遅延として反映される。

各深度での遅延値が既知であるため、ＲＦＳｉｇ（ｄ，ｌ，ｆ１）またはＲＦＳｉｇ（ｄ，ｌ，ｆ２）の何れかにおける各ラインの、（全ての深度に対応する）全てのサンプルを、他のサンプルに適合するように合わせて、遅延させることができる。

これは、ＲＦサンプル遅延部５０１において実行される。

１つの実施形態において、ＲＦサンプル遅延部５０１における遅延は、信号補間を用いて実行される。

別の実施形態では、ＲＦサンプル遅延部５０１における遅延は、分数遅延フィルタを用いて実行される。

図６は、ＲＦ信号ワーピングおよびＲＦサンプル遅延を示す図である。

ＲＦ信号ワーピング処理が、図６に示される。

変位により、ＲＦ信号の遅延が起こることが、先述の図１２により示される。

ＲＦ信号ワーピングの目的は、推定された変位と逆方向に、１つのＲＦ信号を遅延させることによって、この効果を無効にすることである。

この無効化の結果として、ワープされたＲＦ信号が、相互に適合することが期待される。

遅延値計算部５００から計算されたｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）は、１つのＲＦ信号内の各サンプルが、これを達成するために遅延される値を表す。

２つのＲＦ信号６００および６０１（図６）を、この例示のために用いる。

データ６０２は、ある推定ラウンドの後に、遅延値計算部５００から計算されたｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）を示す。

ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）の特定の深度での値が、ＲＦ信号６００における、同じ深度でのサンプルが、遅延値計算部６０１に適合するための遅延に必要な量を表すことが分かる。

しかしながら、ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）における遅延値は、常に整数であるわけではない。

よって、ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）に従って、ＲＦ信号６００における各サンプルを遅延できるようにするため、信号補間または分数遅延フィルタが、２つの選択候補として挙げられる。

変位推定部１０１をＲＦ信号に適用することで、第１ラウンドに用いられる変位と、後続ラウンドに用いられる残留変位を得て、最終的な変位を得るために、それらを組み合わせる。

好ましい実施形態が図９に示されるが、これに限定されるものではない。

図９は、自己相関を用いた変位推定を表す。

ＩＱ復調部９００は、ＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を、ベースバンド信号ＩＱＳｉｇ（ｄ，ｌ，ｆ）に変換する。

自己相関計算部９０１は、以下の式に従って、ａｕｔｏｃｏｒｒ（ｄ，ｌ）を計算する。

ここで、ｃｏｎｊ（）は、複素共役を表す。

最後に、以下の式に従って、変位計算部９０２により、ａｕｔｏｃｏｒｒ（ｄ，ｌ）が、変位に変換される。

ここで、ａｒｇ（）は、複素数の偏角を計算する関数である。

１から開始するラウンド（第１ラウンドを表す）により、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，０）は、ゼロに初期化される。

図２は、ウインドウ計算方法を示す図である。

ウインドウ計算部１０２により、各推定ラウンドにおける、各ラインの各深度に対し、異なる推定ウインドウを決定する。

これは図２に示される。

図３は、ウインドウサイズ計算の具体例を示す図である。

ウインドウサイズ計算部２００（図２）は、各推定ラウンドに対するウインドウサイズの決定をする。

このブロックの出力は、ｗｉｎＳｉｚｅ（ｄ，ｌ）である。

好ましい実施形態において、各ラウンドに対するウインドウサイズは、ラウンドの関数が減少するにつれて、変化する。

これは図３に示される。

ウインドウサイズ計算部２００（図２）のある実施形態において、図３の（ａ）欄に示されるように、ウインドウサイズは、減少関数ｗｉｎＳｉｚｅＲｏｕｎｄ（ｒｏｕｎｄ）により、推定ラウンド数に、直接関係する。

ｃｕｒＲｏｕｎｄによって特定される対象推定ラウンドについて、以下に基づいて、全ての深度に対するウインドウサイズが特定される。

全てのｄおよびｌについて、
ｗｉｎＳｉｚｅ（ｄ，ｌ） ＝ ｗｉｎＳｉｚｅＲｏｕｎｄ（ｃｕｒＲｏｕｎｄ）

「固定ウィンドウイングスキーム」は、変位推定に用いられる全ての深度で推定ウインドウサイズとなるｗｉｎＳｉｚｅ（ｄ，ｌ）の同じ値を選択する処理である。

「固定ウィンドウイングスキーム」の利点は、必要な演算量が多くなく（比較的少なく）、かつ、ＲＦ信号エネルギーが、深度方向に沿って、規則正しく配分されている場合に、適していることである。

ウインドウサイズ計算部２００の別の実施形態において、各深度でのウインドウサイズは、上限であるｍａｘＷｉｎ（ｒｏｕｎｄ）と、下限であるｍｉｎＷｉｎ（ｒｏｕｎｄ）と、ＲＦ信号電力ｓｉｇＰｏｗ（ｄ，ｌ）とに基づいて計算される。

ｍａｘＷｉｎ（ｒｏｕｎｄ）およびｍｉｎＷｉｎ（ｒｏｕｎｄ）は、ラウンドの減少関数である。

これは、「瞬間電力ウィンドウイングスキーム」として、図３の（ｂ）欄に示される。

以下の条件が適用される。

また、図３の（ｂ）欄に示されるように、ｗｉｎＳｉｚｅ（ｄ，ｌ）は、ｓｉｇＰｏｗ（ｄ，ｌ）、つまり、ｓｉｇＰｏｗ ＝ ０に対応する最大ウインドウサイズと、最大ｓｉｇＰｏｗに対応する最小ウインドウサイズとに反比例してもよい。

「瞬間電力ウィンドウイングスキーム」は、「固定ウィンドウイングスキーム」と比較して、必要な演算量が多いが、特に、ＲＦ信号エネルギーが、深度方向に沿って、かなり不均一である場合に、より精度が高い。

信号対ノイズ比の低さを克服するために、低い信号電力を有する領域に対して、より大きな推定ウインドウが指定され、より詳細に変位を分布するために、大きな信号電力を有する領域に対して、より小さな推定ウインドウが指定される。

しかしながら、信号エネルギーが規則正しく配分されているＲＦ信号に関しては、「固定ウィンドウイングスキーム」と「瞬間電力ウィンドウイングスキーム」とは大差がない。

ウインドウサイズ計算部２００の別の実施形態において、全ての深度でのウインドウサイズを、全てのウインドウ内の信号エネルギーが同量に保たれるように計算してもよく、これは、「一定エネルギーウィンドウイングスキーム」と称される。

信号電力計算部２０１（図２）により、フレーム方向に沿う、ＲＦ信号の平均電力が計算される。

信号電力計算部２０１は、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）と、対象ラウンド数ｃｕｒＲｏｕｎｄとを、入力として取得する。

ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｃｕｒＲｏｕｎｄ）内には、複数のフレームがあり、電力が、以下の通り計算される。

ここで、Ｎは、計算に用いるために選択されたフレームの総数である。

図４は、対称ウインドウおよびエネルギー均等ウインドウの例を示す図である。

ウインドウ境界計算部２０２（図２）により、各ラインの各深度でのウインドウ境界が計算される。

これは図４に示される。

ウインドウ境界計算部２０２のある実施形態においては、ウインドウサイズ計算部２００から計算されたｗｉｎＳｉｚｅ（ｄ，ｌ）が入力であり、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）が、以下の式に従って、対応する深度を中心に、ウインドウが対称になるように計算される（「対称ウインドウ」と称される）。

ここで、ＲＯＵＮＤ（）は、最も近い整数に丸める処理を示す。

これは、図４の（ａ）欄に示される。

対称ウインドウの利点は、必要とされる演算量が多くないことである。

これは、ＲＦ信号電力が、深度方向において、サンプル毎に突然変化しない場合に適している。

ウインドウ境界計算部２０２の別の実施形態においては、次の通りである。

つまり、ウインドウサイズ計算部２００から計算されたｗｉｎＳｉｚｅ（ｄ，ｌ）が、ウインドウ境界計算部２０２の入力である。

そして、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）が、以下の式により、各深度の両側での信号エネルギーが等しくなるように（「エネルギー均等ウインドウ」と称される）計算される。

これは、図４の（ｂ）欄に示される。

「エネルギー均等ウインドウ」は、「対称ウインドウ」と比較して、必要な演算量が多いが、従来のウィンドウイング技術においてバイアスの原因となっていた、深度方向に沿う、信号電力の不均一な分布の影響を抑制できる。

図７は、収束に導く方法を示す図である。

図１の収束ガイド部１０３による、本発明において説明した方法で、収束が可能であることが保証される。

収束ガイド部１０３は、図７に示されるように、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を、当該収束ガイド部１０３の入力として取得する。

ＲＦ信号差分計算部７００は、ワープされたＲＦ信号間の差分を計算し、各ラインの各深度について、各推定ラウンドにおける差分を、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）として出力する。

ＲＦ信号ワーピング部１００の主な目的は、ＲＦ信号の１つを、他と適合するように、修正することであるので、この残差（ｄ，ｌ，ｒｏｕｎｄ）は、ＲＦ信号同士が、どの程度適合しているかの尺度になる。

（ｒｏｕｎｄにより特定される）ある特定の推定ラウンドについて、推定された変位ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）の質を評価する目的で、ＲＦ信号ワーピング部１００が、当該ＲＦ信号ワーピング部１００の入力として、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）を取得する。

このことにより、ＲＦ信号がワープされ、適合レベルを評価するために、ワープされた信号間の差分が計算される。

これらのワープされた信号は、後続ラウンドの変位推定のための入力としても用いられ、よってＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ＋１）と表される。

よって、ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ＋１）から、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が計算される。

ＲＦ信号差分計算部７００（図７）で行われる計算の一例として、以下の式（つまり、二乗差分和）が挙げられるが、これに限定されるものではない。

ここで、ｆ１およびｆ２は、推定に用いるために選択された２フレームを表す。

この例において、１ラインの全ての深度に対し同一の信号差分値が指定され、これは、対応するラインのＲＦ信号が、全体的に適合していることを表す。

小さな領域で変位が起こるだけではないときには、「二乗差分和」が適しており、その結果、ライン全体に反復が起こる。

なお、ＲＦ信号差分計算部７００（図７）で行われる計算の別の例として、以下の式（つまり、絶対差分）が挙げられるが、これに限定されるものではない。

ここで、ｆ１およびｆ２は、推定に用いるために選択された２フレームを表す。

この例において、各ラインの各深度は、個別の信号差分値を有し、これは、ＲＦ信号が、局所的に適合していることを表す。

深度方向に小さな領域だけで、変位が起こる場合には、「絶対差分」が適しており、その結果、その領域内での局所的反復が可能になる。

最小残余計算部７０１は、以下の式に基づき、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）における、各深度およびラインについて、残余が最小であるラウンド数を確認する。

これは、最も正確な推定結果を表す。

ここで関数

は、ｆ（ｘ）を最小化するｘの値を返す。

収束ガイド部１０３（図１）の出力は、対象ラウンド数ｃｕｒＲｏｕｎｄと、残余が最小であるラウンド数ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）と、信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とである。

出力スイッチ１０４（図１）は、入力として、ｃｕｒＲｏｕｎｄおよびｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）を取得する。

ｃｕｒＲｏｕｎｄが最大値に達すると、出力スイッチ１０４は、最も精度の高い変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））を、最終的な出力ｄｉｓｐＯｕｔとして選択する。

そうでなければ、対象ラウンドｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｃｕｒＲｏｕｎｄ）が、変位推定される代わりに、次の推定ラウンドのための、ＲＦ信号ワーピング部１００の入力として用いられるｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））が、変位推定される。

ここで、変位推定の多次元（２次元 or ３次元）への拡張について説明する。

つまり、上記の実施の形態では、超音波の深さ方向における１次元の変位を求める方法について説明した。

しかしながら、血管や腫瘍組織などの生体組織の拍動方向は、１次元内に限定されず、また、場所によっても、方向が異なる。

したがって、複数方向から超音波を送受信することで、特定位置における複数方向の変位を求め、組織の２次元あるいは３次元的な変位を測定してもよい。

こうして、多次元で測定することにより、診断精度の更なる向上が見込まれる。

そして、超音波を複数方向から送受信する方法としては、超音波振動子が、１次元に配置されたリニアプローブを用いて、プローブの方向を変化させる、または、超音波振動子が２次元アレイ上に配置されたマトリックスプローブを用いて、ビームフォーミングにより、複数方向での送受信を行うなどの方法が可能である。

ここで、リニアプローブの使用時には、磁気センサや加速度センサ、ジャイロなどの位置センサ、あるいは、カメラにより、プローブの位置を追跡することで、プローブの位置や向きを取得して、送受信時の方向を決定できる。

なお、変位測定は、超音波の深さ方向など、１方向についてのみ行い、この変位ベクトルを、所定の２次元あるいは３次元座標系の各軸方向の成分に分解することで、多次元的に評価してもよい。

この時、超音波の深さ方向と所定の座標系との対応関係が必要となるが、上述の位置センサなどを用いて、プローブの位置と方向を取得することで、両者を対応付ける。

ここで、癌検出、あるいは、動脈硬化（血管の硬さ）の診断などに適用できる歪み計測（エラストグラフィ）との関連について補足説明する。歪みは、変位の時間微分であり、組織の硬さの指標となるが、本発明により、変位推定精度が向上することで、結果として、歪みの測定精度が向上する。エラストグラフィでは微分値を用いるが、その他に、組織の特定部位の変位の時間変化波形なども、組織性状の判定に利用できる。例えば、癌の成長過程においては、癌の周辺および内部に血管新生が起こることが知られている。従って、これら血管の拍動、あるいは、血管の拍動に伴う周囲の組織の拍動などを変位量の時間変化として捉え、変化のパターンや振幅が癌に特有のものであるかどうかを判定することで、癌検出を行ってもよい。

こうして、変位推定装置であって、超音波信号を用いて、変位を反復的に推定する変位推定装置であって、少なくとも１つの超音波信号を、トランスデューサから走査された媒体に送信する手段と、前記走査された媒体から反射した超音波信号を受信する手段と、推定ウインドウサイズを計算する手段と、前記推定されたウインドウサイズに基づいて、前記推定ウインドウの境界を計算する手段と、前記計算された推定ウインドウを用いて、前記超音波信号の深度方向に沿って変位を推定する手段と、前記推定された変位に基づいて、前記超音波信号をワーピングする手段と、前記方法の収束傾向を前記ワープされた超音波信号間の差分として計算することで、収束を導く手段とを備える変位推定装置が開示される。

こうして、開示される方法は、収束を導く方法、精度を限定する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて変位推定の精度を保証する。

つまり、次の問題が解決される。つまり、超音波を用いる変位推定において、先行技術では、推定の精度と解像度を高める方法を提供している。その中には、複数段階からなる推定または反復推定を導入しているものもある。しかし、（１）反復推定の収束を導く方法を組み込んだ包括的な反復手法の欠如、（２）品質評価方法の欠如、および（３）変位の精度を限定する要素を克服する方法の欠如という問題があり、本発明の変位推定方法により、これらの問題が解決される。

なお、本発明は、装置、システム、集積回路などとして実現できるだけでなく、その装置等を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを記録したコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体として実現したり、そのプログラムを示す情報、データ又は信号として実現したりすることもできる。そして、それらプログラム、情報、データ及び信号は、インターネット等の通信ネットワークを介して配信してもよい。

なお、本発明について、実施の形態に基づいて説明したが、本発明は、上記の実施の形態に限定されるものではない。本発明の主旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の範囲内に含まれる。

本技術は、超音波信号を用いて変位を推定する方法を導入するものである。これは、最終的な結果、または、更なる処理のための中間ステップの何れであっても、変位推定が必要な適用例において、用いることができ、医療用および工業用超音波機器に用いることができる。

変位を精度よく推定できる変位推定装置が提供でき、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供可能な変位推定装置が提供できる。

１Ｘ 変位推定装置
１Ｘ１ 送信部
１Ｘ２ 受信部
１Ｘ３ サイズ計算部
１Ｘ４ 境界計算部
１Ｘ５ 推定部
１Ｘ６ ワーピング部
１Ｘ７ 収束制御部
１２０１ｓ 超音波信号
１２０３ｘ 下層構造
１Ｘ３Ｍ ウインドウサイズ
１Ｘ４Ｍ 境界の情報
１ＸｂＭ 変位
１Ｘ７Ｍ 相関値
１００ ＲＦ信号ワーピング部
１０３ 収束ガイド部
１０１ 変位推定部
１０２ ウインドウ計算部
１０３ 収束ガイド部
１０４ 出力スイッチ
６００、１２０４、１２０５ ＲＦ信号
８００、１２００ 超音波トランスデューサ
１２００Ｍ 媒体
１２０８ 遅延時間
１００ａ、１０１ａ、１０２１ａ、１０３ａ、１０４ａ、１０４ｂ、１０４ｃ 情報
２０１ｎ、２０１ｍ、２０１ａ、２００ｍ、２００ａ、２０２ａ 情報
５００ｍ、５００ａ、５０１ｎ、５０１ａ 情報
７００ｍ、７００ａ、７０１ａ、７００ｂ 情報
９００ｍ、９００ａ、９０１ｍ、９０１ａ、９０２ａ 情報

本発明は、超音波信号を用いて変位を推定する方法を導入するものである。これは、最終的な結果、または、更なる処理のための中間ステップの何れであっても、変位推定が必要な適用例において、用いることができ、医療用および工業用超音波機器に用いることができる。

超音波装置は、高周波の機械的音波を、走査対象の媒体に対して送受信することにより、媒体を非侵襲で診断するものである。

そのような装置のトランスデューサは、超音波を、走査対象の媒体に送信する。

超音波は、散乱および反射することにより、その媒体の下層構造と相互作用する。

なお、下層構造は、例えば、媒体の表面よりも、媒体の内部側における構造などをいう。下層構造は、例えば、人間の身体における血管などである。

散乱および反射された超音波は、下層構造の有益な情報を含み、トランスデューサによって受信される。その後、超音波装置によって処理され、例えば画像情報としてユーザに提示される。

図８は、超音波ＲＦ信号を示す図である。

超音波装置が反射波から得る最も基礎的なデータの種類の１つとして、ＲＦ（Radiofrequency Signal）信号が挙げられる（図８）。

これは、受信波を、アナログからデジタルに直接変換したものである。

異なる適用例では、超音波ＲＦ信号から、輝度モード（Ｂモード）画像、ドップラー画像等、他の種類のデータを導出できる。

多くの適用例のうちの１つは、走査された媒体の下層構造の動きを分析することである。

ドップラー効果は、構造の動きの方向および強度を推測する簡易な方法として利用される。

しかし、ドップラー効果の精度は、非常に限られたものである。

医療用超音波における、血流のモニタリングのような、精度が厳密に要求されない適用例には、ドップラー超音波は適している。

しかし、構造が小さく、かつ動きが極めて小さい血管などが対象物であり、高水準の精度が必要とされる場合には、より感度の高い技術が必要である。

最近の新しい適用例として、構造の弾性を推測するために、構造の変位を用いる超音波弾性率計測法がある。

弾性を正確に推定するためには、構造の変位を、受信した超音波信号から、正確に推定することが必要である。

超音波を用いた変位推定には、解像度（つまり、極めて小さな変位を推定する能力）とは別に、高水準の精度も必要である。

多くの先行技術において、構造の変位は、超音波Ｂモード画像から推定される。

しかし、推定された変位の質は、Ｂモード画像の画質と解像度とに大きく依存する。

殆どの超音波装置において、Ｂモード画像の解像度では、マイクロメートル単位の変位を推定できない。

また、受信したＲＦ信号から、直接的に変位を推定することに焦点を当てた技術もある。

最も一般的な技術の１つは、［１］のような相互相関の技術である。

しかしながら、相互相関の技術では、計算量が大きくなり、複数のサンプリング点に対応する変位しか推定できない。

マイクロメートル単位の変位は、通常、サンプリング間隔のひと間隔のごく一部に相当するに過ぎない。

このため、相互相関を用いては推定できない。

そのような状況において、信号補間に基づいて変位を推定する技術がある。

しかし、この方法は、処理時間を増大させ、また、推定品質は、補間方法に依存する。

自己相関は、受信したＲＦ信号の直交復調信号（ベースバンド信号としても知られる）の位相情報に依存する。

この方法は、［２］に記載の、ＲＦ信号のサブサンプルに対応する変位を推定できるという利点を有する。

しかし、この方法では、ノイズが発生しやすく、また、振幅変調効果に影響される。

つまり、推定された変位は、信号電力の高い領域よりに、バイアスがかかっている。

ノイズの影響を克服するために、より大きなセットのサンプルを選択して、自己相関を実行することもできるが、この方法では、より詳細な変位を推定する能力が落ちる。

上記技術から、推定された変位の不正確性を克服する方法がいくつか開発されている。

［３］のように、「ｃｏａｒｓｅ ｔｏ ｆｉｎｅ」アプローチは、異なる推定段階に、異なるウインドウ領域を利用することで、第１段階では、変位を粗く推定し、第２段階では、より細かく推定して精度を向上させる。

しかし、推定品質を評価する方法は提供されていない。

［１］のように、第１段階の結果の精度を向上させるために、第２推定段階を実行し、第１段階の推定結果と組み合わせることで、第１段階の推定結果に基づく信号ワーピングを用いる技術もある。

［２］のように、他の変位修正方法も利用できる。

これらの方法は、発散を発生させるリスクが高いため、段階数を制限しなければならない。

また、収束に導き、かつＲＦ信号電力の不均一な分布によって起こるノイズおよびバイアスの影響を克服する方法はない。

変位推定品質を示すことができ、推定段階数を制限する必要なく、反復推定の収束を導き、ごく小さな変位に対しても精度の高い、より包括的な反復推定方法を提供するため、これらの方法を改善する必要がある。

このように、従来例として、特許文献１〜６の技術が知られる。

米国特許第６２７７０７４号明細書 米国特許第６５０６１５８号明細書（米国特許出願公開第２００１／００３４４８５号明細書） 米国特許出願公開第２００８／００１９６０９号明細書 米国特許出願公開第２００９／０２２１９１６号明細書 米国特許第６２７０４５９号明細書 国際公開第２００８／０３８６１５号

単一段階の変位推定方法に関し、現在の方法は、精度および変位解像度の点で課題を有する。

相互相関に基づく方法は、信号補間なしでは、サブサンプリング間隔に対応する小さな変位を推定することができない。

しかしながら、この方法では、ノイズの影響が少ない。

よって、相互相関ベースの方法では、精度は高いが、変位解像度は低い。

相互相関方法の変位解像度を向上させる目的で、信号補間を適用することができる。

しかし、これは、推定対象の変位が、マイクロメートル単位の場合に、処理電力を大幅に増大させる。

また、この方法の精度は、補間アルゴリズムの質に依存する。

自己相関に基づく方法では、小さな変位を推定することができるが、非常に、ノイズを発生しやすい。

よって、この方法では、変位解像度は高いが、精度は低い。

よって、精度を向上させるため、２段階の推定方法が利用できる。

しかし、これらの方法は、包括的な反復手法の欠如、変位の質を評価する方法の欠如、変位の精度を制限する要素、つまり、ノイズおよび振幅変調効果を解消する方法の欠如という３つの観点で、課題を有する。

包括的な反復手法には、反復推定の収束を導く技術を組み込む必要がある。

本発明の目的は、組織変位を精度よく推定できる変位推定装置（変位推定方法）を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。

なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。

上記の課題を解決するために、本発明の変位推定方法は、超音波信号を用いて、変位を反復的に推定する変位推定方法であって、少なくとも１つの超音波信号を走査して、媒体に送信するステップと、走査された前記媒体から反射した前記超音波信号を受信するステップと、ウインドウサイズを計算するステップと、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む変位推定方法である。

本発明において、発明者は、反復推定手法を用いる超音波を用いて、変位を推定する方法を教示する。

発明者の実験によると、推定ウインドウサイズが同じままか、または、全ての反復推定ラウンドに対して、適切に構成されていない場合には、推定された変位は、収束しない。

よって、この収束に役立つウインドウ計算方法を導入する。

ウインドウ計算方法は、反復推定の収束を導くために、ウインドウサイズを計算することの部分と、信号電力の不均一な分布によって起こるノイズおよびバイアスの影響を克服するために、ウインドウ境界を計算することの部分との２つの部分を含む。

ウインドウサイズの計算により、反復推定ラウンド毎に、異なるウインドウサイズを決定する。

これは、後続のラウンドにおいて、より詳細な推定結果を得るためと、推定を確実に収束させるためという２つの目的を有する。

つまり、不適切なウインドウサイズでの処理がされて、収束が得られず、正しい変位量が得られないことが回避されて、適切なウインドウサイズでの処理がされて、確実に収束が得られ、より正確な変位量が得られる。

ウインドウ境界の計算では、信号電力を利用し、各推定位置での、推定ウインドウの境界を決定する。

従来の推定方法において、ウインドウは通常、その位置を中心に対称である。

これは、不均一な信号電力の分布によって、バイアスが起こる一因である。

信号電力を用いて、ウインドウ境界を決定することにより、そのような影響を克服できる。

さらなる収束に導くため、推定された変位の質の評価を導入し、反復推定方法が、常に最良の結果に収束するよう保証する。

本発明は、収束を導く方法、精度を制限する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて、変位推定の精度を保証するものである。

図１０は、シミュレーションに、バックトラックアルゴリズムを適用した後の、変位推定の改善結果を示す図である。

図１１は、ファントム実験における、同じ結果を示す図である。

図１０は、推定の第１ラウンドにおける、推定された変位（変位１００１）が、シミュレーションされた変位プロファイル（プロファイル１０００）から逸脱していることを示す図である。

しかしながら、収束を導くウィンドウイング方式を用いるバックトラックアルゴリズムの最終ラウンドでは、推定された変位（変位１００２）は、シミュレーションされたプロファイルと、非常によく適合する。

図１１は、グローバル変位が生成されるファントム実験による推定結果を示す図である。

例示のため、単線が選択される（プロファイル１１００を参照）。

グローバル変位が生成されるので、一定の変位プロファイル（プロファイル１１０１）が期待される。

推定の第１ラウンドでは、結果（データ１１０２を参照）は、期待される変位に、あまり適合しない。

しかしながら、最終ラウンドの推定結果（データ１１０３を参照）は、第１ラウンドの推定結果よりも、グローバルな動きを、ずっとよく表している。

また、適切なウインドウサイズでの処理ができる変位推定装置を提供できる。

本発明の目的は、組織変位を精度よく推定できる変位推定装置（変位推定方法）を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。

なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。

図１は、本発明の変位推定方法を示す図である。 図２は、本発明のウインドウ計算方法を示す図である。 図３は、ウインドウサイズ計算の例を示す図である。 図４は、対称ウインドウおよびエネルギー均等ウインドウの例を示す図である。 図５は、本発明のＲＦ信号ワーピング方法を示す図である。 図６は、ＲＦ信号ワーピングおよびＲＦサンプル遅延を示す図である。 図７は、本発明の収束に導く方法を示す図である。 図８は、超音波ＲＦ信号を示す図である。 図９は、自己相関を用いた変位推定を示す図である。 図１０は、シミュレーションに本発明を用いた変位推定結果の例を示す図である。 図１１は、ファントム実験に本発明を用いた変位推定結果の例を示す図である。 図１２は、下層構造の変位、およびＲＦ信号への影響を示す図である。 図１３は、本装置のブロック図である。

以下、図面を参照して詳しく説明される。

実施形態の変位推定方法は、超音波信号（図１３の超音波信号１２０１ｓ）を用いて、変位（図１２の変位１２０７を参照）を反復的に推定する変位推定方法（図１３の変位推定装置１Ｘ参照）である。

つまり、例えば、変位の大きさの特定が複数回行われて、特定の処理が複数回行われることにより、比較的精度のよい、変位の大きさが特定される。

そして、少なくとも１つの超音波信号を走査して、媒体（媒体１２００Ｍ（図１２）、および、媒体１２００Ｍにある、下層構造（変位の測定の対象）１２０３ｘ（図１２、１３）を参照）に送信する（図１３の超音波処理部１Ｘａにある送信部１Ｘ１を参照）。

つまり、例えば、ライン方向（図８の方向８１Ｌ）の複数の位置へと超音波信号が送信されて、それら複数の位置の走査が行われる。

そして、走査された前記媒体から反射した前記超音波信号（超音波信号１２０１ｓ）を受信する（受信部１Ｘ２を参照）。

なお、超音波処理部１Ｘａは、例えば、プローブ（探触子）などである。

そして、ウインドウサイズ（図１３のウインドウサイズ１Ｘ３Ｍ：例えば、図３の（ａ）欄、（ｂ）欄のウインドウサイズ（winSize）など）を計算する（図１３の情報処理部１Ｘｂに含まれるサイズ計算部１Ｘ３を参照）。

つまり、それぞれのラウンドにおけるウインドウサイズが決定される。

なお、サイズ計算部１Ｘ３は、例えば、ウインドウ計算部１０２（図１）を含んでもよいし、ウインドウ計算部１０２の少なくとも一部などでもよい。

すなわち、例えば、第ｎ−１のラウンドでの相関値（図１３の相関値１Ｘ７Ｍ）に基づいて、第ｎのラウンドでのウインドウサイズが決定される（ｎ≧２）。

そして、計算された前記ウインドウサイズに基づいて、ウインドウの境界（境界の情報１Ｘ４Ｍ：図４の（ａ）欄、（ｂ）欄に示される、ウインドウの始端（ｗｉｎＳｔａｒｔ）、および、終端（ｗｉｎＥｎｄ）などを参照）を計算する（境界計算部１Ｘ４参照）。

なお、境界計算部１Ｘ４は、例えば、ウインドウ計算部１０２（図１）の少なくとも一部等でもよい。

つまり、それら始端および終端の間の幅が、計算されたウインドウサイズである始端および終端が特定される。

そして、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位（変位１ＸｂＭ）を推定する（推定部１Ｘ５参照）。

なお、推定部１Ｘ５は、例えば、例えば、変位推定部１０１（図１）の少なくとも一部等でもよい。

つまり、例えば、本装置からの距離が互いに異なる複数の位置における変位が特定される。

そして、推定された前記変位に基づいて、前記超音波信号をワーピングする（ワーピング部１Ｘ６参照）。

なお、ワーピング部１Ｘ６は、例えば、ＲＦ信号ワーピング部１００（図１）の少なくとも一部等でもよい。

そして、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値（相関値１Ｘ７Ｍ）が大きくなるように、この変位推定方法の収束を導く（収束制御部１Ｘ７参照）。

なお、収束制御部１Ｘ７は、例えば、収束ガイド部１０３の少なくとも一部等でもよい。

つまり、ワーピングされた超音波信号から算出される相関値が最大であるラウンドで推定された、変位の大きさが、最も精度の高い大きさとして特定される。

そして、計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される。

つまり、例えば、それぞれのラウンドで計算されるウインドウサイズは、そのラウンドよりも前のラウンドで計算されるウインドウサイズよりも小さいサイズであり、計算されるウインドウサイズが、小さいサイズへと変更される。

そして、前記ウインドウサイズは、全ての（深度の）ウインドウにおける信号エネルギーが等しくなるように計算される。

そして、超音波信号の相関値（第ｎ−１ラウンドでの相関値）が大きくなるほど、より小さい前記ウインドウサイズ（第ｎラウンドでのウインドウサイズ）が計算される。ここで、相関値など各ラウンドにおける収束度合いの評価値とウインドウサイズとの関係は、生体での実験値に基づく所定の関係式、あるいは、両者の対応関係を示すテーブルなどを参照して決定できる。

そして、第ｎ−１ラウンドでの相関値に対して、決定された対応関係が対応させるウインドウサイズが、第ｎラウンドでのウインドウサイズと特定される。

つまり、例えば、本装置において、次の動作が行われる。

すなわち、測定対象物（例えば新生血管を伴う悪性腫瘍など）である下層構造１２０３ｘ（図１２）へと、超音波信号１２０１ｓが送信され、送信された超音波信号１２０１ｓが受信されてもよい。

特定される遅延時間１２０８は、第１のパルスの位置（図１２の（ｂ）欄のＲＦ信号１２０４を参照）と、第２のパルスの位置（ＲＦ信号１２０５を参照）との間における遅延時間である。

ここで、第１のパルスは、受信された超音波信号１２０１ｓにおける、下層構造１２０３ｘが変位１２０７をする前のパルスである。

そして、第２のパルスは、測定対象物が変位１２０７をした後のパルスである。

そして、特定された遅延時間１２０８に基づいて変位１２０７を決定する。

そして、特定された、変位１２０７の大きさでの、下層構造１２０３ｘの動きが、予め定められた動きであるか否かが判定されてもよい。

そして、予め定められた動きと判定される場合に、当該下層構造１２０３ｘが、悪性腫瘍（癌）と判定され、当該動きではないと判定される場合に、良性腫瘍、あるいは、正常組織と判定されてもよい。

なお、このような、癌か否かの判定は、例えば、収束制御部１ｘ７（図１３）により行われてもよい。

このような本装置において、より具体的には、例えば、後で詳しく説明されるように、第ｎ−１のラウンドで特定された相関値が大きいほど、第ｎのラウンドにおいて、より大きいウインドウサイズでの処理がされてもよい。

なお、ここで、ウインドウサイズによる処理は、例えば、第ｎのラウンドにおける相関値が、受信される超音波信号１２０１ｓにおける、当該ウインドウサイズのウインドウの部分のデータから算出される処理でもよい。

本処理により、各ラウンドにおいて超音波の受信信号の信号レベルに応じた適切なウインドウが設定できるため、収束の精度が高まると共に、少ない反復回数で収束が得られ、高精度な変位推定が少ない処理量で実現できる。

以下の実施形態は、単に、様々な発明ステップの原理を説明するものである。

ここに説明する具体例の様々な変形は、当業者には明らかであろう。

したがって、本発明の技術範囲は、本明細書に記載の、具体的かつ例示的な内容によってではなく、特許請求項の範囲によってのみ限定されるものとする。

図１は、変位推定方法を示す図である。

本発明の主な実施形態は、図１に示される。

ＲＦＳｉｇ（ｄ，ｌ，ｆ）は、超音波部から得たＲＦ信号であり、ｄは深度方向（図８の方向８１Ｄ）、ｌはライン方向（方向８１Ｌ）、ｆはフレーム方向（方向８１Ｆ）を表す。

これらに関することは、さらに、図８に示される。

そのようなＲＦ信号の１フレーム内の１ライン（図示される８０１を参照）は、超音波トランスデューサ８００が、当該ラインの位置で走査方向（方向８０２）に向かってパルスを送信するメカニズムによって生成される。

ここで、当該パルスは、反射および散乱しながら、進路に沿って、下層構造と相互作用し、当該超音波トランスデューサ８００は、反射および散乱した信号を受信する。

そして、その信号は、超音波部によって、対応するラインに変換される。

同じメカニズムにより、複数のラインが様々な位置で生成され、複数のフレームが様々なタイムインスタンスに生成される。

図１２の（ａ）欄は、下層構造の変位を示し、（ｂ）欄は、ＲＦ信号への影響を示す。

走査された媒体の下層構造の変位が、ＲＦ信号の遅延として反射されるメカニズムが、図１２により示される。

例示を目的として、１ラインの１パルスが送信されるとし、下層構造には、反射境界が１つだけあるとする。

第１のタイムインスタンス（（ａ）欄の左側部分を参照）において、超音波トランスデューサ１２００は、走査された媒体に向けてパルスを送信する（パルス１２０１を参照）。

このパルスは、下層構造１２０３（１２０３ｘ）の反射境界に到達するまで、媒体を通して伝搬する。

この境界により、パルスはトランスデューサに向かって反射する（反射したパルス１２０２を参照）。

トランスデューサは、このパルスを受信し、ＲＦ信号１２０４（図１２の（ｂ）欄の表の第１行を参照）に変換する。

ＲＦ信号１２０４の全体における、パルスの位置は、送信されたパルス１２０１が、下層構造１２０３に向かって進み、反射して戻ってくるのに必要な時間を示す。

第２のタイムインスタンス（図１２の（ａ）欄の右側を参照）において、ＲＦ信号１２０５（図１２の（ｂ）欄の表の第３行）は、同一メカニズムにより作成されるが、下層構造は、新たな位置に変位している（変位した後の下層構造１２０６を参照）。

変位ｘ（変位１２０７）により、送信されたパルス（パルス１２０１）が、伝搬され、下層構造の境界で反射するのに要する時間が長くなり、その結果、ＲＦ信号の遅延が起こる（遅延時間１２０８を参照）。

なお、この例において、変位ｘの向きは、図１２に示される方向１２０３ｄの向きである。

この遅延時間１２０８の遅延値が分かると、変位（変位１２０７）を推測することができる。

そして、複数の下層構造を有する媒体では、複数の反射パルスが生成される。

それら複数の下層構造が変位する量が互いに異なることが原因で、ＲＦ信号における、反射パルスの遅延は、異なるものになる。

これらのパルスが重複し、かつ相互作用することにより、各下層構造の正確な変位を推定することが難しくなる。

本発明は、反復して推定を行うことにより、下層構造の変位を、正確に推定する方法を提示するものである。

図１に示される主な実施形態によると、本発明は、以下の主なブロックを備える。

つまり、本装置では、ＲＦ信号ワーピング部１００（図５参照）と、変位推定部１０１（図９参照）と、ウインドウ計算部１０２（図２参照）と、収束ガイド部１０３（図７参照）と、出力スイッチ１０４とを備える。

ＲＦ信号ワーピング部１００は、前回の推定ラウンド（第ｎ−１ラウンド）の変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）に基づいて、対象推定ラウンド（第ｎラウンド）においてのワーピングをする。このワーピングでは、対象推定ラウンドについて、ＲＦ信号の、選択されたフレーム内の、選択されたラインをワーピングする。

なお、これは、第１ラウンドでは、ゼロに設定される。

このブロックの出力は、推定ラウンドに用いられる、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）である。

なお、第１の推定ラウンド（つまり、ｒｏｕｎｄ ＝ １）では、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，０）をゼロに初期化したことにより、ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，１）と、超音波部から得たＲＦ信号であるＲＦＳｉｇ（ｄ，ｌ，ｆ） （先述）は、実質的に、互いに同一である。

変位を推定するため、予め定義されたフレームのセットを、このブロックに対する入力として選択することができる。

例示を目的として、２つのフレームが選択され、これらはｆ１およびｆ２で表される。

ＲＦ信号ワーピングの目的は、先行するラウンドで推定された変位に基づいて、１つのＲＦ信号を、他のＲＦ信号と適合するように、修正することである。

そして、その後、修正されたＲＦ信号から、残留変位を推定し、先行するラウンドの変位に追加することで、新たな変位を作成する。

いくつかの数（予め定められた数）の推定ラウンドの後、推定された変位は、修正されたＲＦ信号にほぼ一致するものとなり、残留変位はゼロに収束する。

変位推定部１０１は、ＲＦ信号の時間遅延から導出することで、ＲＦ信号ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）から、変位を推定する。

このタスクの実行に好ましい方法は、各深度における推定ウインドウが、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）によって記述される、自己相関の方法であるが、これに限定されるものではない。

ウインドウ計算部１０２は、各深度での、選択されたラインに用いられる推定ウィンドウパラメータｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）を計算する。

ウインドウ計算部１０２は、ＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を取得し、かつ、対象ラウンド数ｃｕｒＲｏｕｎｄと、ＲＦ信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とを、収束ガイド部１０３からの入力として取得する。

収束ガイド部１０３は、本技術の、収束を導く処理を行う。

つまり、推定された変位の質を決定するために、ワープされたＲＦ信号間の差分が計算される。

この差分が最小である（対象ラウンドを含む）ラウンドが、最も正確な変位推定結果を表していると判定される。

このブロックの出力は、対象推定ラウンド数ｃｕｒＲｏｕｎｄと、ＲＦ信号差分が最小であるラウンド数ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）と、推定された変位の質の測定値としてのＲＦ信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とである。

ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が小さいほど、質がよいことを示し、推定ラウンド全体にわたって、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が減少していると、収束を示す。

なお、ワープされたＲＦ信号間の差分を計算する方法として、「二乗差分和」および「絶対差分」の方法が考えられるが、これらに限定されるものではない。

出力スイッチ１０４は、入力として、ｃｕｒＲｏｕｎｄおよびｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）を取得する。

対象ラウンド数が最大値に到達した場合、最も正確な推定結果が、最終出力として選択される。

そうでなければ、その最も正確な推定結果は、後続の推定ラウンドに用いるために、ＲＦ信号ワーピング部１００の入力として選択され、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）およびｃｕｒＲｏｕｎｄが、入力として変位推定部１０１に送られる。

次の段落に、図１に示される主な実施形態に基づく具体的な実施形態をいくつか示す。

図５は、ＲＦ信号ワーピング方法を示す図である。

本発明の具体的な実施形態として、主な実施形態のＲＦ信号ワーピング部１００は、図５に示される方法によって実現される。

変位の結果、ＲＦ信号は移動され、かつ伸長される。

よって、先行するラウンドの、最も正確な変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））から、ＲＦＳｉｇ（ｄ，ｌ，ｆ１）を、ＲＦＳｉｇ（ｄ，ｌ，ｆ２）に適合させるサンプルにおける遅延値（図１２の遅延時間１２０８を参照）が、遅延値計算部５００において、以下のように算出される。

ここで、ｆｓは、ＲＦ信号のサンプリング周波数であり、ｃは、走査された媒体における音速である。

走査された媒体における変位は、結果として、得られたＲＦ信号における遅延として反映される。

各深度での遅延値が既知であるため、ＲＦＳｉｇ（ｄ，ｌ，ｆ１）またはＲＦＳｉｇ（ｄ，ｌ，ｆ２）の何れかにおける各ラインの、（全ての深度に対応する）全てのサンプルを、他のサンプルに適合するように合わせて、遅延させることができる。

これは、ＲＦサンプル遅延部５０１において実行される。

１つの実施形態において、ＲＦサンプル遅延部５０１における遅延は、信号補間を用いて実行される。

別の実施形態では、ＲＦサンプル遅延部５０１における遅延は、分数遅延フィルタを用いて実行される。

図６は、ＲＦ信号ワーピングおよびＲＦサンプル遅延を示す図である。

ＲＦ信号ワーピング処理が、図６に示される。

変位により、ＲＦ信号の遅延が起こることが、先述の図１２により示される。

ＲＦ信号ワーピングの目的は、推定された変位と逆方向に、１つのＲＦ信号を遅延させることによって、この効果を無効にすることである。

この無効化の結果として、ワープされたＲＦ信号が、相互に適合することが期待される。

遅延値計算部５００から計算されたｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）は、１つのＲＦ信号内の各サンプルが、これを達成するために遅延される値を表す。

２つのＲＦ信号６００および６０１（図６）を、この例示のために用いる。

データ６０２は、ある推定ラウンドの後に、遅延値計算部５００から計算されたｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）を示す。

ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）の特定の深度での値が、ＲＦ信号６００における、同じ深度でのサンプルが、遅延値計算部６０１に適合するための遅延に必要な量を表すことが分かる。

しかしながら、ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）における遅延値は、常に整数であるわけではない。

よって、ｄｅｌａｙＶａｌｕｅｓ（ｄ，ｌ）に従って、ＲＦ信号６００における各サンプルを遅延できるようにするため、信号補間または分数遅延フィルタが、２つの選択候補として挙げられる。

変位推定部１０１をＲＦ信号に適用することで、第１ラウンドに用いられる変位と、後続ラウンドに用いられる残留変位を得て、最終的な変位を得るために、それらを組み合わせる。

好ましい実施形態が図９に示されるが、これに限定されるものではない。

図９は、自己相関を用いた変位推定を表す。

ＩＱ復調部９００は、ＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を、ベースバンド信号ＩＱＳｉｇ（ｄ，ｌ，ｆ）に変換する。

自己相関計算部９０１は、以下の式に従って、ａｕｔｏｃｏｒｒ（ｄ，ｌ）を計算する。

ここで、ｃｏｎｊ（）は、複素共役を表す。

最後に、以下の式に従って、変位計算部９０２により、ａｕｔｏｃｏｒｒ（ｄ，ｌ）が、変位に変換される。

ここで、ａｒｇ（）は、複素数の偏角を計算する関数である。

１から開始するラウンド（第１ラウンドを表す）により、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，０）は、ゼロに初期化される。

図２は、ウインドウ計算方法を示す図である。

ウインドウ計算部１０２により、各推定ラウンドにおける、各ラインの各深度に対し、異なる推定ウインドウを決定する。

これは図２に示される。

図３は、ウインドウサイズ計算の具体例を示す図である。

ウインドウサイズ計算部２００（図２）は、各推定ラウンドに対するウインドウサイズの決定をする。

このブロックの出力は、ｗｉｎＳｉｚｅ（ｄ，ｌ）である。

好ましい実施形態において、各ラウンドに対するウインドウサイズは、ラウンドの関数が減少するにつれて、変化する。

これは図３に示される。

ウインドウサイズ計算部２００（図２）のある実施形態において、図３の（ａ）欄に示されるように、ウインドウサイズは、減少関数ｗｉｎＳｉｚｅＲｏｕｎｄ（ｒｏｕｎｄ）により、推定ラウンド数に、直接関係する。

ｃｕｒＲｏｕｎｄによって特定される対象推定ラウンドについて、以下に基づいて、全ての深度に対するウインドウサイズが特定される。

全てのｄおよびｌについて、
ｗｉｎＳｉｚｅ（ｄ，ｌ） ＝ ｗｉｎＳｉｚｅＲｏｕｎｄ（ｃｕｒＲｏｕｎｄ）

「固定ウィンドウイングスキーム」は、変位推定に用いられる全ての深度で推定ウインドウサイズとなるｗｉｎＳｉｚｅ（ｄ，ｌ）の同じ値を選択する処理である。

「固定ウィンドウイングスキーム」の利点は、必要な演算量が多くなく（比較的少なく）、かつ、ＲＦ信号エネルギーが、深度方向に沿って、規則正しく配分されている場合に、適していることである。

ウインドウサイズ計算部２００の別の実施形態において、各深度でのウインドウサイズは、上限であるｍａｘＷｉｎ（ｒｏｕｎｄ）と、下限であるｍｉｎＷｉｎ（ｒｏｕｎｄ）と、ＲＦ信号電力ｓｉｇＰｏｗ（ｄ，ｌ）とに基づいて計算される。

ｍａｘＷｉｎ（ｒｏｕｎｄ）およびｍｉｎＷｉｎ（ｒｏｕｎｄ）は、ラウンドの減少関数である。

これは、「瞬間電力ウィンドウイングスキーム」として、図３の（ｂ）欄に示される。

以下の条件が適用される。

また、図３の（ｂ）欄に示されるように、ｗｉｎＳｉｚｅ（ｄ，ｌ）は、ｓｉｇＰｏｗ（ｄ，ｌ）、つまり、ｓｉｇＰｏｗ ＝ ０に対応する最大ウインドウサイズと、最大ｓｉｇＰｏｗに対応する最小ウインドウサイズとに反比例してもよい。

「瞬間電力ウィンドウイングスキーム」は、「固定ウィンドウイングスキーム」と比較して、必要な演算量が多いが、特に、ＲＦ信号エネルギーが、深度方向に沿って、かなり不均一である場合に、より精度が高い。

信号対ノイズ比の低さを克服するために、低い信号電力を有する領域に対して、より大きな推定ウインドウが指定され、より詳細に変位を分布するために、大きな信号電力を有する領域に対して、より小さな推定ウインドウが指定される。

しかしながら、信号エネルギーが規則正しく配分されているＲＦ信号に関しては、「固定ウィンドウイングスキーム」と「瞬間電力ウィンドウイングスキーム」とは大差がない。

ウインドウサイズ計算部２００の別の実施形態において、全ての深度でのウインドウサイズを、全てのウインドウ内の信号エネルギーが同量に保たれるように計算してもよく、これは、「一定エネルギーウィンドウイングスキーム」と称される。

信号電力計算部２０１（図２）により、フレーム方向に沿う、ＲＦ信号の平均電力が計算される。

信号電力計算部２０１は、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）と、対象ラウンド数ｃｕｒＲｏｕｎｄとを、入力として取得する。

ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｃｕｒＲｏｕｎｄ）内には、複数のフレームがあり、電力が、以下の通り計算される。

ここで、Ｎは、計算に用いるために選択されたフレームの総数である。

図４は、対称ウインドウおよびエネルギー均等ウインドウの例を示す図である。

ウインドウ境界計算部２０２（図２）により、各ラインの各深度でのウインドウ境界が計算される。

これは図４に示される。

ウインドウ境界計算部２０２のある実施形態においては、ウインドウサイズ計算部２００から計算されたｗｉｎＳｉｚｅ（ｄ，ｌ）が入力であり、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）が、以下の式に従って、対応する深度を中心に、ウインドウが対称になるように計算される（「対称ウインドウ」と称される）。

ここで、ＲＯＵＮＤ（）は、最も近い整数に丸める処理を示す。

これは、図４の（ａ）欄に示される。

対称ウインドウの利点は、必要とされる演算量が多くないことである。

これは、ＲＦ信号電力が、深度方向において、サンプル毎に突然変化しない場合に適している。

ウインドウ境界計算部２０２の別の実施形態においては、次の通りである。

つまり、ウインドウサイズ計算部２００から計算されたｗｉｎＳｉｚｅ（ｄ，ｌ）が、ウインドウ境界計算部２０２の入力である。

そして、ｗｉｎＳｔａｒｔ（ｄ，ｌ）およびｗｉｎＥｎｄ（ｄ，ｌ）が、以下の式により、各深度の両側での信号エネルギーが等しくなるように（「エネルギー均等ウインドウ」と称される）計算される。

これは、図４の（ｂ）欄に示される。

「エネルギー均等ウインドウ」は、「対称ウインドウ」と比較して、必要な演算量が多いが、従来のウィンドウイング技術においてバイアスの原因となっていた、深度方向に沿う、信号電力の不均一な分布の影響を抑制できる。

図７は、収束に導く方法を示す図である。

図１の収束ガイド部１０３による、本発明において説明した方法で、収束が可能であることが保証される。

収束ガイド部１０３は、図７に示されるように、ワープされたＲＦ信号 ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ）を、当該収束ガイド部１０３の入力として取得する。

ＲＦ信号差分計算部７００は、ワープされたＲＦ信号間の差分を計算し、各ラインの各深度について、各推定ラウンドにおける差分を、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）として出力する。

ＲＦ信号ワーピング部１００の主な目的は、ＲＦ信号の１つを、他と適合するように、修正することであるので、この残差（ｄ，ｌ，ｒｏｕｎｄ）は、ＲＦ信号同士が、どの程度適合しているかの尺度になる。

（ｒｏｕｎｄにより特定される）ある特定の推定ラウンドについて、推定された変位ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）の質を評価する目的で、ＲＦ信号ワーピング部１００が、当該ＲＦ信号ワーピング部１００の入力として、ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｒｏｕｎｄ）を取得する。

このことにより、ＲＦ信号がワープされ、適合レベルを評価するために、ワープされた信号間の差分が計算される。

これらのワープされた信号は、後続ラウンドの変位推定のための入力としても用いられ、よってＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ＋１）と表される。

よって、ＲＦＳｉｇＷａｒｐ（ｄ，ｌ，ｆ，ｒｏｕｎｄ＋１）から、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）が計算される。

ＲＦ信号差分計算部７００（図７）で行われる計算の一例として、以下の式（つまり、二乗差分和）が挙げられるが、これに限定されるものではない。

ここで、ｆ１およびｆ２は、推定に用いるために選択された２フレームを表す。

この例において、１ラインの全ての深度に対し同一の信号差分値が指定され、これは、対応するラインのＲＦ信号が、全体的に適合していることを表す。

小さな領域で変位が起こるだけではないときには、「二乗差分和」が適しており、その結果、ライン全体に反復が起こる。

なお、ＲＦ信号差分計算部７００（図７）で行われる計算の別の例として、以下の式（つまり、絶対差分）が挙げられるが、これに限定されるものではない。

ここで、ｆ１およびｆ２は、推定に用いるために選択された２フレームを表す。

この例において、各ラインの各深度は、個別の信号差分値を有し、これは、ＲＦ信号が、局所的に適合していることを表す。

深度方向に小さな領域だけで、変位が起こる場合には、「絶対差分」が適しており、その結果、その領域内での局所的反復が可能になる。

最小残余計算部７０１は、以下の式に基づき、ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）における、各深度およびラインについて、残余が最小であるラウンド数を確認する。

これは、最も正確な推定結果を表す。

ここで関数

は、ｆ（ｘ）を最小化するｘの値を返す。

収束ガイド部１０３（図１）の出力は、対象ラウンド数ｃｕｒＲｏｕｎｄと、残余が最小であるラウンド数ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）と、信号差分ｒｅｓｉｄｕｅ（ｄ，ｌ，ｒｏｕｎｄ）とである。

出力スイッチ１０４（図１）は、入力として、ｃｕｒＲｏｕｎｄおよびｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ）を取得する。

ｃｕｒＲｏｕｎｄが最大値に達すると、出力スイッチ１０４は、最も精度の高い変位推定結果ｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））を、最終的な出力ｄｉｓｐＯｕｔとして選択する。

そうでなければ、対象ラウンドｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｃｕｒＲｏｕｎｄ）が、変位推定される代わりに、次の推定ラウンドのための、ＲＦ信号ワーピング部１００の入力として用いられるｄｉｓｐＲｏｕｎｄ（ｄ，ｌ，ｍｉｎＲｅｓｉｄｕｅＲｏｕｎｄ（ｄ，ｌ））が、変位推定される。

ここで、変位推定の多次元（２次元 or ３次元）への拡張について説明する。

つまり、上記の実施の形態では、超音波の深さ方向における１次元の変位を求める方法について説明した。

しかしながら、血管や腫瘍組織などの生体組織の拍動方向は、１次元内に限定されず、また、場所によっても、方向が異なる。

したがって、複数方向から超音波を送受信することで、特定位置における複数方向の変位を求め、組織の２次元あるいは３次元的な変位を測定してもよい。

こうして、多次元で測定することにより、診断精度の更なる向上が見込まれる。

そして、超音波を複数方向から送受信する方法としては、超音波振動子が、１次元に配置されたリニアプローブを用いて、プローブの方向を変化させる、または、超音波振動子が２次元アレイ上に配置されたマトリックスプローブを用いて、ビームフォーミングにより、複数方向での送受信を行うなどの方法が可能である。

ここで、リニアプローブの使用時には、磁気センサや加速度センサ、ジャイロなどの位置センサ、あるいは、カメラにより、プローブの位置を追跡することで、プローブの位置や向きを取得して、送受信時の方向を決定できる。

なお、変位測定は、超音波の深さ方向など、１方向についてのみ行い、この変位ベクトルを、所定の２次元あるいは３次元座標系の各軸方向の成分に分解することで、多次元的に評価してもよい。

この時、超音波の深さ方向と所定の座標系との対応関係が必要となるが、上述の位置センサなどを用いて、プローブの位置と方向を取得することで、両者を対応付ける。

ここで、癌検出、あるいは、動脈硬化（血管の硬さ）の診断などに適用できる歪み計測（エラストグラフィ）との関連について補足説明する。歪みは、変位の時間微分であり、組織の硬さの指標となるが、本発明により、変位推定精度が向上することで、結果として、歪みの測定精度が向上する。エラストグラフィでは微分値を用いるが、その他に、組織の特定部位の変位の時間変化波形なども、組織性状の判定に利用できる。例えば、癌の成長過程においては、癌の周辺および内部に血管新生が起こることが知られている。従って、これら血管の拍動、あるいは、血管の拍動に伴う周囲の組織の拍動などを変位量の時間変化として捉え、変化のパターンや振幅が癌に特有のものであるかどうかを判定することで、癌検出を行ってもよい。

こうして、変位推定装置であって、超音波信号を用いて、変位を反復的に推定する変位推定装置であって、少なくとも１つの超音波信号を、トランスデューサから走査された媒体に送信する手段と、前記走査された媒体から反射した超音波信号を受信する手段と、推定ウインドウサイズを計算する手段と、前記推定されたウインドウサイズに基づいて、前記推定ウインドウの境界を計算する手段と、前記計算された推定ウインドウを用いて、前記超音波信号の深度方向に沿って変位を推定する手段と、前記推定された変位に基づいて、前記超音波信号をワーピングする手段と、前記方法の収束傾向を前記ワープされた超音波信号間の差分として計算することで、収束を導く手段とを備える変位推定装置が開示される。

こうして、開示される方法は、収束を導く方法、精度を限定する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて変位推定の精度を保証する。

つまり、次の問題が解決される。つまり、超音波を用いる変位推定において、先行技術では、推定の精度と解像度を高める方法を提供している。その中には、複数段階からなる推定または反復推定を導入しているものもある。しかし、（１）反復推定の収束を導く方法を組み込んだ包括的な反復手法の欠如、（２）品質評価方法の欠如、および（３）変位の精度を限定する要素を克服する方法の欠如という問題があり、本発明の変位推定方法により、これらの問題が解決される。

なお、本発明は、装置、システム、集積回路などとして実現できるだけでなく、その装置等を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを記録したコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体として実現したり、そのプログラムを示す情報、データ又は信号として実現したりすることもできる。そして、それらプログラム、情報、データ及び信号は、インターネット等の通信ネットワークを介して配信してもよい。

なお、本発明について、実施の形態に基づいて説明したが、本発明は、上記の実施の形態に限定されるものではない。本発明の主旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の範囲内に含まれる。

本技術は、超音波信号を用いて変位を推定する方法を導入するものである。これは、最終的な結果、または、更なる処理のための中間ステップの何れであっても、変位推定が必要な適用例において、用いることができ、医療用および工業用超音波機器に用いることができる。

変位を精度よく推定できる変位推定装置が提供でき、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供可能な変位推定装置が提供できる。

１Ｘ 変位推定装置
１Ｘ１ 送信部
１Ｘ２ 受信部
１Ｘ３ サイズ計算部
１Ｘ４ 境界計算部
１Ｘ５ 推定部
１Ｘ６ ワーピング部
１Ｘ７ 収束制御部
１２０１ｓ 超音波信号
１２０３ｘ 下層構造
１Ｘ３Ｍ ウインドウサイズ
１Ｘ４Ｍ 境界の情報
１ＸｂＭ 変位
１Ｘ７Ｍ 相関値
１００ ＲＦ信号ワーピング部
１０３ 収束ガイド部
１０１ 変位推定部
１０２ ウインドウ計算部
１０４ 出力スイッチ
６００、１２０４、１２０５ ＲＦ信号
８００、１２００ 超音波トランスデューサ
１２００Ｍ 媒体
１２０８ 遅延時間
１００ａ、１０１ａ、１０２１ａ、１０３ａ、１０４ａ、１０４ｂ、１０４ｃ 情報
２０１ｎ、２０１ｍ、２０１ａ、２００ｍ、２００ａ、２０２ａ 情報
５００ｍ、５００ａ、５０１ｎ、５０１ａ 情報
７００ｍ、７００ａ、７０１ａ、７００ｂ 情報
９００ｍ、９００ａ、９０１ｍ、９０１ａ、９０２ａ 情報

Claims (10)

1. 超音波信号を用いて、変位を反復的に推定する変位推定方法であって、
   少なくとも１つの超音波信号を走査して、媒体に送信するステップと、
   走査された前記媒体から反射した前記超音波信号を受信するステップと、
   ウインドウサイズを計算するステップと、
   計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、
   計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、
   推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、
   ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む
   変位推定方法。
2. 計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される
   請求項１記載の変位推定方法。
3. 連続する複数の前記ラウンドの前記ウインドウサイズの上限と下限とのそれぞれを、徐々に減少するように変化させるステップと、
   各深度での信号電力に基づいて、それぞれの深度における前記ウインドウサイズとして、その深度における上限と下限との間のウインドウサイズを計算するステップとを含む
   請求項１記載の変位推定方法。
4. 前記ウインドウサイズは、全てのウインドウにおける信号エネルギーが等しくなるように計算される
   請求項１記載の変位推定方法。
5. 超音波信号の相関値が大きくなるほど、より小さい前記ウインドウサイズが計算される
   請求項１記載の変位推定方法。
6. 前記境界は、当該境界による前記ウインドウが、対応する深度を中心に、対称のウインドウになるように、当該深度の両側に延びる
   請求項１記載の変位推定方法。
7. 前記境界は、対応する深度での、当該境界による前記ウインドウの両側における信号エネルギーが互いに等しくなるように、対応する前記深度の両側に延びる
   請求項１記載の変位推定方法。
8. 前記超音波信号を、計算された遅延値により、前記超音波信号内の各サンプルを遅延させることにより、ワーピングするステップを含む
   請求項１記載の変位推定方法。
9. 各ラウンドの後に、ワーピングされたＲＦ信号間の信号差分を計算するステップと、
   計算された前記信号差分が最小になるラウンドを決定するステップとを含む
   請求項１記載の変位推定方法。
10. 超音波信号を用いて、変位を反復的に推定する変位推定装置であって、
    少なくとも１つの超音波信号を走査して、媒体に送信する送信部と、
    走査された前記媒体から反射した前記超音波信号を受信する受信部と、
    ウインドウサイズを計算するサイズ計算部と、
    計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算する境界計算部と、
    計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定する推定部と、
    推定された前記変位に基づいて、前記超音波信号をワーピングするワーピング部と、
    ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、当該変位推定装置が実行する変位推定方法の収束を導く収束制御部とを含む
    変位推定装置。

Images