"深さ方向に広がる超音波ビーム走査空間を順次形成することによって取得される第1ビームデータ群と第2ビームデータ群との間で注目点の変位を演算する超音波データ処理装置において、前記第1ビームデータ群と前記第2ビームデータ群との間における第1マッチング処理の結果に基づいて、前記注目点の変位を表す整数部として、深さ方向の整数値及び走査方向の整数値を演算する整数値演算手段と、前記第1マッチング処理の結果に対して第1サブピクセル処理を適用することによって、前記注目点の変位を表す小数部として、少なくとも走査方向の第1小数値を演算する第1小数値演算手段と、前記第1ビームデータ群と前記第2ビームデータ群との間における第2マッチング処理の結果に基づいて、前記注目点の変位を表す小数部として、少なくとも走査方向の第2小数値を演算する第2小数値演算手段と、前記注目点の変位を表す高精度小数部として、前記第1小数値及び前記第2小数値に基づく高精度小数値を演算する高精度小数値演算手段と、前記注目点の変位として、前記整数部及び前記高精度小数部に基づく実数情報を出力し得る出力部と、を含むことを特徴とする超音波データ処理装置。請求項1記載の装置において、前記第1ビームデータ群及び前記第2ビームデータ群は第1受信フレーム及び第2受信フレームを構成し、前記第1受信フレーム及び前記第2受信フレームに対して補間処理を伴う座標変換処理を実行して第1表示フレーム及び第2表示フレームを生成する変換手段を含み、前記変位演算手段は前記変換手段の前段に設けられた、ことを特徴とする超音波データ処理装置。請求項1又は2記載の装置において、前記高精度小数値演算手段は、前記第1小数値及び前記第2小数値の平均値として前記高精度小数値を演算する、ことを特徴とする超音波データ処理装置。請求項1乃至3のいずれか1項に記載の装置において、前記第2小数値演算手段は、前記第1ビームデータ群を半ピクセルシフトさせた補正ビームデータ群を生成する手段と、前記補正ビームデータ群と前記第2ビームデータ群との間において、前記第1マッチング処理と同じ処理としての前記第2マッチング処理を実行することにより、少なくとも走査方向の第2小数値を演算する手段と、を含むことを特徴とする超音波データ処理装置。請求項1乃至4のいずれか1項に記載の装置において、前記出力部は、前記注目点の変位として、前記整数部及び前記高精度小数部に基づく実数情報、及び、前記整数部及び前記第1小数部に基づく実数情報、を選択的に出力し得る、ことを特徴とする超音波データ処理装置。請求項5記載の装置において、前記出力部は、更に、前記注目点の変位として、前記整数部を選択的に出力し得る、ことを特徴とする超音波データ処理装置。請求項5又は6記載の装置において、前記出力部は、少なくとも前記注目点の深さに応じて出力対象を切り替える、ことを特徴とする超音波データ処理装置。請求項5又は6記載の装置において、前記出力部は、少なくともビーム間隔に応じて出力対象を切り替える、ことを特徴とする超音波データ処理装置。本発明の目的は、末広がりのビーム走査空間が形成され、座標変換前のデータに対して変位演算を行う場合に、特に深い部位における変位演算の精度を高めることにある。あるいは、本発明の目的は、状況に応じた相応しい変位推定法が適用されるようにすることにある。本発明は、深さ方向に広がる超音波ビーム走査空間を順次形成することによって取得される第1ビームデータ群と第2ビームデータ群との間で注目点の変位を演算する超音波データ処理装置において、前記第1ビームデータ群と前記第2ビームデータ群との間における第1マッチング処理の結果に基づいて、前記注目点の変位を表す整数部として、深さ方向の整数値及び走査方向の整数値を演算する整数値演算手段と、前記第1マッチング処理の結果に対して第1サブピクセル処理を適用することによって、前記注目点の変位を表す小数部として、少なくとも走査方向の第1小数値を演算する第1小数値演算手段と、前記第1ビームデータ群と前記第2ビームデータ群との間における第2マッチング処理の結果に基づいて、前記注目点の変位を表す小数部として、少なくとも走査方向の第2小数値を演算する第2小数値演算手段と、前記注目点の変位を表す高精度小数部として、前記第1小数値及び前記第2小数値に基づく高精度小数値を演算する高精度小数値演算手段と、前記注目点の変位として、前記整数部及び前記高精度小数部に基づく実数情報を出力し得る出力部と、を含むことを特徴とする。上記構成によれば、第1ビームデータ群と第2ビームデータ群との間で注目点の変位が演算される。各ビームデータ群は望ましくは二次元の受信フレーム又は三次元のボリュームデータを構成するものである。変位の演算に際しては、まず第1マッチング処理が実行され、これにより整数部として、深さ方向及び走査方向についての整数値が演算される。この場合、通常、パターンマッチング法が利用され、一方のフレーム上に注目点を含むテンプレートが設定され、他方のフレーム上において検索ブロックテンプレート相当のブロックがその位置を変えながら順次設定され、両者間でのパターンマッチング結果が順次演算される。最良のマッチング結果を生じさせる位置として、注目点の変位の整数部が特定される。しかし、整数部だけでは特に深い部位においてビーム間隔が広がっているところで変位推定精度が大きく低下してしまう。そこで、整数部の演算結果に対してサブピクセル処理が適用されて、整数部に付加される小数部が演算される。この場合には、少なくとも走査方向について小数値が演算されるが、更に深さ方向について小数値が演算されるようにしてもよい。以上によって整数部と小数部の組合せとして注目点の変位が実数として推定されることになる。サブピクセル処理においては誤差の変動が認められており、それを低減することによってより高精度に変位を推定することが望まれる。そこで、本発明では、第1ビームデータ群と第2ビームデータ群との間において第1マッチング処理とは異なる第2マッチング処理を適用し、更にサブピクセル処理を適用して第2小数値を演算するようにしている。第1小数値と第2小数値とを併せて考慮すれば小数値演算上の誤差を低減できることが確認されており、そのような性質を利用して、注目点の変位をより高精度に推定するものである。第2小数値の演算は、少なくとも走査方向について行われるのが望ましいが、更に深さ方向について行われてもよい。望ましくは、前記第1ビームデータ群及び前記第2ビームデータ群は第1受信フレーム及び第2受信フレームを構成し、前記第1受信フレーム及び前記第2受信フレームに対して補間処理を伴う座標変換処理を実行して第1表示フレーム及び第2表示フレームを生成する変換手段を含み、前記変位演算手段は前記変換手段の前段に設けられる。座標変換処理における補間処理によって誤差が生じるため、変位の演算は座標変換処理前に行うのが望ましい。望ましくは、前記高精度小数値演算手段は、前記第1小数値及び前記第2小数値の平均値として前記高精度小数値を演算する。このため、第2小数値の演算に際しては、第1小数値との比較、統合によって高精度小数値が得られるように、その演算条件が定められる。望ましくはその際に半ピクセル補正処理が適用され、その場合にシフトさせる方向が第1小数値によって定められる。望ましくは、前記第2小数値演算手段は、前記第1ビームデータ群を半ピクセルシフトさせた補正ビームデータ群を生成する手段と、前記補正ビームデータ群と前記第2ビームデータ群との間において、前記第1マッチング処理と同じ処理としての前記第2マッチング処理を実行することにより、少なくとも走査方向の第2小数値を演算する手段と、を含む。第2小数値の演算に際しては半ピクセルシフトに相当する半ピクセル加減算処理が適用される。望ましくは、前記出力部は、前記注目点の変位として、前記整数部及び前記高精度小数部に基づく実数情報、及び、前記整数部及び前記第1小数部に基づく実数情報、を選択的に出力し得る。望ましくは、前記出力部は、更に、前記注目点の変位として、前記整数部を選択的に出力し得る。望ましくは、前記出力部は、少なくとも前記注目点の深さに応じて出力対象を切り替える。望ましくは、少なくともビーム間隔に応じて出力対象を切り替える。演算量と推定精度との兼ね合いにおいて演算方式を動的に切り替えるのが望ましい。本発明によれば、末広がりのビーム走査空間が形成される場合において、特に深い部位における注目点の変位推定の演算を高められる。以下、本発明の好適な実施形態を図面に基づいて説明する。図1には本発明に係る超音波データ処理装置の好適な実施形態が示されている。この超音波データ処理装置は超音波診断装置として構成されている。プローブ10は、被検体の体表面上に当接して用いられ超音波の送受波を行うものである。プローブ10は1Dアレイ振動子を備えている。1Dアレイ振動子は複数の振動素子で構成され、それによって超音波ビームが形成され、それが電子的に走査される。本実施形態では電子セクタ走査が適用されている。その他にコンベックス走査等が適用されてもよい。それらの走査方式においては、扇状に広がる複数のビームラインが構成され、深い部位ほどビーム間隔が増大することになる。図1においてはセクタ走査面が示され、rは深さ方向を示し、θは電子走査方向を示している。走査面は繰り返し形成され、その都度、受信フレーム複数ビームデータの配列が構成される。送受信部12は送信ビームフォーマー及び受信ビームフォーマーとして機能する。送信時において、複数の送信信号がアレイ振動子に供給され、これによって送信ビームが形成される。受信時において、生体内からの反射波がアレイ振動子にて受波され、アレイ振動子から複数の受信信号が出力される。送受信部12では複数の受信信号に対する整相加算処理によってビームデータを生成する。そのビームデータは画像形成部18側へ出力される。変位推定部14は隣接フレーム間において注目点の変位二次元移動ベクトルを演算するモジュールであり、その機能は例えばソフトウエア機能によって実現される。変位推定部14の具体的内容については後に詳述するが、推定された変位の情報は必要に応じて変位補正部16、画像形成部18及び表示処理部20へ出力される。変位補正部16は、推定された変位を利用してフレームに対して補正処理を適用するものである。各画素が注目画素とされる場合、各注目画素の変位を利用して各注目画素の移動を補正によってキャンセルすることが可能となる。変位補正部16は必要に応じて設けられるものである。画像形成部18は、デジタルスキャンコンバータDSCを含み、各受信フレームに対して座標変換処理や補間処理を適用し、表示フレームを順次生成する。この画像形成部18によってBモード画像が形成される。表示フレームは表示処理部20へ送られる。なお、画像形成部18においてフレーム相関処理が行われる場合、上記の変位の情報を利用して変位分をキャンセルしつつフレーム相関処理を行うようにしてもよい。表示処理部20は、画像合成機能を備え、変位に基づいてグラフその他の画像を形成する機能も備えている。表示部22には超音波画像が表示される。制御部24は、CPU及び動作プログラムによって構成され、図1に示される各構成の動作を制御している。入力部26は操作パネルによって構成される。図2には、図1に示された変位推定部14の構成例が示されている。この変位推定部14は、ブロック抽出部30、マッチング実行部32、整数推定部34、第1小数推定部36、第2小数推定部38、出力部40、サイズ決定部42及び動作制御部44を有している。サイズ決定部42は、ライン密度D1の情報に基づいてマッチング処理で参照するテンプレート及び検索ブロックのサイズを適応的に設定するものである。但し、ライン密度や深さ方向によらずにそのサイズを一定にするようにしてもよい。動作制御部44は、ブロック深さR1の情報及びライン密度D1の情報に基づいて、変位の推定方式を切り替えるものであり、具体的には、第1小数推定部36及び第2小数推定部38の動作をオンオフ制御している。出力部40においても動作の有無が参照されている。深さによって段階的に推定精度を高めるのが望ましい。浅い領域において簡易な推定方法を選択するならば演算量を削減できるから効率的である。図2に示すブロック抽出部30、マッチング実行部32及び整数推定部34は以下に説明する整数値演算法を実行するものである。第1小数推定部36は以下に説明する第1のサブピクセル推定法を実行するものである。第2小数推定部38は以下に説明する第2のサブピクセル推定法を実行するものである。それは、実際のところ、ブロック抽出部30及びマッチング実行部32に相当する機能を備えているが、図1においてはそれらについて図示省略されている。変位前のデータ及び変位後のデータが第2小数推定部38にも入力されている図1に示す概念図においてはそれらのデータは整数推定部34及び第1小数推定部36を経由して第2小数推定部38へ送られている。1整数値演算法:ブロックマッチングBlockMatching法では、変位前の基準画像実際には受信フレーム上において注目点を含む第1テンプレートブロックRθ,rが設定され、その第1テンプレートブロックRθ,rが変位後の画像実際には受信フレーム上において検索エリア中のどこに移動したかの総当り探索が実行される。図3のAには変位前の受信フレームの一部分が示されている。図3のBには変位後の受信フレームの一部分検索エリアが示されている。同図においては各格子点がデータ座標を示している。変位前の受信フレーム上においては、注目点を中心として所定形状をもったテンプレートブロック46が設定される。一方、後の受信フレームにおいては注目点に対応する対応点を中心として大きな検索エリアが設定され、その範囲内において検索ブロック48を各方向に1ピクセルずつ順次移動させながら例えばラスタースキャンの実行、マッチング処理が実行される。これにより各位置において相関指標が計算される。相関指標としては、相互相関と、差分の絶対和が広く知られている。ここでは以下の式に示すSADを説明する。SADの値が小さいほど、相関は高いと評価される。このため、図4に示すように、SAD1Δθ,Δrが最小値となる検索ブロックをマッチングブロック最適位置にある検索ブロックとし、テンプレートブロックとマッチングブロックの座標差を整数単位での推定変位Δθ0,Δr0とする図4には方位方向についてのSAD関数が示されている。なお、相互相関を採用する場合、それは値が高いほど、相関が高いと評価されるため、値が最大値となる検索ブロックをマッチングブロックとすることになる。以降の説明ではθ,rは、各ブロック中心点V0θ,rV1θ,rの座標とする。上記において、ブロックのサイズ=MBMBであり、検索エリアのサイズ=MAMAであり、SAD1のサイズ=MCNC=MA-MB1NA-NB1であり、Δθの範囲=-MC/2~MC/2-1であり、Δrの範囲=-NC/2~NC/2-1である。Δθ0,Δr0=SAD1Δθ,Δrである。図2において、ブロック抽出部30は、変位前の受信フレームV0θ,rに対してテンプレートブロックを設定し、変位後の受信フレームV1θ,rに対して検索エリアSθ,r及び検索ブロックを設定するものである。マッチング実行部32は、テンプレートブロックと各位置に設定された検索ブロックとの間で相関指標としてSADを演算するものである。整数推定部34は、最良値をとるSAD1Δθ,Δrに基づいて整数部として推定変位Δθ0,Δr0を求める。2第1のサブピクセル推定法:上記の整数値演算法で求めたSAD1Δθ,Δrと整数単位の推定変位Δθ0,Δr0から、変位の小数部として、サブピクセルSubPixeldθ1,dr1を推定する。なお、dθ1,dr1は整数成分を含まない整数未満の成分のみとする。つまり、|dθ1|≦0.5、かつ、|dr1|≦0.5である。サブピクセル推定法としては、キュービックスプラインCubicSpline補間法やパラボラフィッティングParabolicFitting法等、多くの方式が提案されている。ここでは、最も一般的なパラフィッティング法を説明する。同法は、整数値演算法で求めたSADの分布を2次曲線と仮定することで、図5に示すように、以下の式および式でSubPixelを算出する方法である。なお、方向毎のパラボラフィッティング法に使用するSAD1Δθ,Δrは、整数単位の変位Δθ0,Δr0前後の3個のみとなる。それらが図6及び図7に図解されている。図2に示した第1小数推定部は以上の演算を行うものである。推定結果として、変位の小数部dθ1,dr1が得られる。3第2のサブピクセル推定法高精度サブピクセル推定法、半ピクセル補正法:上記の伝統的なサブピクセル推定では何れの方式でも、整数単位もしくは±12に偏る誤差が含まれる。これに関し、図8に示すように、第1テンプレートブロックでの推定誤差と、正負に05ピクセル並進した第2テンプレートブロックgm,nによるサブピクセル推定誤差は、ほぼ正負対称となる。半ピクセル補正は、この特性を利用し、整数値演算法及び上記サブピクセル推定法の適用による第1テンプレートブロックによるサブピクセル推定値dθ1と、第2テンプレートブロックによるサブピクセル推定値dθ2とを平均し、サブピクセル推定誤差を低減するものである。手順は以下の通りである。まず、±05並進した第2テンプレートブロックgm,nを作成する。それが図9に示されている。第1小数部dθ1が負なら、正側へ0.5並進させた第2テンプレートブロックを線形補間処理により生成する図9のA。第1小数部dθ1が正なら、負側へ―0.5並進させた第2テンプレートブロックを線形補間処理により生成する図9のB。図9のCには変位後の受信フレームが示されている。上記処理を数式で表現すれば以下の4式となる。次に、第2テンプレートブロックgm,nを用いて、相関指標を計算する。その場合、整数値演算法で利用した相関指標と同じ相関指標を用いる。上記の例ではSADが利用されていたのでここではSADが演算される。すなわち、以下の5式が演算される。ちなみに、演算量の削減の観点からは、距離方向深さ方向については整数値演算法で求められたΔr0をそのまま用い、方位方向Δθ0の近傍のみで探索を行うのが望ましい。パラボラフィッティング法によるSAD2ΔθはΔθ=-1,0,1の3個のみでよいので、その意味で半ピクセル補正法の適用にする処理時間の増加はあまり大きくはない。続いて、第2テンプレートブロックgm,nのサブピクセルdθ2を計算する。この場合には上記のサブピクセル推定法と同じ推定法が用いられる。すなわち、本例ではパラボラフィッティング法が利用され、以下の6式が演算される。次に、以下の7式に示すように、±05並進した分を差し引き、|dθ2|≦0.5とする。この段階で、dθ1とdθ2の間の誤差は、ほぼ正負対称となる図8参照。以上の4式から7式までを実行するのが図2に示した第2小数推定部38である。図2に示した出力部40は、図10に示される条件に従って出力を合成あるいは選択する。すなわち、第1の場合には、整数推定部のみが機能し、それによって推定された整数部が最終的な出力となる。例えば、浅い領域、ビーム密度が大きい領域についてこの条件が適用される。第2の場合には、整数推定部及び第1小数推定部が機能し、それによって推定された整数部及び第1小数部が合成された上で出力される。例えば、中間領域、ビーム密度が中程度の領域についてこの条件が適用される。図10においては、方位方向及び深さ方向の両方向について整数値に小数値が付加されている。但し、一方方向だけについて小数値の演算、付加を行うようにしてもよい。第3の場合には、整数推定部、第1小数推定部及び第2小数推定部の3つが機能する。特に、第1小数値と第2小数値の平均値が演算され、それは高精度小数値と言いうるものである。整数値に高精度小数値が付加され、それが出力とされる。例えば、深い領域、ビーム密度が小さい領域についてこの条件が適用される。図10に示す例では方位方向について高精度小数値が演算されているが、両向についてそのような演算を行うようにしてもよい図2においては両方向について第2小数値が演算されている。上記実施形態は二次元データに対する処理であったが三次元データに対して上記同様の処理を適用することも可能である。"