JP2012152601A - 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム - Google Patents

Abstract

【課題】特定の組織の運動情報をリアルタイムに提供することが可能な超音波診断装置を提供することを目的とする。
【解決手段】輪郭検出部１１は所定の心時相に取得された断層像データに基づいて、心臓の内膜（外膜）の輪郭を検出する。輪郭追跡部１２は、各心時相に取得された断層像データにおける内膜（外膜）を構成する各点の位置を時相ごとにパターンマッチングによって求める。演算部２０は、各心時相の輪郭を構成する各点の位置に基づいて、各心時相における内膜（外膜）の運動状態を表す運動情報を求める。表示制御部７は、各心時相に取得された断層像データに基づく断層像を時相ごとに表示部８１に表示させ、さらに、各心時相における運動情報を時相ごとにその都度、表示部８１に表示させる。
【選択図】図１

Description

この発明は、超音波によって被検体の超音波画像を取得し、その超音波画像を利用して被検体の運動状態を解析する超音波診断装置に関する。また、超音波画像を利用して被検体の運動状態を解析する超音波画像処理装置、及び超音波画像処理プログラムに関する。

心臓の心筋などの生体組織について、その機能を客観的かつ定量的に評価することは、その生体組織の診断にとって非常に重要である。例えば、超音波診断装置によって心臓の画像データを取得し、その画像データに基づく定量的な評価方法が提案されている。

例えば、画像のスペックルパターンを追跡することで、組織の変位や歪みなどの運動情報を求める手法が提案されている（例えば特許文献１）。この手法は、画像のスペックルパターンを用いてパターンマッチングを行っており、スペックルトラッキング（Ｓｐｅｃｋｌｅ Ｔｒａｃｋｉｎｇ：ＳＴ）と称されている。以下、ＳＴ法と称する場合がある。

具体例として心臓の心筋の機能を評価する場合、心臓に対して超音波を送信することで、取得された時間が連続する複数の断層像データを取得する。そして、ＳＴ法によって心臓の内膜のパターンマッチングを行うことで、内膜の変位や歪みなどの壁運動情報を求めることができる。

また、画像の輝度情報に基づいて組織領域と血液領域との境界を検出し、リアルタイムにその境界の位置を検出する手法が提案されている（例えば非特許文献１）。この手法は、Ａｃｕｓｔｉｃ Ｑｕａｎｔｉｆｉｃａｔｉｏｎ（ＡＱ）法と称されている。以下、ＡＱ法と称する場合がある。

特開２００３−１７５０４１号公報

Ｕｌｔｒａｓｏｎｉｃ Ｉｍａｇｉｎｇ ５，３００−３０７（１９８３）

しかしながら、非特許文献１に記載の手法は、組織領域の輪郭を検出しているだけであり、例えば内膜上の各点の動きを追跡していないため、内膜上の各点の動きを表すベクトル情報は得られない。そのため、心筋の壁運動を解析することができない。

また、ＳＴ法によれば、例えば内膜上の各点の動きを表すベクトル情報を取得することができるため、心筋の壁運動を解析することが可能である。しかしながら、所定の時相において、追跡対象となる関心領域の境界位置を操作者が指定する必要があるため、リアルタイムに関心領域の運動を解析することが困難であった。例えば、内膜の輪郭を追跡する場合、所定の時相に取得された断層像（静止画）を表示装置に表示し、操作者はその断層像を観察して内膜の境界位置を指定する必要があった。そのため、内膜の壁運動をリアルタイムに解析することが困難であった。

また、特許文献１に記載の手法においても、ＳＴ法を用いて関心領域の運動情報をリアルタイムに表示するための具体的な手法については開示されていない。

以上のように、従来技術によると、ベクトル情報を用いて関心領域の運動をリアルタイムに解析して表示することが困難であった。

この発明は上記の問題点を解決するものであり、特定の組織の運動情報をリアルタイムに提供することが可能な超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラムを提供することを目的とする。

請求項１に記載の発明は、被検体を超音波で走査することで、時相ごとに超音波画像データを取得する画像取得手段と、所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡手段と、前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算手段と、前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示手段に表示させ、さらに、前記演算手段によって求められた各時相における運動情報を時相ごとにその都度、前記表示手段に表示させる表示制御手段と、を有することを特徴とする超音波診断装置である。
また、請求項８に記載の発明は、被検体を超音波で走査することで時相ごとに取得された超音波画像データを受け付け、所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡手段と、前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算手段と、前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示手段に表示させ、さらに、前記演算手段によって求められた各時相における運動情報を時相ごとにその都度、前記表示手段に表示させる表示制御手段と、を有することを特徴とする超音波画像処理装置である。
また、請求項９に記載の発明は、コンピュータに、被検体を超音波で走査することで時相ごとに取得された超音波画像データを受け付け、所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡機能と、前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算機能と、前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示装置に表示させ、さらに、前記演算機能によって求められた各時相における運動情報を時相ごとにその都度、前記表示装置に表示させる表示制御機能と、を実行させることを特徴とする超音波画像処理プログラムである。

この発明の第１実施形態に係る超音波診断装置を示すブロック図である。 この発明の第１実施形態に係る超音波診断装置で取得された画像と、心機能評価の結果を示す図である。 この発明の第１実施形態に係る超音波診断装置による一連の動作を示すフローチャートである。 この発明の第１実施形態に係る超音波診断装置の変形例によって取得された画像を示す図である。 この発明の第２実施形態に係る超音波診断装置を示すブロック図である。 この発明の第２実施形態に係る超音波診断装置で取得された画像を示す図である。

［第１の実施の形態］
（構成）
この発明の第１実施形態に係る超音波診断装置の構成について、図１を参照して説明する。図１は、この発明の第１実施形態に係る超音波診断装置を示すブロック図である。

超音波プローブ２には、複数の超音波振動子が所定方向（走査方向）に１列に配列された１次元アレイプローブ、又は、複数の超音波振動子が２次元的に配置された２次元アレイプローブが用いられる。また、超音波振動子が所定方向（走査方向）に配列され、超音波振動子を走査方向に直交する方向（揺動方向）に機械的に揺動可能な１次元アレイプローブを用いても良い。

送受信部３は送信部と受信部とを備え、超音波プローブ２に電気信号を供給して超音波を発生させ、超音波プローブ２が受信したエコー信号を受信する。

送受信部３の送信部は、図示しないクロック発生回路、送信遅延回路、及びパルサ回路を備えている。クロック発生回路は、超音波信号の送信タイミングや送信周波数を決めるクロック信号を発生する回路である。送信遅延回路は、超音波の送信時に遅延を掛けて送信フォーカスを実施する回路である。パルサ回路は、各超音波振動子に対応した個別経路（チャンネル）の数分のパルサを内蔵し、遅延が掛けられた送信タイミングで駆動パルスを発生し、超音波プローブ２の各超音波振動子に供給するようになっている。

送受信部３の受信部は、プリアンプ回路、Ａ／Ｄ変換回路、及び受信遅延・加算回路を備えている。プリアンプ回路は、超音波プローブ２の各超音波振動子から出力されるエコー信号を受信チャンネルごとに増幅する。Ａ／Ｄ変換回路は、増幅されたエコー信号をＡ／Ｄ変換する。受信遅延・加算回路は、Ａ／Ｄ変換後のエコー信号に対して受信指向性を決定するのに必要な遅延時間を与え、加算する。その加算により、受信指向性に応じた方向からの反射成分が強調される。なお、この送受信部３によって加算処理された信号を「ＲＦデータ（生データ）」と称する場合がある。送受信部３は、ＲＦデータを信号処理部４に出力する。

なお、超音波プローブ２と送受信部３がこの発明の「画像取得手段」の１例に相当する。

信号処理部４は、Ｂモード処理部やＣＦＭ処理部などを備えて構成されている。Ｂモード処理部は、エコーの振幅情報の映像化を行う。具体的には、Ｂモード処理部は送受信部３から出力された受信信号に対してバンドパスフィルタ処理を行い、その後、出力信号の包絡線を検波し、検波されたデータに対して対数変換による圧縮処理を施す。また、ＣＦＭ処理部は、動いている血流情報の映像化を行う。血流情報には、速度、分散、パワー等の情報があり、血流情報は２値化情報として得られる。

画像生成部５は、走査線信号列で表される信号処理後のデータを、空間座標に基づいた座標系のデータに変換する（デジタルスキャンコンバージョン）。例えば、画像生成部５は、Ｂモード処理部から出力された信号処理後のデータに対してスキャンコンバージョン処理を施すことで、被検体の組織形状を表すＢモード画像データ（以下、「断層像データ」と称する場合がある）を生成する。そして、画像生成部５は、断層像データなどの超音波画像データを記憶部６に出力する。

また、ボリュームスキャンが行なわれている場合、画像生成部５は、信号処理部４からボリュームデータを受け、そのボリュームデータにボリュームレンダリングを行うことで３次元画像データを生成するようにしても良い。さらに、画像生成部５は、ボリュームデータにＭＰＲ処理（Ｍｕｌｔｉ Ｐｌａｎｎａｒ Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ）を施すことにより、ＭＰＲ画像データ（任意断面の画像データ）を生成するようにしても良い。そして、画像生成部５は、３次元画像データやＭＰＲ画像データなどの超音波画像データを記憶部６に出力する。

画像生成部５によって生成された断層像データや３次元画像データなどの超音波画像データは、記憶部６に記憶される。また、被検体のＥＣＧ信号が取得されている場合、制御部９は、ＥＣＧ信号（心電波形）を超音波診断装置１の外部から受け付け、超音波画像データにその超音波画像データが生成されたタイミングで受け付けた心時相を対応付けて記憶部６に記憶させる。

第１実施形態に係る超音波診断装置１は、被検体の心臓を超音波で走査することで、心臓を表す断層像データを心時相ごとに取得する。すなわち、超音波診断装置１は、心臓の動画像データを取得する。例えば、１心周期以上に亘って被検体の心臓を超音波で走査することで、１心周期以上に亘って複数の断層像データ（心臓の動画像データ）を取得する。また、ＥＣＧ信号が取得されている場合、制御部９は、各断層像データにその断層像データが生成されたタイミングで受け付けた心時相を対応付けて記憶部６に記憶させる。これにより、複数の断層像データのそれぞれに、断層像データが生成された心時相が対応付けられて記憶部６に記憶される。

表示制御部７は、記憶部６から断層像データを読み込み、その断層像データに基づく断層像を表示部８１に表示させる。例えば、操作者が操作部８２を用いて任意の時相を指定すると、表示制御部７は指定された時相が対応付けられた断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。

画像処理部１０は、輪郭検出部１１、輪郭追跡部１２、マーカ生成部１３、及び補正部１４を備えている。画像処理部１０は、心臓を表す断層像に基づいて心臓の内膜の輪郭を検出し、取得された心時相が異なる２つの断層像間のパターンマッチングによって、各心時相における内膜の輪郭の位置を求める。

輪郭検出部１１は記憶部６から断層像データを読み込み、その断層像データを構成する各画素の輝度値に基づいて、組織領域と血液領域の境界を検出する。第１実施形態では１例として、輪郭検出部１１は、断層像データに基づき、上述したＡＱ法を用いて心臓の内膜の２次元的な輪郭を検出する。輪郭検出部１１は、予め設定された心時相にて取得された断層像データを記憶部６から読み込み、その断層像データの輝度分布に基づいて内膜の２次元的な輪郭を検出する。この予め設定された心時相は、操作者によって任意の心時相に変更することが可能である。例えば、輪郭検出部１１は、拡張末期（Ｒ波が検出された心時相）に取得された断層像データ、又は、収縮末期（Ｒ波が検出された心時相から所定時間経過後の心時相）に取得された断層像データを記憶部６から読み込み、内膜の２次元的な輪郭を検出する。断層像データには、その断層像データが生成された心時相が対応付けられて記憶部６に記憶されているため、輪郭検出部１１は、拡張末期や収縮末期などの心時相に取得された断層像データを記憶部６から読み込んで、その心時相における内膜の２次元的な輪郭を検出する。そして、輪郭検出部１１は、検出した内膜の輪郭の座標情報を輪郭追跡部１２に出力する。輪郭検出部１１にて検出された内膜の２次元的な輪郭が、輪郭追跡部１２において追跡対象となる内膜の初期輪郭に設定される。例えば、Ｒ波が検出された心時相における内膜の２次元的な輪郭が初期輪郭に設定される。

さらに、輪郭検出部１１は、検出した内膜の各位置における法線ベクトルを求め、内膜の位置からその法線ベクトル方向へ一定距離外側の位置を、心臓の外膜の２次元的な輪郭と定義しても良い。例えば、輪郭検出部１１は、内膜の位置から８ｍｍ外側の位置を外膜の輪郭と定義する。この一定距離は、操作者によって任意の値に変えることが可能である。このように外膜の２次元的な輪郭が検出されると、輪郭検出部１１は、内膜の輪郭の座標情報と、外膜の輪郭の座標情報を輪郭追跡部１２に出力する。輪郭検出部１１にて検出された外膜の２次元的な輪郭が、輪郭追跡部１２において追跡対象となる外膜の初期輪郭に設定される。例えば、Ｒ波が検出された心時相における外膜の輪郭が初期輪郭に設定される。

また、輪郭検出部１１は、ＡＱ法に代えて、ＡＳＭ法（Ａｃｔｉｖｅ Ｓｈａｐｅ Ｍｏｄｅｌ）を用いて内膜の２次元的な輪郭や外膜の２次元的な輪郭を検出しても良い。ＡＳＭ法による輪郭検出の方法については、例えば、「Ｃｏｏｌｅｓ ｅｔ ａｌ，“Ａｃｔｉｖｅ ｓｈａｐｅ ｍｏｄｅｌｓ Ｔｈｅｉｒ ｔｒａｉｎｉｎｇ ａｎｄ ａｐｐｌｉｃａｔｉｏｎ” Ｃｏｍｐｕｔ．Ｖｉｓ．Ｉｍａｇｅ Ｕｎｄｅｓｔａｎｄ．，ｖｏｌ．６１，ｎｏ．１，ｐｐ．３８−５９，Ｊａｎ．１９９５」に記載されている。ＡＳＭ法においては、所望の組織の一般的な形状を予め形状辞書に登録しておき、輪郭検出部１１は、断層像データの輝度情報と、その形状辞書に登録されている形状とに基づいて、所望の組織の輪郭を検出する。例えば、心臓の内膜や外膜の一般的な形状を形状辞書に予め登録しておき、輪郭検出部１１は、断層像データの輝度情報と、内膜（外膜）の一般的な形状とに基づいて、内膜（外膜）の輪郭を検出する。

ＡＱ法では、処理が比較的容易であるという利点があるが、輪郭検出のゲイン設定依存が大きく、検出結果がゲインの影響を受けやすい。また、ＡＱ法では、心臓の内膜だけでなく、乳頭筋や腱索などの組織の位置を内膜の輪郭位置として検出してしまうおそれがある。一方、ＡＳＭ法によると、ＡＱ法よりも処理が複雑になるが、ＡＱ法よりもゲイン設定の影響を受けにくく、また、乳頭筋や腱索を除外して内膜の輪郭を検出することができるという利点がある。第１実施形態においては、ＡＱ法とＡＳＭ法のうちいずれかの方法を用いて輪郭を検出すれば良い。

以上のように、輪郭検出部１１によって所定の心時相における内膜の２次元的な輪郭（内膜の初期輪郭）が検出されると、輪郭追跡部１２は、スペックルパターンを用いた２画像間のパターンマッチングによって、各心時相で生成された断層像データごとに、初期輪郭に設定された内膜の２次元的な輪郭を構成する各点の位置を求める。そして、輪郭追跡部１２は、内膜の２次元的な輪郭を構成する各点を時間的に追跡（トラッキング）する。この実施形態では１例として、輪郭追跡部１２は、上述したＳＴ法を用いて内膜の２次元的な輪郭を構成する各点を時間的に追跡する。

例えば、輪郭追跡部１２は、初期輪郭に設定された内膜の輪郭を構成する各点の座標情報を受け、さらに、その初期輪郭が検出された断層像データ（以下、「断層像データＡ」と称する場合がある）の次の心時相に生成された断層像データ（以下、「断層像データＢ」と称する場合がある）を記憶部６から読み込む。そして、輪郭追跡部１２は、スペックルパターンを用いた時間的に連続する２つの断層像間（断層像Ａと次の断層像Ｂ）のパターンマッチングを行うことにより、初期輪郭に設定された内膜の輪郭を構成する各点の移動ベクトルを求める。この移動ベクトルは、輪郭を構成する各点の変位と、各点が変位した移動方向を表している。つまり、輪郭追跡部１２は、２つの断層像間のパターンマッチングを行い、スペックルの移動量を算出することで、輪郭を構成する各点の移動ベクトルを求める。このように輪郭を構成する各点の移動ベクトルを求めることで、初期輪郭が設定された断層像データＡの次の断層像データＢが生成された心時相における内膜の輪郭を構成する各点の位置が求められる。

さらに輪郭追跡部１２は、断層像データＢの次の心時相に生成された断層像データ（以下、「断層像データＣ」と称する場合がある）を記憶部６から読み込み、スペックルパターンを用いた時間的に連続する２つの断層像間（断層像Ｂと断層像Ｃ）のパターンマッチングによって、内膜の輪郭を構成する各点の移動ベクトルを求める。これにより、断層像データＣが生成された心時相における内膜の輪郭を構成する各点の位置が求められる。

以上のようにして、輪郭追跡部１２は、スペックルパターンを用いたパターンマッチングによって、初期輪郭に設定された内膜の輪郭を構成する各点における移動ベクトルを各断層像データが生成された心時相ごとに求め、内膜の輪郭を構成する各点における移動ベクトルを時間的に追跡する。その結果、内膜の２次元的な輪郭を構成する各点を時間的に追跡することが可能となる。例えば、輪郭追跡部１２は、１心周期に亘って取得された全ての断層像データについて、各心時相における内膜の２次元的な輪郭を構成する各点の位置を求める。これにより、１心周期に亘って、各心時相における内膜の２次元的な輪郭を構成する各点の位置が求められる。

また、輪郭検出部１１によって所定の心時相における外膜の２次元的な輪郭（外膜の初期輪郭）が検出されると、輪郭追跡部１２は、内膜の追跡と同様に、スペックルパターンを用いた２画像間のパターンマッチングによって、各心時相で生成された断層像データごとに、外膜の２次元的な輪郭を構成する各点の位置を求める。そして、輪郭追跡部１２は、外膜の２次元的な輪郭を構成する各点を時間的に追跡する。

以上のように、各心時相の内膜の２次元的な輪郭を構成する各点が追跡されると、演算部２０は、輪郭検出部１１及び輪郭追跡部１２から各心時相の内膜の２次元的な輪郭を構成する各点の座標情報を受けて、内膜の壁運動情報を求める。例えば演算部２０は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報に基づいて、各心時相における内膜の円周方向（輪郭の接線方向）の歪み（変位の変化率）を求める。例えば、演算部２０は、Ｒ波が検出された心時相を基準時相とし、その基準時相における内膜の２次元的な輪郭と、その他の心時相における内膜の２次元的な輪郭とを比べて、各心時相における歪みを求める。そして、演算部２０は、心拍の経過に伴って、解析対象となっている心拍のＲ波の心時相を新たな基準時相に設定して、各心時相における歪みを求める。

例えば、輪郭検出部１１は、Ｒ波が検出された心時相における内膜の２次元的な輪郭（初期輪郭）を検出し、その内膜の２次元的な輪郭の座標情報を演算部２０に出力する。また、輪郭追跡部１２は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報を演算部２０に出力する。演算部２０は、輪郭検出部１１からＲ波が検出された心時相における内膜の２次元的な輪郭（初期輪郭）の座標情報を受け、さらに、輪郭追跡部１２から各心時相における内膜の２次元的な輪郭を構成する各点の座標情報を受けて、Ｒ波が検出された心時相を基準時相として、各心時相における内膜の円周方向の歪みを求める。

また、輪郭検出部１１及び輪郭追跡部１２によって各心時相における外膜の２次元的な輪郭が求められている場合は、演算部２０は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報と外膜の２次元的な輪郭を構成する各点の座標情報とに基づいて、各心時相における壁厚方向の歪み（内膜と外膜との間の厚さ方向の歪み）を求めても良い。

そして、演算部２０は、内膜の歪みなどの壁運動情報を表示制御部７に出力する。表示制御部７は演算部２０から壁運動情報を受けて、その壁運動情報をその都度、表示部８１に表示させる。演算部２０によって各心時相における壁運動情報が求められると、表示制御部７は、各心時相の壁運動情報を時相ごとにその都度、表示部８１に表示させる。また、心拍間での結果の安定性を向上させるために、演算部２０は、現在の心拍の演算結果を含む過去の数心拍の演算結果を重み付け平均し、表示制御部７はその平均値を表示部８１に表示させても良い。

また、演算部２０は、歪み以外にも、内膜（外膜）の変位や回転角を求めても良い。さらに、演算部２０は、歪み、変位、又は回転角などの時間微分を行い、内膜（外膜）の歪み率、速度、又は回転率などの運動情報を求めても良い。例えば、演算部２０は、各心時相における歪み、変位、又は回転角の差分を求めることで、歪み率、速度、又は回転率などを求める。そして、演算部２０は、歪み率などを表示制御部７に出力する。表示制御部７は、演算部２０から各心時相における歪み率などの壁運動情報を受けて、各心時相の壁運動情報を時相ごとにその都度、表示部８１に表示させる。

補正部１４は、複数の心拍の経過に伴うＳＴ法による追跡誤差が蓄積されるのを軽減するため、追跡対象となる内膜の輪郭の位置を所定のタイミングで補正する。例えば、（１）補正部１４は、内膜の輪郭の位置を定期的に補正する。例えば、補正部１４は、Ｎ心拍ごと（Ｎは整数）に内膜の輪郭の位置を補正する。Ｎ＝１の場合、補正部１４は、毎心拍、輪郭の位置を補正することになる。すなわち、補正部１４は、１心拍以上の時間の間隔で内膜の輪郭の位置を補正する。また、（２）補正部１４は、操作者が操作部８２によって指示を与えたタイミングで、内膜の輪郭の位置を補正しても良い。

（１）定期的に内膜の輪郭の位置を補正する場合、補正を行うタイミングを制御部９に予め設定しておく。制御部９は、設定されたタイミングで補正の指示を画像処理部１０に与える。補正部１４はその指示に従って、内膜の輪郭の位置を補正する。例えば、制御部９が所定数のＲ波を受信するたびに、補正の指示を画像処理部１０に与える。また、（２）操作者が操作部８２を用いて補正の指示を与えると、その指示が制御部９に出力され、制御部９は、補正の指示を画像処理部１０に与える。例えば、操作部８２に補正実行ボタンなどの入力装置を設けておき、操作者がそのボタンを押下すると、その押下に応じた信号が制御部９に出力される。制御部９はその信号を受けると、補正の指示を画像処理部１０に与える。

制御部９によって補正の指示が与えられると、まず、輪郭検出部１１がＡＱ法を実行することにより、内膜の２次元的な輪郭を検出する。輪郭検出部１１によって検出された２次元的な輪郭の位置（座標）を、位置Ｐａｑとする。そして、その時点において、輪郭追跡部１２がＳＴ法によって求めた内膜の２次元的な輪郭を構成する各点の位置（座標）を、位置Ｐｓｔとする。補正部１４は、輪郭検出部１１から位置Ｐａｑの座標情報を受け、さらに、輪郭追跡部１２から位置Ｐｓｔの座標情報を受けて、位置Ｐａｑと位置Ｐｓｔに基づいて、輪郭追跡部１２に設定する内膜の新たな初期輪郭の位置を求める。例えば、補正部１４は、以下の式（１）に従って、輪郭追跡部１２に新たに設定する内膜の初期輪郭の位置Ｐを求める。
（式１）
初期輪郭の位置Ｐ＝α×Ｐａｑ＋（１−α）×Ｐｓｔ
（但し、０≦α≦１）

上記式（１）に示すように、補正部１４は、輪郭検出部１１によって検出された位置Ｐａｑと、輪郭追跡部１２によって求められた位置Ｐｓｔについて、重み付け加算を行うことで、輪郭追跡部１２に新たに設定する初期輪郭の位置Ｐを求める。例えば、α＝１とすることで、輪郭追跡部１２に新たに設定する位置Ｐを、輪郭検出部１１が新たに検出した位置Ｐｓｔに置き換えることになる。

補正部１４は、式（１）に従って求めた内膜の初期輪郭の位置Ｐを輪郭追跡部１２に出力する。輪郭追跡部１２は、新たに求められた初期輪郭の位置Ｐに従って、ＳＴ法を実行することで、内膜の２次元的な輪郭を追跡する。また、補正部１４は、外膜の初期輪郭の位置についても、内膜と同様の処理によって補正しても良い。

以上のように、補正部１４によって、定期的に又は任意のタイミングで、追跡対象となる内膜（外膜）の初期輪郭の位置を補正することで、操作者による追跡対象の再設定の手間が省け、壁運動情報をリアルタイムかつ高精度に求めることが可能となる。

また、内膜の２次元的な輪郭（及び外膜の２次元的な輪郭）を表すマーカを表示部８１に表示しても良い。マーカ生成部１３は、輪郭検出部１１から内膜の初期輪郭の座標情報を受けると、その初期輪郭を表すマーカを生成する。そして、表示制御部７は、輪郭検出部１１にて初期輪郭が求められた断層像データに基づく断層像を表示部８１に表示させ、さらに、その断層像に初期輪郭を表すマーカを重ねて表示部８１に表示させる。

また、マーカ生成部１３は、輪郭追跡部１２から各心時相における内膜の２次元的な輪郭を構成する各点の座標情報を受けると、各心時相における内膜の２次元的な輪郭を表すマーカを生成する。そして、表示制御部７は、各心時相に生成された断層像データに基づく断層像を心時相ごとに、順次、表示部８１に表示させる。さらに、表示制御部７は、各心時相における内膜の輪郭を表すマーカを、各心時相に生成された断層像に重ねて、順次、表示部８１に表示させる。そして、表示制御部７は、順次、断層像とマーカを更新して表示部８１に表示させる。

また、輪郭検出部１１によって外膜の２次元的な輪郭が検出され、輪郭追跡部１２によって各心時相の外膜の２次元的な輪郭が求められている場合、マーカ生成部１３は、内膜の２次元的な輪郭を表すマーカと同様に、外膜の２次元的な輪郭を表すマーカを生成する。そして、表示制御部７は、内膜の輪郭を表すマーカと外膜の輪郭を表すマーカを、断層像に重ねて表示部８１に表示させる。

以上のようにして求められた内膜（外膜）の輪郭、運動情報、及び断層像の１例を図２に示す。図２は、この発明の第１実施形態に係る超音波診断装置によって取得された画像と、壁運動情報を示す図である。

表示制御部７は、記憶部６から断層像データを読み込んでその断層像データに基づく断層像１００を表示部８１に表示させる。図２に示す例では、断層像１００は心臓の短軸像を表している。そして、表示制御部７は、その断層像データが生成された心時相に対応する内膜の輪郭を表すマーカ１０１を断層像１００に重ねて表示部８１に表示させる。さらに、外膜の輪郭が追跡されて輪郭のマーカが生成されている場合は、表示制御部７は、外膜の輪郭を表すマーカ１０２をマーカ１０１とともに断層像１００に重ねて表示部８１に表示させる。そして、表示制御部７は、各心時相に取得された断層像１００、内膜の輪郭を表すマーカ１０１、及び外膜の輪郭を表すマーカ１０２を、時相ごとに順次更新して表示部８１に表示させる。なお、外膜の輪郭が追跡されていない場合は、表示制御部７は、内膜の輪郭を表すマーカ１０１のみを断層像１００に重ねて表示部８１に表示させても良い。

さらに、表示制御部７は、演算部２０によって求められた各心時相の壁運動情報を心時相ごとに表示部８１に表示させる。例えば演算部２０が、各心時相における壁厚方向（内膜と外膜との間の厚さ方向）の歪みを局所的に求め、図２に示すように、表示制御部７は、その局所的な歪みの大きさに応じた色を内膜（マーカ１０１）と外膜（マーカ１０２）の間の領域に割り当てて、断層像１００に重ねて表示部８１に表示させる。すなわち、演算部２０は、予め設定された所定箇所における壁厚方向の歪みを求め、表示制御部７は、その所定箇所における歪みの大きさに応じた色をマーカ１０１とマーカ１０２との間の領域に割り当てて、表示部８１に表示させる。そして、表示制御部７は、各心時相に取得された断層像１００、内膜の輪郭を表すマーカ１０１、外膜の輪郭を表すマーカ１０２、及び壁運動情報を、時相ごとに順次更新して表示部８１に表示させる。

また、図２に示すように、演算部２０は、心臓の短軸像に重心位置を設定し、その重心位置から放射状に直線を引くことにより、ＡＳＥ（米国心エコー図学会）が推奨する領域ごとに短軸像を６分割しても良い。この場合、演算部２０は、各領域に含まれる複数箇所の歪みを求め、それら複数箇所における歪みの平均値を求めることで、領域ごとに歪みの平均値を求める。図２に示す例では、演算部２０は、短軸像に相当する断層像１００を、ｉｎｔ（下壁）：Ｄ１、ｐｓｔ（後壁）：Ｄ２、ｌａｔ（側壁）：Ｄ３、ａｎｔ（前壁）：Ｄ４、ａｓｐ（前壁中隔）：Ｄ５、ｓｐ（後壁中隔）：Ｄ６の６つの領域に分割している。そして、演算部２０は、Ｄ１〜Ｄ６の領域ごとに、壁厚方向（内膜と外膜との間の厚さ方向）の歪みの大きさの平均値を求める。

例えば、演算部２０は、１心拍中における歪みのピーク値（図２中のＭａｘ１）や、基準時相であるＲ波が検出された心時相からそのピーク値が得られる心時相までの時間（図２中のＭａｘ１Ｔ）を領域ごとに求め、それらの値を表示制御部７に出力する。表示制御部７は、歪みのピーク値（Ｍａｘ１）とピーク値が得られる時間（Ｍａｘ１Ｔ）を表示部８１に表示させる。

また、演算部２０は、各領域における歪みの平均値の時間変化をグラフで表しても良い。表示制御部７は演算部２０から歪みの平均値の時間変化を表すグラフデータを受けて、そのグラフを表示部８１に表示させる。例えば図２に示すように、表示制御部７は、演算部２０によって作成された、各領域における歪みの平均値の時間変化を表すグラフ１０４を表示部８１に表示させる。図２中のグラフ１０４において、横軸は時間［ｍｓ］を示し、縦軸が歪み（Ｒａｄｉａｌ Ｓｔｒａｉｎ）を示している。

また、演算部２０は、全領域の歪みの平均値や、その平均値の時間変化を表すグラフデータを生成しても良い。表示制御部７は、全領域の歪みの平均値やグラフを表示部８１に表示させる。例えば図２中のｇｌｏｂａｌが、全領域の歪みの平均値を表している。

さらに、表示制御部７は、ＥＣＧ波形を表示部８１に表示させ、表示部８１に表示している断層像データが取得された心時相を表すバー１０３を、そのＥＣＧ波形上に重ねて表示部８１に表示させても良い。そして、表示制御部７は、バー１０３をＥＣＧ波形上で移動させて表示させ、バー１０３が示す心時相における断層像１００、マーカ１０１、１０２、グラフ１０４を表示部８１に表示させる。このように、心時相に応じて断層像１００やグラフ１０４を更新して表示部８１に表示させる。

ユーザインターフェース（ＵＩ）８は表示部８１と操作部８２を備えて構成されている。表示部８１は、ＣＲＴや液晶ディスプレイなどのモニタで構成されており、画面上に断層像や３次元画像などが表示される。操作部８２は、キーボード、マウス、トラックボール、又はＴＣＳ（Ｔｏｕｃｈ Ｃｏｍｍａｎｄ Ｓｃｒｅｅｎ）などで構成されており、操作者の操作によって各種の指示が与えられる。例えば、内膜（外膜）の初期輪郭の補正のタイミングが操作部８２によって行われる。

制御部９は、超音波診断装置１の各部に接続されて、各部の動作を制御する。

また、画像処理部１０は、ＣＰＵと、ＲＯＭ、ＲＡＭ、ＨＤＤなどの記憶部（図示しない）を備えて構成されている。記憶部には、画像処理部１０の各部の機能を実行するための画像処理プログラムが記憶されている。その画像処理プログラムには、輪郭検出部１１の機能を実行するための輪郭検出プログラム、輪郭追跡部１２の機能を実行するための輪郭追跡プログラム、マーカ生成部１３の機能を実行するためのマーカ生成プログラム、及び補正部１４の機能を実行するための補正プログラムが含まれている。そして、ＣＰＵが、輪郭検出プログラムを実行することにより、断層像データから心臓の内膜（外膜）の輪郭を検出する。また、ＣＰＵが、輪郭追跡プログラムを実行することにより、各心時相における内膜（外膜）の輪郭を求める。また、ＣＰＵが、マーカ生成プログラムを実行することにより、内膜（外膜）の輪郭を表すマーカを生成する。また、ＣＰＵが、補正プログラムを実行することにより、予め設定されたタイミング、又は操作者からの指示を受けたタイミングに、内膜（外膜）の輪郭の位置を補正する。

また、演算部２０は、ＣＰＵと、ＲＯＭ、ＲＡＭ、ＨＤＤなどの記憶部（図示しない）を備えて構成されている。記憶部には、演算部２０の機能を実行するための演算プログラムが記憶されている。そして、ＣＰＵがその演算プログラムを実行することにより、内膜（外膜）の歪みや変位などの壁運動情報を求める。

また、表示制御部７は、ＣＰＵと、ＲＯＭ、ＲＡＭ、ＨＤＤなどの記憶部（図示しない）を備えて構成されている。記憶部には、表示制御部７の機能を実行するための表示制御プログラムが記憶されている。そして、ＣＰＵがその表示制御プログラムを実行することにより、断層像、マーカ、及び壁運動情報を表示部８１に表示させる。

なお、この実施形態における画像処理プログラム、演算プログラム、及び表示制御プログラムが、この発明の「超音波画像処理プログラム」の１例に相当する。

以上のように、内膜（外膜）の初期輪郭を検出し、パターンマッチングによって内膜（外膜）の輪郭を時間的に追跡することで、輪郭を追跡するための初期輪郭の設定が自動化されて、内膜（外膜）の壁運動情報をリアルタイムに求めて表示することが可能となる。

（超音波画像処理装置）
また、内膜（外膜）の輪郭を追跡して壁運動情報を求める超音波画像処理装置を、超音波診断装置の外部に設けても良い。この超音波画像処理装置は、上述した記憶部６、表示制御部７、ユーザインターフェース８、画像処理部１０、及び演算部２０を備えている。そして、超音波画像処理装置は、取得された時間が連続する複数の断層像データを外部の超音波診断装置から取得し、それら複数の断層像データに基づいて、内膜（外膜）の輪郭を追跡して壁運動情報を求める。

超音波画像処理装置の外部に設置された超音波診断装置によって心臓を超音波で走査することで、心時相ごとに断層像データを取得する。そして、超音波画像処理装置は、超音波診断装置によって取得された複数の断層像データを受け付けて、それら複数の断層像データを記憶部６に記憶させる。そして、超音波画像処理装置の画像処理部１０は、例えばＡＱ法によって内膜（外膜）の２次元的な輪郭を検出し、さらに、ＳＴ法によって各心時相における内膜（外膜）の２次元的な輪郭を構成する各点の位置を求めることで、内膜（外膜）の輪郭を追跡する。そして、超音波画像処理装置の演算部２０は、画像処理部１０によって追跡された内膜（外膜）の２次元的な輪郭を構成する各点の位置に基づいて、内膜（外膜）の歪みや変位などの壁運動情報を求める。

以上のように、超音波診断装置の外部に設けられた超音波画像処理装置によって各心時相における内膜（外膜）の輪郭を検出して追跡することで、上述した超音波診断装置１と同様に、内膜（外膜）の壁運動情報をリアルタイムに求めて表示することが可能となる。

（動作）
次に、この発明の第１実施形態に係る超音波診断装置の動作について、図３を参照して説明する。図３は、この発明の第１実施形態に係る超音波診断装置による一連の動作を説明するためのフローチャートである。この実施形態では、心臓を診断部位として、取得された時間が連続する複数の断層像データ（動画像データ）を取得し、その断層像データに基づいて、心機能の評価に供される心臓の内膜と外膜の２次元的な輪郭を検出する。

（ステップＳ０１）
まず、超音波プローブ２を被検体に当てて、心臓に対して超音波を送信し、画像生成部５によって、各心時相の断層像データ（心臓の動画像データ）を生成する。制御部９は、超音波診断装置１の外部からＥＣＧ信号を受け付けて、生成された断層像データにその断層像データが生成された心時相を対応付けて記憶部６に記憶させる。

（ステップＳ０２）
輪郭検出部１１は、予め設定された心時相（例えばＲ波が検出された心時相）に取得された断層像データを記憶部６から読み込み、その断層像データの輝度分布に基づいて心臓の内膜と外膜の２次元的な輪郭を検出する。輪郭検出部１１にて検出された内膜と外膜の２次元的な輪郭が、輪郭追跡部１２において内膜と外膜の初期輪郭に設定される。

（ステップＳ０３）
輪郭追跡部１２は、初期輪郭が求められた断層像データの次の心時相にて取得された断層像データを記憶部６から読み込み、スペックルパターンを用いた２画像間のパターンマッチングによって、その心時相における内膜と外膜の２次元的な輪郭を構成する各点の位置を求める。

（ステップＳ０４）
演算部２０は、その心時相における内膜と外膜の２次元的な輪郭を構成する各点の座標情報を受けて、その心時相における歪みや変位などの壁運動情報を求める。

（ステップＳ０５）
マーク生成部１３は、その心時相における内膜と外膜の２次元的な輪郭を構成する各点の座標情報を受けて、内膜と外膜の輪郭を表すマーカを生成する。

（ステップＳ０６）
表示制御部７は、その心時相に生成された断層像データに基づく断層像を表示部８１に表示させ、内膜と外膜の２次元的な輪郭を表すマーカをその断層像に重ねて表示部８１に表示させる。さらに、表示制御部７は、演算部２０によって求められた壁運動情報を表示部８１に表示させる。

（ステップＳ０７、Ｓ０８）
そして、初期輪郭の位置を補正するタイミングになると（ステップＳ０７、Ｙｅｓ）、制御部９は、画像処理部１０に補正の指示を与える。例えば、初期輪郭の位置を補正するタイミングを制御部９に設定しておき、制御部９は、そのタイミングで補正の指示を画像処理部１０に与える。また、操作者が操作部８２を用いて、補正の指示を与えても良い。補正部１４は、上述した式（１）に従って、内膜と外膜の新たな初期輪郭の位置を求める（ステップＳ０８）。補正部１４によって求められた初期輪郭は、輪郭追跡部１２に設定され、輪郭追跡部１２は、その初期輪郭が示す輪郭を追跡する。そして、ステップＳ０１、及びステップＳ０３〜ステップＳ０６の処理を繰り返すことで、内膜と外膜の２次元的な輪郭を追跡し、各心時相における壁運動情報を求める。

一方、初期輪郭の位置を補正するタイミングにならない場合（ステップＳ０７、Ｎｏ）、ステップＳ０１、及びステップＳ０３〜ステップＳ０６の処理を繰り返すことで、内膜と外膜の２次元的な輪郭を追跡し、各心時相における壁運動情報を求める。

（変形例）
次に、上述した第１実施形態の変形例について図４を参照して説明する。図４は、この発明の第１実施形態に係る超音波診断装置の変形例によって取得された画像を示す図である。上述した第１実施形態においては、輪郭検出部１１が断層像データに基づいて内膜（外膜）の初期輪郭を検出する。この変形例においては、内膜（外膜）の初期輪郭を操作者が指定する。すなわち、輪郭追跡部１２は、所定の心時相に取得された断層像に対する内膜（外膜）の２次元的な輪郭の指定を受け付け、各心時相における内膜（外膜）の２次元的な輪郭を構成する各点の位置を心時相ごとにパターンマッチングによって求める。

変形例においては、表示制御部７は、初期輪郭を設定するための断層像と、運動情報を表示するための断層像とを並べて表示部８１に表示させる。例えば図４に示すように、表示制御部７は、予め設定された心時相に取得された断層像データを記憶部６から読み込んで、その断層像データに基づく断層像２００を表示部８１に表示させる。図４に示す例では、断層像２００は心臓の短軸像を表している。例えば、表示制御部７は、拡張末期（Ｒ波が検出された心時相）に取得された断層像データ、又は、収縮末期に取得された断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。さらに、表示制御部７は、予め設定された周期で、断層像２００を更新して表示部８１に表示させる。例えば、表示制御部７は、１心拍以上の時間の間隔をおいて、断層像２００を更新して表示部８１に表示させる。１例として、表示制御部７は、Ｒ波が検出された心時相に取得された断層像データを毎心拍ごとに記憶部６から読み込み、断層像２００を更新して表示部８１に表示させる。この断層像２００は、内膜（外膜）の初期輪郭を設定するために用いられる画像であり、この発明の「初期輪郭設定用の断層像」の１例に相当する。

操作者は、断層像２００を観察し、内膜（外膜）の初期輪郭を指定する。例えば、操作者は操作部８２を用いて、断層像２００に表されている内膜（外膜）の２次元的な輪郭をなぞることで、内膜（外膜）の初期輪郭を指定する。このように内膜（外膜）の２次元的な輪郭が指定されると、ユーザインターフェース（ＵＩ）８から内膜（外膜）の２次元的な輪郭の座標情報が制御部９を介して画像処理部１０に出力される。マーカ生成部１３は、操作者によって指定された内膜（外膜）の２次元的な輪郭の座標情報に基づいて、内膜（外膜）の初期輪郭を表すマーカを生成する。表示制御部７は、その初期輪郭を表すマーカを受けると、断層像２００にそのマーカを重ねて表示部８１に表示させる。図４に示す例では、表示制御部７は、内膜の初期輪郭を表すマーカ２０１と、外膜の初期輪郭を表すマーカ２０２を断層像２００に重ねて表示部８１に表示させる。このように、操作者によって指定された内膜（外膜）の輪郭が、内膜（外膜）の初期輪郭として、輪郭追跡部１２に設定される。

また、表示制御部７は、初期輪郭を設定するための輪郭設定用マーカを断層像２００に重ねて表示部８１に表示させても良い。例えば、表示制御部７は、予め設定された大きさの円形状のマーカや、楕円形状のマーカを断層像２００に重ねて表示部８１に表示させる。そして、操作者が操作部８２を用いて輪郭設定用マーカを画面上で所望の位置に移動させることで、内膜（外膜）の初期輪郭を指定する。例えば、輪郭設定用マーカの形状とほぼ合致する位置に輪郭設定用マーカを移動させて、内膜（外膜）の初期輪郭を指定する。

また、心臓の心尖像を取得して表示部８１に表示している場合、操作者が操作部８２を用いて心尖の位置と両弁輪の位置を指定し、輪郭追跡部１２は、指定された３点に基づいて、内膜（外膜）の２次元的な輪郭の位置を推定しても良い。そして、輪郭追跡部１２は、推定した内膜（外膜）の輪郭をパターンマッチングによって時間的に追跡する。さらに、弁輪の形状を予め形状辞書に登録しておき、パターンマッチングを用いた自己探索型の処理によって、内膜（外膜）の２次元的な輪郭を推定しても良い。この場合、操作者は心尖の位置のみを指定すれば良いため、操作者の負担を軽減することができる。

以上のように操作者によって内膜（外膜）の初期輪郭が指定されると、輪郭追跡部１２は、上述した第１実施形態と同様に、ＳＴ法を用いて、内膜（外膜）の２次元的な輪郭を追跡する。そして、マーカ生成部１３は、各心時相における内膜（外膜）の２次元的な輪郭を表すマーカを生成する。表示制御部７は、各心時相に取得された断層像データを記憶部６から読み込み、各心時相の断層像データに基づく断層像を心時相ごとに順次、表示部８１に表示させる。さらに、表示制御部７は、各心時相における内膜（外膜）の２次元的な輪郭を表すマーカをマーカ生成部１３から受けて、各心時相のマーカを心時相ごとに順次、断層像に重ねて表示部８１に表示させる。例えば、図４に示すように、表示制御部７は、各心時相の断層像２１０を心時相ごとに順次、表示部８１に表示させ、さらに、各心時相の内膜の輪郭を表すマーカ２１１を心時相ごとに順次、断層像２１０に重ねて表示部８１に表示させる。図４に示す例では、表示制御部７は、初期輪郭を設定するための断層像２００を予め設定された周期で表示部８１に表示させ、更に、断層像２００と並べて、断層像２１０を心時相ごとに順次、表示部８１に表示させる。

さらに、表示制御部７は、ＥＣＧ波形を表示部８１に表示させ、表示部８１に表示している断層像データが取得された心時相を表すバー２１３を、そのＥＣＧ波形上に重ねて表示部８１に表示させても良い。そして、表示制御部７は、バー２１３をＥＣＧ波形上で移動させて表示させ、バー２１３が示す心時相における断層像２１０、マーカ２１１、２１２を表示部８１に表示させる。

また、演算部２０は、上述した第１実施形態と同様に、各心時相における内膜（外膜）の２次元的な輪郭を構成する各点の座標情報に基づいて、各心時相における歪みや変位などの壁運動情報を求め、表示制御部７は各心時相の壁運動情報を心時相ごとに表示部８１に表示させる。

さらに、任意の時に、内膜（外膜）の初期輪郭を再度設定しても良い。上述した初期輪郭の指定方法と同様に、任意の時に、操作者は操作部８２を用いて、表示部８１に表示されている初期輪郭設定用の断層像２００上で内膜（外膜）の輪郭を指定する。このように再度、内膜（外膜）の初期輪郭が指定されると、ユーザインターフェース（ＵＩ）８から新たに指定された初期輪郭の座標情報が画像処理部１０に出力される。そして、輪郭追跡部１２は、新たに指定された内膜（外膜）の輪郭を追跡する。

また、初期輪郭を再設定する場合においても、表示制御部７は、内膜と外膜の初期輪郭を指定するための輪郭設定用マーカを断層像２００に重ねて表示部８１に表示させても良い。このとき、表示制御部７は、輪郭追跡部１２によって求められた所定の心時相における内膜（外膜）の輪郭を表すマーカを輪郭設定用マーカとして、断層像２００に重ねて表示部８１に表示させても良い。この輪郭設定用マーカは、マーカ生成部１３によって生成される。このように、ＳＴ法によって求められた内膜（外膜）の輪郭に基づいて輪郭設定用マーカを生成することで、輪郭設定用マーカの形状を実際の心筋の位置や輪郭形状により適合させることが可能となる。

また、表示制御部７は、初期輪郭設定用の断層像２００を常時、表示部８１に表示させても良いし、初期輪郭を設定する時のみ、断層像２００を表示部８１に表示させても良い。例えば、操作者が操作部８２を用いて初期輪郭の設定指示を与えると、表示制御部７は、予め設定された心時相（例えばＲ波が検出された心時相）に取得された断層像データを記憶部６から読み込み、その断層像データに基づく断層像２００を初期輪郭設定用の断層像として表示部８１に表示させる。そして、操作者によって初期輪郭が指定されると、表示制御部７は、壁運動情報を表示するための断層像２１０のみを表示部８１に表示させる。

また、表示制御部７は、各心時相における断層像２１０を心時相ごとに表示部８１に表示させ、所定のタイミングで、その断層像２１０に代えて、予め設定された心時相に取得された初期輪郭設定用の断層像２００を、１心拍以上の時間の間隔で更新して表示部８１に表示させても良い。そして、操作者は操作部８２を用いて、表示部８１に表示されている初期輪郭設定用の断層像２００上で内膜（外膜）の輪郭を指定する。操作者によって初期輪郭が指定されると、表示制御部７は、断層像２００に代えて、各心時相における断層像２１０を心時相ごとに表示部８１に表示させる。

なお、この変形例においては、操作者が内膜（外膜）の初期輪郭を指定するため、画像処理部１０に輪郭検出部１１を設けなくても良い。

また、上述した第１実施形態と同様に、内膜（外膜）の初期輪郭の指定を受け付け、内膜（外膜）の輪郭を追跡して壁運動情報を求める超音波画像処理装置を、超音波診断装置の外部に設けても良い。

［第２の実施の形態］
次に、この発明の第２実施形態に係る超音波診断装置について図５及び図６を参照して説明する。図５は、この発明の第２実施形態に係る超音波診断装置を示すブロック図である。図６は、この発明の第２実施形態に係る超音波診断装置で取得された画像を示す図である。

上述した第１実施形態においては、２次元画像としての断層像データに基づいて、２次元平面上における内膜（外膜）の輪郭を追跡し、壁運動情報を求めている。第２実施形態においては、３次元画像としてのボリュームデータに基づいて、内膜（外膜）の３次元的な輪郭を追跡し、壁運動情報を求める。

第２実施形態に係る超音波診断装置１Ａは、第１実施形態に係る超音波診断装置１の画像処理部１０の代わりに、画像処理部１０Ａを備えている。画像処理部１０Ａは、輪郭検出部１１Ａ、輪郭追跡部１２、マーカ生成部１３、及び補正部１４を備えている。

第２実施形態においては、超音波プローブ２に２次元アレイプローブを用いてボリュームスキャンを実行することで、心時相ごとにボリュームデータを取得する。そして、画像生成部５は、予め設定された心時相（例えばＲ波が検出された心時相）に取得されたボリュームデータにＭＰＲ処理を施すことで、任意断面におけるＭＰＲ画像データを生成する。例えば、画像生成部５は、ボリュームデータにＭＰＲ処理を施すことで、複数の異なる断面に沿ったＭＰＲ画像データをそれぞれの断面ごとに求める。１例として、画像生成部５は、互いに直交する断面に沿った２つのＭＰＲ画像データを生成する。例えば図６（ａ）に示すように、画像生成部５は、心臓の長軸方向に沿った心尖四腔像（４Ｃ）３００と、心尖二腔像（２Ｃ）３１０を生成する。断面は予め画像生成部５に設定されていても良いし、操作者が操作部８２を用いることで指定しても良い。画像生成部５によって生成されたＭＰＲ画像データは記憶部６に記憶される。また、表示制御部７は、記憶部６に記憶されているＭＰＲ画像データを読み込んで、そのＭＰＲ画像データに基づくＭＰＲ画像を表示部８１に表示させても良い。例えば、表示制御部７は、心尖四腔像（４Ｃ）３００と、心尖二腔像（２Ｃ）３１０を表示部８１に表示させても良い。

輪郭検出部１１Ａは、２次元輪郭検出部１５と補間部１６を備えている。２次元輪郭検出部１５は、予め設定された心時相（例えばＲ波が検出された心時相）に取得された複数のＭＰＲ画像データを記憶部６から読み込み、ＡＱ法又はＡＳＭ法を用いて、各ＭＰＲ画像データから内膜（外膜）の２次元的な輪郭を検出する。図６（ａ）に示す例では、２次元輪郭検出部１５は、心尖四腔像（４Ｃ）３００から内膜の２次元的な輪郭３０１と外膜の２次元的な輪郭３０２を検出する。また、２次元輪郭検出部１５は、心尖二腔像（２Ｃ）３１０から内膜の２次元的な輪郭３１１と外膜の２次元的な輪郭３１２を検出する。このように、第２実施形態においては１例として、心臓の長軸方向に沿ったＭＰＲ画像データに基づいて、長軸方向における内膜（外膜）の２次元的な輪郭を検出する。そして、輪郭検出部１１は、内膜（外膜）の２次元的な輪郭の座標情報を補間部１６に出力する。

補間部１６は、心尖四腔像（４Ｃ）３００に基づいて求められた内膜の２次元的な輪郭３０１と、心尖二腔像（２Ｃ）３１０に基づいて求められた内膜の２次元的な輪郭３１１とに基づいて、スプライン関数などを用いて空間的に内膜の輪郭を補間することで、内膜の３次元的な輪郭の位置を求める。また、補間部１６は、心尖四腔像（４Ｃ）３００に基づいて求められた外膜の２次元的な輪郭３０２と、心尖二腔像（２Ｃ）３１０に基づいて求められた外膜の２次元的な輪郭３１２とに基づいて、スプライン関数などを用いて空間的に外膜の輪郭を補間することで、外膜の３次元的な輪郭の位置を求める。すなわち、補間部１６は、長軸方向に直交する方向（短軸方向）に沿った各断面（各Ｃ面）における内膜（外膜）の輪郭の位置を求めることで、内膜（外膜）の３次元的な輪郭の位置を求める。

以上のように、補間部１６は、長軸方向に沿った断面で定義された輪郭３０１、３０２、３１１、３１２に基づいて、長軸方向の各深さにおける内膜（外膜）の輪郭の位置を求めることで、内膜（外膜）の３次元的な輪郭の位置を求める。補間部１６によって求められた内膜（外膜）の３次元的な輪郭の座標情報は輪郭追跡部１２に出力され、輪郭追跡部１２に内膜（外膜）の初期輪郭として設定される。

また、画像生成部５は、長軸方向に沿った断面のＭＰＲ画像データを生成し、さらに、長軸方向に直交する方向（短軸方向）に沿った断面（Ｃ面）のＭＰＲ画像データを生成しても良い。例えば図６（ａ）に示すように、画像生成部５は、心臓の長軸方向に沿った心尖四腔像（４Ｃ）３００と、心尖二腔像（２Ｃ）３１０を生成する。さらに、図６（ｂ）に示すように、画像生成部５は、長軸方向の所定深さにおけるラインＡ、ラインＭ、及びラインＢに沿った断面におけるＭＰＲ画像３２０、３３０、３４０を生成する。画像生成部５によって生成されたＭＰＲ画像データは記憶部６に記憶される。ラインＡ、ラインＭ、及びラインＢによって指定される断面は、予め画像生成部５に設定されていても良いし、操作者が操作部８２を用いて断面を指定しても良い。

２次元輪郭検出部１５は、予め設定された心時相（例えばＲ波が検出された心時相）に取得された各断面のＭＰＲ画像データを記憶部６から読み込んで、各断面における内膜（外膜）の２次元的な輪郭を検出する。例えば、２次元輪郭検出部１５は、長軸方向に沿った心尖四腔像（４Ｃ）３００に基づいて、内膜の２次元的な輪郭３０１と、外膜の２次元的な輪郭３０２を検出する。また、２次元輪郭検出部１５は、長軸方向に沿った心尖二腔像（２Ｃ）３１０に基づいて、内膜の２次元的な輪郭３１１と、外膜の２次元的な輪郭３１２を検出する。さらに、２次元輪郭検出部１５は、ラインＡに沿った断面のＭＰＲ画像３２０に基づいて、内膜の２次元的な輪郭３２１と、外膜の２次元的な輪郭３２２を検出する。また、２次元輪郭検出部１５は、ラインＭに沿った断面のＭＰＲ画像３３０に基づいて、内膜の２次元的な輪郭３３１と、外膜の２次元的な輪郭３３２を検出する。さらに、２次元輪郭検出部１５は、ラインＢに沿った断面のＭＰＲ画像３４０に基づいて、内膜の２次元的な輪郭３４１と、外膜の２次元的な輪郭３４２を検出する。２次元輪郭検出部１５は、各断面に沿ったＭＰＲ画像から検出した内膜（外膜）の２次元的な輪郭の座標情報を補間部１６に出力する。

補間部１６は、２次元輪郭検出部１５によって検出された複数の２次元的な輪郭の座標情報に基づいて、スプライン関数などを用いて空間的に内膜（外膜）の輪郭を補間することで、内膜（外膜）の３次元的な輪郭の位置を求める。上述の例では、補間部１６は、心尖四腔像（４Ｃ）３００、心尖二腔像（２Ｃ）３１０、ＭＰＲ画像３２０、３３０、及び３４０から検出された内膜（外膜）の２次元的な輪郭に基づいて、スプライン関数などを用いて空間的に内膜（外膜）の輪郭を補間することで、内膜（外膜）の３次元的な輪郭の位置を求める。補間部１６によって求められた内膜（外膜）の３次元的な輪郭の座標情報は輪郭追跡部１２に出力され、輪郭追跡部１２に内膜（外膜）の初期輪郭として設定される。

輪郭追跡部１２は、補間部１６によって求められた内膜（外膜）の初期輪郭に基づき、ＳＴ法を用いて、各心時相において取得されたボリュームデータ間のパターンマッチングを行うことで、各心時相における内膜（外膜）の３次元的な輪郭を構成する各点の位置を求める。そして、輪郭追跡部１２は、内膜（外膜）の３次元的な輪郭を追跡する。

演算部２０は、第１実施形態と同様に、各心時相における内膜（外膜）の３次元的な輪郭を構成する各点の座標情報に基づいて、各心時相における歪みや変位などの壁運動情報を求める。そして、表示制御部７は、各心時相における壁運動情報を心時相ごとに表示部８１に表示させる。

以上のように、各心時相における内膜（外膜）の３次元的な輪郭を検出して追跡することで、上述した第１実施形態と同様に、内膜（外膜）の壁運動情報をリアルタイムに求めて表示することが可能となる。

また、表示制御部７は、各心時相におけるＭＰＲ画像を心時相ごとに表示部８１に表示させても良い。例えば図６（ａ）に示すように、表示制御部７は、各心時相における心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０を心時相ごとに表示部８１に表示させる。また、図６（ｂ）に示すように、表示制御部７は、短軸方向の断面に沿ったＭＰＲ画像３２０、３３０、３４０を心時相ごとに表示部８１に表示させても良い。

また、輪郭追跡部１２は、各ＭＰＲ画像上で設定された内膜（外膜）の輪郭を追跡しても良い。そして、マーカ生成部１３は、各心時相における内膜（外膜）の輪郭を表すマーカを生成する。表示制御部７は、各心時相におけるＭＰＲ画像を心時相ごとに表示部８１に表示させ、さらに、各心時相における内膜（外膜）の輪郭を表すマーカを、各心時相に生成されたＭＰＲ画像に重ねて、順次、表示部８１に表示させても良い。

例えば、輪郭追跡部１２は、図６（ｂ）に示すＭＰＲ画像３２０上で設定された内膜（外膜）の２次元的な輪郭を追跡し、マーカ生成部１３は、各心時相における内膜（外膜）の２次元的な輪郭を表すマーカを生成する。そして、表示制御部７は、各心時相におけるＭＰＲ画像３２０を心時相ごとに表示部８１に表示させ、さらに、表示制御部７は、各心時相における内膜のマーカと外膜のマーカを、ＭＰＲ画像３２０に重ねて表示部８１に表示させる。ここでは、ＭＰＲ画像３２０を例にして説明したが、他の断面に沿ったＭＰＲ画像とマーカを重ねて表示部８１に表示させても良い。

さらに、表示制御部７は、上述した第１実施形態と同様に、演算部２０によって求められた歪みの大きさに応じた色を内膜と外膜の間の領域に割り当てて、ＭＰＲ画像に重ねて表示部８１に表示させても良い。

また、演算部２０は、第１実施形態と同様に、ＡＳＥが推奨する領域ごとに歪みなどの壁運動情報を求めても良い。そして、演算部２０は、各領域の壁運動情報のピーク値や、壁運動情報をグラフ化して表示部８１に表示させても良い。

また、第１実施形態と同様に、複数の心拍の経過に伴うＳＴ法による追跡誤差が蓄積されるのを軽減するため、追跡対象となる内膜（外膜）の３次元的な輪郭の位置を所定のタイミングで補正しても良い。例えば、補正部１４は、Ｎ心拍（Ｎは整数）ごとに内膜（外膜）の３次元的な輪郭の位置を補正する。また、補正部１４は、操作者によって補正の指示が与えられたタイミングで、内膜（外膜）の３次元的な輪郭の位置を補正しても良い。例えば、制御部９が所定数のＲ波を受信するたびに、補正の指示を画像処理部１０Ａに与える。また、操作者が操作部８２を用いて補正の指示を与えると、その指示が制御部９に出力され、制御部９は、補正の指示を画像処理部１０Ａに与える。

制御部９によって補正の指示が与えられると、２次元輪郭検出部１５がＡＱ法を実行することにより、ＭＰＲ画像データに基づいて内膜（外膜）の２次元的な輪郭を検出する。例えば、２次元輪郭検出部１５は、心尖四腔像（４Ｃ）３００から内膜（外膜）の２次元的な輪郭を検出し、さらに、心尖二腔像（２Ｃ）３１０から内膜（外膜）の２次元的な輪郭を検出する。補間部１６は、２つのＭＰＲ画像から求められた内膜（外膜）の２次元的な輪郭に基づいて、スプライン関数を用いることで、内膜（外膜）の３次元的な輪郭の位置Ｐａｑを求める。そして、その時点において、輪郭追跡部１２がＳＴ法によって求めた内膜（外膜）の３次元的な輪郭の位置を、位置Ｐｓｔとする。そして、補正部１４は、上述した式（１）に従って、輪郭追跡部１２に新たに設定する内膜（外膜）の初期輪郭の位置Ｐを求める。そして、輪郭追跡部１２は、式（１）に従って求められた内膜（外膜）の初期輪郭を追跡する。

また、第１実施形態の変形例と同様に、操作者が内膜（外膜）の初期輪郭を指定し、補間部１６は、指定された初期輪郭に基づいて、３次元空間における内膜と外膜の輪郭を補間して求めても良い。例えば、表示制御部７は、予め設定された心時相に取得された心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０を表示部８１に表示させる。例えば、表示制御部７は、拡張末期（Ｒ波が検出された心時相）に取得された心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０を表示部８１に表示させる。そして、表示制御部７は、１心拍以上の時間の間隔をおいて、心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０を更新して表示部８１に表示させる。心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０は、初期輪郭を設定するための画像であり、この発明の「初期輪郭設定用の断層像」の１例に相当する。

操作者は操作部８２を用いて、心尖四腔像（４Ｃ）３００と心尖二腔像（２Ｃ）３１０上に内膜（外膜）の２次元的な輪郭を指定する。ユーザインターフェース８から内膜（外膜）の輪郭の座標情報が画像処理部１０Ａに出力される。補間部１６は、指定された内膜（外膜）の輪郭の座標情報に基づいて、スプライン関数などを用いることで、３次元空間における内膜（外膜）の初期輪郭を求める。このように操作者によって内膜（外膜）の初期輪郭が指定されると、輪郭追跡部１２は、ＳＴ法を用いて、内膜（外膜）の３次元的な輪郭を追跡する。演算部２０は、各心時相における内膜（外膜）の３次元的な輪郭に基づいて、各心時相における歪みや変位などの壁運動情報を求め、表示制御部７は各心時相の壁運動情報を心時相ごとに表示部８１に表示させる。

また、表示制御部７は、各心時相のＭＰＲ画像３２０を心時相ごとに表示部８１に表示させ、さらに、ＭＰＲ画像３２０と並べて、初期輪郭設定用のＭＰＲ画像を表示部８１に表示させても良い。このとき、表示制御部７は、初期輪郭設定用のＭＰＲ画像を１心拍以上の時間の間隔で更新して表示部８１に表示させる。また、表示制御部７は、初期輪郭設定用のＭＰＲ画像を常時、表示部８１に表示させても良いし、初期輪郭を設定する時のみ、初期輪郭設定用のＭＰＲ画像を表示部８１に表示させても良い。ここでは、ＭＰＲ画像３２０を例にして説明したが、他の断面に沿ったＭＰＲ画像を初期輪郭設定用のＭＰＲ画像と並べて表示部８１に表示させても良い。

また、表示制御部７は、各心時相のＭＰＲ画像３２０を心時相ごとに表示部８１に表示させ、所定のタイミングで、そのＭＰＲ画像３２０に代えて、予め設定された心時相に取得された初期輪郭設定用のＭＰＲ画像を表示部８１に表示させても良い。そして、操作者は操作部８２を用いて、表示部８１に表示されている初期輪郭設定用のＭＰＲ画像上で内膜（外膜）の輪郭を指定する。操作者によって初期輪郭が指定されると、表示制御部７は、初期輪郭設定用のＭＰＲ画像に代えて、各心時相におけるＭＰＲ画像３２０を心時相ごとに表示部８１に表示させる。

なお、操作者が内膜（外膜）の初期輪郭を指定する場合は、輪郭検出部１１Ａに２次元輪郭検出部１５を設けなくても良い。

また、内膜（外膜）の３次元的な輪郭を追跡して壁運動情報を求める超音波画像処理装置を、超音波診断装置の外部に設けても良い。この超音波画像処理装置は、上述した記憶部６、表示制御部７、ユーザインターフェース８、画像処理部１０Ａ、及び演算部２０を備えている。そして、超音波画像処理装置は、取得された時間が連続する複数のボリュームデータを外部の超音波診断装置から取得し、それら複数のボリュームデータに基づいて、内膜（外膜）の３次元的な輪郭を追跡して壁運動情報を求める。

なお、画像処理部１０Ａは、ＣＰＵと、ＲＯＭ、ＲＡＭ、ＨＤＤなどの記憶部（図示しない）を備えて構成されている。記憶部には、画像処理部１０Ａの各部の機能を実行するための画像処理プログラムが記憶されている。その画像処理プログラムには、２次元輪郭検出部１５の機能を実行するための輪郭検出プログラム、補間部１６の機能を実行するための補間プログラム、輪郭追跡部１２の機能を実行するための輪郭追跡プログラム、マーカ生成部１３の機能を実行するためのマーカ生成プログラム、及び補正部１４の機能を実行するための補正プログラムが含まれている。この実施形態における画像処理プログラム、演算部２０の機能を実行するための演算プログラム、及び表示制御部７の機能を実行するための表示制御プログラムが、この発明の「超音波画像処理プログラム」の１例に相当する。

１、１Ａ 超音波診断装置
２ 超音波プローブ
３ 送受信部
４ 信号処理部
５ 画像生成部
６ 記憶部
７ 表示制御部
８ ユーザインターフェース
９ 制御部
１０、１０Ａ 画像処理部
１１、１１Ａ 輪郭検出部
１２ 輪郭追跡部
１３ マーカ生成部
１４ 補正部
１５ ２次元輪郭検出部
１６ 補間部
２０ 演算部
８１ 表示部
８２ 操作部

Claims (9)

1. 被検体を超音波で走査することで、時相ごとに超音波画像データを取得する画像取得手段と、
   所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡手段と、
   前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算手段と、
   前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示手段に表示させ、さらに、前記演算手段によって求められた各時相における運動情報を時相ごとにその都度、前記表示手段に表示させる表示制御手段と、
   を有することを特徴とする超音波診断装置。
2. 前記画像取得手段は、前記被検体を超音波で走査することで、前記超音波画像データとして、時相ごとに断層像データを取得し、
   前記輪郭追跡手段は、前記所定時相に取得された断層像データに基づく断層像に対する前記特定の組織の２次元的な輪郭の指定を受け付け、前記各時相に取得された前記断層像データにおける前記指定された２次元的な輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求め、
   前記演算手段は、前記各時相の２次元的な輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求めることを特徴とする請求項１に記載の超音波診断装置。
3. 前記表示制御手段は、前記被検体の心拍情報を受けて、前記所定時相に取得された断層像データに基づく初期輪郭設定用の断層像を、１心拍以上の時間の間隔で更新して前記表示手段に表示させ、
   前記初期輪郭設定用の断層像上で前記特定の組織の２次元的な輪郭が指定されると、前記輪郭追跡手段は、その指定を受け付け、前記各時相に取得された前記断層像データにおける前記指定された２次元的な輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求めることを特徴とする請求項２に記載の超音波診断装置。
4. 画像生成手段と補間手段を更に有し、
   前記画像取得手段は、前記被検体を超音波で走査することで、前記超音波画像データとして、時相ごとにボリュームデータを取得し、
   前記画像生成手段は、前記所定時相に取得されたボリュームデータに基づいて、複数の異なる断面に沿った断層像データをそれぞれの断面ごとに求め、
   前記補間手段は、前記所定時相における複数の断層像データに基づく断層像のそれぞれに対する前記特定の組織の２次元的な輪郭の指定を受け付け、各断面における２次元的な輪郭の位置の間を補間することで、前記組織の３次元的な輪郭を求め、
   前記輪郭追跡手段は、前記各時相に取得された前記ボリュームデータにおける前記補間によって求められた３次元的な輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求め、
   前記演算手段は、前記各時相の３次元的な輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求めることを特徴とする請求項１に記載の超音波診断装置。
5. 前記表示制御手段は、前記被検体の心拍情報を受けて、前記所定時相に取得された断層像データに基づく初期輪郭設定用の断層像を、１心拍以上の時間の間隔で更新して前記表示手段に表示させ、
   前記更新された複数の初期輪郭設定用の断層像上で前記特定の組織の２次元的な輪郭が指定されると、前記補間手段は、その指定を受け付け、前記各断面における２次元的な輪郭の位置の間を補間することで、前記組織の３次元的な輪郭を求めることを特徴とする請求項４に記載の超音波診断装置。
6. 前記表示制御手段は、前記各時相に取得された断層像データに基づく断層像を時相ごとに前記表示手段に表示させ、さらにその断層像と並べて、前記所定時相に取得された断層像データに基づく初期輪郭設定用の断層像を、前記１心拍以上の時間の間隔で更新して前記表示手段に表示させることを特徴とする請求項３又は請求項５のいずれかに記載の超音波診断装置。
7. 前記表示制御手段は、前記各時相に取得された断層像データに基づく断層像を時相ごとに前記表示手段に表示させ、所定のタイミングで、その断層像に代えて、前記所定時相に取得された断層像データに基づく初期輪郭設定用の断層像を、前記１心拍以上の時間の間隔で更新して前記表示手段に表示させることを特徴とする請求項３又は請求項５のいずれかに記載の超音波診断装置。
8. 被検体を超音波で走査することで時相ごとに取得された超音波画像データを受け付け、所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡手段と、
   前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算手段と、
   前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示手段に表示させ、さらに、前記演算手段によって求められた各時相における運動情報を時相ごとにその都度、前記表示手段に表示させる表示制御手段と、
   を有することを特徴とする超音波画像処理装置。
9. コンピュータに、
   被検体を超音波で走査することで時相ごとに取得された超音波画像データを受け付け、所定時相に取得された超音波画像データに基づく超音波画像に対する特定の組織の輪郭の指定を受け付け、各時相に取得された前記超音波画像データにおける前記指定された輪郭を構成する各点の位置を時相ごとにパターンマッチングによって求める輪郭追跡機能と、
   前記各時相の輪郭を構成する各点の位置に基づいて、前記各時相における前記特定の組織の運動状態を表す運動情報を求める演算機能と、
   前記各時相に取得された超音波画像データに基づく超音波画像を時相ごとに表示装置に表示させ、さらに、前記演算機能によって求められた各時相における運動情報を時相ごとにその都度、前記表示装置に表示させる表示制御機能と、
   を実行させることを特徴とする超音波画像処理プログラム。