"心臓を含む三次元空間に対して超音波を送受するプローブと、前記プローブを制御することにより三次元空間から受信信号を得る送受信部と、前記受信信号に基づいて形成される三次元空間のボリュームデータ内において、前記心臓の心筋表面を抽出する表面抽出部と、抽出された前記心筋表面に複数の参照点を設定する参照点設定部と、前記ボリュームデータ内において複数の時相に亘って前記各参照点の位置を三次元的に追跡する追跡処理部と、前記複数の参照点に関する追跡の結果を反映させた拍動解析画像を形成する画像形成部と、を有する、ことを特徴とする超音波診断装置。請求項1に記載の超音波診断装置において、前記画像形成部は、前記心筋表面を立体的に映し出した心筋画像上において前記各参照点の位置にマーキングを付した三次元画像を形成する、ことを特徴とする超音波診断装置。請求項1または2に記載の超音波診断装置において、前記画像形成部は、前記三次元空間に対応した三次元座標系内において前記各参照点の位置にマーキングを付した座標系画像を形成する、ことを特徴とする超音波診断装置。請求項3に記載の超音波診断装置において、前記画像形成部は、前記三次元座標系内において前記各参照点に関する複数の時相における位置にマーキングを付した座標系画像を形成する、ことを特徴とする超音波診断装置。請求項1から4のいずれか1項に記載の超音波診断装置において、前記画像形成部は、前記各参照点ごとに時相間における変位量を示した変位量グラフの画像を形成する、ことを特徴とする超音波診断装置。上述した背景事情において、本願の発明者は、超音波を利用した心筋の診断について研究開発を重ねてきた。特に、心筋の運動を評価する技術に注目した。本発明は、その研究開発の過程において成されたものであり、その目的は、心筋の運動を評価するための拍動解析画像を提供することにある。上記目的にかなう好適な超音波診断装置は、心臓を含む三次元空間に対して超音波を送受するプローブと、前記プローブを制御することにより三次元空間から受信信号を得る送受信部と、前記受信信号に基づいて形成される三次元空間のボリュームデータ内において前記心臓の心筋表面を抽出する表面抽出部と、抽出された前記心筋表面に複数の参照点を設定する参照点設定部と、前記ボリュームデータ内において複数の時相に亘って前記各参照点の位置を三次元的に追跡する追跡処理部と、前記複数の参照点に関する追跡の結果を反映させた拍動解析画像を形成する画像形成部と、を有することを特徴とする。望ましい具体例において、前記画像形成部は、前記心筋表面を立体的に映し出した心筋画像上において前記各参照点の位置にマーキングを付した三次元画像を形成する、ことを特徴とする。望ましい具体例において、前記画像形成部は、前記三次元空間に対応した三次元座標系内において前記各参照点の位置にマーキングを付した座標系画像を形成する、ことを特徴とする。望ましい具体例において、前記画像形成部は、前記三次元座標系内において前記各参照点に関する複数の時相における位置にマーキングを付した座標系画像を形成する、ことを特徴とする。望ましい具体例において、前記画像形成部は、前記各参照点ごとに時相間における変位量を示した変位量グラフの画像を形成する、ことを特徴とする。本発明により、心筋の運動を評価するための拍動解析画像が提供される。例えば、本発明の好適な実施形態によれば、各時相ごとに心筋の運動を評価するための拍動解析画像が提供される。以下に本発明の好適な実施形態を説明する。図1は、本発明の実施において好適な超音波診断装置の全体構成を示す図である。プローブ10は、対象組織を含む領域に対して超音波を送受波する超音波プローブである。プローブ10は、超音波を送受する複数の振動素子を備えており、複数の振動素子が送受信部12によって送信制御されて送信ビームが形成される。また、複数の振動素子が対象組織を含む領域内から得られる超音波を受波し、これにより得られた信号が送受信部12へ出力され、送受信部12が受信ビームを形成して受信ビームに沿ってエコーデータが収集される。プローブ10は、超音波ビーム送信ビームと受信ビームを三次元空間内において走査して立体的にエコーデータを収集する三次元プローブが好適である。例えば、一次元的に配列された複数の振動素子1Dアレイ振動子によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが立体的に走査される。また、二次元的に配列された複数の振動素子2Dアレイ振動子を電子的に制御して超音波ビームを立体的に走査してもよい。三次元空間内において超音波ビームが走査されてエコーデータが収集されると、その三次元空間に対応した三次元データ空間を構成する複数のボクセルについてのエコーデータボクセルデータが図示しないメモリなどに記憶される。複数のボクセルデータで構成された三次元データ空間はボリュームデータと呼ばれる。そして、そのボリュームデータに対して、心筋表面抽出部20以降の各部において各種の処理が実行される。本発明における好適な対象組織は心臓であり、本実施形態においては、その心臓の心筋の運動に関する伝播を評価するための拍動解析画像が形成される。そこで、その拍動解析画像が形成されるまでの各種の処理について説明する。なお、図1に示した部分構成については、以下の説明において図1の符号を利用する。心筋表面抽出部20は、複数のボクセルデータで構成されるボリュームデータ内において心臓の心筋表面を抽出する。心筋表面抽出部20は、例えば、複数のボクセルデータを閾値により二値化処理するなどして、心臓の心筋表面に対応した複数のボクセルデータを識別する。なお、心臓内の心腔に対応したボクセルは、ボクセル値エコー値が比較的小さいため、閾値を利用して心腔に対応した複数のボクセルを識別し、心筋表面として心腔の表面を抽出するようにしてもよい。図2は、ボリュームデータ内において抽出された心筋100を示す図である。心臓は、左房と左室と右房と右室を備えており、さらに、心筋100の表面には左冠動脈や右冠動脈などの血管が存在する。心筋100の表面が抽出されると、参照点設定部30は、抽出された心筋100の表面に複数の参照点を設定する。図2において、心筋100の表面に示される数字を付された複数の丸印が参照点マーキングポジションである。複数の参照点は、心筋100の表面において特徴的な複数箇所に設定される。例えば、左冠動脈の分岐点、右冠動脈の分岐点、動脈静脈などの血管の接合部分などに複数の参照点が設定される。なお、複数の参照点を設定するにあたっては、例えば、図2に示す心筋100の表示画像を表示部60に表示させ、タッチパネルやマウスなどの操作デバイスを介して、ユーザが表示画像上に各参照点の位置を指定し、その指定に応じて参照点設定部30が各参照点を設定するようにしてもよい。複数の参照点が設定されると、追跡処理部40は、ボリュームデータ内において複数の時相に亘って各参照点の位置を三次元的に追跡する。追跡処理部40は、例えば、ボリュームデータに対応した三次元座標系内において、パターンマッチングなどを利用して、複数の時相に亘って各参照点の座標値を追跡する。例えば、参照点1について、設定された時点である時刻0における座標値が確認され、その参照点1の近傍の複数のボクセルを含むようにテンプレートが設定される。そして、次の時相である時刻1のボリュームデータ内において、テンプレート内の複数のボクセルと最も相関関係が強いボリューム部分が特定され、その特定された位置に参照点1が移動されて時刻1における参照点1の座標値が確認される。図3は、複数の参照点に関する追跡の結果を説明するための図である。図3には、追跡処理部40における追跡処理の結果として得られたデータが図示されている。図3Aには、各参照点の座標データが示されている。つまり、各参照点nについて時刻0から順に各時刻t時相ごとに確認された座標値xnynznが示されている。また、複数の参照点のうちの一つの参照点を基準参照点とし、その基準参照点との間の差分の座標値を算出してもよい。図3Bには、各参照点の差分座標データが示されている。つまり、参照点1を基準参照点として常に原点000とした場合に、他の各参照点nについて時刻0から順に各時刻t時相ごとに確認された、参照点1との間の差分の座標値dxndyndznが示されている。追跡処理部40における追跡処理の結果として例えば図3に示すデータが得られると、画像形成部50は、複数の参照点に関する追跡の結果を反映させた拍動解析画像を形成する。画像形成部50は、拍動解析画像として、例えば、心筋表面を立体的に映し出した心筋画像上において各参照点の位置にマーキングを付した三次元画像を形成する。その心筋画像の形成においては、例えば公知のサーフェスレンダリングが利用される。つまり、複数のボクセルデータで構成されるボリュームデータに対して、視線となる複数のレイが設定され、各レイごとにレンダリング演算が行われて心筋表面が検出され、複数のレイから得られるレンダリング演算の結果に基づいて、例えば図2に示すような心筋100の表面を立体的に映し出した心筋画像が形成される。そして、その心筋画像が複数の時相に亘って各時相ごとに形成される。さらに、例えば図3ABのデータに基づいて、各時相ごとに各時刻ごとに心筋表面の画像上の各参照点の位置にマーキングが付される。これにより、例えば、図2に示す心筋100の画像と、数字を付された複数の参照点マーキングポジションを示した三次元の動画像が得られる。なお、心筋100の動画像の形成においては、特許第3537594号公報などに記載された公知の再構成処理再構築処理などを利用して、心筋の高速な動きに対応するようにしてもよい。複数の時相に亘って得られた三次元画像を表示する際には、通常の速度で再生する動画再生の他、スロー再生や早送り再生などの特殊再生が実現されてもよい。また、ユーザによって指定された時相の静止画像を表示するようにしてもよい。各種の再生機能を利用することにより、ユーザ検査者は、例えば図2の心筋100の表面に付された複数の参照点のうち、特異な動きをする参照点などを発見して、その参照点に対応した部位の異常な動作などを診断することが可能になる。さらに、画像形成部50は、拍動解析画像として、三次元空間に対応した三次元座標系内において各参照点の位置にマーキングを付した座標系画像を形成してもよい。図4は、座標系画像の具体例を示す図である。画像形成部50は、例えば図3Bのデータに基づいて、図4の座標系画像を形成する。つまり、ボリュームデータが得られた三次元空間に対応したXYZ直交座標系内に、図3Bの差分座標データに基づいて、各参照点nの位置に丸印などの表示を施して図4の座標系画像を形成する。座標系画像は各時刻ごとにその時刻のみの各参照点nを示してもよいし、図4に示すように、ある時刻における各参照点nの位置を実線丸印で示しつつ、その後の時刻に得られる各参照点nの位置を破線丸印などで示してもよい。また、基準となる参照点1からの相対的な位置を示す相対位置ベクトル実線矢印が表示されてもよいし、各参照点nごとにその参照点nについての時相間の変位を示す変位ベクトル破線矢印が表示されてもよい。変位ベクトルを表示させることにより、ユーザ検査者は、複数の参照点のうち、極端に動きの鈍い参照点などを発見して、その参照点に対応した部位の異常な動作などを診断することが可能になる。さらに、画像形成部50は、拍動解析画像として、各参照点ごとに時相間における変位量を示した変位量グラフの画像を形成してもよい。図5は、変位量グラフの具体例を示す図である。画像形成部50は、例えば図3Bのデータに基づいて図5の変位量グラフを形成する。例えば、図3Bのデータから各参照点nごとに、時刻0と時刻tとの間における座標間の距離移動距離が算出され、その距離が各参照点nの変位量とされ、図5に示すような横軸を参照点nとして縦軸に変位量を示したグラフが形成される。変位量グラフから、ユーザ検査者は、複数の参照点のうち、極端に動きの鈍い参照点図5の例では参照点4,5などを発見して、その参照点に対応した部位の異常狭窄などを診断することが可能になる。なお、例えば、ユーザが任意の2つの参照点を指定して、各時刻ごとにその2つの参照点間の距離が算出されてグラフなどで表示されてもよい。画像形成部50において形成された拍動解析画像は、表示部60に表示される。例えば心筋表面の画像上に各参照点の位置を示した三次元画像図2と座標系画像図4と変位量グラフの画像図5のうち、ユーザによって選択された画像が表示部60に表示される。また、三次元画像と座標系画像と変位量グラフの画像のうちのいずれか又は全てを組み合わせて並べて表示した表示画像が形成されてもよい。以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。"