(57)【特許請求の範囲】
【請求項１】
超音波プローブからの出力に基づいて、心腔内の流体のドプラ情報を生成するドプラ情報生成部と、
前記超音波プローブからの出力に基づいて、心筋領域の各位置における３次元運動ベクトルを生成する運動情報生成部と、
前記３次元運動ベクトルを、互いに直交する３つの方向に分解し、前記３つの方向のうちの少なくとも２つの方向にそれぞれ沿った複数の断面各々において、前記分解された３次元運動ベクトルと前記ドプラ情報とに基づき前記各々の断面内の各位置の２次元血流ベクトルを計算し、前記各位置において前記計算された複数の２次元血流ベクトルをそれぞれ合成することにより、前記心腔内の３次元領域の各位置における前記流体の速さに関する３次元のベクトル情報を生成する処理部と、
を具備する超音波診断装置。
【請求項２】
前記ベクトル情報は、前記３次元領域の各位置における前記流体の移動方向と移動速度とを有する、請求項１に記載の超音波診断装置。
【請求項３】
前記処理部は、前記ベクトル情報に対応するマーカを前記各位置に対応付けて表示可能なベクトル場画像を生成する、請求項２に記載の超音波診断装置。
【請求項４】
前記処理部は、前記マーカの大きさを、前記移動速度に応じて決定する、請求項３に記載の超音波診断装置。
【請求項５】
前記処理部は、前記マーカの色相と、明度と、彩度とのうち少なくとも１つを、前記移動速度に応じて決定する、請求項３又は４に記載の超音波診断装置。
【請求項６】
前記３次元運動ベクトルと心臓に関するボリュームデータとに基づいて前記心臓の運動情報を有するレンダリング画像を生成する画像生成部と、
前記レンダリング画像と前記ベクトル場画像とを並列表示または重畳表示する表示部と、
をさらに具備する請求項３に記載の超音波診断装置。
【請求項７】
前記処理部は、前記レンダリング画像の生成に関するレンダリング処理を、前記ベクトル情報に適用することにより、前記ベクトル場画像を生成する請求項６に記載の超音波診断装置。
【請求項８】
前記表示部は、前記レンダリング画像における前記運動情報の表示態様と、前記ベクトル場画像における前記心腔の壁部分に隣接する前記ベクトル情報の表示態様とを、同一の表示態様として表示する請求項６または７に記載の超音波診断装置。
【請求項９】
前記表示態様は、色相と明度と彩度とのうち少なくとも一つである請求項８に記載の超音波診断装置。
【請求項１０】
前記処理部は、
前記３次元運動ベクトルを、前記３つの方向に対応する、走査線方向と、前記走査線方向に直交しかつ互いに直交する第１方向と第２方向とにそれぞれ投影することにより、前記３次元運動ベクトルを分解し、
前記走査線方向と前記第１方向とにより規定され、前記第２方向に沿った複数の第１断面各々において、前記分解された３次元運動ベクトルと前記ドプラ情報とを用いて、前記第１断面内の各位置における前記流体の速さに関する第１の２次元血流ベクトルを計算し、
前記走査線方向と前記第２方向とにより規定され、前記第１方向に沿った複数の第２断面各々において、前記分解された３次元運動ベクトルと前記ドプラ情報とを用いて、前記第２断面内の各位置における前記流体の速さに関する第２の２次元血流ベクトルを計算し、
前記第１の２次元血流ベクトルと前記第２の２次元血流ベクトルとを合成することにより、前記３次元のベクトル情報を生成する請求項１乃至９のうちいずれか一項に記載の超音波診断装置。
【請求項１１】
超音波プローブからの出力を基に生成された心筋領域の各位置における３次元運動ベクトルを、互いに直交する３つの方向に分解し、前記３つの方向のうちの少なくとも２つの方向にそれぞれ沿った複数の断面各々において、前記分解された３次元運動ベクトルと前記超音波プローブからの出力を基に生成された心腔内の流体のドプラ情報とに基づき前記各々の断面内の各位置の２次元血流ベクトルを計算し、前記各位置において前記計算された複数の２次元血流ベクトルをそれぞれ合成することにより、前記心腔内の３次元領域の各位置における前記流体の速さに関する３次元のベクトル情報を生成する処理部を具備する医用画像処理装置。
【請求項１２】
コンピュータに、
超音波プローブからの出力を基に生成された心筋領域の各位置における３次元運動ベクトルを、互いに直交する３つの方向に分解し、前記３つの方向のうちの少なくとも２つの方向にそれぞれ沿った複数の断面各々において、前記分解された３次元運動ベクトルと前記超音波プローブからの出力を基に生成された心腔内の流体のドプラ情報とに基づき前記各々の断面内の各位置の２次元血流ベクトルを計算し、前記各位置において前記計算された複数の２次元血流ベクトルをそれぞれ合成することにより、前記心腔内の３次元領域の各位置における前記流体の速さに関する３次元のベクトル情報を生成すること、
を実現させる医用画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、血流および組織の速度情報を取得可能な超音波診断装置と医用画像処理装置と医用画像処理プログラムとに関する。
【背景技術】
【０００２】
従来、超音波診断装置では、血流や組織からドプラ情報（速度情報）を取得し、それらを２次元マッピングすることにより２次元血流画像を構成するカラードプラ法がよく知られている。一般に、カラードプラ法で表示できる速度情報は、超音波ビーム方向の成分である。また、２次元領域内において、血流などの流れの情報として得られる血流などの速度ベクトルを利用して、血流や組織などの運動状態を反映させた物理量を求める技術がある。
【０００３】
カラードプラ情報と壁運動情報とを用いて、局所領域の質量保存の法則から心腔内の血流速度をベクトル表示する技術がある。血流ベクトルを表示する上記手法において、壁運動情報として、実際の心臓はねじれ運動をしているため、心壁運動情報には３次元的なねじれ運動の情報を用いる必要がある。しかしながら、従来の方法では２次元の平面上での血流ベクトルの計算しか考慮していないため、心臓本来のねじれ運動情報を利用できない問題がある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２４４８３５号公報
【非特許文献】
【０００５】
【非特許文献１】Ｉｔａｔａｎｉ．Ｋ  ｅｔ  ａｌ．  “Ｉｎｔｒａｖｅｎｔｒｉｃｕｌａｒ  Ｆｌｏｗ  Ｖｅｌｏｃｉｔｙ  Ｖｅｃｔｏｒ  Ｖｉｓｕａｌｉｚａｔｉｏｎ  Ｂａｓｅｄ  ｏｎ  ｔｈｅ  Ｃｏｎｔｉｎｕｉｔｙ  Ｅｑｕａｔｉｏｎ  ａｎｄ  Ｍｅｓｕｒｅｍｅｎｔｓ  ｏｆ  Ｖｏｒｔｉｃｉｔｙ  ａｎｄ  Ｗａｌｌ  Ｓｈｅａｒ  Ｓｔｒｅｓｓ”  Ｐｒｏｃｅｅｄｉｎｇｓ  ｏｆ  Ｓｙｍｐｏｓｉｕｍ  ｏｎ  Ｕｌｔｒａｓｏｎｉｃ  Ｅｌｅｃｔｒｏｎｉｃｓ，  Ｖｏｌ．３３（２０１２）  ｐｐ．１４５−１４６  １３−１５  Ｎｏｖｅｍｂｅｒ，２０１２．
【発明の概要】
【発明が解決しようとする課題】
【０００６】
目的は、心臓の３次元的な運動を考慮して、実際の血流に近い血流ベクトルを求めることが可能な超音波診断装置、医用画像処理装置および医用画像処理プログラムを提供することにある。
【課題を解決するための手段】
【０００７】
本実施形態に係る超音波診断装置は、超音波プローブからの出力に基づいて、心腔内の流体のドプラ情報を生成するドプラ情報生成部と、前記超音波プローブからの出力に基づいて、心筋領域の各位置における３次元運動ベクトルを生成する運動情報生成部と、前記３次元運動ベクトルを、互いに直交する３つの方向に分解し、前記３つの方向のうちの少なくとも２つの方向にそれぞれ沿った複数の断面各々において、前記分解された３次元運動ベクトルと前記ドプラ情報とに基づき前記各々の断面内の各位置の２次元血流ベクトルを計算し、前記各位置において前記計算された複数の２次元血流ベクトルをそれぞれ合成することにより、前記心腔内の３次元領域の各位置における前記流体の速さに関する３次元のベクトル情報を生成する処理部と、を具備する。
【図面の簡単な説明】
【０００８】
【図１】図１は、本実施形態に係る超音波診断装置１の構成を示す模式図である。
【図２】図２は、本実施形態に係るボリュームデータの生成に関する撮像ボリューム（被走査領域）の一例を示す図である。
【図３】図３は、本実施形態に係る撮像ボリュームのｘｙ断面における３次元運動ベクトルの分解の一例を示す図である。
【図４】図４は、本実施形態に係る断面の複数の部分領域（微小領域）各々における２次元血流ベクトルの計算過程の一例を示す図である。
【図５】図５は、本実施形態に係る一つの第１断面内の各位置における第１の２次元血流ベクトルの一例を示す図である。
【図６】図６は、本実施形態に係る一つの第２断面内の各位置における第２の２次元血流ベクトルの一例を示す図である。
【図７】図７は、本実施形態に係り、レンダリング処理により生成されたベクトル場画像をレンダリング画像に重畳した第１重畳画像の一例を示す図である。
【図８】図８は、本実施形態に係り、レンダリング処理により生成されたベクトル場画像の一例を示す図である。
【図９】図９は、本実施形態に係り、３つのＭＰＲ画像に３つのＭＰＲ画像をそれぞれ重畳した３つの第２重畳画像を示す図である。
【図１０】図１０は、本実施形態に係り、ベクトル場画像生成表示処理の処理手順の一例を示すフローチャートである。
【発明を実施するための形態】
【０００９】
以下、図面を参照しながら本実施形態に係わる超音波診断装置を説明する。なお、以下の説明において、略同一の構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１０】
図１は、本実施形態に係る超音波診断装置１の構成を示す模式図である。同図に示すように、超音波診断装置１は、超音波プローブ１１と、装置本体１３と、装置本体１３に接続され操作者からの各種指示・命令・情報を装置本体１３に取り込むための入力インタフェース回路（入力部）１５と、表示回路（表示部）１７とを有する。
【００１１】
超音波プローブ１１は、複数の圧電振動子と、整合層と、複数の圧電振動子の背面側に設けられるバッキング材とを有する。複数の圧電振動子は、圧電セラミックス等の音響／電気可逆的変換素子である。複数の圧電振動子は並列され、超音波プローブ１１の先端に装備される。以下、説明を簡単にするために一つの圧電振動子が一チャンネルを構成するものとして説明する。複数の圧電振動子各々は、後述する送受信回路２３から供給される駆動信号に応答して超音波を発生する。
【００１２】
以下、超音波プローブ１１は、アジマス（Ａｚｉｍｕｔｈ）方向とエレベーション（Ｅｌｅｖａｔｉｏｎ）方向とに複数の圧電振動子を配列した２次元アレイプローブとして説明する。アジマス方向とは例えば、圧電振動子の配列方向に沿った電子走査方向に対応する。エレベーション方向とは、超音波走査面の揺動方向に対応する。なお、本実施形態に係る超音波プローブ１１は、３次元エコー法を実行可能な任意の超音波プローブであってもよい。
【００１３】
なお、超音波プローブ１１は、２次元アレイプローブに限定されず、メカニカル４次元プローブであってもよい。メカニカル４次元プローブは、１列に配列された複数の圧電振動子（以下、１次元アレイ振動子と呼ぶ）と、複数の圧電振動子の配列方向（アジマス方向）に直交する方向（エレベーション方向）に１次元アレイ振動子を揺動させる揺動機構とを有する。揺動機構は、例えば、アジマス方向を回転軸として１次元アレイ振動子を所定の角度範囲で回転させる回転機構と、回転機構に伝達する動力を発生するモータとを有する。
【００１４】
回転機構は、例えば、１次元アレイ振動子の端点を支持する支持フレームと、ベアリング等を介して支持フレームを回転させる回転フレームと、モータにより発生された動力を伝達する伝達部材とを有する。伝達部材は、例えば、ラックアンドピニオン、チェーンおよびスプロケット、ベルトおよびプーリーなどである。
【００１５】
２次元アレイプローブおよびメカニカル４次元プローブは、３次元的なエコー信号を取得することができる超音波プローブ１１である。また、超音波プローブ１１は、１次元アレイプローブであってもよい。このとき、エレベーション方向に超音波プローブ１１を揺動させる操作者の操作により、３次元的なエコー信号が取得される。
【００１６】
超音波プローブ１１を介して被検体Ｐに超音波が送信されると、送信された超音波（以下、送信超音波と呼ぶ）は、被検体Ｐ内の生体組織における音響インピーダンスの不連続面で反射される。圧電振動子は、反射された超音波を受信し、エコー信号を発生する。エコー信号の振幅は、超音波の反射に関する不連続面を境界とする音響インピーダンスの差に依存する。また、送信超音波が移動している血流で反射された場合のエコー信号の周波数は、ドプラ効果により、移動体（血流等）の超音波送信方向の速度成分に依存して偏移する。
【００１７】
整合層は、被検体Ｐに対する超音波の送受信を効率よくするために、複数の圧電振動子の超音波放射面側に設けられる。バッキング材は、圧電振動子の後方への超音波の伝搬を防止する。
【００１８】
装置本体１３は、バス（データバス：ｄａｔａ  ｂｕｓ）２１と、送受信回路（送受信部）２３と、Ｂモード回路（Ｂモード情報生成部）２５と、ドプラモード回路（ドプラ情報生成部）２７と、画像生成回路（画像生成部）２９と、処理回路（処理部）３１と、通信インタフェース回路３３と、記憶回路（記憶部）３５と、制御回路（制御部）３７とを有する。
【００１９】
バス２１は、例えば、入力インタフェース回路１５と、送受信回路２３と、Ｂモード回路２５と、ドプラモード回路２７と、画像生成回路２９と、処理回路３１と、通信インタフェース回路３３と、記憶回路３５と、制御回路３７とを電気的に接続する信号路である。
【００２０】
送受信回路２３は、制御回路３７による制御のもとで、超音波プローブ１１における複数の圧電振動子各々に駆動信号を供給する。送受信回路２３は、各圧電振動子によって発生された受信エコー信号に基づいて、受信信号を発生する。
【００２１】
具体的には、送受信回路２３は、いずれも図示していないパルス発生器と、送信遅延回路と、パルサ回路と、プリアンプと、アナログディジタル（Ａｎａｌｏｇ  ｔｏ  ｄｉｇｉｔａｌ（以下、Ａ／Ｄと呼ぶ））変換器と、受信遅延回路と、加算器とを有する。
【００２２】
パルス発生器は、所定のレート周波数ｆｒＨｚ（周期：１／ｆｒ秒）で、送信超音波を形成するためのレートパルスを繰り返し発生する。発生されたレートパルスは、チャンネル数に分配され、送信遅延回路に送られる。
【００２３】
送信遅延回路は、複数のチャンネルごとに、送信超音波をビーム状に収束し、かつ送信指向性を決定するために必要な遅延時間（以下、送信遅延時間と呼ぶ）を、各レートパルスに与える。送信超音波の送信方向または送信遅延時間（以下、送信遅延パターンと呼ぶ）は、記憶回路３５に記憶される。記憶回路３５に記憶された送信遅延パターンは、制御回路３７により超音波の送信時に参照される。
【００２４】
パルサ回路は、このレートパルスに基づくタイミングで、超音波プローブ１１の圧電振動子ごとに電圧パルス（駆動信号）を印加する。これにより、超音波ビームが被検体Ｐに送信される。プリアンプは、超音波プローブ１１を介して取り込まれた被検体Ｐからのエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器は、増幅された受信エコー信号をディジタル信号に変換する。
【００２５】
受信遅延回路は、ディジタル信号に変換された受信エコー信号に、受信指向性を決定するために必要な遅延時間（以下、受信遅延時間と呼ぶ）を与える。エコー信号の受信方向または受信遅延時間（以下、受信遅延パターンと呼ぶ）は、記憶回路３５に記憶される。記憶回路３５に記憶された受信遅延パターンは、制御回路３７により超音波の受信時に参照される。
【００２６】
加算器は、遅延時間が与えられた複数のエコー信号を加算する。この加算により、送受信回路２３は、受信指向性に応じた方向からの反射成分を強調した受信信号（ＲＦ（ｒａｄｉｏｆｒｅｑｕｅｎｃｙ）信号ともいう）を生成する。この送信指向性と受信指向性とにより超音波送受信の総合的な指向性が決定される。この総合的な指向性により、超音波ビーム（いわゆる「超音波走査線」）が決まる。
【００２７】
Ｂモード回路２５は、送受信回路２３から出力された受信信号に基づいて、Ｂモードデータを生成する。Ｂモード回路２５は、Ｂモードデータに基づいて、Ｂモードに関するボリュームデータを生成する。Ｂモード回路２５は、いずれも図示していない包絡線検波器、対数変換器などを有する。
【００２８】
包絡線検波器は、送受信回路２３から出力された受信信号に対して包絡線検波を実行する。包絡線検波器は、包絡線検波された信号を、対数変換器に出力する。対数変換器は、包絡線検波された信号に対して対数変換を行い、包絡線検波された信号における弱い信号を相対的に強調する。Ｂモード回路２５は、対数変換器により強調された信号に基づいて、各走査線における深さごとの信号値を生成する。この信号値は、Ｂモードデータと称される。
【００２９】
Ｂモード回路２５は、被走査領域におけるアジマス方向、エレベーション方向、深さ方向（以下レンジ（Ｒａｎｇｅ）方向と呼ぶ）にそれぞれ対応付けて配列された複数の信号値に基づいて、ボリュームデータを生成する。レンジ方向とは、走査線上の深さ方向（超音波ビーム方向）である。なお、ボリュームデータは、複数の画素値または複数の輝度値などを、走査線に沿って、アジマス方向、エレベーション方向、レンジ方向にそれぞれ対応付けて配列させたデータであってもよい。なお、Ｂモード回路２５は、複数の走査面各々に対応するＢモードデータを用いた所定の補間処理により、ボリュームデータを生成してもよい。
【００３０】
図２は、ボリュームデータの生成に関する撮像ボリューム（被走査領域）ＩＶの一例を示す図である。被走査領域は、例えば、対象組織である心臓を含む空間である。図２における座標系において、ｚ方向は、超音波プローブ１１における超音波送受信面の中心を通るレンジ方向に対応する。図２におけるｘ方向とｙ方向とは、ｚ方向に直交し、互いに直交する２方向である。図２におけるｘ方向とｙ方向とは、例えば、アジマス方向とエレベーション方向とにそれぞれ対応する。撮像ボリュームＩＶは、被検体Ｐにおける心臓に対応する領域であるものとする。なお、撮像ボリュームＩＶは、心臓に限定されず、任意の臓器であってもよい。
【００３１】
ドプラモード回路２７は、送受信回路２３から出力された受信信号に基づいて、ドプラデータを生成する。ドプラモード回路２７は、いずれも図示していないミキサー、低域通過フィルタ（Ｌｏｗ  Ｐａｓｓ  Ｆｉｌｔｅｒ：以下、ＬＰＦと呼ぶ）、速度／分散／Ｐｏｗｅｒ演算デバイス等を有する。
【００３２】
ミキサーは、送受信回路２３から出力された受信信号に、送信周波数と同じ周波数ｆ０を有する基準信号を掛け合わせる。ミキサーは、この掛け合わせにより、ドプラ偏移周波数ｆｄの成分の信号と、（２ｆ０＋ｆｄ）の周波数成分を有する信号とを生成する。ＬＰＦは、ミキサーからの２種の周波数成分を有する信号のうち、高い周波数成分（２ｆ０＋ｆｄ）の信号を取り除く。ドプラモード回路２７は、高い周波数成分（２ｆ０＋ｆｄ）の信号を取り除くことにより、ドプラ偏移周波数ｆｄの成分を有するドプラ信号を生成する。
【００３３】
なお、ドプラモード回路２７は、ドプラ信号を生成するために、直交検波方式を用いてもよい。このとき、受信信号（ＲＦ信号）は、直交検波されＩＱ信号に変換される。ドプラモード回路２７は、ＩＱ信号を複素フーリエ変換することにより、ドプラ偏移周波数ｆｄの成分を有するドプラ信号を生成する。ドプラ信号は、例えば、血流、組織、造影剤によるドプラ成分である。
【００３４】
速度／分散／Ｐｏｗｅｒ演算デバイスは、いずれも図示していないＭＴＩ（Ｍｏｖｉｎｇ  Ｔａｒｇｅｔ  Ｉｎｄｉｃａｔｏｒ）フィルタ、ＬＰＦフィルタ、自己相関演算器等を有する。なお、自己相関演算器の代わりに相互相関演算器を有していてもよい。ＭＴＩフィルタは、生成されたドプラ信号に対して、臓器の呼吸性移動や拍動性移動などに起因するドプラ成分（クラッタ成分）を除去する。
【００３５】
ＭＴＩフィルタは、ドプラ信号から血流に関するドプラ成分（以下、血流ドプラ成分と呼ぶ）を抽出するために用いられる。ＬＰＦは、ドプラ信号から組織の移動に関するドプラ成分（以下、組織ドプラ成分と呼ぶ）を抽出するために用いられる。
【００３６】
自己相関演算器は、血流ドプラ成分および組織ドプラ成分に対して自己相関値を計算する。自己相関演算器は、計算された自己相関値に基づいて、血流および組織の平均速度値、分散値、ドプラ信号の反射強度（パワー）等を算出する。速度／分散／Ｐｏｗｅｒ演算デバイスは、複数のドプラ信号に基づく血流および組織の平均速度値、分散値、ドプラ信号の反射強度等に基づいて、被走査領域の各位置におけるカラードプラデータを生成する。
【００３７】
以下、ドプラ信号とカラードプラデータとをまとめて、ドプラデータと呼ぶ。ドプラデータは、処理回路３１等の各回路に出力される。ドプラデータとボリュームデータとは、例えば、心電波形における心位相の時相に基づいて対応付けられる。
【００３８】
撮像ボリュームＩＶに対応する被走査領域が心臓である場合、血流ドプラ成分は、心腔内の血液（流体）のドプラ情報（以下、血流ドプラ情報と呼ぶ）に対応する。すなわち、ドプラモード回路２７は、超音波プローブ１１からの出力に基づいて、心腔内の血液に関する血流ドプラ情報を生成する。血流ドプラ情報は、心腔内において走査線方向（超音波ビーム方向）に沿った複数の位置における血流の速さに相当する。血流の速さとは、血液の移動方向および血液の移動速度に対応する。血流ドプラ情報は、処理回路３１に出力される。
【００３９】
画像生成回路２９は、いずれも図示していないプロセッサ（ＣＰＵ（ｃｅｎｔｒａｌ  ｐｒｏｃｅｓｓｉｎｇ  ｕｎｉｔ））と、メモリと、ディジタルスキャンコンバータ（Ｄｉｇｉｔａｌ  Ｓｃａｎ  Ｃｏｎｖｅｒｔｅｒ：以下ＤＳＣと呼ぶ）と、画像メモリとを有する。画像生成回路２９は、画像生成に関する各種プログラムを記憶回路３５から読み出す。画像生成回路２９は、読み出した各種プログラムを自身のメモリに展開し、実行することで、各種プログラムに対応する機能を実現する。換言すると、各種プログラムを読み出した状態の画像生成回路２９は、図１に示すように、画像生成機能２９１、運動情報生成機能２９３等の各機能を有する。
【００４０】
具体的には、画像生成回路２９は、画像生成機能２９１に関する画像生成プログラムを記憶回路３５から読み出す。画像生成回路２９は、読み出した画像生成プログラムを自身のメモリに展開し、実行することで、画像生成機能２９１を実現する。このとき、画像生成回路２９は、画像生成部として機能する。
【００４１】
画像生成機能２９１は、ＤＳＣと、画像メモリとを制御する。画像生成機能２９１は、ＤＳＣに対して、座標変換処理（リサンプリング）等を実行する。座標変換処理とは、例えば、Ｂモードデータおよびドプラデータからなる超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換する処理である。画像生成機能２９１は、座標変換処理により、表示画像としての超音波画像を生成する。
【００４２】
具体的には、画像生成機能２９１は、Ｂモードデータに基づいて、Ｂモード画像（以下、２次元Ｂモード画像と呼ぶ）を生成する。画像生成機能２９１は、カラードプラデータに基づいて、平均速度画像、分散画像などのカラー血流画像と、ドプラ信号の反射強度を示すパワー血流画像などを生成する。画像生成機能２９１は、ボリュームデータに基づいて、３次元的なＢモード画像（以下、３次元Ｂモード画像と呼ぶ）を生成する。
【００４３】
画像生成機能２９１は、時系列に沿った一連のボリュームデータの生成に応じて、時系列に沿った一連の３次元Ｂモード画像を生成する。画像生成機能２９１は、例えば、少なくとも１心拍において、時系列に沿った一連の３次元Ｂモード画像を生成する。画像生成機能２９１は、１心拍において、１心拍の時間にフレームレートを乗じた数に対応する複数の３次元Ｂモード画像を生成する。
【００４４】
画像生成機能２９１は、ボリュームデータに対してレンダリング処理を実行することにより、レンダリング画像を生成する。レンダリング処理に関する種々のパラメータ（レンダリングパラメータ：例えば、視点、視線（レンダリング方向）、不透明度など）は、入力インタフェース回路１５を介して入力される。すなわち、画像生成機能２９１は、入力されたレンダリングパラメータを用いて、ボリュームデータに対してレンダリング処理を実行する。
【００４５】
画像生成機能２９１は、断面変換（ｍｕｌｔｉ−ｐｌａｎａｒ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ：以下、ＭＰＲと呼ぶ）処理により、少なくとも一つのＭＰＲ画像を生成する。ＭＰＲ処理に関する種々のパラメータ（ＭＰＲパラメータ：例えば、断面位置、断面方向など）は、入力インタフェース回路１５を介して入力される。すなわち、画像生成機能２９１は、入力されたＭＰＲパラメータを用いて、ボリュームデータに対してＭＰＲ処理を実行する。
【００４６】
撮像ボリュームＩＶが心臓である場合、ＭＰＲ画像は、例えば、４つの内腔（右心房、右心室、左心房、左心室：以下、４ｃｈ（ｃｈａｍｂｅｒ）と呼ぶ）を含む断面（以下、４ｃｈビューと呼ぶ）と、４ｃｈのうち２つの内腔（以下、２ｃｈ（ｃｈａｍｂｅｒ）と呼ぶ）を含む断面（以下、２ｃｈビューと呼ぶ）と、４ｃｈビューと２ｃｈビューとに垂直な断面（以下、垂直断面と呼ぶ）とにそれぞれ対応する３断面画像（以下、直交３断面画像と呼ぶ）である。なお、ＭＰＲ画像は、上記直交３断面画像に限定されず、任意の断面に対応する断面画像でもよい。
【００４７】
画像生成機能２９１は、処理回路３１により生成されたベクトル場画像を、レンダリング画像に、それぞれの心位相の時相と位置とを略一致させて重畳する。この重畳により、画像生成機能２９１は、ベクトル場画像をレンダリング画像に重畳した第１重畳画像を生成する。このとき、ベクトル場画像は、後述する３次元のベクトル情報に対応する３次元ベクトルデータをレンダリングした画像に対応する。ベクトル場画像については、処理回路３１において詳述する。
【００４８】
画像生成機能２９１は、ＭＰＲ画像の断面に対応するベクトル場画像をＭＰＲ画像にそれぞれの心位相の時相と断面位置とを略一致させて重畳する。この重畳により、画像生成機能２９１は、ベクトル場画像をＭＰＲ画像に重畳した第２重畳画像を生成する。ＭＰＲ画像の断面に対応するベクトル場画像の生成については、処理回路３１において詳述する。
【００４９】
画像生成機能２９１は、Ｂモード画像にカラー血流画像を重畳させた重畳画像を生成する。画像生成機能２９１は、Ｂモード画像にパワー血流画像を重畳させた重畳画像を生成する。Ｂモード画像にカラー血流画像を重畳させた重畳画像と、Ｂモード画像にパワー血流画像を重畳させた重畳画像とをまとめてカラードプラ画像と呼ぶ。
【００５０】
画像生成機能２９１は、ドプラ信号に基づいて、ドプラ波形を示すドプラ波形画像を生成する。画像生成機能２９１は、カラードプラ画像とドプラ波形画像とを有する超音波画像を生成する。なお、超音波画像は、Ｂモード画像とドプラ波形画像とを有していてもよい。
【００５１】
画像生成機能２９１は、Ｂモード画像、カラードプラ画像、ドプラ波形画像、第１重畳画像、第２重畳画像、ベクトル場画像などに、種々のパラメータ、文字情報、目盛等を合成する。画像生成機能２９１は、種々のパラメータ、文字情報、目盛等を合成したＢモード画像、ドプラ画像、カラードプラ画像、ドプラ波形画像、第１重畳画像、第２重畳画像、ベクトル場画像などを表示回路１７に出力する。
【００５２】
画像メモリは、生成された各種画像に対応するデータ（以下、画像データと呼ぶ）を記憶する。画像メモリに記憶された画像データは、入力インタフェース回路１５を介した操作者の指示により、読み出される。画像メモリは、例えば、フリーズする直前の一連のフレームに対応する各種画像を保存するメモリである。画像メモリに記憶されている画像を、表示回路１７に連続表示（シネ表示）させることで、各種画像の動画像が表示回路１７に表示される。
【００５３】
画像生成回路２９は、運動情報生成機能２９３に関する運動情報生成プログラムを記憶回路３５から読み出す。画像生成回路２９は、読み出した運動情報生成プログラムを自身のメモリに展開し、実行することで、運動情報生成機能２９３を実現する。このとき、画像生成回路２９は、運動情報生成部として機能する。
【００５４】
運動情報生成機能２９３は、超音波プローブ１１からの出力に基づいて、心腔に関する心臓の運動情報を生成する。運動情報生成機能２９３は、生成した運動情報を処理回路３１に出力する。心臓の運動情報とは、例えば、心臓における心筋（または心壁）の各位置における３次元的な速さを示す情報である。具体的には、心臓の運動情報は、心臓壁を含む心筋に対応する３次元的な領域（以下、心筋領域と呼ぶ）の各位置において、心筋の運動に伴う３次元的な心筋の速度ベクトル（以下、３次元運動ベクトルと呼ぶ）に対応する。
【００５５】
すなわち、３次元運動ベクトルは、心筋の運動（ねじれ運動）を示す３次元ベクトルである。心臓の運動情報は、３次元的な心筋領域内の複数の位置にそれぞれ対応する複数の３次元運動ベクトルと、心位相（隣接する２つのＲ波間における位置）とを有する。以下、運動情報の生成について詳述する。
【００５６】
運動情報生成機能２９３は、一連の３次元Ｂモード画像各々におけるスペックルパターン（ｓｐｅｃｋｌｅ  ｐａｔｔｅｒｎ）を追跡すること（トラッキング：ｔｒａｃｋｉｎｇ）により、運動情報を生成する。具体的には、運動情報生成機能２９３は、一連の３次元Ｂモード画像のうち時間的に隣接する２つの３次元Ｂモード画像（換言すれば、互いに隣接する２フレームにそれぞれ対応する２つの３次元Ｂモード画像）各々におけるスペックルパターンに対して、パターンマッチングを実行する。
【００５７】
このパターンマッチングにより、フレーム間において、スペックルパターンがトラッキングされ、心筋の動きが追跡される。運動情報生成機能２９３は、上記処理をフレームごとに繰り返すことで、時間とともに変化する心筋（局所組織）の位置、すなわち心筋の移動を示す運動情報を追跡することができる。
【００５８】
なお、運動情報生成機能２９３は、時系列に沿った一連の２次元的なＢモード画像におけるスペックルパターンを用いて、運動情報を生成してもよい。また、運動情報生成機能２９３は、心筋に関する組織ドプラ成分を有するドプラデータ（以下、組織ドプラ情報）を用いて、運動情報を生成してもよい。
【００５９】
上記説明においては、単一の画像生成回路２９において画像生成機能及び運動情報生成機能等の各種機能が実行されるものとして説明したが、複数の独立したプロセッサを組み合わせて画像生成回路２９を構成し、各プロセッサがプログラムを実行することにより各種機能を実現するものとしても構わない。また、画像生成機能２９１と運動情報生成機能２９３とは、それぞれ異なる回路で実現されてもよいし、処理回路３１、制御回路３７等に組み込まれてもよい。
【００６０】
処理回路（Ｐｒｏｃｅｓｓｉｎｇ  ｃｉｒｃｕｉｔｒｙ）３１は、いずれも図示していないＣＰＵとメモリを備える。処理回路３１は、心腔内の３次元領域の各位置における流体（血液）の速さに関する３次元のベクトル情報の生成に関する各種プログラムを自身のメモリに展開し、実行することで、各種プログラムに対応する機能を実現する。換言すると、各種プログラムを読み出した状態の処理回路３１は、図１に示すように、運動ベクトル分解機能３１１と、血流ベクトル計算機能３１３と、ベクトル場画像生成機能３１５等の各機能を有する。
【００６１】
処理回路３１は、上記読み出した各種プログラムを実行することにより、ドプラ情報と運動情報とに基づいて、心腔内の３次元領域の各位置における流体（血液）の速さに関する３次元のベクトル情報（３次元ベクトルデータ）を生成する。生成したベクトル情報は、撮像ボリュームＩＶにおける心腔領域に対応する３次元領域の各位置において、流体の移動方向と移動速度とを有する。以下、３次元ベクトルデータの生成のための処理内容について詳述する。
【００６２】
処理回路３１は、運動ベクトル分解機能３１１に関する運動ベクトル分解プログラムを記憶回路３５から読み出す。処理回路３１は、読み出した運動ベクトル分解プログラムを自身のメモリに展開し、実行することで、運動ベクトル分解機能３１１を実現する。このとき、処理回路３１は、運動ベクトル分解部として機能する。
【００６３】
運動ベクトル分解機能３１１は、撮像ボリュームＩＶに対応するボリュームデータにおいて、走査線方向（超音波ビーム方向）に直交し、かつ互いに直交する第１方向と第２方向とを特定する。第１方向とは、例えば、走査線方向に直交し、エレべーション方向に沿った方向である。第２方向とは、走査線方向に直交し、アジマス方向に沿った方向である。
【００６４】
例えば、超音波プローブ１１が１次元アレイプローブである場合、第１方向は、走査線方向に直交し、かつボリュームデータの収集時における複数の走査面各々を規定する方向である。すなわち、撮像ボリュームＩＶにおける複数の走査面各々は、走査線方向と第１方向とにより規定される。第２方向は、走査線方向に直交し、かつボリュームデータの収集時において走査面を揺動させる方向である。
【００６５】
運動ベクトル分解機能３１１は、撮像ボリュームＩＶに対応する３次元領域の各位置における３次元運動ベクトルを、第１方向と第２方向と走査線方向とにそれぞれ投影する。この投影により、運動ベクトル分解機能３１１は、３次元領域の各位置において、３次元運動ベクトルを、第１方向に沿った第１成分と、第２方向に沿った第２成分と、走査線方向に沿った成分とに分解する。
【００６６】
図３は、撮像ボリュームＩＶのｘｙ断面における３次元運動ベクトルの分解の一例を示す図である。図３に示すＭＹは、心筋領域を示している。図に示す矢印ＭＶ１、ＭＶ２は、ｘｙ平面に投影された３次元運動ベクトルを示している。図３に示すように、ｘｙ平面に投影された３次元運動ベクトルＭＶ２は、ｘ方向（第１方向）成分ＭＶｘと、ｙ方向（第２方向）成分ＭＶｙとに分解される。
【００６７】
なお、第１方向および第２方向は、超音波プローブ１１における超音波送受信面の中心と撮像ボリュームＩＶの下面の中心とを結ぶ方向（以下、中心軸方向と呼ぶ）に直交し、かつ互いに直交する２方向であってもよい。例えば、被検体Ｐの体表面に超音波プローブ１１が鉛直方向に平行に沿って当接されている場合、中心軸方向は、鉛直方向に平行となる。
【００６８】
また、第１方向と第２方向とが、中心軸方向に対して直交し、かつ互いに直交する場合、中心軸方向と、第１方向と、第２方向とは、直交３軸を構成する。このとき、血流ベクトル計算機能３１３で用いられる血流ドプラ情報には、心腔内において中心軸方向に沿った複数の位置における血流の速さに相当する。すなわち、血流ドプラ情報は、心腔に対応する領域内の複数の位置各々において、中心軸方向に平行な走査線方向を除いて、ドプラモード回路２７において生成された血流ドプラ成分を、中心軸方向に投影したドプラデータとなる。
【００６９】
処理回路３１は、血流ベクトル計算機能３１３に関する血流ベクトル計算プログラムを記憶回路３５から読み出す。処理回路３１は、読み出した血流ベクトル計算プログラムを自身のメモリに展開し、実行することで、血流ベクトル計算機能３１３を実現する。このとき、処理回路３１は、血流ベクトル計算部として機能する。以下、血流ベクトル計算機能３１３は、心腔内の血液の流れを計算する機能として説明する。なお、血流ベクトル計算機能３１３における計算対象である血液の流れは、心腔内に限定されず、他の臓器における血液の流れであってもよい。
【００７０】
血流ベクトル計算機能３１３は、第１方向と走査線方向とにより規定され第２方向に沿った複数の第１断面を特定する。血流ベクトル計算機能３１３は、第２方向と走査線方向とにより規定され第１方向に沿った複数の第２断面を特定する。例えば、超音波プローブ１１が１次元アレイプローブである場合、複数の第１断面は、走査面の揺動方向に直交する複数の走査面に対応する。また、このとき、複数の第２断面は、揺動方向に平行であって、複数の走査面に直交する面となる。なお、第１断面と第２断面とは、中心軸方向に垂直な面に直交し、かつ互いに直交する面であってもよい。
【００７１】
複数の第１断面各々は、後述する第１の２次元血流ベクトルが計算される複数の部分領域（微小領域）に分割されているものとする。加えて、複数の第２断面各々は、後述する第２の２次元血流ベクトルが計算される複数の部分領域に分割されているものとする。複数の部分領域各々において、血液の流量に関して質量保存則が成り立つものと仮定する。
【００７２】
より詳細には、複数の部分領域各々において、部分領域に流入および流出する血液（流量）に関して連続の式が成り立っているものとする。なお、部分領域内において、血流（ベクトルｕ）の沸き出しおよび吸い込みはないもの（ｄｉｖ  ｕ＝０）とする。なお、質量保存則の代わりに、運動量保存則が用いられてもよい。
【００７３】
具体的には、血流ベクトル計算機能３１３は、被検体Ｐの心腔領域内の複数の第１断面各々を、例えば、Ｍ行Ｎ列（Ｍ、Ｎは任意の自然数）に格子状に分割する。血流ベクトル計算機能３１３は、被検体Ｐの心腔領域内の複数の第２断面各々を、例えば、Ｍ行Ｌ列（Ｌは任意の自然数）に格子状に分割する。なお、心腔領域を分割した部分領域の形状は、任意の多角形であってもよい。
【００７４】
血流ベクトル計算機能３１３は、複数の第１断面（ｘ−ｚ平面）各々において、分解された３次元運動ベクトルの第１成分とドプラ情報とを用いて、第１断面内の各位置における流体（血液）の速さに関する第１の２次元血流ベクトルを計算する。すなわち、血流ベクトル計算機能３１３は、第１方向に沿った心筋の運動を示すベクトル（第１成分）と走査線方向に沿ったドプラ情報とを用いて、複数の第１断面各々における各位置における血液の速度と移動方向とを示す２次元血流ベクトルを計算する。以下、２次元血流ベクトルの計算について、具体的に説明する。
【００７５】
図４は、第１断面における複数の部分領域各々における第１の２次元血流ベクトルの計算過程の一例を示す図である。図４におけるＭＹは、被検体Ｐの心筋領域を示している。図４におけるＥｘ０１、Ｅｘ０２、Ｅｘ０３、Ｅｘ１１、Ｅｘ１２、Ｅｘ１３、Ｅｘ２１、Ｅｘ２２、Ｅｘ２３は、第１断面における複数の部分領域の一部分を示している。以下、Ｅｘ１１、Ｅｘ１２各々において、第１方向に沿った血流ベクトルの計算について、図４を参照して説明する。
【００７６】
以下説明を簡単にするために、第１方向（図４のｘ方向）に平行な向きの符号と、走査線方向（図４のｚ方向）に平行な向きの符号とを正（プラス）として、説明する。なお、第１方向（図４のｘ方向）に反平行な向きの符号と、走査線方向（図４のｚ方向）に反平行な向きの符号とは、負（マイナス）となる。
【００７７】
血流ベクトル計算機能３１３は、部分領域Ｅｘ１１において、運動ベクトル分解機能３１１によって分解された３次元運動ベクトルの第１成分Ｅｘ１１ｘｍ（心壁ＭＹの速度）と、走査線方向に沿った２つの血流ドプラ成分Ｅｘ１１ｚｕ、Ｅｘ１１ｚｄと、質量保存則とを用いて、第１方向に沿ったベクトルＥｘ１１ｘｃを計算する。
【００７８】
具体的には、図４における部分領域Ｅｘ１１と心壁ＭＹとの境界から、部分領域Ｅｘ１１には、Ｅｘ１１ｘｍに関する血流が流入する。すなわち、部分領域Ｅｘ１１において、心壁の動きＭＹＭにより心壁ＭＹを介して、第１成分Ｅｘ１１ｘｍに対応する運動量が、部分領域Ｅｘ１１に流入する。
【００７９】
また、Ｅｘ１１ｚｄ（走査線に沿った血流速度）に関する血流が、図４における部分領域Ｅｘ１１の下面から流入する。すなわち、血流成分Ｅｘ１１ｚｄに対応する運動量が、部分領域Ｅｘ１１の下面から部分領域Ｅｘ１１に流入する。一方、Ｅｘ１１ｚｕ（走査線に沿った血流速度）に関する血流が、図４における部分領域Ｅｘ１１の上面から部分領域Ｅｘ０１に流出する。すなわち、血流成分Ｅｘ１１ｚｕに対応する運動量が、部分領域Ｅｘ１１の上面から部分領域Ｅｘ０１に流出する。
【００８０】
部分領域Ｅｘ１１において、第１成分Ｅｘ１１ｘｍ（心壁ＭＹの速度）と、走査線方向に沿った２つの血流ドプラ成分Ｅｘ１１ｚｕ、Ｅｘ１１ｚｄとは既知である。このため、血流ベクトル計算機能３１３は、部分領域Ｅｘ１１から部分領域Ｅｘ１２へ流出する血流Ｅｘ１１ｘ１を、以下に示す質量保存則（連続の式）により計算する。  
Ｅｘ１１ｘｍ＋Ｅｘ１１ｚｄ＝Ｅｘ１１ｚｕ＋Ｅｘ１１ｘ１  
Ｅｘ１１ｘ１＝Ｅｘ１１ｘｍ＋Ｅｘ１１ｚｄ−Ｅｘ１１ｚｕ
【００８１】
次いで、血流ベクトル計算機能３１３は、部分領域Ｅｘ１１に隣接する部分領域Ｅｘ１２について、部分領域Ｅｘ１２から部分領域Ｅｘ１３へ流出する血流Ｅｘ１１ｘ２を、以下に示す質量保存則（連続の式）により計算する。  
Ｅｘ１１ｘ１＋Ｅｘ１２ｚｄ＝Ｅｘ１２ｚｕ＋Ｅｘ１１ｘ２  
Ｅｘ１１ｘ２＝Ｅｘ１１ｘ１＋Ｅｘ１２ｚｄ−Ｅｘ１２ｚｕ
【００８２】
同様にして、血流ベクトル計算機能３１３は、複数の第１断面各々の各部分領域（各位置）において、心筋領域に隣接する部分領域から第１方向に沿って心腔内部の領域に向けて、第１方向に沿った血流量および血流の向きを、逐次的に計算する。血流ベクトル計算機能３１３は、複数の第１断面各々の各部分領域（各位置）から流出する血流量および血流の向きを用いて、複数の第１断面各々の各部分領域（各位置）における第１の２次元血流ベクトルを計算する。第１の２次元血流ベクトルは、第１断面内の各位置における血流の速度および方向を２次元的に示すベクトルである。
【００８３】
図５は、一つの第１断面内の各位置における第１の２次元血流ベクトルの一例を示す図である。図５において、複数の部分領域から構成される立方体は、心腔領域に対応する。図５におけるｘｚ平面に垂直な矢印Ｘ１、Ｘ２は、複数の３次元運動ベクトルをそれぞれ分解した複数の第１成分を模式的に示している。図５におけるｘｙ平面に垂直な矢印Ｚ１、Ｚ２は、走査線方向に沿った複数の血流ドプラ成分を模式的に示している。図５におけるｘｚ平面における矢印Ａ１は、血流ベクトル計算機能３１３により計算された第１の２次元血流ベクトルを示している。
【００８４】
血流ベクトル計算機能３１３は、複数の第２断面各々の各部分領域（各位置）において、心筋領域に隣接する部分領域から第２方向に沿って心腔内部の領域に向けて、第２方向に沿った血流量および血流の向きを、分解された３次元運動ベクトルの第２成分とドプラ情報とを用いて、逐次的に計算する。血流ベクトル計算機能３１３は、複数の第２断面各々の各部分領域（各位置）から流出する血流量および血流の向きを用いて、複数の第２断面各々の各部分領域（各位置）における第２の２次元血流ベクトルを計算する。第２の２次元血流ベクトルは、第２断面内の各位置における血流の速度および方向を２次元的に示すベクトルである。
【００８５】
具体的な計算手法は、第１の２次元血流ベクトルの導出に関する計算と実質的に同様であり、例えば、図４に関する説明において、「第１方向」、「第１断面」、「ｘ」を、「第２方向」、「第２断面」、「ｙ」にそれぞれ置換することにより、説明される。
【００８６】
図６は、一つの第２断面内の各位置における第２の２次元血流ベクトルの一例を示す図である。図６において、複数の部分領域から構成される立方体は、心腔領域に対応する。図６におけるｙｚ平面に垂直な矢印Ｙ１、Ｙ２は、心壁における複数の３次元運動ベクトルをそれぞれ分解した複数の第２成分を模式的に示している。図６におけるｘｙ平面に垂直な矢印Ｚ１、Ｚ２は、走査線方向に沿った複数の血流ドプラ成分を模式的に示している。図６におけるｙｚ平面における矢印Ｂ１は、計算された第２の２次元血流ベクトルを示している。
【００８７】
なお、血流ベクトル計算機能３１３は、心腔領域の上端及び下端（心尖部近傍）に隣接する部分領域に関して、走査線方向に沿って分解された３次元運動ベクトルの成分（以下、走査線方向成分とよぶ）または組織ドプラ成分を用いて、第１の２次元血流ベクトルおよび第２の２次元血流ベクトルを計算してもよい。
【００８８】
また、第１の２次元血流ベクトルの計算の開始箇所における心壁に略対向する心壁（第１方向に沿った計算の終了位置の部分領域に隣接する心壁）に関する第１成分（図５のＸ１またはＸ２）は、第１の２次元血流ベクトルの計算における整合性および計算精度を向上させるために、例えば、丸め誤差を低減させるために利用されてもよい。
【００８９】
また、第２の２次元血流ベクトルの計算の開始時点における心壁に略対向する心壁（第２方向に沿った計算の終了位置の部分領域に隣接する心壁）に関する第２成分（図６のＹ１またはＹ２）は、第２の２次元血流ベクトルの計算における整合性および計算精度を向上させるために、例えば、丸め誤差を低減させるために利用されてもよい。
【００９０】
血流ベクトル計算機能３１３は、第１の２次元血流ベクトルと第２の２次元血流ベクトルとを合成することにより、心腔内の３次元領域の各位置における流体（血液）の速さに関する３次元のベクトル情報（３次元ベクトルデータ）を生成する。具体的には、血流ベクトル計算機能３１３は、隣接する２つの第１断面において略対向する２つの部分領域と、これら２つの部分領域に接続する２つの第２断面における２つの部分領域とを、心腔内の３次元領域において特定する。
【００９１】
すなわち、血流ベクトル計算機能３１３は、１ボクセルの周囲の４つの部分領域（以下、同一局所領域と呼ぶ）を特定する。同一局所領域のうち対向する２つの部分領域は、隣接する２つの第１断面に属する。同一局所領域のうち対向する他方の２つの部分領域は、隣接する２つの第２断面に属する。
【００９２】
血流ベクトル計算機能３１３は、第１断面に属する２つの部分領域における第１の２次元血流ベクトルと、第２断面に属する２つの部分領域における第２の２次元血流ベクトルとを合成する。例えば、血流ベクトル計算機能３１３は、同一局所領域に属する４つのベクトル（隣接する２つの第１の２次元血流ベクトルと隣接する２つの第２の２次元血流ベクトル）各々の成分（第１方向に沿った成分、第２方向に沿った成分、走査線方向に沿った成分）ごとに平均化することにより、４つのベクトルを合成してもよい。
【００９３】
以上の処理により、血流ベクトル計算機能３１３は、互いに対向しかつ隣接する４つの部分領域により規定される重心位置において、３次元血流ベクトルを生成する。なお、４つの部分領域のうち少なくとも一つの部分領域が心壁に対応する場合、血流ベクトル計算機能３１３は、心壁に対応する部分領域を除いた複数の部分領域に対応する第１の２次元血流ベクトルと第２の２次元血流ベクトルとを合成することにより、３次元血流ベクトルを生成する。
【００９４】
第１断面における部分領域の数がＭ×Ｎであって、第２断面における部分領域の数がＭ×Ｌであって、心腔領域が直方体である場合、重心位置の総数はＬ×Ｍ×Ｎ＋２×（Ｌ×Ｍ＋Ｍ×Ｎ＋Ｎ×Ｌ）となる。すなわち、血流ベクトル計算機能３１３は、複数の重心位置各々において、３次元血流ベクトルを生成する。これにより、血流ベクトル計算機能３１３は、心腔内の３次元領域の各位置における血液（流体）の速さに関する３次元のベクトル情報を生成する。３次元のベクトル情報は、心腔内の３次元領域の各位置における血液（流体）の移動方向と移動速度とを有する。
【００９５】
血流ベクトル計算機能３１３は、複数の重心位置にそれぞれ対応する複数の３次元血流ベクトルを、３次元のベクトル情報として、ベクトル場画像生成機能３１５に出力する。血流ベクトル計算機能３１３は、３次元血流ベクトルの生成に関連する各種データの時相（心位相の時相）を、３次元のベクトル情報に付帯させる。なお、３次元のベクトル情報は、３次元血流ベクトルの生成に関連する各種データの取得に応じて、リアルタイムに生成されてもよい。
【００９６】
処理回路３１は、ベクトル場画像生成機能３１５に関するベクトル場画像生成プログラムを記憶回路３５から読み出す。処理回路３１は、読み出したベクトル場画像生成プログラムを自身のメモリに展開し、実行することで、ベクトル場画像生成機能３１５を実現する。このとき、処理回路３１は、ベクトル場画像生成部として機能する。
【００９７】
ベクトル場画像生成機能３１５は、３次元のベクトル情報の表示に関する種々のパラメータ（以下、ベクトル表示パラメータと呼ぶ）を、３次元のベクトル情報に基づいて決定する。ベクトル表示パラメータとは、３次元血流ベクトルの表示に関するパラメータである。３次元血流ベクトルの表示とは、例えば、マーカである。マーカとは例えば、矢印である。
【００９８】
具体的には、ベクトル場画像生成機能３１５は、３次元のベクトル情報における複数の３次元血流ベクトル各々に対応するマーカの大きさ（長さ）を、血液の移動速度（３次元血流ベクトルの大きさ）に応じて決定する。このとき、マーカの大きさは、血液の移動速度に対応する。ベクトル場画像生成機能３１５は、マーカの色相と、明度と、彩度とのうち少なくとも１つを、血液の移動速度に応じて決定する。このとき、マーカの色相と、明度と、彩度は、血液の移動速度に対応する。
【００９９】
例えば、ベクトル場画像生成機能３１５は、血液の移動速度に対するベクトル表示パラメータの対応表を、記憶回路３５から読み出す。ベクトル場画像生成機能３１５は、読み出した対応表と、３次元のベクトル情報における血液の移動速度の大きさとを用いて、腔内の３次元領域の各位置におけるベクトル表示パラメータを決定する。
【０１００】
ベクトル場画像生成機能３１５は、決定されたベクトル表示パラメータと３次元のベクトル情報とに対応するマーカを心腔内の３次元領域の各位置に対応付けて表示可能なベクトル場画像を生成する。ベクトル場画像生成機能３１５は、生成したベクトル場画像を、画像生成回路２９に出力する。以下、ベクトル場画像の生成について詳述する。
【０１０１】
ベクトル場画像生成機能３１５は、レンダリングパラメータを用いて、３次元のベクトル情報における３次元血流ベクトルに対してレンダリング処理を実行する。ベクトル場画像生成機能３１５は、３次元血流ベクトルに対するレンダリング処理の結果と、ベクトル表示パラメータとを用いて、レンダリングベクトル場画像を生成する。レンダリングベクトル場画像とは、ベクトル表示パラメータに従って３次元的に血流ベクトルを表示可能なレンダリング画像である。
【０１０２】
なお、レンダリングパラメータは、レンダリング画像の生成に関するレンダリングパラメータと略同一であってもよい。このとき、ベクトル場画像生成機能３１５は、レンダリング画像の生成に関するレンダリング処理を３次元血流ベクトルに適用することにより、レンダリング処理を実行する。
【０１０３】
ベクトル場画像生成機能３１５は、ＭＰＲパラメータを用いて、３次元のベクトル情報における３次元血流ベクトルに対してＭＰＲ処理を実行する。具体的には、ベクトル場画像生成機能３１５は、ＭＰＲ処理に関する断面に応じて、３次元血流ベクトルの成分を計算する。ベクトル場画像生成機能３１５は、３次元血流ベクトルに対するＭＰＲ処理の結果と、ベクトル表示パラメータとを用いて、ＭＰＲベクトル場画像を生成する。３次元血流ベクトルに対するＭＰＲ処理の結果は、ＭＰＲ処理に関する断面上の各位置における３次元血流ベクトルの成分である。
【０１０４】
なお、ＭＰＲパラメータは、ＭＰＲ画像の生成に関するＭＰＲパラメータと略同一であってもよい。このとき、ベクトル場画像生成機能３１５は、ＭＰＲ画像の生成に関するＭＰＲ処理を３次元血流ベクトルに適用することにより、ＭＰＲ処理を実行する。
【０１０５】
上記説明においては、単一の処理回路３１において運動ベクトル分解機能３１１、血流ベクトル計算機能３１３、及びベクトル場画像生成機能３１５等の各種機能が実行されるものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路３１を構成し、各プロセッサがプログラムを実行することにより各種機能を実現するものとしても構わない。また、運動ベクトル分解機能３１１と、血流ベクトル計算機能３１３と、ベクトル場画像生成機能３１５とは、それぞれ異なる回路で実現されてもよいし、制御回路３７等に組み込まれてもよい。
【０１０６】
通信インタフェース回路３３は、ネットワークを介して他の医用画像診断装置および医用画像保管装置等の各種装置に接続される。医用画像診断装置は、例えば、Ｘ線コンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ：以下、ＣＴと呼ぶ）、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ：以下ＭＲＩと呼ぶ）装置、Ｘ線診断装置などのボリュームデータを生成可能な各種モダリティである。
【０１０７】
通信インタフェース回路３３は、ネットワークを介して被検体Ｐのボリュームデータ等の各種データを、医用画像診断装置または医用画像保管装置から受信する。通信インタフェース回路３３は、受信したボリュームデータ等の各種データを、記憶回路３５に出力する。処理回路３１により処理された各種データは、通信インタフェース回路３３およびネットワークを介して、各種装置へ転送可能である。加えて、通信インタフェース回路３３には、心電計、心音計、脈波計などに代表される不図示の生体信号計測センサが、ネットワークを介して接続されてもよい。
【０１０８】
通信インタフェース回路３３は、心臓の運動情報の生成に関する各種データ、３次元のベクトル情報における３次元血流ベクトルの生成に関する各種データ、ベクトル場画像の生成に関する各種データを、ネットワークを介して他の医用画像診断装置から受信してもよい。このとき、画像生成機能２９１、運動情報生成機能２９３、運動ベクトル分解機能３１１、血流ベクトル計算機能３１３、ベクトル場画像生成機能３１５等の各機能は、受信したデータを用いて各種処理を実行する。
【０１０９】
記憶回路３５は、各種メモリ、ＨＤＤ（ハードディスクドライブ）、ＳＳＤ（ソリッドステートドライブ）、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどにより構成される。記憶回路３５は、通信インタフェース回路３３およびネットワークを介して受信した各種データを記憶する。記憶回路３５は、本超音波診断装置１を統括的に制御する装置制御プログラム、ディスプレイに表示される各種画像の表示制御に関する表示制御プログラム等の各種プログラムを記憶する。
【０１１０】
記憶回路３５は、フォーカス深度の異なる複数の受信遅延パターンおよび複数の送信遅延パターン、診断プロトコル、送受信条件等の各種データ群を記憶する。記憶回路３５は、通信インタフェース回路３３から出力された各種データを記憶する。記憶回路３５は、Ｂモード回路２５により生成されたＢモードデータ、ボリュームデータ等の各種データを記憶する。記憶回路３５は、ドプラモード回路２７により生成されたドプラデータ、血流ドプラ情報、組織ドプラ情報等を記憶する。
【０１１１】
記憶回路３５は、画像生成回路２９において実行される画像生成プログラム（レンダリング処理に関するプログラム、ＭＰＲ処理に関するプログラム等）、運動機能生成プログラム（例えば、スペックルパターンのトラッキングに関するパターンマッチング等のプログラム）を記憶する。記憶回路３５は、レンダリングパラメータ、ＭＰＲパラメータ等の種々のパラメータ、文字情報、目盛等を記憶する。記憶回路３５は、画像生成回路２９により生成された種々の画像（Ｂモード画像、ドプラ画像、カラードプラ画像、ドプラ波形画像、第１重畳画像、第２重畳画像等）のデータ、心臓（心筋）の運動情報（３次元運動ベクトル）を記憶する。
【０１１２】
記憶回路３５は、処理回路２９において実行される運動ベクトル分解プログラム、血流ベクトル計算プログラム、及びベクトル場画像生成プログラムなどの各種プログラムと、ベクトル表示パラメータと、血液の移動速度に対するベクトル表示パラメータの対応表等とを記憶する。記憶回路３５は、処理回路３１により生成される３次元のベクトル情報（３次元血流ベクトル）、ベクトル場画像（レンダリングベクトル場画像、ＭＰＲベクトル場画像）等を記憶する。
【０１１３】
記憶回路３５は、運動情報生成機能２９３により生成された心臓の運動情報において、心臓（心筋）の運動量に応じたカラーマップを記憶する。心臓の運動量とは、例えば、心筋領域の各位置における３次元運動ベクトルの大きさに対応する。
【０１１４】
入力インタフェース回路１５は、操作者からの各種指示・命令・情報・選択・設定を本超音波診断装置１に取り込む。例えば、入力インタフェース回路１５は、トラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、および表示画面とタッチパッドとが一体化されたタッチパネルでディスプレイ等によって実現される。
【０１１５】
入力インタフェース回路１５は、表示画面上に表示されるカーソルの座標を検出し、検出した座標を制御回路３７に出力する。入力インタフェース回路１５は、電磁誘導式、電磁歪式、感圧式等の座標読み取り原理でタッチ指示された座標を検出し、検出した座標を制御回路３７に出力する。
【０１１６】
入力インタフェース回路１５は、バス２１を介して処理回路３１および制御回路３７に接続される。入力インタフェース回路１５は、操作者から受け取った入力操作を電気信号に変換する。入力インタフェース回路１５は、変換した電気信号を処理回路３１および制御回路３７へ出力する。なお、本明細書において入力インタフェース回路１５は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限らない。例えば、本超音波診断装置１とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、受け取った電気信号を処理回路３１および制御回路３７へ出力するような電気信号の処理回路も入力インタフェース回路１５の例に含まれる。
【０１１７】
入力インタフェース回路１５は、レンダリングパラメータ、ＭＰＲパラメータ等の種々のパラメータを入力する。入力インタフェース回路１５は、レンダリング画像およびＭＰＲ画像を任意に回転させる操作を入力してもよい。また、入力インタフェース回路１５は、フリーズ操作、シネ表示の開始・終了の指示を入力する。
【０１１８】
制御回路３７は、本超音波診断装置１における各回路等を制御するプロセッサ（Ｐｒｏｃｅｓｓｉｎｇ  ｃｉｒｃｕｉｔｒｙ）である。制御回路３７は、いずれも図示していないＣＰＵとメモリを備える。制御回路３７は、本超音波診断装置１における各回路等を制御するための各種プログラム（装置制御プログラム、表示制御プログラム等）を記憶回路３５から読み出す。制御回路３７は、読み出した各種プログラムを自身のメモリに展開し、実行することで、各種プログラムに対応する機能を実現する。換言すると、記憶回路３５からプログラムを読み出した状態の制御回路３７は、システム制御機能、表示制御機能等の各機能を有する。
【０１１９】
具体的には、制御回路３７は、装置制御プログラムを記憶回路３５から読み出す。制御回路３７は、読み出した装置制御プログラムを自身のメモリに展開し、実行することで、各回路に対する装置制御機能を実現する。このとき、制御回路３７は、システム制御部として機能する。
【０１２０】
装置制御機能は、操作者により入力インタフェース回路１５を介して入力されたフレームレート、被走査深度、送信開始・終了等に基づいて、記憶回路３５に記憶された送信遅延パターン、受信遅延パターンと装置制御プログラムとを読み出し、これらに従って装置本体１３の送受信回路２３を制御する。操作者による入力インタフェース回路１５における終了ボタンの押下に応じて、装置制御機能は、超音波の送受信を終了するために送受信回路を制御する。このとき、装置本体１３は一時停止状態となる。
【０１２１】
制御回路３７は、表示制御プログラムを記憶回路３５から読み出し、表示回路１７を制御する。このとき、制御回路３７は、表示制御部として機能する。表示制御機能は、例えば、入力インタフェース回路１５を介した操作者の指示により、Ｂモード画像、ドプラ画像、カラードプラ画像、ドプラ波形画像、第１重畳画像、第２重畳画像等の各種画像を切り換えて表示するために、表示回路１７を制御する。表示制御機能は、動画表示された上記各種画像に対するフリーズ操作の入力に応答して、動画表示されている画像をフリーズさせるために、表示回路１７を制御する。
【０１２２】
表示制御機能は、第１重畳画像および第２重畳画像各々において、心筋に対応する領域の各位置における３次元運動ベクトルの大きさに対応するカラーマップを重畳させるように、表示回路１７を制御する。例えば、表示制御機能は、第１重畳画像および第２重畳画像各々において、心筋領域に隣接する３次元血流ベクトルを示すマーカの色相を、３次元血流ベクトルに最近接する心筋領域のカラーマップに対応する色相で表示するために、表示回路１７を制御する。
【０１２３】
これにより、レンダリング画像における運動情報（３次元運動ベクトル）と、ベクトル場画像における心腔の壁部分に隣接する３次元のベクトル情報（３次元血流ベクトル）とが、同一の表示態様として表示される。表示態様とは、例えば、色相と明度と彩度とのうち少なくとも一つである。
【０１２４】
上記説明においては、単一の制御回路３７において装置制御機能及び表示制御機能等の各種機能が実行されるものとして説明したが、複数の独立したプロセッサを組み合わせて制御回路３７を構成し、各プロセッサがプログラムを実行することにより各種機能を実現するものとしても構わない。また、装置制御機能と表示制御機能とは、それぞれ異なる回路で実現されてもよいし、処理回路３１等に組み込まれてもよい。
【０１２５】
また、制御回路３７は、画像生成回路２９における画像生成機能２９１と運動情報生成機能２９３とを有していてもよい。また、制御回路３７は、処理回路３１における運動ベクトル分解機能３１１と血流ベクトル計算機能３１３とベクトル場画像生成機能３１５とを有していてもよい。
【０１２６】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ、ＧＰＵ（ｇｒａｐｈｉｃｓ  ｐｒｏｃｅｓｓｉｎｇ  ｕｎｉｔ）、或いは、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｌｅｘ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。
【０１２７】
表示回路１７は、生成された各種画像（Ｂモード画像、ドプラ画像、カラードプラ画像、ドプラ波形画像、第１重畳画像、第２重畳画像）および生体信号波形等を、ディスプレイに表示する。表示回路１７は、表示された画像に対して、ブライトネス、コントラスト、ダイナミックレンジ、γ補正などの調整および、カラーマップの割り当てを実行して表示する。
【０１２８】
表示回路１７は、カラードプラ画像をドプラ波画像とともに示すトリプレックス表示を実行してもよい。表示回路１７は、所定の範囲の時間幅に含まれる同一種類の複数の画像を連続表示（シネ表示）させることで、動画像を表示することも可能である。表示回路１７は、入力インタフェース回路１５を開始入力された回転操作に従って、レンダリング画像およびＭＰＲ画像を任意に回転させて、ディスプレイに表示する。
【０１２９】
表示回路１７は、レンダリングベクトル場画像をレンダリング画像に心位相の時相および位置を略一致させて重畳した第１重畳画像を、ディスプレイに表示する。なお、表示回路１７は、第１重畳画像における心筋領域を、心臓の運動量に応じたカラーマップでディスプレイに表示してもよい。また、表示回路１７は、心筋領域に隣接する３次元血流ベクトルの色相を、３次元血流ベクトルに最近接する心筋領域のカラーマップに対応する色相で表示してもよい。また、表示回路１７は、レンダリングベクトル場画像とレンダリング画像とを、心位相の時相を略一致させて並列表示してもよい。
【０１３０】
図７は、レンダリングベクトル場画像ＲＶＦ１をレンダリング画像ＲＩに重畳した第１重畳画像の一例を示す図である。図７に示すように、レンダリング画像ＲＩにおける心筋領域ＭＹには、心筋の運動量に応じたカラーマップが割り当てられる。図７において、心筋領域に重畳されるカラーマップは、心筋の運動量に応じてドットの密度が異なるように示されている。
【０１３１】
なお、図７に記載のレンダリングベクトル場画像ＲＶＦ１における３次元血流ベクトルは説明の便宜上２次元状に記載されているが、第１重畳画像における３次元血流ベクトルは、奥行きの位置に応じて不透明度が異なるように、３次元的（立体的に）に表示される。また、図７において、心壁に近い複数の３次元血流ベクトルＮＶ１、ＮＶ２は、最も近い心壁のカラーマップに対応する色相で表示されてもよい。
【０１３２】
図８は、レンダリングベクトル場画像の一例を示す図である。図８に示すように、３次元血流ベクトルは、第１重畳画像において、立体的に表示される。なお、図８における３次元血流ベクトルは、奥行きの位置に応じて不透明度が異なるように、ディスプレイに表示されてもよい。
【０１３３】
表示回路１７は、ＭＰＲベクトル場画像をＭＰＲ画像に心位相の時相および断面位置を略一致させて重畳した第２重畳画像を、ディスプレイに表示する。なお、表示回路１７は、第２重畳画像における心筋領域に、心臓の運動量に応じたカラーマップを重畳させて、ディスプレイに表示してもよい。また、表示回路１７は、心筋領域に隣接する血流ベクトルの色相を、血流ベクトルに最近接する心筋領域のカラーマップに対応する色相で表示してもよい。
【０１３４】
図９は、３つのＭＰＲ画像に３つのＭＰＲベクトル場画像をそれぞれ重畳した３つの第２重畳画像を示す図である。図９における２ｃｈは、２ｃｈビューに対応するＭＰＲ画像に、対応するＭＰＲベクトル場画像をそれぞれの心位相の時相と断面位置とを略一致させて重畳した第２重畳画像（以下、２ｃｈ重畳画像と呼ぶ）を示している。
【０１３５】
図９における４ｃｈは、４ｃｈビューに対応するＭＰＲ画像に、対応するＭＰＲベクトル場画像をそれぞれの心位相の時相と断面位置とを略一致させて重畳した第２重畳画像（以下、４ｃｈ重畳画像と呼ぶ）を示している。図９におけるＯｒｔｈは、２ｃｈビューおよび４ｃｈビューに直交する垂直断面に対応するＭＰＲ画像に、対応するＭＰＲベクトル場画像をそれぞれの心位相の時相と断面位置とを略一致させて重畳した第２重畳画像（以下、垂直重畳画像と呼ぶ）を示している。
【０１３６】
図９に示すように、２ｃｈ重畳画像、４ｃｈ重畳画像、垂直重畳画像各々において、心筋領域ＭＹには、心筋の運動量に応じたカラーマップが割り当てられる。図９において、カラーマップは、心筋の運動量に応じてドットの密度が異なるように示されている。好適には、３つのＭＰＲ画像（直交３断面画像）、および直交３断面画像にそれぞれ対応する３つのＭＰＲベクトル場画像は、それぞれ同一の心位相である。また、図９において、心壁に近い複数の３次元血流ベクトルＮＶ３、ＮＶ４、ＮＶ５は、最も近い心壁のカラーマップに対応する色相で表示されてもよい。
【０１３７】
図７乃至図９において、矢印は、３次元血流ベクトルに対応するマーカである。図７乃至図９におけるマーカの大きさ（矢印の長さ）は、血液の移動速度（３次元血流ベクトルの大きさ）に対応する。図７および図９におけるマーカの向き（矢印の向き）は、血液の移動方向（３次元血流ベクトルの向き）に対応する。
【０１３８】
（ベクトル場画像生成表示処理）  
ベクトル場画像生成表示処理とは、運動情報生成機能２９３と、運動ベクトル分解機能３１１と、血流ベクトル計算機能３１３と、ベクトル場画像生成機能３１５とにより実行される一連の処理である。
【０１３９】
図１０は、ベクトル場画像生成表示処理の処理手順の一例を示すフローチャートである。以下、図１０に記載のフローチャートを参照して、ベクトル場画像生成表示処理を説明する。
【０１４０】
心臓を被走査領域とした超音波の送受信により、時系列に沿った一連のボリュームデータが生成される（ステップＳａ１）。運動情報生成機能２９３により、一連のボリュームデータに基づいて、心筋の運動に関する３次元運動ベクトル（心臓の運動情報）が生成される（ステップＳａ２）。運動ベクトル分解機能３１１により、走査線方向（超音波ビーム方向）と、走査線方向に直交しかつ互いに直交する第１方向と第２方向とに、３次元運動ベクトルが投影される。この投影により、３次元運動ベクトルが、第１方向と第２方向と走査線方向とに分解される（ステップＳａ３）。
【０１４１】
心腔を被走査領域とした超音波の送受信により、走査線方向に沿ったドプラ情報が生成される（ステップＳａ４）。ステップＳａ４の処理は、ステップＳａ１の処理の前に実行されてもよい。走査線方向と第１方向とにより規定され第２方向に沿った複数の第１断面が、血流ベクトル計算機能３１３により特定される（ステップＳａ５）。走査線方向と第２方向とにより規定され第１方向に沿った複数の第２断面が、血流ベクトル計算機能３１３により特定される（ステップＳａ６）。
【０１４２】
第１断面各々において、第１断面内の各位置における血液の速さを示す第１の２次元血流ベクトルが、分解された３次元運動ベクトルの第１成分とドプラ情報とを用いて計算される（ステップＳａ７）。第１の２次元血流ベクトルは、走査線方向と第１方向とのうち少なくとも一つの方向をベクトルの成分として有する２次元ベクトルである。ステップＳａ７の処理により、第１断面各々には、第１断面内の複数の部分領域にそれぞれ対応する複数の第１の２次元血流ベクトルが生成される。すなわち、第１断面各々には、第１の２次元血流ベクトルの分布が生成される。
【０１４３】
第２断面各々において、第２断面内の各位置における血液の速さを示す第２の２次元血流ベクトルが、分解された３次元運動ベクトルの第２成分とドプラ情報とを用いて計算される（ステップＳａ８）。第２の２次元血流ベクトルは、走査線方向と第２方向とのうち少なくとも一つの方向をベクトルの成分として有する２次元ベクトルである。ステップＳａ８の処理により、第２断面各々には、第２断面内の複数の部分領域にそれぞれ対応する複数の第２の２次元血流ベクトルが生成される。すなわち、第２断面各々には、第２の２次元血流ベクトルの分布が生成される。なお、ステップＳａ７の処理とステップＳａ８の処理とは、処理手順が逆であってもよい。
【０１４４】
同一局所領域に属する第１の２次元血流ベクトルと第２の２次元血流ベクトルを合成することにより、心腔領域内の各位置における血流の速さを示す３次元血流ベクトルが生成される（ステップＳａ９）。３次元血流ベクトルの分布を示すベクトル場画像が、ディスプレイに表示される（ステップＳａ１０）。ベクトル場画像とは、心腔領域内の各位置において、血液の移動方向と血液の移動速度とに対応するマーカ（矢印）を用いて、３次元血流ベクトルの分布を示す画像である。このとき、マーカの大きさは、３次元血流ベクトルの大きさ、すなわち血液の移動速度に対応する。また、マーカの向きは、３次元血流ベクトルの向き、すなわち血液の移動方向に対応する。
【０１４５】
例えば、レンダリング画像内の心腔領域の各微小領域の位置における３次元血流ベクトルの向きおよび大きさに応じた表示処理の後に、レンダリング処理後のベクトル場画像（レンダリングベクトル場画像）が、レンダリング画像に重畳表示される。このとき、ベクトル場画像における複数の３次元血流ベクトルは、レンダリング方向に沿った所定の不透明度で表示される。また、ＭＰＲ処理後のベクトル場画像（ＭＰＲベクトル場画像）が、ＭＰＲ画像に重畳表示される。
【０１４６】
上記重畳表示において、３次元運動ベクトルの大きさに対応するカラーマップが心筋領域に重畳されてもよい。また、上記重畳表示において、心筋領域に隣接する３次元血流ベクトルを示すマーカは、最近接する心筋領域の運動情報（３次元運動ベクトル）に応じた色相が割り当てられて表示されてもよい。また、図１０に記載のフローチャートを繰り返すことで、第１重畳画像または第２重畳画像の生成に応じて、第１重畳画像または第２重畳画像は、動画としてリアルタイムにディスプレイに表示されてもよい。
【０１４７】
以上に述べた構成によれば、以下の効果を得ることができる。  
本実施形態における超音波診断装置１によれば、心臓本来の３次元的なねじれ運動を考慮して、心腔内の３次元領域の各位置における血液（流体）の速さに関する３次元のベクトル情報（３次元血流ベクトル）を生成することができる。次いで、本超音波診断装置１は、心腔領域内の各位置における３次元血流ベクトルを、３次元血流ベクトルの大きさおよび方向を表示可能なマーカとして、心筋の運動に伴う運動情報とともに、表示することができる。さらに、心筋領域の運動情報は、３次元運動ベクトルの大きさに応じたカラーマップとして表示することができる。
【０１４８】
加えて、本実施形態における超音波診断装置１によれば、心筋領域に隣接する３次元血流ベクトルを示すマーカの色相を、３次元血流ベクトルに最近接する心筋領域のカラーマップに対応する色相で表示することができる。
【０１４９】
本実施形態における超音波診断装置１によれば、３次元ベクトルデータ（３次元血流ベクトル）に対するレンダリング処理の結果として、３次元血流ベクトルを３次元的（立体的）に表示することができる。このとき、視線方向（レンダリング方向）に沿った奥行きの位置に応じて異なる不透明度が、３次元血流ベクトルに割り当てられる。
【０１５０】
以上のことから、本超音波診断装置１によれば、より実際の血流に近い血流ベクトルを生成して、３次元的に表示することができる。これにより、操作者は、心腔内における血流の状態および心筋の運動状態を、容易にかつ正確に把握および観察することができる。
【０１５１】
また、本実施形態の変形例として、本超音波診断装置１の技術的思想を医用画像処理装置１９で実現する場合には、例えば図１の構成図における破線内の構成要素を有するものとなる。このとき、本実施形態に係る画像生成機能２９１、運動情報生成機能２９３、運動ベクトル分解機能３１１、血流ベクトル計算機能３１３、ベクトル場画像生成機能３１５等の各種機能は、医用画像処理装置１９においても本実施形態と同様にして実行される。例えば、本医用画像処理装置１９は、超音波プローブ１１からの出力を基に生成された心腔内の流体のドプラ情報と、超音波プローブ１１からの出力を基に生成された前記心腔に関する心臓の運動情報とに基づいて、心腔内の３次元領域の各位置における流体の速さに関する３次元のベクトル情報を生成する処理回路（処理部）３１を有する。
【０１５２】
医用画像処理装置１９では、図１０に示すフローチャートにおけるステップＳａ１の処理は、「記憶回路３５から一連のボリュームデータを読み出し」となる。なお、ステップＳａ１の処理は、「通信インタフェース回路３３を介して、超音波診断装置または他の医用画像診断装置から一連のボリュームデータを取得」という処理であってもよい。
【０１５３】
また、医用画像処理装置１９では、図１０に示すフローチャートにおけるステップＳａ４の処理は、「記憶回路３５からドプラ情報を読み出し」となる。なお、ステップＳａ４の処理は、「通信インタフェース回路３３を介して、超音波診断装置からドプラ情報をを取得」という処理であってもよい。
【０１５４】
加えて、本実施形態に係る画像生成機能２９１、運動情報生成機能２９３、運動ベクトル分解機能３１１、血流ベクトル計算機能３１３、ベクトル場画像生成機能３１５等の各種機能は、当該各種機能およびベクトル場画像生成表示処理を実行するプログラム（医用画像処理プログラム）をワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、医用画像処理プログラムは、コンピュータに、超音波プローブ１１からの出力を基に生成された心腔内の流体のドプラ情報と、超音波プローブ１１からの出力を基に生成された心腔に関する心臓の運動情報とに基づいて、心腔内の３次元領域の各位置における前記流体の速さに関する３次元のベクトル情報を生成することを実現させる。また、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク、光ディスク、半導体メモリなどの各種可搬型記憶媒体に格納して頒布することも可能である。
【０１５５】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１５６】
１…超音波診断装置、１１…超音波プローブ、１３…装置本体、１５…入力インタフェース回路、１７…表示回路、２１…バス、２３…送受信回路、２５…Ｂモード回路、２７…ドプラモード回路、２９…画像生成回路、３１…処理回路、３３…通信インタフェース回路、３５…記憶回路、３７…制御回路、２９１…画像生成機能、２９３…運動情報生成機能、３１１…運動ベクトル分解機能、３１３…血流ベクトル計算機能、３１５…ベクトル場画像生成機能。
【図１】
0006571461.tif 000002
【図２】
0006571461.tif 000003
【図３】
0006571461.tif 000004
【図４】
0006571461.tif 000005
【図５】
0006571461.tif 000006
【図６】
0006571461.tif 000007
【図７】
0006571461.tif 000008
【図８】
0006571461.tif 000009
【図９】
0006571461.tif 000010
【図１０】
0006571461.tif 000011