(57)【要約】【課題】超音波を利用して心臓の診断情報を得る改良技術を提供する。【解決手段】画像形成処理部２０は、超音波の受信信号に基づいて生体内の画像データを形成する。ドプラ処理部３０は、超音波の受信信号に基づいて生体内のドプラ情報を生成する。速度ベクトル演算部４０は、生体内のドプラ情報に基づいて血流の速度情報（速度ベクトル）を生成する。関心領域設定部５０は、画像データ内において心臓内腔に対応した関心領域を設定する。エネルギー演算部７０は、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出する。【選択図】図１
2015208522.tif 000002
【特許請求の範囲】
【請求項１】
超音波を送受するプローブと、
プローブを送信制御して生体内から超音波の受信信号を得る送受信部と、
超音波の受信信号に基づいて生体内のドプラ情報を得るドプラ処理部と、
生体内のドプラ情報に基づいて血流の速度情報を得る速度情報生成部と、
超音波の受信信号に基づいて生体内の画像データを得る画像形成処理部と、
画像データ内において心臓内腔に対応した関心領域を設定する関心領域設定部と、
心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出するエネルギー演算部と、
算出されたエネルギーの損失量を示す表示画像を形成する表示処理部と、
を有する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記エネルギー演算部は、心臓内腔に対応した関心領域内の複数サンプル点について、各サンプル点ごとに前記損失量を算出する、
ことを特徴とする超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
前記表示処理部は、生体内の画像データに対して、関心領域内の前記各サンプル点ごとに当該サンプル点の前記損失量に応じた表示処理を施すことにより、心臓内腔における前記損失量の分布状態を示した表示画像を形成する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項２または３に記載の超音波診断装置において、
前記エネルギー演算部は、心臓内腔に対応した関心領域内の複数サンプル点における前記損失量に基づいて、当該関心領域内における前記損失量の総和を算出し、
前記表示処理部は、各時相ごとに算出される前記損失量の総和を複数時相に亘って示したグラフを形成する、
ことを特徴とする超音波診断装置。
【請求項５】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記関心領域設定部は、画像データ内において心臓内腔の外縁として設定された内腔ラインにより囲まれた領域を前記関心領域とする、
ことを特徴とする超音波診断装置。
【請求項６】
請求項５に記載の超音波診断装置において、
前記関心領域設定部は、画像データ内において、心臓内腔へ流入する血流の流路に設定された流入ラインと、心臓内腔から流出する血流の流路に設定された流出ラインと、心臓内腔の外縁として設定された前記内腔ラインと、により囲まれた領域を前記関心領域とする、
ことを特徴とする超音波診断装置。
【請求項７】
請求項１から６のいずれか１項に記載の超音波診断装置において、
前記エネルギー演算部は、心臓内腔へ流入する血流の流路に設定された流入ラインに対応した血流の速度情報に基づいて、心臓内腔へ流入する血流の運動エネルギーである流入エネルギーを算出し、心臓内腔から流出する血流の流路に設定された流出ラインに対応した血流の速度情報に基づいて、心臓内腔から流出する血流の運動エネルギーである流出エネルギーを算出し、
前記表示処理部は、各時相ごとに算出される流入エネルギーと流出エネルギーの少なくとも一方を複数時相に亘って示したグラフを形成する、
ことを特徴とする超音波診断装置。
【請求項８】
請求項１から７のいずれか１項に記載の超音波診断装置において、
前記速度情報生成部は、生体内のドプラ情報と画像データの少なくとも一方に基づいて、心臓内腔の外縁として設定された内腔ライン上における血流と心筋の少なくとも一方の速度情報を生成し、
前記エネルギー演算部は、前記内腔ライン上における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーを算出する、
ことを特徴とする超音波診断装置。
【請求項９】
請求項８に記載の超音波診断装置において、
前記速度情報生成部は、前記内腔ライン上の複数サンプル点について、各サンプル点ごとに血流または心筋の速度情報を生成し、
前記エネルギー演算部は、前記内腔ライン上の複数サンプル点における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーの総量を算出する、
ことを特徴とする超音波診断装置。
【請求項１０】
請求項９に記載の超音波診断装置において、
前記エネルギー演算部は、前記内腔ライン上の複数サンプル点における速度情報から、心臓内腔の内側に向かう速度情報に基づく運動エネルギーの総量と、心臓内腔の外側に向かう速度情報に基づく運動エネルギーの総量を算出する、
ことを特徴とする超音波診断装置。
【請求項１１】
超音波を利用して得られた生体内のドプラ情報に基づいて、血流の速度情報を得る速度情報生成部と、
超音波を利用して得られた生体に係る画像データ内において、心臓内腔に対応した関心領域を設定する関心領域設定部と、
心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出するエネルギー演算部と、
を有する、
ことを特徴とする情報処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、心臓の診断情報を得る技術に関する。
【背景技術】
【０００２】
血流などの流体に超音波を送受して得られる受信信号から流体に係る診断情報を得る技術が知られている。例えば、特許文献１には、生体内の流体に対して超音波を送受して得られる受信信号（エコーデータ）に基づいて、観測面内の複数点において流体に関する２次元の速度ベクトルを得る技術が記載されている。観測面内の複数点における２次元の速度ベクトルの分布から、流体の流れを示す流線などの診断情報を得ることが可能になり、例えば、心臓等の診断への応用が期待される。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−１９２６４３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上述した背景技術に鑑み、本願の発明者は、超音波を利用して心臓の診断情報を得る技術について研究開発を重ねてきた。
【０００５】
本発明は、その研究開発の過程において成されたものであり、その目的は、超音波を利用して心臓の診断情報を得る改良技術を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的にかなう好適な超音波診断装置は、超音波を送受するプローブと、プローブを送信制御して生体内から超音波の受信信号を得る送受信部と、超音波の受信信号に基づいて生体内のドプラ情報を得るドプラ処理部と、生体内のドプラ情報に基づいて血流の速度情報を得る速度情報生成部と、超音波の受信信号に基づいて生体内の画像データを得る画像形成処理部と、画像データ内において心臓内腔に対応した関心領域を設定する関心領域設定部と、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出するエネルギー演算部と、算出されたエネルギーの損失量を示す表示画像を形成する表示処理部と、を有することを特徴とする。
【０００７】
上記装置において、血流の速度情報とは、血流についての部分的な又は全体的な運動に関する情報であり、具体的には、血流内の各箇所における速度と方向を示す速度ベクトルや、各箇所における移動量と方向を示す移動ベクトルなどが含まれる。なお、血流の速度ベクトルや移動ベクトルの導出には、例えば、特許文献１に記載された技術（２次元の速度ベクトルの分布）を利用することが望ましいものの、他の公知の技術を利用して速度ベクトル等を得るようにしてもよい。
【０００８】
また、血流内で失われたエネルギーの損失量とは、例えば、血流内における摩擦等により熱エネルギー等に変換されて失われた運動エネルギー等の損失量であり、心臓や血管の機能を評価する際の一つの評価値となることが期待される。例えば、心臓内腔の血流内で失われるエネルギーの損失量だけ、心臓内腔における血流の運動エネルギーが低下してしまうため、血流の流量を維持するために、心臓から血流に大きな運動エネルギーを与えることとなり、心臓の負担（心負荷）が大きいことが予想できる。また、心臓内腔に流入する血流の運動エネルギーや心臓内腔から流出する血流の運動エネルギーが小さいと、血流の流量が維持できていないことが予想される。
【０００９】
上記装置によれば、心臓内腔の血流内で失われるエネルギーの損失量が算出され、心臓や血管の機能を評価する際の一つの評価値とすることができる。そのため、エネルギーの損失量に基づいて心臓や血管の機能を評価することが可能になり、さらに好適には、心臓内腔に流入する血流の運動エネルギーや心臓内腔から流出する血流の運動エネルギー等を評価値に加えることにより、心臓や血管の機能について、複数の評価値に基づいた総合的な診断が可能になる。
【００１０】
望ましい具体例において、前記エネルギー演算部は、心臓内腔に対応した関心領域内の複数サンプル点について、各サンプル点ごとに前記損失量を算出することを特徴とする。
【００１１】
望ましい具体例において、前記表示処理部は、生体内の画像データに対して、関心領域内の前記各サンプル点ごとに当該サンプル点の前記損失量に応じた表示処理を施すことにより、心臓内腔における前記損失量の分布状態を示した表示画像を形成する、ことを特徴とする。
【００１２】
望ましい具体例において、前記エネルギー演算部は、心臓内腔に対応した関心領域内の複数サンプル点における前記損失量に基づいて、当該関心領域内における前記損失量の総和を算出し、前記表示処理部は、各時相ごとに算出される前記損失量の総和を複数時相に亘って示したグラフを形成する、ことを特徴とする。
【００１３】
望ましい具体例において、前記関心領域設定部は、画像データ内において心臓内腔の外縁として設定された内腔ラインにより囲まれた領域を前記関心領域とする、ことを特徴とする。
【００１４】
望ましい具体例において、前記関心領域設定部は、画像データ内において、心臓内腔へ流入する血流の流路に設定された流入ラインと、心臓内腔から流出する血流の流路に設定された流出ラインと、心臓内腔の外縁として設定された前記内腔ラインと、により囲まれた領域を前記関心領域とする、ことを特徴とする。
【００１５】
望ましい具体例において、前記エネルギー演算部は、心臓内腔へ流入する血流の流路に設定された流入ラインに対応した血流の速度情報に基づいて、心臓内腔へ流入する血流の運動エネルギーである流入エネルギーを算出し、心臓内腔から流出する血流の流路に設定された流出ラインに対応した血流の速度情報に基づいて、心臓内腔から流出する血流の運動エネルギーである流出エネルギーを算出し、前記表示処理部は、各時相ごとに算出される流入エネルギーと流出エネルギーの少なくとも一方を複数時相に亘って示したグラフを形成する、ことを特徴とする。
【００１６】
望ましい具体例において、前記速度情報生成部は、生体内のドプラ情報と画像データの少なくとも一方に基づいて、心臓内腔の外縁として設定された内腔ライン上における血流と心筋の少なくとも一方の速度情報を生成し、前記エネルギー演算部は、前記内腔ライン上における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーを算出する、ことを特徴とする。
【００１７】
望ましい具体例において、前記速度情報生成部は、前記内腔ライン上の複数サンプル点について、各サンプル点ごとに血流または心筋の速度情報を生成し、前記エネルギー演算部は、前記内腔ライン上の複数サンプル点における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーの総量を算出する、ことを特徴とする。
【００１８】
望ましい具体例において、前記エネルギー演算部は、前記内腔ライン上の複数サンプル点における速度情報から、心臓内腔の内側に向かう速度情報に基づく運動エネルギーの総量と、心臓内腔の外側に向かう速度情報に基づく運動エネルギーの総量を算出する、ことを特徴とする。
【００１９】
また、上記目的にかなう好適な情報処理装置は、超音波を利用して得られた生体内のドプラ情報に基づいて、血流の速度情報を得る速度情報生成部と、超音波を利用して得られた生体に係る画像データ内において、心臓内腔に対応した関心領域を設定する関心領域設定部と、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出するエネルギー演算部と、を有することを特徴とする。
【００２０】
上記情報処理装置は、コンピュータ等により実現することができる。例えば、超音波を利用して得られた生体内のドプラ情報に基づいて、血流の速度情報を得る速度情報生成機能と、超音波を利用して得られた生体に係る画像データ内において、心臓内腔に対応した関心領域を設定する関心領域設定機能と、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出するエネルギー演算機能と、をコンピュータに実現させるプログラムにより、コンピュータを上記情報処理装置として機能させることができる。なお、そのプログラムは、例えば、ディスクやメモリなどのコンピュータが読み取り可能な記憶媒体に記憶され、その記憶媒体を介してコンピュータに提供されてもよいし、インターネットなどの電気通信回線を介してコンピュータに提供されてもよい。
【発明の効果】
【００２１】
本発明により、超音波を利用して心臓の診断情報を得る改良技術が提供される。例えば本発明の好適な態様によれば、心臓内腔の血流内で失われるエネルギーの損失量が算出され、心臓や血管の機能を評価する際の一つの評価値とすることができる。そのため、エネルギーの損失量に基づいて心臓や血管の機能を評価することが可能になり、さらに好適には、心臓内腔に流入する血流の運動エネルギーや心臓内腔から流出する血流の運動エネルギー等を評価値に加えることにより、心臓や血管の機能について、複数の評価値に基づいた総合的な診断が可能になる。
【図面の簡単な説明】
【００２２】
【図１】本発明の実施において好適な超音波診断装置の全体構成図である。
【図２】心臓内腔に対応した関心領域の具体例を示す図である。
【図３】流入エネルギーと流出エネルギーの算出例を説明するための図である。
【図４】エネルギー損失等の時間的な変化を示すグラフの具体例を示す図である。
【図５】エネルギー損失の分布に係る表示画像の具体例を示す図である。
【図６】内腔ラインから運動エネルギーを算出する具体例を示す図である。
【図７】運動エネルギー総量の時間的な変化を示すグラフの具体例を示す図である。
【発明を実施するための形態】
【００２３】
図１は、本発明の実施において好適な超音波診断装置の全体構成図である。図１の超音波診断装置は、心臓に係る診断情報を生成する機能を備えており、特に、心臓内腔の血流内で失われるエネルギーの損失量等を算出することができる。
【００２４】
プローブ１０は、生体内の心臓を含む領域に超音波を送受波する超音波探触子である。プローブ１０は、複数の振動素子を備えており、複数の振動素子が電子的に走査制御されて、心臓を含む空間内で超音波ビームが走査される。プローブ１０は、例えば、医師等のユーザ（検査者）に把持されて被検者の体表面上に当接して用いられる。なお、プローブ１０は、被検者の体腔内に挿入して用いられるものであってもよいし、電子的な走査と機械的な走査とを組み合わせた探触子であってもよい。
【００２５】
送受信部１２は、送信ビームフォーマーおよび受信ビームフォーマーとしての機能を備えている。つまり、送受信部１２は、プローブ１０が備える複数の振動素子の各々に対して送信信号を出力することにより送信ビームを形成し、さらに、複数の振動素子から得られる複数の受波信号に対して整相加算処理などを施して受信ビームを形成する。これにより、超音波ビーム（送信ビームと受信ビーム）が走査面内において走査され、超音波ビームに沿って受信信号が形成される。なお、超音波の受信信号を得るにあたって、超音波ビームが三次空間内で立体的に走査されてもよいし、送信開口合成等の技術が利用されてもよい。
【００２６】
画像形成処理部２０は、走査面内から得られる超音波の受信信号に基づいて、超音波画像用のデータ（画像データ）を形成する。画像形成処理部２０は、例えば、超音波の受信信号に対して、検波処理やフィルタ処理やＡＤ変換処理等を施すことにより、Ｂモード画像用のフレームデータを形成する。もちろん、Ｂモード画像以外の公知の超音波画像に係る画像データが形成されてもよい。
【００２７】
ドプラ処理部３０は、超音波ビームに沿って得られる受信信号に含まれるドプラシフト量を計測する。ドプラ処理部３０は、例えば公知のドプラ処理により、血流によって超音波の受信信号内に生じるドプラシフトを計測し、血流についての超音波ビーム方向の速度情報を得る。
【００２８】
速度ベクトル演算部４０は、血流についての超音波ビーム方向の速度情報から、走査面内における２次元の速度ベクトルの分布を形成する。速度ベクトル演算部４０は、例えば特許文献１（特開２０１３−１９２６４３号公報）に説明されるように、血流についての超音波ビーム方向の速度情報に加えて、心臓壁の運動情報を利用して、走査面内の各位置における血流の２次元速度ベクトルを得る。
【００２９】
なお、超音波ビーム方向に沿った１次元の速度情報から、走査面内における２次元の速度ベクトルの分布を形成するにあたっては、公知の様々な手法を利用することができる。もちろん、互いに方向が異なる２本の超音波ビームを形成して、２本の超音波ビームの各々から速度情報を得て、２次元の速度ベクトルを形成するようにしてもよい。
【００３０】
速度ベクトル演算部４０は、超音波が送受される空間に対応した演算用座標系において複数のサンプル点について、各サンプル点ごとに速度ベクトルを得る。例えば、演算用座標系をｘｙｚ直交座標系で表し、超音波の走査面に対応したｘｙ平面内で、各サンプル点ごとに速度ベクトルを得て２次元の速度ベクトルの分布を形成する。なお、超音波の走査に対応した走査座標系、例えばビームの深さ方向ｒとビームの走査方向θによるｒθ座標系において２次元の速度ベクトルの分布が形成されてもよい。
【００３１】
関心領域設定部５０は、画像形成処理部２０による処理で得られた画像データ内において、心臓内腔に対応した関心領域を設定する。関心領域設定部５０は、画像データ内において、心臓内腔の外縁として設定された内腔ラインと、心臓内腔へ流入する血流の流路に設定された流入ラインと、心臓内腔から流出する血流の流路に設定された流出ラインにより囲まれた領域を関心領域とする。関心領域設定部５０における処理と関心領域の具体例については後に詳述する。
【００３２】
内腔ライン速度演算部６０は、画像形成処理部２０において形成された画像データに基づいて、内腔ライン上における心筋（心臓壁）の速度情報を生成する。内腔ライン速度演算部６０は、内腔ライン上の複数サンプル点について、各サンプル点ごとに心筋の速度情報を生成する。
【００３３】
内腔ライン速度演算部６０は、例えば、複数フレームに亘って得られる画像データのフレーム間において、画像データの画素値（輝度値等）に基づく相関演算等を利用したパターンマッチングにより、内腔ライン上の各サンプル点ごとに、複数フレームに亘ってそのサンプル点の移動位置を２次元平面内で追跡する。これにより、各サンプル点ごとに、２次元的な移動情報が得られ、例えば、フレーム間における移動量とフレーム間の時間に基づいて、２次元の速度ベクトルが算出される。画像データがｘｙ直交座標系に対応したデータであれば、ｘｙ直交座標系における速度ベクトルが算出され、画像データがｒθ座標系に対応したデータであれば、ｒθ座標系における速度ベクトルが算出される。
【００３４】
なお、速度ベクトル演算部４０において、特許文献１に説明される手法により血流の２次元速度ベクトルを得る場合には、心臓壁の運動情報を利用するが、この場合には、心臓壁の運動情報として、内腔ライン速度演算部６０において算出された内腔ライン上の各サンプル点における速度ベクトルを利用してもよい。
【００３５】
このように、速度ベクトル演算部４０と内腔ライン速度演算部６０が速度情報生成部として機能し、心臓心腔内における血流の速度情報（速度ベクトル）や心筋（心臓壁）の速度情報（速度ベクトル）が生成される。
【００３６】
エネルギー演算部７０は、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出する。また、エネルギー演算部７０は、心臓内腔へ流入する血流の流路に設定された流入ラインに対応した血流の速度情報に基づいて、心臓内腔へ流入する血流の運動エネルギーである流入エネルギーを算出し、心臓内腔から流出する血流の流路に設定された流出ラインに対応した血流の速度情報に基づいて、心臓内腔から流出する血流の運動エネルギーである流出エネルギーを算出する。さらに、エネルギー演算部７０は、内腔ライン上における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーを算出する。エネルギー演算部７０における処理の具体例については後に詳述する。
【００３７】
表示処理部８０は、画像形成処理部２０から得られる超音波画像の画像データと、エネルギー演算部７０から得られるエネルギーに係る算出値などに基づいて、表示画像を形成する。表示処理部８０において形成された表示画像は表示部８２に表示される。
【００３８】
制御部９０は、図１に示す超音波診断装置内を全体的に制御する。なお、図１の超音波診断装置は、例えば、マウス、キーボード、トラックボール、タッチパネル、ジョイスティック等の操作デバイスを備えていることが望ましい。そして、制御部９０による全体的な制御には、操作デバイス等を介してユーザから受け付けた指示も反映される。
【００３９】
図１に示す構成（符号を付した各部）のうち、送受信部１２，画像形成処理部２０，ドプラ処理部３０，速度ベクトル演算部４０，関心領域設定部５０，内腔ライン速度演算部６０，エネルギー演算部７０，表示処理部８０の各部は、例えば電気電子回路やプロセッサ等のハードウェアを利用して実現することができ、その実現において必要に応じてメモリ等のデバイスが利用されてもよい。また、上記各部に対応した機能が、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現されてもよい。表示部８２の好適な具体例は液晶ディスプレイ等である。制御部９０は、例えば、ＣＰＵやプロセッサやメモリ等のハードウェアと、ＣＰＵやプロセッサの動作を規定するソフトウェア（プログラム）との協働により実現することができる。
【００４０】
図１の超音波診断装置の概要は以上のとおりである。次に、図１の超音波診断装置により実現される機能の具体例について詳述する。なお、図１に示した構成（符号を付した各部）については、以下の説明において図１の符号を利用する。
【００４１】
図２は、心臓内腔に対応した関心領域の具体例を示す図である。図２には、画像形成処理部２０において得られる画像データ２２の具体例が図示されており、図２の画像データ２２内には、周囲を心筋や弁によって囲まれた心臓内腔が含まれている。
【００４２】
関心領域設定部５０は、画像データ２２内において、心臓内腔の外縁として設定された内腔ライン５２と、心臓内腔へ流入する血流の流路に設定された流入ライン５４と、心臓内腔から流出する血流の流路に設定された流出ライン５６により囲まれた領域を関心領域とする。
【００４３】
内腔ライン５２は、心臓内腔の外縁に対応した複数のトレースポイントに基づいて形成される。例えば、画像データ２２に対応した表示画像が表示部８２に表示され、医師等のユーザがその表示画像を見ながら操作デバイス等を利用して、心臓内腔と心筋の境界上または境界付近に複数のトレースポイントを設定する。そして、ユーザによって設定された複数のトレースポイントに基づいて、例えば補間処理等により複数のトレースポイントを結ぶように、内腔ライン５２が形成される。なお、画像データ２２に対する二値化処理等の画像処理により、心臓内腔と心筋の境界が特定され、その境界に沿って内腔ライン５２が形成されてもよい。
【００４４】
流入ライン５４と流出ライン５６は、医師等のユーザにより設定される。例えば、画像データ２２に対応した表示画像を見ながら、医師等のユーザが、始点Ｓと終点Ｅの位置を指定することにより、流入ライン５４と流出ライン５６を設定する。
【００４５】
ユーザにより流入ライン５４と流出ライン５６が初期設定されると、関心領域設定部５０は、内腔ライン５２と流入ライン５４を接続するように流入ライン５４を修正し、内腔ライン５２と流出ライン５６を接続するように流出ライン５６を修正し、さらに、流入ライン５４と流出ライン５６とを接続する。
【００４６】
例えば、関心領域設定部５０は、流入ライン５４の始点Ｓを、その始点Ｓに最も近い内腔ライン５２上のサンプル点の位置に移動する。また、関心領域設定部５０は、流出ライン５６の終点Ｅを、その終点Ｅに最も近い内腔ライン５２上のサンプル点の位置に移動する。さらに、関心領域設定部５０は、流入ライン５４の終点Ｅと流出ライン５６の始点Ｓとを結ぶ直線または曲線を形成する。
【００４７】
こうして、関心領域設定部５０は、内腔ライン５２と流入ライン５４と流出ライン５６によって囲まれた領域を形成し、その領域を関心領域として設定する。
【００４８】
心臓内腔に対応した関心領域が設定されると、エネルギー演算部７０は、心臓内腔に対応した関心領域内における血流の速度情報に基づいて、心臓内腔の血流内で失われたエネルギーの損失量を算出する。
【００４９】
エネルギー演算部７０は、関心領域内の複数サンプル点について、次式に基づいて、各サンプル点ごとにエネルギー損失を算出し、さらに、関心領域内の複数サンプル点におけるエネルギー損失の総和を算出する。
【００５０】
【数１】
2015208522.tif 000003
【００５１】
ｘｙ座標系を利用した演算であれば、数１式において、例えばｉ方向がｘ方向となりｊ方向がｙ方向となる。また、ｒθ座標系を利用した演算であれば、例えばｉ方向がｒ方向となりｊ方向がθ方向となる。なお、ｒθ座標系を利用した演算においては、数１式がｒθ座標系に適した数式に変更されてもよい。
【００５２】
エネルギー演算部７０は、例えば数１式に基づいて、複数時相に亘って各時相ごとにエネルギー損失の総和を算出する。
【００５３】
また、エネルギー演算部７０は、流入ライン５４に対応した血流の速度情報に基づいて心臓内腔へ流入する血流の運動エネルギーである流入エネルギーを算出し、流出ライン５６に対応した血流の速度情報に基づいて、心臓内腔から流出する血流の運動エネルギーである流出エネルギーを算出する。
【００５４】
図３は、流入エネルギーと流出エネルギーの算出例を説明するための図である。流入エネルギーを算出するにあたり、エネルギー演算部７０は、まず、流入ライン５４を複数のサンプル直線（長さｄＬ）に分割する。各サンプル直線の長さｄＬは、例えば、画像データの１ピクセル程度の大きさが望ましい。
【００５５】
次に、エネルギー演算部７０は、各サンプル直線ごとに、そのサンプル直線上の速度ベクトルＶを得る。例えば、各サンプル直線上における１つのサンプル点、例えばサンプル直線の中心点における速度ベクトルＶが、そのサンプル直線の速度ベクトルＶとされる。もちろん、各サンプル直線上における複数サンプル点に基づいて、例えば複数サンプル点における速度ベクトルＶのベクトル平均が、そのサンプル直線の速度ベクトルＶとされてもよい。
【００５６】
さらに、エネルギー演算部７０は、各サンプル直線ごとに、そのサンプル直線の速度ベクトルＶに基づいて、次式により、運動エネルギーを算出する。
【００５７】
【数２】
2015208522.tif 000004
【００５８】
そして、エネルギー演算部７０は、流入ライン５４を構成する複数のサンプル直線についての運動エネルギーの総和を算出し、その総和を流入エネルギーとする。エネルギー演算部７０は、複数時相に亘って各時相ごとに流入エネルギーを算出する。
【００５９】
また、エネルギー演算部７０は、流出エネルギーを算出するにあたり、流出ライン５６を複数のサンプル直線（長さｄＬ）に分割し、流入エネルギーの場合と同様な処理により流出エネルギーを算出する。つまり、エネルギー演算部７０は、流出ライン５６の各サンプル直線ごとに、数２式により運動エネルギーを算出し、流出ライン５６を構成する複数のサンプル直線についての運動エネルギーの総和を算出し、その総和を流出エネルギーとする。エネルギー演算部７０は、複数時相に亘って各時相ごとに流出エネルギーを算出する。
【００６０】
なお、流入ライン５４のみが設定されて流入エネルギーのみが算出されてもよいし、流出ライン５６のみが設定されて流出エネルギーのみが算出されてもよい。
【００６１】
エネルギー演算部７０において、エネルギー損失と流入エネルギーと流出エネルギーの算出値が得られると、表示処理部８０は、それらの算出値を示す表示画像を形成する。表示処理部８０は、表示画像として、例えばエネルギー損失等の時間的な変化を示すグラフを形成する。
【００６２】
図４は、エネルギー損失等の時間的な変化を示すグラフの具体例を示す図である。図４には、横軸を時間軸として縦軸に算出値を示したグラフが図示されている。
【００６３】
表示処理部８０は、各時相ごとに算出されるエネルギー損失の算出値（関心領域内における総和）を複数時相に亘って示したグラフを形成する。これにより、図４に実線で示すエネルギー損失の波形が得られる。
【００６４】
また、表示処理部８０は、各時相ごとに算出される流入エネルギーの算出値（流入ライン上における総和）を複数時相に亘って示したグラフと、各時相ごとに算出される流出エネルギーの算出値（流出ライン上における総和）を複数時相に亘って示したグラフを形成する。これにより、図４に一点鎖線で示す流入エネルギーの波形と破線で示す流出エネルギーの波形が得られる。
【００６５】
なお、流入エネルギーと流出エネルギーとエネルギー損失の波形は、互いに時間軸を揃えて、つまり互いに同じ時相に対応した算出値を同じ時相に対応付けて表示されることが望ましい。もちろん、流入エネルギーと流出エネルギーとエネルギー損失の波形が、それぞれ個別に表示されてもよい。
【００６６】
また、表示処理部８０は、心臓内腔に関する超音波画像の画像データに対して、関心領域内の各サンプル点ごとに、そのサンプル点のエネルギー損失に応じた表示処理を施すことにより、心臓内腔におけるエネルギー損失の分布状態を示した表示画像を形成してもよい。
【００６７】
図５は、エネルギー損失の分布に係る表示画像の具体例を示す図である。図５＜Ａ＞には、関心領域が設定された心臓内腔に関する超音波画像（Ｂモード画像）に対して、関心領域内の複数サンプル点について、各サンプル点ごとにそのサンプル点におけるエネルギー損失値に応じた色相により色付け処理を施した表示画像が図示されている。
【００６８】
表示処理部８０は、例えば、エネルギー損失値が比較的大きな各サンプル点の画像部分に黄色を基調とした色づけ処理を施し、エネルギー損失値が比較的小さな各サンプル点の画像部分に赤色を基調とした色づけ処理を施すことにより、図５＜Ａ＞の表示画像を形成する。もちろん、黄色と赤色とは異なる色を基調としてエネルギー損失値の大きさが表現されてもよいし、輝度の大きさや模様の変化等によりエネルギー損失値の大きさが表現されてもよい。
【００６９】
また、表示処理部８０は、関心領域内の複数サンプル点における速度ベクトルに基づいて、血流の流れを示す流線を形成し、図５＜Ａ＞に示す画像上に、流線を重ねて表示することにより、図５＜Ｂ＞に示す表示画像を形成してもよい。
【００７０】
なお、流線は、例えば、複数の開始点について、各開始点ごとにその開始点を起点として２次元の速度ベクトルの分布に従って血流の流れを追跡することにより得ることができる。例えば、各開始点から、その開始点の位置における血流の速度ベクトルの方向に進んで追跡点が探索され、さらに、その追跡点の位置における血流の速度ベクトルが参照されて、その速度ベクトルの方向に進んで次の追跡点が探索される。こうして、各開始点ごとに血流の速度ベクトルの分布に従って次々に追跡点ＴＰが探索され、血流の流線が形成される。これにより、図５＜Ｂ＞に示すように複数の開始点に対応した複数の流線を形成することができる。
【００７１】
次に、心臓内腔に及ぼされる運動エネルギーについて説明する。エネルギー演算部７０は、内腔ライン上における速度情報に基づいて、心臓内腔に及ぼされる運動エネルギーを算出する。
【００７２】
図６は、内腔ライン５２から運動エネルギーを算出する具体例を示す図である。運動エネルギーの算出には内腔ライン５２上の複数のトレースポイントＰが利用される。また、互いに隣接する２つのトレースポイント（Ｐ１，Ｐ２）を結ぶ直線が演算ラインとされ、演算ラインの中点（xh，yh）を始点として、演算ラインに対して直交する法線ベクトル（xd，yd）が設定される。そして、各演算ラインごとに、次式により、運動エネルギーが算出される。
【００７３】
【数３】
2015208522.tif 000005
【００７４】
数３式では、演算ライン上における複数サンプル点について、各サンプル点ごとに速度ベクトル（Ｕ，Ｖ）に基づいて運動エネルギーが算出され、演算ライン上における複数サンプル点についての運動エネルギーの総和をサンプル総数で割り算し、さらに面積変化Ｑを乗算している。面積変化Ｑは、演算ライン上の中点における流速ベクトル（wＵ，wＶ）と法線ベクトル（xd，yd）の内積（（xd×wＵ）＋（yd×wＶ））により得られる。
【００７５】
なお、演算ライン上における各サンプル点の速度ベクトル（Ｕ，Ｖ）は、演算ライン上（内腔ライン５２上）の画像データを利用したパターンマッチングにより算出されてもよいし、演算ライン上の近傍における血流の速度ベクトルに基づいて算出されてもよい。
【００７６】
エネルギー演算部７０は、内腔ライン５２内における複数の演算ラインについて、各演算ラインごとに、数３式に基づいて運動エネルギーを算出する。そして、例えば複数の演算ライン（望ましくは全ての演算ライン）から得られる運動エネルギーの総和を算出し、その総和を心臓内腔に及ぼされる運動エネルギーの総量とする。エネルギー演算部７０は複数時相に亘って各時相ごとに運動エネルギーの総量を算出する。
【００７７】
図７は、運動エネルギー総量の時間的な変化を示すグラフの具体例を示す図である。図７には、横軸を時間軸として縦軸に算出値を示したグラフが図示されている。
【００７８】
表示処理部８０は、各時相ごとに算出される運動エネルギーの総量を複数時相に亘って示したグラフを形成する。これにより、図７に破線で示す運動エネルギーの波形が得られる。なお、表示処理部８０は、各時相ごとに算出されるエネルギー損失の算出値を複数時相に亘って示したグラフを形成し、図７に示すように、運動エネルギーの波形（破線）と共にエネルギー損失の波形（実線）を表示させてもよい。もちろん、運動エネルギーの波形のみを個別に表示させてもよい。また、エネルギー損失と運動エネルギーの比率（エネルギー損失／運動エネルギー）が算出され、心臓の負荷状態を評価する指標値として表示されてもよい
【００７９】
図６に戻り、エネルギー演算部７０は、運動エネルギーの総量を算出するにあたり、心臓内腔の内側に向かう速度情報に基づく運動エネルギーの総量と、心臓内腔の外側に向かう速度情報に基づく運動エネルギーの総量を算出してもよい。
【００８０】
例えば、数３式において、面積変化Ｑが正（プラス）となる複数の演算ラインについての運動エネルギーの総和と、面積変化Ｑが負（マイナス）となる複数の演算ラインについての運動エネルギーの総和が、別々に算出されてもよい。また、運動エネルギーの波形を形成する際に、各時相ごとに、面積変化Ｑが正となる運動エネルギーの総和と、面積変化Ｑが負となる運動エネルギーの総和のうち、絶対値の大きい方を選択して、複数時相に亘る運動エネルギー波形を形成してもよい。もちろん、各時相ごとに、面積変化Ｑが正となる運動エネルギーの総和と、面積変化Ｑが負となる運動エネルギーの総和について、２つの総和の合計値や差分値を算出し、合計値や差分値の波形を形成してもよい。
【００８１】
図１の超音波診断装置を利用することにより、例えば、流入エネルギーと流出エネルギーに基づいて血流の流量が十分か否かを診断し、さらに、エネルギー損失に基づいて心臓の負担（心負荷）が大きいか否かを診断すること等が可能になる。
【００８２】
例えば、流入エネルギーや流出エネルギーが比較的大きければ、例えば基準値よりも大きければ、被検体の体内において十分な血流の流量を維持できていることが予想される。また、十分な血流の流量が維持できているうえに、心臓内腔の血流内におけるエネルギー損失が比較的小さければ、例えば基準値よりも小さければ、心臓が正常である可能性が高い。しかし、十分な血流の流量が維持できていたとしても、心臓内腔の血流内におけるエネルギー損失が大きいと、例えば基準値よりも大きいと、十分な血流の流量を維持するために、心臓が血流に大きな運動エネルギーを与えることとなり、心臓の負担（心負荷）が大きいことが予想される。
【００８３】
なお、上記診断例はあくまでも一例に過ぎない。また、図１の超音波診断装置を利用した診断は、医師等の専門家による指導に基づいて、慎重に行われるべきことは言うまでもない。
【００８４】
以上、本発明の実施において好適な超音波診断装置について説明したが、例えば、図１に示した速度ベクトル演算部４０と関心領域設定部５０と内腔ライン速度演算部６０とエネルギー演算部７０と表示処理部８０のうちの少なくとも一つをコンピュータにより実現し、そのコンピュータを本発明に係る情報処理装置として機能させてもよい。
【００８５】
なお、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【符号の説明】
【００８６】
１０  プローブ、１２  送受信部、２０  画像形成処理部、３０  ドプラ処理部、４０  速度ベクトル演算部、５０  関心領域設定部、６０  内腔ライン速度演算部、７０  エネルギー演算部、８０  表示処理部、８２  表示部、９０  制御部。
【図６】
2015208522.tif 000006
【図７】
2015208522.tif 000007
【図１】
2015208522.tif 000008
【図２】
2015208522.tif 000009
【図３】
2015208522.tif 000010
【図４】
2015208522.tif 000011
【図５】
2015208522.tif 000012