"超音波信号を用いて、変位を反復的に推定する変位推定方法であって、少なくとも1つの超音波信号を走査して、媒体に送信するステップと、走査された前記媒体から反射した前記超音波信号を受信するステップと、ウインドウサイズを計算するステップと、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む変位推定方法。計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される請求項1記載の変位推定方法。連続する複数の前記ラウンドの前記ウインドウサイズの上限と下限とのそれぞれを、徐々に減少するように変化させるステップと、各深度での信号電力に基づいて、それぞれの深度における前記ウインドウサイズとして、その深度における上限と下限との間のウインドウサイズを計算するステップとを含む請求項1記載の変位推定方法。前記ウインドウサイズは、全てのウインドウにおける信号エネルギーが等しくなるように計算される請求項1記載の変位推定方法。超音波信号の相関値が大きくなるほど、より小さい前記ウインドウサイズが計算される請求項1記載の変位推定方法。前記境界は、当該境界による前記ウインドウが、対応する深度を中心に、対称のウインドウになるように、当該深度の両側に延びる請求項1記載の変位推定方法。前記境界は、対応する深度での、当該境界による前記ウインドウの両側における信号エネルギーが互いに等しくなるように、対応する前記深度の両側に延びる請求項1記載の変位推定方法。前記超音波信号を、計算された遅延値により、前記超音波信号内の各サンプルを遅延させることにより、ワーピングするステップを含む請求項1記載の変位推定方法。各ラウンドの後に、ワーピングされたRF信号間の信号差分を計算するステップと、計算された前記信号差分が最小になるラウンドを決定するステップとを含む請求項1記載の変位推定方法。超音波信号を用いて、変位を反復的に推定する変位推定装置であって、少なくとも1つの超音波信号を走査して、媒体に送信する送信部と、走査された前記媒体から反射した前記超音波信号を受信する受信部と、ウインドウサイズを計算するサイズ計算部と、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算する境界計算部と、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定する推定部と、推定された前記変位に基づいて、前記超音波信号をワーピングするワーピング部と、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、当該変位推定装置が実行する変位推定方法の収束を導く収束制御部とを含む変位推定装置。単一段階の変位推定方法に関し、現在の方法は、精度および変位解像度の点で課題を有する。相互相関に基づく方法は、信号補間なしでは、サブサンプリング間隔に対応する小さな変位を推定することができない。しかしながら、この方法では、ノイズの影響が少ない。よって、相互相関ベースの方法では、精度は高いが、変位解像度は低い。相互相関方法の変位解像度を向上させる目的で、信号補間を適用することができる。しかし、これは、推定対象の変位が、マイクロメートル単位の場合に、処理電力を大幅に増大させる。また、この方法の精度は、補間アルゴリズムの質に依存する。自己相関に基づく方法では、小さな変位を推定することができるが、非常に、ノイズを発生しやすい。よって、この方法では、変位解像度は高いが、精度は低い。よって、精度を向上させるため、2段階の推定方法が利用できる。しかし、これらの方法は、包括的な反復手法の欠如、変位の質を評価する方法の欠如、変位の精度を制限する要素、つまり、ノイズおよび振幅変調効果を解消する方法の欠如という3つの観点で、課題を有する。包括的な反復手法には、反復推定の収束を導く技術を組み込む必要がある。本発明の目的は、組織変位を精度よく推定できる変位推定装置変位推定方法を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。上記の課題を解決するために、本発明の変位推定方法は、超音波信号を用いて、変位を反復的に推定する変位推定方法であって、少なくとも1つの超音波信号を走査して、媒体に送信するステップと、走査された前記媒体から反射した前記超音波信号を受信するステップと、ウインドウサイズを計算するステップと、計算された前記ウインドウサイズに基づいて、ウインドウの境界を計算するステップと、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位を推定するステップと、推定された前記変位に基づいて、前記超音波信号をワーピングするステップと、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値が大きくなるように、この変位推定方法の収束を導くステップとを含む変位推定方法である。本発明において、発明者は、反復推定手法を用いる超音波を用いて、変位を推定する方法を教示する。発明者の実験によると、推定ウインドウサイズが同じままか、または、全ての反復推定ラウンドに対して、適切に構成されていない場合には、推定された変位は、収束しない。よって、この収束に役立つウインドウ計算方法を導入する。ウインドウ計算方法は、反復推定の収束を導くために、ウインドウサイズを計算することの部分と、信号電力の不均一な分布によって起こるノイズおよびバイアスの影響を克服するために、ウインドウ境界を計算することの部分との2つの部分を含む。ウインドウサイズの計算により、反復推定ラウンド毎に、異なるウインドウサイズを決定する。これは、後続のラウンドにおいて、より詳細な推定結果を得るためと、推定を確実に収束させるためという2つの目的を有する。つまり、不適切なウインドウサイズでの処理がされて、収束が得られず、正しい変位量が得られないことが回避されて、適切なウインドウサイズでの処理がされて、確実に収束が得られ、より正確な変位量が得られる。ウインドウ境界の計算では、信号電力を利用し、各推定位置での、推定ウインドウの境界を決定する。従来の推定方法において、ウインドウは通常、その位置を中心に対称である。これは、不均一な信号電力の分布によって、バイアスが起こる一因である。信号電力を用いて、ウインドウ境界を決定することにより、そのような影響を克服できる。さらなる収束に導くため、推定された変位の質の評価を導入し、反復推定方法が、常に最良の結果に収束するよう保証する。本発明は、収束を導く方法、精度を制限する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて、変位推定の精度を保証するものである。図10は、シミュレーションに、バックトラックアルゴリズムを適用した後の、変位推定の改善結果を示す図である。図11は、ファントム実験における、同じ結果を示す図である。図10は、推定の第1ラウンドにおける、推定された変位変位1001が、シミュレーションされた変位プロファイルプロファイル1000から逸脱していることを示す図である。しかしながら、収束を導くウィンドウイング方式を用いるバックトラックアルゴリズムの最終ラウンドでは、推定された変位変位1002は、シミュレーションされたプロファイルと、非常によく適合する。図11は、グローバル変位が生成されるファントム実験による推定結果を示す図である。例示のため、単線が選択されるプロファイル1100を参照。グローバル変位が生成されるので、一定の変位プロファイルプロファイル1101が期待される。推定の第1ラウンドでは、結果データ1102を参照は、期待される変位に、あまり適合しない。しかしながら、最終ラウンドの推定結果データ1103を参照は、第1ラウンドの推定結果よりも、グローバルな動きを、ずっとよく表している。また、適切なウインドウサイズでの処理ができる変位推定装置を提供できる。本発明の目的は、組織変位を精度よく推定できる変位推定装置変位推定方法を提供することを含み、ひいては、推定された変位に基づいて、悪性腫瘍と良性腫瘍、あるいは正常組織を判別するために適した情報を提供する変位推定装置を提供することを含む。なお、本発明の他の目的は、反復的な変位推定処理を、少ない反復回数で高精度に収束させることができる変位推定装置を提供することを含む。以下、図面を参照して詳しく説明される。実施形態の変位推定方法は、超音波信号図13の超音波信号1201sを用いて、変位図12の変位1207を参照を反復的に推定する変位推定方法図13の変位推定装置1X参照である。つまり、例えば、変位の大きさの特定が複数回行われて、特定の処理が複数回行われることにより、比較的精度のよい、変位の大きさが特定される。そして、少なくとも1つの超音波信号を走査して、媒体媒体1200M図12、および、媒体1200Mにある、下層構造変位の測定の対象1203x図12、13を参照に送信する図13の超音波処理部1Xaにある送信部1X1を参照。つまり、例えば、ライン方向図8の方向81Lの複数の位置へと超音波信号が送信されて、それら複数の位置の走査が行われる。そして、走査された前記媒体から反射した前記超音波信号超音波信号1201sを受信する受信部1X2を参照。なお、超音波処理部1Xaは、例えば、プローブ探触子などである。そして、ウインドウサイズ図13のウインドウサイズ1X3M:例えば、図3のa欄、b欄のウインドウサイズwinSizeなどを計算する図13の情報処理部1Xbに含まれるサイズ計算部1X3を参照。つまり、それぞれのラウンドにおけるウインドウサイズが決定される。なお、サイズ計算部1X3は、例えば、ウインドウ計算部102図1を含んでもよいし、ウインドウ計算部102の少なくとも一部などでもよい。すなわち、例えば、第n-1のラウンドでの相関値図13の相関値1X7Mに基づいて、第nのラウンドでのウインドウサイズが決定されるn≧2。そして、計算された前記ウインドウサイズに基づいて、ウインドウの境界境界の情報1X4M:図4のa欄、b欄に示される、ウインドウの始端winStart、および、終端winEndなどを参照を計算する境界計算部1X4参照。なお、境界計算部1X4は、例えば、ウインドウ計算部102図1の少なくとも一部等でもよい。つまり、それら始端および終端の間の幅が、計算されたウインドウサイズである始端および終端が特定される。そして、計算された前記境界による前記ウインドウを用いて、前記超音波信号のそれぞれの深度における変位変位1XbMを推定する推定部1X5参照。なお、推定部1X5は、例えば、例えば、変位推定部101図1の少なくとも一部等でもよい。つまり、例えば、本装置からの距離が互いに異なる複数の位置における変位が特定される。そして、推定された前記変位に基づいて、前記超音波信号をワーピングするワーピング部1X6参照。なお、ワーピング部1X6は、例えば、RF信号ワーピング部100図1の少なくとも一部等でもよい。そして、ワーピングされた前記超音波信号を用いて、前記超音波信号の相関値相関値1X7Mが大きくなるように、この変位推定方法の収束を導く収束制御部1X7参照。なお、収束制御部1X7は、例えば、収束ガイド部103の少なくとも一部等でもよい。つまり、ワーピングされた超音波信号から算出される相関値が最大であるラウンドで推定された、変位の大きさが、最も精度の高い大きさとして特定される。そして、計算される前記ウインドウサイズは、連続する複数のラウンドに用いられる前記ウインドウサイズが、徐々に減少するように、変化される。つまり、例えば、それぞれのラウンドで計算されるウインドウサイズは、そのラウンドよりも前のラウンドで計算されるウインドウサイズよりも小さいサイズであり、計算されるウインドウサイズが、小さいサイズへと変更される。そして、前記ウインドウサイズは、全ての深度のウインドウにおける信号エネルギーが等しくなるように計算される。そして、超音波信号の相関値第n-1ラウンドでの相関値が大きくなるほど、より小さい前記ウインドウサイズ第nラウンドでのウインドウサイズが計算される。ここで、相関値など各ラウンドにおける収束度合いの評価値とウインドウサイズとの関係は、生体での実験値に基づく所定の関係式、あるいは、両者の対応関係を示すテーブルなどを参照して決定できる。そして、第n-1ラウンドでの相関値に対して、決定された対応関係が対応させるウインドウサイズが、第nラウンドでのウインドウサイズと特定される。つまり、例えば、本装置において、次の動作が行われる。すなわち、測定対象物例えば新生血管を伴う悪性腫瘍などである下層構造1203x図12へと、超音波信号1201sが送信され、送信された超音波信号1201sが受信されてもよい。特定される遅延時間1208は、第1のパルスの位置図12のb欄のRF信号1204を参照と、第2のパルスの位置RF信号1205を参照との間における遅延時間である。ここで、第1のパルスは、受信された超音波信号1201sにおける、下層構造1203xが変位1207をする前のパルスである。そして、第2のパルスは、測定対象物が変位1207をした後のパルスである。そして、特定された遅延時間1208に基づいて変位1207を決定する。そして、特定された、変位1207の大きさでの、下層構造1203xの動きが、予め定められた動きであるか否かが判定されてもよい。そして、予め定められた動きと判定される場合に、当該下層構造1203xが、悪性腫瘍癌と判定され、当該動きではないと判定される場合に、良性腫瘍、あるいは、正常組織と判定されてもよい。なお、このような、癌か否かの判定は、例えば、収束制御部1x7図13により行われてもよい。このような本装置において、より具体的には、例えば、後で詳しく説明されるように、第n-1のラウンドで特定された相関値が大きいほど、第nのラウンドにおいて、より大きいウインドウサイズでの処理がされてもよい。なお、ここで、ウインドウサイズによる処理は、例えば、第nのラウンドにおける相関値が、受信される超音波信号1201sにおける、当該ウインドウサイズのウインドウの部分のデータから算出される処理でもよい。本処理により、各ラウンドにおいて超音波の受信信号の信号レベルに応じた適切なウインドウが設定できるため、収束の精度が高まると共に、少ない反復回数で収束が得られ、高精度な変位推定が少ない処理量で実現できる。以下の実施形態は、単に、様々な発明ステップの原理を説明するものである。ここに説明する具体例の様々な変形は、当業者には明らかであろう。したがって、本発明の技術範囲は、本明細書に記載の、具体的かつ例示的な内容によってではなく、特許請求項の範囲によってのみ限定されるものとする。図1は、変位推定方法を示す図である。本発明の主な実施形態は、図1に示される。RFSigd,l,fは、超音波部から得たRF信号であり、dは深度方向図8の方向81D、lはライン方向方向81L、fはフレーム方向方向81Fを表す。これらに関することは、さらに、図8に示される。そのようなRF信号の1フレーム内の1ライン図示される801を参照は、超音波トランスデューサ800が、当該ラインの位置で走査方向方向802に向かってパルスを送信するメカニズムによって生成される。ここで、当該パルスは、反射および散乱しながら、進路に沿って、下層構造と相互作用し、当該超音波トランスデューサ800は、反射および散乱した信号を受信する。そして、その信号は、超音波部によって、対応するラインに変換される。同じメカニズムにより、複数のラインが様々な位置で生成され、複数のフレームが様々なタイムインスタンスに生成される。図12のa欄は、下層構造の変位を示し、b欄は、RF信号への影響を示す。走査された媒体の下層構造の変位が、RF信号の遅延として反射されるメカニズムが、図12により示される。例示を目的として、1ラインの1パルスが送信されるとし、下層構造には、反射境界が1つだけあるとする。第1のタイムインスタンスa欄の左側部分を参照において、超音波トランスデューサ1200は、走査された媒体に向けてパルスを送信するパルス1201を参照。このパルスは、下層構造12031203xの反射境界に到達するまで、媒体を通して伝搬する。この境界により、パルスはトランスデューサに向かって反射する反射したパルス1202を参照。トランスデューサは、このパルスを受信し、RF信号1204図12のb欄の表の第1行を参照に変換する。RF信号1204の全体における、パルスの位置は、送信されたパルス1201が、下層構造1203に向かって進み、反射して戻ってくるのに必要な時間を示す。第2のタイムインスタンス図12のa欄の右側を参照において、RF信号1205図12のb欄の表の第3行は、同一メカニズムにより作成されるが、下層構造は、新たな位置に変位している変位した後の下層構造1206を参照。変位x変位1207により、送信されたパルスパルス1201が、伝搬され、下層構造の境界で反射するのに要する時間が長くなり、その結果、RF信号の遅延が起こる遅延時間1208を参照。なお、この例において、変位xの向きは、図12に示される方向1203dの向きである。この遅延時間1208の遅延値が分かると、変位変位1207を推測することができる。そして、複数の下層構造を有する媒体では、複数の反射パルスが生成される。それら複数の下層構造が変位する量が互いに異なることが原因で、RF信号における、反射パルスの遅延は、異なるものになる。これらのパルスが重複し、かつ相互作用することにより、各下層構造の正確な変位を推定することが難しくなる。本発明は、反復して推定を行うことにより、下層構造の変位を、正確に推定する方法を提示するものである。図1に示される主な実施形態によると、本発明は、以下の主なブロックを備える。つまり、本装置では、RF信号ワーピング部100図5参照と、変位推定部101図9参照と、ウインドウ計算部102図2参照と、収束ガイド部103図7参照と、出力スイッチ104とを備える。RF信号ワーピング部100は、前回の推定ラウンド第n-1ラウンドの変位推定結果dispRoundd,l,roundに基づいて、対象推定ラウンド第nラウンドにおいてのワーピングをする。このワーピングでは、対象推定ラウンドについて、RF信号の、選択されたフレーム内の、選択されたラインをワーピングする。なお、これは、第1ラウンドでは、ゼロに設定される。このブロックの出力は、推定ラウンドに用いられる、ワープされたRF信号RFSigWarpd,l,f,roundである。なお、第1の推定ラウンドつまり、round=1では、dispRoundd,l,0をゼロに初期化したことにより、RFSigWarpd,l,f,1と、超音波部から得たRF信号であるRFSigd,l,f先述は、実質的に、互いに同一である。変位を推定するため、予め定義されたフレームのセットを、このブロックに対する入力として選択することができる。例示を目的として、2つのフレームが選択され、これらはf1およびf2で表される。RF信号ワーピングの目的は、先行するラウンドで推定された変位に基づいて、1つのRF信号を、他のRF信号と適合するように、修正することである。そして、その後、修正されたRF信号から、残留変位を推定し、先行するラウンドの変位に追加することで、新たな変位を作成する。いくつかの数予め定められた数の推定ラウンドの後、推定された変位は、修正されたRF信号にほぼ一致するものとなり、残留変位はゼロに収束する。変位推定部101は、RF信号の時間遅延から導出することで、RF信号RFSigWarpd,l,f,roundから、変位を推定する。このタスクの実行に好ましい方法は、各深度における推定ウインドウが、winStartd,lおよびwinEndd,lによって記述される、自己相関の方法であるが、これに限定されるものではない。ウインドウ計算部102は、各深度での、選択されたラインに用いられる推定ウィンドウパラメータwinStartd,lおよびwinEndd,lを計算する。ウインドウ計算部102は、RF信号RFSigWarpd,l,f,roundを取得し、かつ、対象ラウンド数curRoundと、RF信号差分residued,l,roundとを、収束ガイド部103からの入力として取得する。収束ガイド部103は、本技術の、収束を導く処理を行う。つまり、推定された変位の質を決定するために、ワープされたRF信号間の差分が計算される。この差分が最小である対象ラウンドを含むラウンドが、最も正確な変位推定結果を表していると判定される。このブロックの出力は、対象推定ラウンド数curRoundと、RF信号差分が最小であるラウンド数minResidueRoundd,lと、推定された変位の質の測定値としてのRF信号差分residued,l,roundとである。residued,l,roundが小さいほど、質がよいことを示し、推定ラウンド全体にわたって、residued,l,roundが減少していると、収束を示す。なお、ワープされたRF信号間の差分を計算する方法として、二乗差分和および絶対差分の方法が考えられるが、これらに限定されるものではない。出力スイッチ104は、入力として、curRoundおよびminResidueRoundd,lを取得する。対象ラウンド数が最大値に到達した場合、最も正確な推定結果が、最終出力として選択される。そうでなければ、その最も正確な推定結果は、後続の推定ラウンドに用いるために、RF信号ワーピング部100の入力として選択され、residued,l,roundおよびcurRoundが、入力として変位推定部101に送られる。次の段落に、図1に示される主な実施形態に基づく具体的な実施形態をいくつか示す。図5は、RF信号ワーピング方法を示す図である。本発明の具体的な実施形態として、主な実施形態のRF信号ワーピング部100は、図5に示される方法によって実現される。変位の結果、RF信号は移動され、かつ伸長される。よって、先行するラウンドの、最も正確な変位推定結果dispRoundd,l,minResidueRoundd,lから、RFSigd,l,f1を、RFSigd,l,f2に適合させるサンプルにおける遅延値図12の遅延時間1208を参照が、遅延値計算部500において、以下のように算出される。ここで、fsは、RF信号のサンプリング周波数であり、cは、走査された媒体における音速である。走査された媒体における変位は、結果として、得られたRF信号における遅延として反映される。各深度での遅延値が既知であるため、RFSigd,l,f1またはRFSigd,l,f2の何れかにおける各ラインの、全ての深度に対応する全てのサンプルを、他のサンプルに適合するように合わせて、遅延させることができる。これは、RFサンプル遅延部501において実行される。1つの実施形態において、RFサンプル遅延部501における遅延は、信号補間を用いて実行される。別の実施形態では、RFサンプル遅延部501における遅延は、分数遅延フィルタを用いて実行される。図6は、RF信号ワーピングおよびRFサンプル遅延を示す図である。RF信号ワーピング処理が、図6に示される。変位により、RF信号の遅延が起こることが、先述の図12により示される。RF信号ワーピングの目的は、推定された変位と逆方向に、1つのRF信号を遅延させることによって、この効果を無効にすることである。この無効化の結果として、ワープされたRF信号が、相互に適合することが期待される。遅延値計算部500から計算されたdelayValuesd,lは、1つのRF信号内の各サンプルが、これを達成するために遅延される値を表す。2つのRF信号600および601図6を、この例示のために用いる。データ602は、ある推定ラウンドの後に、遅延値計算部500から計算されたdelayValuesd,lを示す。delayValuesd,lの特定の深度での値が、RF信号600における、同じ深度でのサンプルが、遅延値計算部601に適合するための遅延に必要な量を表すことが分かる。しかしながら、delayValuesd,lにおける遅延値は、常に整数であるわけではない。よって、delayValuesd,lに従って、RF信号600における各サンプルを遅延できるようにするため、信号補間または分数遅延フィルタが、2つの選択候補として挙げられる。変位推定部101をRF信号に適用することで、第1ラウンドに用いられる変位と、後続ラウンドに用いられる残留変位を得て、最終的な変位を得るために、それらを組み合わせる。好ましい実施形態が図9に示されるが、これに限定されるものではない。図9は、自己相関を用いた変位推定を表す。IQ復調部900は、RF信号RFSigWarpd,l,f,roundを、ベースバンド信号IQSigd,l,fに変換する。自己相関計算部901は、以下の式に従って、autocorrd,lを計算する。ここで、conjは、複素共役を表す。最後に、以下の式に従って、変位計算部902により、autocorrd,lが、変位に変換される。ここで、argは、複素数の偏角を計算する関数である。1から開始するラウンド第1ラウンドを表すにより、dispRoundd,l,0は、ゼロに初期化される。図2は、ウインドウ計算方法を示す図である。ウインドウ計算部102により、各推定ラウンドにおける、各ラインの各深度に対し、異なる推定ウインドウを決定する。これは図2に示される。図3は、ウインドウサイズ計算の具体例を示す図である。ウインドウサイズ計算部200図2は、各推定ラウンドに対するウインドウサイズの決定をする。このブロックの出力は、winSized,lである。好ましい実施形態において、各ラウンドに対するウインドウサイズは、ラウンドの関数が減少するにつれて、変化する。これは図3に示される。ウインドウサイズ計算部200図2のある実施形態において、図3のa欄に示されるように、ウインドウサイズは、減少関数winSizeRoundroundにより、推定ラウンド数に、直接関係する。curRoundによって特定される対象推定ラウンドについて、以下に基づいて、全ての深度に対するウインドウサイズが特定される。全てのdおよびlについて、winSized,l=winSizeRoundcurRound固定ウィンドウイングスキームは、変位推定に用いられる全ての深度で推定ウインドウサイズとなるwinSized,lの同じ値を選択する処理である。固定ウィンドウイングスキームの利点は、必要な演算量が多くなく比較的少なく、かつ、RF信号エネルギーが、深度方向に沿って、規則正しく配分されている場合に、適していることである。ウインドウサイズ計算部200の別の実施形態において、各深度でのウインドウサイズは、上限であるmaxWinroundと、下限であるminWinroundと、RF信号電力sigPowd,lとに基づいて計算される。maxWinroundおよびminWinroundは、ラウンドの減少関数である。これは、瞬間電力ウィンドウイングスキームとして、図3のb欄に示される。以下の条件が適用される。また、図3のb欄に示されるように、winSized,lは、sigPowd,l、つまり、sigPow=0に対応する最大ウインドウサイズと、最大sigPowに対応する最小ウインドウサイズとに反比例してもよい。瞬間電力ウィンドウイングスキームは、固定ウィンドウイングスキームと比較して、必要な演算量が多いが、特に、RF信号エネルギーが、深度方向に沿って、かなり不均一である場合に、より精度が高い。信号対ノイズ比の低さを克服するために、低い信号電力を有する領域に対して、より大きな推定ウインドウが指定され、より詳細に変位を分布するために、大きな信号電力を有する領域に対して、より小さな推定ウインドウが指定される。しかしながら、信号エネルギーが規則正しく配分されているRF信号に関しては、固定ウィンドウイングスキームと瞬間電力ウィンドウイングスキームとは大差がない。ウインドウサイズ計算部200の別の実施形態において、全ての深度でのウインドウサイズを、全てのウインドウ内の信号エネルギーが同量に保たれるように計算してもよく、これは、一定エネルギーウィンドウイングスキームと称される。信号電力計算部201図2により、フレーム方向に沿う、RF信号の平均電力が計算される。信号電力計算部201は、ワープされたRF信号RFSigWarpd,l,f,roundと、対象ラウンド数curRoundとを、入力として取得する。RFSigWarpd,l,f,curRound内には、複数のフレームがあり、電力が、以下の通り計算される。ここで、Nは、計算に用いるために選択されたフレームの総数である。図4は、対称ウインドウおよびエネルギー均等ウインドウの例を示す図である。ウインドウ境界計算部202図2により、各ラインの各深度でのウインドウ境界が計算される。これは図4に示される。ウインドウ境界計算部202のある実施形態においては、ウインドウサイズ計算部200から計算されたwinSized,lが入力であり、winStartd,lおよびwinEndd,lが、以下の式に従って、対応する深度を中心に、ウインドウが対称になるように計算される対称ウインドウと称される。ここで、ROUNDは、最も近い整数に丸める処理を示す。これは、図4のa欄に示される。対称ウインドウの利点は、必要とされる演算量が多くないことである。これは、RF信号電力が、深度方向において、サンプル毎に突然変化しない場合に適している。ウインドウ境界計算部202の別の実施形態においては、次の通りである。つまり、ウインドウサイズ計算部200から計算されたwinSized,lが、ウインドウ境界計算部202の入力である。そして、winStartd,lおよびwinEndd,lが、以下の式により、各深度の両側での信号エネルギーが等しくなるようにエネルギー均等ウインドウと称される計算される。これは、図4のb欄に示される。エネルギー均等ウインドウは、対称ウインドウと比較して、必要な演算量が多いが、従来のウィンドウイング技術においてバイアスの原因となっていた、深度方向に沿う、信号電力の不均一な分布の影響を抑制できる。図7は、収束に導く方法を示す図である。図1の収束ガイド部103による、本発明において説明した方法で、収束が可能であることが保証される。収束ガイド部103は、図7に示されるように、ワープされたRF信号RFSigWarpd,l,f,roundを、当該収束ガイド部103の入力として取得する。RF信号差分計算部700は、ワープされたRF信号間の差分を計算し、各ラインの各深度について、各推定ラウンドにおける差分を、residued,l,roundとして出力する。RF信号ワーピング部100の主な目的は、RF信号の1つを、他と適合するように、修正することであるので、この残差d,l,roundは、RF信号同士が、どの程度適合しているかの尺度になる。roundにより特定されるある特定の推定ラウンドについて、推定された変位dispRoundd,l,roundの質を評価する目的で、RF信号ワーピング部100が、当該RF信号ワーピング部100の入力として、dispRoundd,l,roundを取得する。このことにより、RF信号がワープされ、適合レベルを評価するために、ワープされた信号間の差分が計算される。これらのワープされた信号は、後続ラウンドの変位推定のための入力としても用いられ、よってRFSigWarpd,l,f,round1と表される。よって、RFSigWarpd,l,f,round1から、residued,l,roundが計算される。RF信号差分計算部700図7で行われる計算の一例として、以下の式つまり、二乗差分和が挙げられるが、これに限定されるものではない。ここで、f1およびf2は、推定に用いるために選択された2フレームを表す。この例において、1ラインの全ての深度に対し同一の信号差分値が指定され、これは、対応するラインのRF信号が、全体的に適合していることを表す。小さな領域で変位が起こるだけではないときには、二乗差分和が適しており、その結果、ライン全体に反復が起こる。なお、RF信号差分計算部700図7で行われる計算の別の例として、以下の式つまり、絶対差分が挙げられるが、これに限定されるものではない。ここで、f1およびf2は、推定に用いるために選択された2フレームを表す。この例において、各ラインの各深度は、個別の信号差分値を有し、これは、RF信号が、局所的に適合していることを表す。深度方向に小さな領域だけで、変位が起こる場合には、絶対差分が適しており、その結果、その領域内での局所的反復が可能になる。最小残余計算部701は、以下の式に基づき、residued,l,roundにおける、各深度およびラインについて、残余が最小であるラウンド数を確認する。これは、最も正確な推定結果を表す。ここで関数は、fxを最小化するxの値を返す。収束ガイド部103図1の出力は、対象ラウンド数curRoundと、残余が最小であるラウンド数minResidueRoundd,lと、信号差分residued,l,roundとである。出力スイッチ104図1は、入力として、curRoundおよびminResidueRoundd,lを取得する。curRoundが最大値に達すると、出力スイッチ104は、最も精度の高い変位推定結果dispRoundd,l,minResidueRoundd,lを、最終的な出力dispOutとして選択する。そうでなければ、対象ラウンドdispRoundd,l,curRoundが、変位推定される代わりに、次の推定ラウンドのための、RF信号ワーピング部100の入力として用いられるdispRoundd,l,minResidueRoundd,lが、変位推定される。ここで、変位推定の多次元2次元or3次元への拡張について説明する。つまり、上記の実施の形態では、超音波の深さ方向における1次元の変位を求める方法について説明した。しかしながら、血管や腫瘍組織などの生体組織の拍動方向は、1次元内に限定されず、また、場所によっても、方向が異なる。したがって、複数方向から超音波を送受信することで、特定位置における複数方向の変位を求め、組織の2次元あるいは3次元的な変位を測定してもよい。こうして、多次元で測定することにより、診断精度の更なる向上が見込まれる。そして、超音波を複数方向から送受信する方法としては、超音波振動子が、1次元に配置されたリニアプローブを用いて、プローブの方向を変化させる、または、超音波振動子が2次元アレイ上に配置されたマトリックスプローブを用いて、ビームフォーミングにより、複数方向での送受信を行うなどの方法が可能である。ここで、リニアプローブの使用時には、磁気センサや加速度センサ、ジャイロなどの位置センサ、あるいは、カメラにより、プローブの位置を追跡することで、プローブの位置や向きを取得して、送受信時の方向を決定できる。なお、変位測定は、超音波の深さ方向など、1方向についてのみ行い、この変位ベクトルを、所定の2次元あるいは3次元座標系の各軸方向の成分に分解することで、多次元的に評価してもよい。この時、超音波の深さ方向と所定の座標系との対応関係が必要となるが、上述の位置センサなどを用いて、プローブの位置と方向を取得することで、両者を対応付ける。ここで、癌検出、あるいは、動脈硬化血管の硬さの診断などに適用できる歪み計測エラストグラフィとの関連について補足説明する。歪みは、変位の時間微分であり、組織の硬さの指標となるが、本発明により、変位推定精度が向上することで、結果として、歪みの測定精度が向上する。エラストグラフィでは微分値を用いるが、その他に、組織の特定部位の変位の時間変化波形なども、組織性状の判定に利用できる。例えば、癌の成長過程においては、癌の周辺および内部に血管新生が起こることが知られている。従って、これら血管の拍動、あるいは、血管の拍動に伴う周囲の組織の拍動などを変位量の時間変化として捉え、変化のパターンや振幅が癌に特有のものであるかどうかを判定することで、癌検出を行ってもよい。こうして、変位推定装置であって、超音波信号を用いて、変位を反復的に推定する変位推定装置であって、少なくとも1つの超音波信号を、トランスデューサから走査された媒体に送信する手段と、前記走査された媒体から反射した超音波信号を受信する手段と、推定ウインドウサイズを計算する手段と、前記推定されたウインドウサイズに基づいて、前記推定ウインドウの境界を計算する手段と、前記計算された推定ウインドウを用いて、前記超音波信号の深度方向に沿って変位を推定する手段と、前記推定された変位に基づいて、前記超音波信号をワーピングする手段と、前記方法の収束傾向を前記ワープされた超音波信号間の差分として計算することで、収束を導く手段とを備える変位推定装置が開示される。こうして、開示される方法は、収束を導く方法、精度を限定する要素を克服する方法、および結果の質を評価する方法を含み、包括的な反復手法を用いて変位推定の精度を保証する。つまり、次の問題が解決される。つまり、超音波を用いる変位推定において、先行技術では、推定の精度と解像度を高める方法を提供している。その中には、複数段階からなる推定または反復推定を導入しているものもある。しかし、1反復推定の収束を導く方法を組み込んだ包括的な反復手法の欠如、2品質評価方法の欠如、および3変位の精度を限定する要素を克服する方法の欠如という問題があり、本発明の変位推定方法により、これらの問題が解決される。なお、本発明は、装置、システム、集積回路などとして実現できるだけでなく、その装置等を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを記録したコンピュータ読み取り可能なCD-ROMなどの記録媒体として実現したり、そのプログラムを示す情報、データ又は信号として実現したりすることもできる。そして、それらプログラム、情報、データ及び信号は、インターネット等の通信ネットワークを介して配信してもよい。なお、本発明について、実施の形態に基づいて説明したが、本発明は、上記の実施の形態に限定されるものではない。本発明の主旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の範囲内に含まれる。"