JP2007117611A - 心機能解析装置及びその方法 - Google Patents

Abstract

【課題】心機能の異常を明確に反映する評価値を取得する心機能解析装置１０を提供する。
【解決手段】心臓を被写体とした時系列の画像データからの運動を解析する心機能解析装置１０であって、時系列の画像データから心筋の厚みにおける時系列のストレイン値を得るストレイン値取得部１６と、時系列ストレイン値を正規化する正規化ストレイン値算出部１８と、正規化ストレイン値の時系列データを出力する出力部２２とを有する。
【選択図】 図１

Description

本発明は、心機能解析装置及び心機能解析方法に係わり、特に診断に適した形でストレイン値の時系列データを出力できるようにした心機能解析装置及び心機能解析方法に関する。

心筋梗塞や狭心症といった心臓病の詳細な診断を目的として、心臓を被写体とする時系列の画像データから心臓の動きを測定し、解析する心機能解析装置について開発がなされている。

特許文献１には、心臓の左心室内腔の大きさの時間的な変化に基づき心機能の評価値を算出するシステムについての記載がある。左心室内腔の大きさは断面像においては内腔断面積であり、内腔断面積は心内膜境界（心筋壁の内側境界）を求め、その心内膜境界の形状に基づき算出されることになる。つまり、左心室内腔の大きさの時間的な変化とは、心内膜境界位置の時間的な変化から算出されるものである。

しかし内膜境界の変化に基づく方法では、心筋の異常を明確に反映しない場合があった。例えば、心筋の拡張／収縮の機能が局部的に低下した場合等において、周辺の正常な心筋の動きに引っ張られることによって異常部分についても内膜境界の動きは明確には低下せずに、その結果、内膜境界位置の変化から算出される評価値に明確な違いが現れないといった場合があった。

また、心筋を局所毎に分割して評価する場合には、分割された内腔断面積を算出するための中心点が必要となる。中心点を固定したままであると、心臓全体が平行移動するような場合に、分割内腔断面積は心臓全体の移動の影響を受けてしまうため、影響を受けないように心臓全体の移動に応じて中心点を定める必要性が生じていた。
特開２００３−３２５５２１公報

上述したように、従来技術には、心機能の評価値が心筋の異常を明確に反映しない場合があるという問題点があった。

また、心筋を局所毎に分割して評価する場合には、心臓全体の移動等の動きの影響を受けないように中心点を定める必要があるという問題点があった。

そこで、本発明は、上記問題点を解決するためになされたものであって、心機能の異常を明確に反映する評価値を得ることを可能とし、かつ心筋を局所毎に分割して評価する際の中心点の設定を必要としない心機能解析装置及びその方法を提供する。

本発明は、心臓を被写体とした時系列の画像データを用いて前記心臓の運動を解析する心機能解析装置であって、前記時系列の画像データに基づいて前記心臓上の部位における時系列のストレイン値を取得するストレイン値取得部と、前記時系列のストレイン値を正規化して時系列の正規化ストレイン値を算出する正規化部と、前記時系列の正規化ストレイン値を出力する出力部と、を備えることを特徴とする心機能解析装置である。

本発明によれば、正規化ストレイン値の時間的変化を算出することにより心機能の異常を明確に反映することが可能であり、また、心筋を局所毎に分割して評価する際の中心点の設定を必要としない。

以下、本発明の一実施形態の心機能解析装置１０について図面に基づいて説明する。

（１）心機能解析装置１０の構成
図１は、本実施形態に係わる心機能解析装置１０を示すブロック図である。

心機能解析装置１０は、時系列の画像データを入力する画像入力部１２と、画像データを記憶する画像バッファ１４と、画像データから心臓上のストレイン値の時系列データを得るストレイン値取得部１６と、時系列ストレイン値データを正規化する正規化ストレイン値算出部１８と、ストレイン値や正規化ストレイン値のデータを記憶するメモリ２０と、正規化ストレイン値の時系列データを出力する出力部２２とを備えている。

心機能解析装置１０は、例えば、汎用のコンピュータ装置１０を基本ハードウェアとして用いることでも実現することが可能である。すなわち、画像入力部１２、ストレイン値取得部１６、正規化ストレイン値算出部１８及び出力部２２は、上記のコンピュータに搭載されたプロセッサにプログラムを実行させることにより実現することができる。このとき、心機能解析装置１０は、上記のプログラムをコンピュータ装置１０に予めインストールすることで実現してもよいし、ＣＤ−ＲＯＭなどの記憶媒体に記憶して、あるいはネットワークを介して上記のプログラムを配布して、このプログラムをコンピュータに適宜インストールすることで実現してもよい。また、画像バッファ１４及びメモリ２０は、上記のコンピュータに内蔵あるいは外付けされたメモリ２０、ハードディスクもしくはＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＡＭ、ＤＶＤ−Ｒなどの記憶媒体などを適宜利用して実現することができる。

（２）心機能解析装置１０の動作
次に、図２のフローチャートを用いて、心機能解析装置１０の動作について説明する。

（２−１）画像入力部１２
まず、ステップ１において、画像入力部１２では時系列画像データを入力する。時系列画像データは超音波診断装置、Ｘ線ＣＴ装置、ＭＲＩ装置等の医用画像診断装置から入力しても良いし、画像サーバー等に保存された画像データを入力としても良い。入力された画像データは画像バッファ１４に格納される。

（２−２）ストレイン値取得部１６
ステップ２において、ストレイン値取得部１６が、時系列画像データから心臓の時系列のストレイン値を例えば下記の方法で取得する。

（２−２−１）第１の取得方法
時系列のストレイン値の取得の第１の方法について図３〜図６に基づいて説明する。

第１の方法は、次のような画像解析の手順によって行う。

ストレイン値取得のための画像解析手順は、最初に画像中に追跡点を配置する。追跡点は、輝度値のピーク位置やコーナー位置等の画像パターンが特徴的な位置を選択すると良い。

次に、追跡点の位置を基準にしてその周辺の画像パターンをテンプレート画像とし、次フレーム画像との間でテンプレートマッチングを行う。テンプレートマッチングではテンプレート画像と最も類似した位置を探す。したがって、現フレームにおける各追跡点が次フレームにおいてどこに移動したかが検出される（図３参照）。

これらの追跡点の配置とテンプレートマッチングの処理を、フレーム１とフレーム２、フレーム２とフレーム３というように隣接フレーム間で行うことにより、時系列画像における各追跡点上の動きが検出されることになる。

次に、ストレイン値Ｓを式（１）で計算する。まずストレイン値の定義について説明する。「ストレイン値」は計測対象の長さの変化を基準時間における長さで割った値として定義される。

Ｓ＝（Ｌ−Ｌ_ｒｅｆ）／Ｌ_ｒｅｆ （１）

心臓の壁厚を計測対象とした場合には、基準時間における壁厚をＬｒｅｆ、時間ｔにおける壁厚をＬ_ｔと表すと、ストレイン値ＷＴは式（２）のように書くことができる。

ＷＴ＝（Ｌ_ｔ−Ｌ_ｒｅｆ）／Ｌ_ｒｅｆ （２）

この場合には、このストレイン値ＷＴは壁厚増加率、または、壁厚変化率と呼ばれることもある。

また、これらのストレイン値を１００倍し、百分率で表現することもある。

このストレイン値を計測するために、本実施形態ではまず基準となるフレームを指定し、そのフレームで計測対象を決める。計測対象の決定は、計測対象の両端点の位置を指定することで行うと良い（図４参照）。両端点の位置の指定はマウス等のポインティングデバイスを用いて行うと良い。あるいは、後述するような方法で自動的または半自動的に両端点を設定するようにしても良い。

両端点の位置が指定されたら、その位置が次フレームでどの位置に移動するかを推定する。移動先の推定は、テンプレートマッチングで得られた動きの情報を用いて行う。例えば端点位置に最も近い追跡点の動きによって端点の動きを近似する方法をとることができる（図５参照）。あるいは、端点に近い複数追跡点の動きを補間して端点の動きを推定しても良い。補間方法としては、線形補間や２次元的なスプライン補間等の方法を用いると良い。

このようにして計測対象の両端点の位置を時系列な各フレームにわたって追跡することができ、その結果、両端点間の距離を各フレームについて算出することができる。算出された両端点間の距離から、式（１）または式（２）によってストレイン値を算出する。

計測対象が複数個所ある場合には複数の両端点を設定し、上述のストレイン値の算出手順によってそれぞれの計測対象について処理を行う。心臓を対象にする場合には、心筋上に複数の計測対象を設定すると心筋上のストレイン値の分布を把握するのに好都合である。例えば、図６のように心筋の内側と外側に両端点を複数配置することで、心筋全体にわたってストレイン値を算出することができる。図６のような両端点の配置を簡便に行うためには、例えば心筋の内側と外側の輪郭を抽出または指定し、その輪郭上に端点を配置すると良い。

（２−２−２）第２の取得方法
第２のストレイン値取得の方法としては、超音波のドップラモードで得られる組織の速度情報を用いて、計測対象の両端点の次フレームでの移動先を算出しストレイン値を計算しても良い。

この方法では、超音波パルスを複数回発射し、反射パルスの遅延時間の差から対象物の移動速度を検出する。対象物が移動している場合には、複数の反射パルスに時間的ずれが生じるため速度を検出できる。検出される速度情報は超音波ビーム方向の速度成分であるので、別途移動方向を仮定または手動で教示することで実際の速度を推定することができる。このようにして得られた速度情報を用いて、心筋の境界位置を追跡してストレイン値を計算することができる。

（２−２−３）第３の取得方法
第３の方法としては、ストレイン値の算出は別途算出しておき、算出されたストレイン値のデータのみを入力する構成としても良い。

例えば、心筋の内側と外側の境界位置を手動で指示することで追跡し、心筋の厚さの変化からストレイン値を計算することもできる。

以上のような手順によって、１つまたは複数箇所のストレイン値がフレーム毎に算出される。

算出されたストレイン値はメモリ２０に格納される。

（２−３）正規化ストレイン値算出部１８
次に、ステップ３において、正規化ストレイン値算出部１８ではフレーム毎に算出された時系列ストレイン値をそれぞれ正規化する。

まず、心筋の厚さ方向のストレイン値の通常の変化の仕方について述べる。心室が最も拡張した時点から心室が収縮するにつれ心筋の厚さは増大するためストレイン値は増加する。収縮の後、拡張期に入ると心筋の厚さは減少するのでストレイン値は減少し、最初のストレイン値におおよそ戻る。

以下、その正規化方法について説明する。

（２−３−１）第１の正規化方法
第１の正規化方法としては、評価対象とする時間区間を別途定め、式（３）に示すようにその計測対象におけるストレイン値の最大変動幅と、現在のストレイン値と最小ストレイン値の差との比で正規化ストレイン値ＮＳ_１を定義する方法がある。

ＮＳ_１=（Ｓ_ｔ−Ｓ_ｍｉｎ）／（Ｓ_ｍａｘ−Ｓ_ｍｉｎ） （３）

ここで、Ｓ_ｔは時間ｔにおけるストレイン値、Ｓ_ｍａｘ，Ｓ_ｍｉｎはそれぞれその計測対象におけるストレイン値の最大値、最小値である。最もストレイン値が小さい時相を基準にした場合等、ストレイン値の最小値が０となる場合には、ストレイン値ＮＳ_２は式（４）のようになる。

ＮＳ_２=Ｓ_ｔ／Ｓ_ｍａｘ （４）

このような正規化ストレイン値の定義によれば、ストレイン値の値そのものではなく最大変動幅に対する割合として表現されることになる。

正常の動きをする心臓においても、心筋の部位毎に最大ストレイン値は異なることが想定されるため、複数の異なる心筋部分を同時に評価する場合に値の大きさに影響されずに同様の尺度で評価できることになる。とりわけストレイン値の時間的変化の急緩を比較評価したい場合に効果を発揮する。

図７はこのことを模式的に示している。図７（ａ）は、ストレイン値そのものの時間変化をグラフ表示した例であり、急峻に変化しているのかゆっくり変化しているのかを比較することは容易ではない。これに対して図７（ｂ）は正規化ストレイン値の時間変化をグラフ表示した例であり、ストレイン値の最大値の大きさによらないため、時間変化が緩やかなデータを容易に判別できる。

（２−３−２）第２の正規化方法
第２の正規化方法としては、評価対象とする時間区間を別途定め、式（５）のストレイン値ＮＳ_３のように時間区間の開始時と終了時におけるストレイン値の差によって正規化する方法としても良い。

ＮＳ_３=（Ｓ_ｔ−Ｓ_{ｓｔａｒｔ}）／（Ｓ_ｅｎｄ−Ｓ_{ｓｔａｒｔ}） （５）

ここで、Ｓ_{ｓｔａｒｔ}，Ｓ_ｅｎｄ，Ｓ_ｔはそれぞれ時間区間の開始時、終了時、時間ｔにおけるストレイン値を表している。

このようにすることで指定した時間区間全体でのストレイン値の変化幅を１００％とし、その途中時間までにおいて何％変化したかという情報を得ることができる。その結果、心臓の収縮期または拡張期あるいは任意の時間区間におけるストレイン値の時間的変化の急緩等を理解することに対して効果的である。

このようにして算出された正規化ストレイン値ＮＳはメモリ２０に格納される。

（２−４）出力部２２
最後に、ステップ４において出力部２２では、正規化ストレイン値を出力する。

例えば、正規化ストレイン値の変化のグラフ表示、設定された時間区間中の指定経過時間後の正規化ストレイン値の表示、正規化ストレイン値に応じた色の画像上計測対象位置への重畳表示の少なくとも１つを表示する構成とすると良い。

グラフ表示は、例えば図８のように横軸を時間に縦軸を正規化ストレイン値にとり、正規化ストレイン値の時間変化を示す。さらに、設定された時間区間中の指定経過時間後（図８の例では「４０％」の時間経過後が示されている）の正規化ストレイン値（図８の例では「８０％」）を表示すると、決められた時間での値を比較する場合などに効率が良い。

また、図９のように設定された正規化ストレイン値に到達する時間に応じて、その部位をカラー表示すると、正規化ストレイン値の変化と心臓部位の位置関係が判別しやすい。このカラー表示は上述の指定経過時間後の正規化ストレイン値の値に応じて行うようにしても良い。なお、図９ではカラー表示を異なる模様で表示している。

以上のように出力部２２で表示することで、心筋の正規化ストレイン値の部分毎の違い、時間的変化の仕方の違いを理解することが容易となる。

（３）変更例
本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。

（３−１）変更例１
複数の領域毎に分けて正規化ストレイン値の代表値を算出することもできる。例えば、部分領域内で計測された複数の正規化ストレイン値の平均や中央値を代表値とすることで、領域毎のスコアリング等の用途に対して効率の良い表示法が可能となる。図９はその表示法の一例でもある。

（３−２）変更例２
正規化ストレイン値の標準的な変化パターンとの違いを数値的に算出して表示しても良い。

例えば、図１０に示すように標準的な変化パターンと計測された変化パターンとのグラフ上での間の面積や、差の総和等によって違いを数値化し表示することができる。このようにすることで、正規化ストレイン値が急峻に変化した場合と、緩やかに変化した場合の違いを定量的に得ることができる。

（３−３）変更例３
時間区間の設定の仕方は任意であるが、心臓の収縮期のみを設定区間にしたり、拡張期のみを設定区間にしたり、あるいは１心拍を設定区間にしたりすることができる。

収縮期を設定時間区間とした場合には、収縮過程におけるストレイン値の変化割合、変化タイミングを判断するのに好適である。

また、拡張期を時間区間とした場合には、拡張過程における変化割合、変化タイミングを判断するのに好適である。心臓病の診断において、収縮時のストレイン値の増加の遅れや、拡張時の拡張タイミングの遅れが診断の指標となりうると言われており、正規化ストレイン値を用いることで、心臓の収縮や拡張の変化の様子の判別が容易となり、診断上に役立つ情報が可能となる。

（３−４）変更例４
上記実施形態では、心臓の壁厚（心筋の厚み方向の距離）をストレイン値の計測対象としたが、これに限らず心筋の長さ方向の距離をストレイン値の計測対象としてもよい。

（４）本実施形態の効果
正規化ストレイン値を心機能の評価に用いることで、従来心筋内側輪郭から求められる心室内の面積による評価の際に生じていた問題が解決できる。例えば、心筋の拡張／収縮の機能が局部的に低下した場合等において、周辺の正常な心筋の動きに引っ張られることによって異常部分についても内膜境界の動きは明確には低下せずに、その結果、内膜境界位置の変化から算出される評価値に明確な違いが現れないといった問題を回避でき、心機能の異常をより明確に反映する評価値を得ることが可能となる。

また、心筋を局所毎に分割して評価する場合には、従来の心室内の面積による評価では分割された内腔断面積を算出するための中心点が必要となり、中心点を固定したままであると心臓全体が平行移動するような状況において分割内腔断面積は心臓全体の移動の影響を受けてしまうため、影響を受けないように心臓全体の移動に応じて中心点を定める必要があるという問題も解決できる。すなわち、正規化ストレイン値を用いる場合には、例えば心筋に沿った長さによって分割することで局所毎に分けることができるため、そのような中心点を定める必要は無くなる。

以上のように本実施形態によれば、正規化ストレイン値の時間的変化を算出することにより心機能の異常を明確に反映する評価値を得ることが可能となり、かつ心筋を局所毎に分割して評価する際の中心点の設定を必要としない心機能解析装置１０を提供することが可能となる。

本発明の一実施形態に係わる心機能解析装置１０の構成を示すブロック図である。 本実施形態の動作を示すフローチャートである。 追跡点位置におけるフレーム間の移動検出を説明する図である。 計測対象の指定方法を説明する図である。 計測対象の端点の追跡を説明する図である。 心筋全体のストレイン値を取得する場合の計測対象の配置を説明する模式図である。 図７（ａ）は、ストレイン値そのものの時間変化をグラフ表示した例であり、図７（ｂ）は正規化ストレイン値の時間変化をグラフ表示した例である。 正規化ストレイン値の時間変化を示すグラフである。 正規化ストレイン値に到達する時間に応じて、その部位をカラー表示した例である。 標準的な変化パターンと計測された変化パターンとのグラフである。

符号の説明

１０ 心機能解析装置
１２ 画像入力部
１４ 画像バッファ
１６ ストレイン値取得部
１８ 正規化ストレイン値算出部
２０ メモリ
２２ 出力部

Claims (30)

1. 心臓を被写体とした時系列の画像データを用いて前記心臓の運動を解析する心機能解析装置であって、
   前記時系列の画像データに基づいて前記心臓上の部位における時系列のストレイン値を取得するストレイン値取得部と、
   前記時系列のストレイン値を正規化して時系列の正規化ストレイン値を算出する正規化部と、
   前記時系列の正規化ストレイン値を出力する出力部と、
   を備える
   ことを特徴とする心機能解析装置。
2. 前記ストレイン値取得部は、前記ストレイン値を前記心臓の心筋の厚み方向、または、心筋の長さ方向の距離から求める
   ことを特徴とする請求項１記載の心機能解析装置。
3. 前記時系列の画像データのうちの基準となる画像データにおいて、前記心臓上の部位に２つの追跡点を設定する追跡点設定部と、
   前記時系列の画像データの各々における、前記部位の２つの追跡点の位置を求める追跡部と、
   をさらに備え、
   前記ストレイン値取得部は、前記心臓上の部位における時系列のストレイン値を前記時系列の画像データにおける前記２つの追跡点間の距離に基づいて求める、
   ことを特徴とする請求項１に記載の心機能解析装置。
4. 前記追跡点設定部は、前記心臓の心筋の内側の輪郭上と外側の輪郭上とに１点ずつ前記追跡点を設定する、
   ことを特徴とする請求項３に記載の心機能解析装置。
5. 前記追跡点設定部は、前記心臓の心筋に沿って２つの前記追跡点を設定する
   ことを特徴とする請求項３または４に記載の心機能解析装置。
6. 前記正規化部は、ストレイン値を複数の部位間で共通な尺度の正規化ストレイン値に正規化する
   ことを特徴とする請求項１または２記載の心機能解析装置。
7. 前記正規化部は、前記正規化ストレイン値を、前記部位における解析時間区間内のストレイン値の最大変動幅と、現在のストレイン値と最小ストレイン値の差との比で定義する
   ことを特徴とする請求項１から６のいずれか一項に記載の心機能解析装置。
8. 前記正規化部は、前記正規化ストレイン値を、前記部位における解析時間区間の開始時ストレイン値と終了ストレイン値の差と、現在のストレイン値と前記開始時ストレイン値の差との比で定義する
   ことを特徴とする請求項１から６のいずれか一項に記載の心機能解析装置。
9. 前記正規化ストレイン値の時間的変化パターンと、規定の標準変化パターンとの違いを示す評価値を算出する評価値算出部をさらに備える
   ことを特徴とする請求項１から８のいずれか一項に記載の心機能解析装置。
10. 前記正規化ストレイン値が所定の値に達する時間を算出する到達時間算出部をさらに備える
    ことを特徴とする請求項１から８のいずれか一項に記載の心機能解析装置。
11. 前記出力部が、前記解析時間区間内の所定の時間における正規化ストレイン値に応じた色を、前記心臓の画像に重畳表示する
    ことを特徴とする請求項７または８記載の心機能解析装置。
12. 前記正規化部が、複数の領域に分割された心筋領域毎に前記正規化ストレイン値の代表値を算出する
    ことを特徴とする請求項１から８のいずれか一項に記載の心機能解析装置。
13. 心臓を被写体とした時系列の画像データを用いて前記心臓の運動を解析する心機能解析方法であって、
    前記時系列の画像データに基づいて前記心臓上の部位における時系列のストレイン値を取得し、
    前記時系列のストレイン値を正規化して時系列の正規化ストレイン値を算出し、
    前記時系列の正規化ストレイン値を出力する
    ことを特徴とする心機能解析方法。
14. 前記ストレイン値を、前記心臓の心筋の厚み方向、または、心筋の長さ方向の距離から求める
    ことを特徴とする請求項１３記載の心機能解析方法。
15. ストレイン値を複数の部位間で共通な尺度の正規化ストレイン値に正規化する
    ことを特徴とする請求項１３または１４記載の心機能解析方法。
16. 前記正規化ストレイン値を、前記部位における解析時間区間内のストレイン値の最大変動幅と、現在のストレイン値と最小ストレイン値の差との比で定義する
    ことを特徴とする請求項１３から１５のいずれか一項に記載の心機能解析方法。
17. 前記正規化ストレイン値を、前記部位における解析時間区間の開始時ストレイン値と終了ストレイン値の差と、現在のストレイン値と前記開始時ストレイン値の差との比で定義する
    ことを特徴とする請求項１３から１５のいずれか一項に記載の心機能解析方法。
18. 前記正規化ストレイン値の時間的変化パターンと、規定の標準変化パターンとの違いを示す評価値をさらに算出する
    ことを特徴とする請求項１３から１７のいずれか一項に記載の心機能解析方法。
19. 前記正規化ストレイン値が所定の値に達する時間をさらに算出する
    ことを特徴とする請求項１３から１７のいずれか一項に記載の心機能解析方法。
20. 前記解析時間区間内の所定の時間における正規化ストレイン値に応じた色を、前記心臓の画像に重畳表示する
    ことを特徴とする請求項１６または１７記載の心機能解析方法。
21. 複数の領域に分割された心筋領域毎に前記正規化ストレイン値の代表値を算出する
    ことを特徴とする請求項１３から１７のいずれか一項に記載の心機能解析方法。
22. 心臓を被写体とした時系列の画像データを用いて前記心臓の運動をコンピュータによって解析するための心機能解析プログラムであって、
    前記時系列の画像データに基づいて前記心臓上の部位における時系列のストレイン値を取得するストレイン値取得機能と、
    前記時系列のストレイン値を正規化して時系列の正規化ストレイン値を算出する正規化機能と、
    前記時系列の正規化ストレイン値を出力する出力機能と、
    を実行する
    ことを特徴とする心機能解析プログラム。
23. 前記正規化機能は、前記ストレイン値を、前記心臓の心筋の厚み方向、または、心筋の長さ方向の距離から求める
    ことを特徴とする請求項２２記載の心機能解析プログラム。
24. ストレイン値を複数の部位間で共通な尺度の正規化ストレイン値に正規化する
    ことを特徴とする請求項２２または２３記載の心機能解析方法。
25. 前記正規化機能は、前記正規化ストレイン値を、前記部位における解析時間区間内のストレイン値の最大変動幅と、現在のストレイン値と最小ストレイン値の差との比で定義する
    ことを特徴とする請求項２２から２４のいずれか一項に記載の心機能解析プログラム。
26. 前記正規化機能は、前記正規化ストレイン値を、前記部位における解析時間区間の開始時ストレイン値と終了ストレイン値の差と、現在のストレイン値と前記開始時ストレイン値の差との比で定義する
    ことを特徴とする請求項２２から２４のいずれか一項に記載の心機能解析プログラム。
27. 前記正規化ストレイン値の時間的変化パターンと、規定の標準変化パターンとの違いを示す評価値を算出する評価値算出機能をさらに実行する
    ことを特徴とする請求項２２から２６のいずれか一項に記載の心機能解析プログラム。
28. 前記正規化ストレイン値が所定の値に達する時間を算出する到達時間算出機能をさらに実行する
    ことを特徴とする請求項２２から２６のいずれか一項に記載の心機能解析プログラム。
29. 前記出力機能が、前記解析時間区間内の所定の時間における正規化ストレイン値に応じた色を、前記心臓の画像に重畳表示する
    ことを特徴とする請求項２５または２６記載の心機能解析プログラム。
30. 前記正規化機能が、複数の領域に分割された心筋領域毎に前記正規化ストレイン値の代表値を算出する
    ことを特徴とする請求項２２から２６のいずれか一項に記載の心機能解析プログラム。
    

Images