(57)【特許請求の範囲】
【請求項１】
投与した放射性医薬品の挙動を解析する心筋バイアビリティの評価方法において、カウント評価の統一基準として患者ごとの拡張末期の最大カウントを基準とするように設定し、この統一基準に基づいて任意の心時相の心筋断層データのカウントを評価するようにしたことを特徴とする心筋バイアビリティの評価方法。
【請求項２】
拡張末期の最大カウントを１００％に設定したことを特徴とする請求項１記載の心筋バイアビリティの評価方法。
【請求項３】
投与した放射性医薬品の挙動を解析する心筋バイアビリティの評価方法において、カウント評価の統一基準として患者ごとの拡張末期の最大カウントを基準とするように設定し、この基準に基づいて任意の断層で各心筋断層像を順次並べ、これらを時相ごとに配列したことを特徴とする心筋バイアビリティの評価方法。
【請求項４】
左室全体を網羅できるように各心筋断層像を順次並べて配列したことを特徴とする請求項３記載の心筋バイアビリティの評価方法。
【請求項５】
投与した放射性医薬品の挙動を解析する心筋バイアビリティの評価方法において、任意の心時相の短軸断層データを設定する場合、任意の一短軸断層面を基準とし、収縮または拡張に伴う他の短軸断層面における回転ずれを回転軸を中心として局所が一致するように補正するようにしたことを特徴とする請求項１，２記載の心筋バイアビリティの評価方法。
【請求項６】
予め記憶されたカウント評価の統一基準として患者ごとに設定される拡張末期の最大カウントに基づいて評価する統一基準データ設定手段と、このカウントに基づいて任意の心時相に対応する心筋断層データを設定する心時相データ設定手段と、上記各データに基づいて予め決められた分布に従って任意の各心時相の各心筋イメージを作成する心筋イメージ作成手段とを備える心筋断層画像処理装置。
【請求項７】
統一基準データ設定手段が拡張末期の最大カウントを基準として評価する拡張末期データ設定手段であって、心筋イメージ作成手段が拡張末期と任意の心時相の各心筋イメージを作成する請求項６記載の心筋断層画像処理装置。
【請求項８】
拡張末期データ設定手段は、拡張末期における最大カウントを１００％に設定することを特徴とする請求項７記載の心筋断層画像処理装置。
【請求項９】
心筋イメージ作成手段は、予め決められた分布の範囲内で任意の範囲に設定されることを特徴とする請求項６記載の心筋断層画像処理装置。
【請求項１０】
任意の心時相の心筋イメージの分布の範囲に拡張末期の心筋イメージの分布の範囲よりもレベルの高い分布を追加することを特徴とする請求項９記載の心筋断層画像処理装置。
【請求項１１】
心時相データ設定手段は、心筋の心基部から心尖部までの任意の心時相の短軸断層データを設定する場合、任意の一短軸断層面を基準として収縮または拡張に伴う他の短軸断層面の回転を補正するようにしたことを特徴とする請求項６記載の心筋断層画像処理装置。
【請求項１２】
心時相データ設定手段は、任意の断層を設定する場合、複数の時相から時相間における局所の最高カウントを合成するようにしたことを特徴とする請求項６記載の心筋断層画像処理装置。
【請求項１３】
任意の断層として短軸断層から極座標または展開図が設定されることを特徴とする請求項１２記載の心筋断層画像処理装置。
【請求項１４】
収縮末期の極座標または展開図を設定する場合、複数の時相から時相間における局所の最高カウントを合成するようにしたことを特徴とする請求項１３記載の心筋断層画像処理装置。
【請求項１５】
予め記憶されたカウント評価の統一基準として患者ごとに設定される拡張末期の最大カウントに基づいて評価する統一基準データ設定手段と、このカウントに基づいて任意の心時相に対応する心筋断層データを設定する心時相データ設定手段と、任意の心筋断層データとその他の任意の心時相の心筋断層データとの間のカウントの変化率に基づいて局所的な心機能データを設定する心機能データ設定手段と、上記各データに基づいて予め決められた分布に従って任意の心時相とその他の任意の心時相の各心筋イメージを作成するとともに心機能イメージを作成するイメージ作成手段とを備える心筋断層画像処理装置。
【請求項１６】
統一基準データ設定手段が拡張末期の最大カウントを基準として評価する拡張末期データ設定手段であって、心機能データ設定手段が任意の一つの心筋断層データとして拡張末期の心筋断層データを用い、イメージ作成手段が任意の一つの心時相として拡張末期を用いる請求項１５記載の心筋断層画像処理装置。
【請求項１７】
イメージ作成手段は、予め決められた分布の範囲内で任意の範囲に設定されることを特徴とする請求項１５記載の心筋断層画像処理装置。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、投与した放射性医薬品の挙動を解析する心筋バイアビリティの評価方法及びこの評価方法を用いた心筋断層画像処理装置に関する。本発明において心筋バイアビリティの評価には心機能の評価を含むものである。
【０００２】
【従来の技術】
近年、冠状動脈インターベンションの発達により、ますます正確な心筋バイアビリティの診断が切望されている。代表的な診断方法として、ＲＩ（放射性同位元素）を人体に投与してその挙動を動態解析するという心筋バイアビリティの評価手法がある。
【０００３】
そこでＲＩの挙動を動態解析するには、例えば、 201Ｔｌ（塩化タリウム）を用いて検査、読影を行うという 201Ｔｌ（以下単にＴｌと称する）心筋シンチグラフィがある。このＴｌシンチグラフィには、運動負荷や薬物負荷を施行する方法と安静時に行う方法とに大別される。
【０００４】
運動負荷や薬物負荷を施行するＴｌ心筋シンチグラフィは、負荷をすることで安静時のＴｌ心筋シンチグラフィでは得られない虚血を誘発させる方法である。これによって、負荷直後と負荷後３〜４時間経過させた安静後の各像を比較すると、その差が大きくなるために心筋バイアビリティの診断がしやすいという利点がある。さらに、診断能力を向上させるために、再静注法や遅延再分布法等の工夫もなされている。
【０００５】
また、安静時のＴｌ心筋シンチグラフィでは、梗塞部心筋のバイアビリティ評価について、安静時に撮像されることで心筋に負担をかけずに患者にとって好条件のシンチグラフィとなり、負荷による過度の虚血を起こさせるような心配はないという利点がある。
【０００６】
以上のＴｌ心筋シンチグラフィを支援するものが心筋断層画像処理装置である。この心筋断層画像処理装置は、ガンマシンチカメラによって負荷直後と負荷後３〜４時間経過させた安静後とに分けて撮影した心筋スペクト像（以下「心筋断層像」という）を入力し、心時相を考慮せずに上記負荷直後像と上記安静後像の各断層データの最高カウントをそれぞれ求め、各最高カウントを基準にしてそれぞれのレベルに応じて上記負荷直後像と上記安静後像をカラー表示したり、あるいは短軸断層像を極座標（同心円）や展開図でカラー表示するものである。
【０００７】
【発明が解決しようとする課題】
しかしながら、上述した従来の心筋バイアビリティの評価方法では、運動負荷を施行するＴｌ心筋シンチグラフィの場合、運動負荷量が医師の主観に頼るためにその負荷量によって心筋の虚血状態に差が生じ、再現性が乏しいという問題があった。また安静時のＴｌ心筋シンチグラフィの場合、負荷直後像との比較ができないために判別し難いとともにシンチグラフィが症例の最高値に対する相対評価しかできないという問題があった。このように、いずれのＴｌ心筋シンチグラフィでもバイアブルな心筋を見落としてしまうという要素があり、正確に心筋バイアビリティを評価することができないという大きな問題があった。
【０００８】
また従来の心筋バイアビリティの評価方法では、心拍や心時相を何ら考慮せずに診断していたので、時相ごとに起こる変化を全く把握できないものであった。例えば局所心筋の機能の遅れや進み具合、あるいはあらゆる方向での心機能等は何ら解明できないままであった。
【０００９】
また、負荷直後と負荷後３〜４時間の安静時とでは、Ｔｌが時間とともに減衰したり人体から排出されるため、カウントを同じ条件で計測することは不可能である。これにもまして、前述した従来の心筋断層画像処理装置では、負荷直後、安静時について、それぞれの計測の際に得られた最高カウントのレベルを基準にして、各レベルに応じた分布により心筋イメージが表示される。したがって、負荷直後と安静時の各心筋イメージとを比較しても、それぞれの分布の基準が相違することから、正確に心筋バイアビリティを評価するための資料としては不十分であるという問題があった。
【００１０】
さらに、特開平２ー２２６０８９号公報には、正常者に係る標準デ−タを格納する心筋スペクトエキスパートシステムを用いて、心筋断層像を極座標で整理して標準デ−タと個々の患者デ−タとを比較してデータの解釈を自動化しようとするものが開示されている。しかしながら、正常者に係る標準デ−タがそれぞれ異なる患者の心筋疾患に十分に対応できないばかりか、従来の心筋断層画像処理装置で得られる断層データに基づいているために正確な解釈を求められるものではなく、実用化には遠いものと言わざるを得なかった。
【００１１】
請求項１の発明は、上記のような問題点を解消するために研究されたもので、心拍同期法の心筋シンチグラフィによりカウント評価の統一基準として患者ごとの拡張末期の最大カウントを基準とするように設定し、この統一基準に基づいて心筋バイアビリティを正確に評価できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１２】
削除
【００１３】
請求項２の発明は、請求項１の発明の目的に加え、拡張末期と他の心時相を同一の基準で評価できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１４】
請求項３の発明は、患者ごとの拡張末期の最大カウントを基準として、任意の断層の各心筋断層データを時相ごとに順次配列して心筋血流のみならず、立体的なあらゆる方向での心機能を把握できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１５】
請求項４の発明は、請求項３の目的に加え、左室全体を網羅して時相ごとに順次表示して心機能を把握できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１６】
請求項５の発明は、収縮または拡張に伴う心筋の回転を考慮して短軸断層面の任意の各心時相の局所を一致させることができるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１７】
本発明の関連発明１は、局所心筋の機能の遅れや進みを把握できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１８】
本発明の関連発明２は、短軸断層における極座標または展開図から局所心筋の機能の遅れや進みを把握できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００１９】
本発明の関連発明３は、収縮末期における局所心筋の収縮の遅れを見落としせずに把握できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００２０】
本発明の関連発明４は、任意の局所の心機能を反映させて心筋バイアビリティを診断するための有効な資料を作成できるようにした評価方法を提供することを目的とする。
【００２１】
本発明の関連発明５は、関連発明４の目的に加え、拡張末期を基準として任意の局所の心機能を反映させて心筋バイアビリティを診断するための有効な資料を作成できるようにした評価方法を提供することを目的とする。
【００２２】
本発明の関連発明６は、関連発明４，５の発明の目的に加え、任意の領域の僅かなカウントの変化を明瞭に評価できるようにした心筋バイアビリティの評価方法を提供することを目的とする。
【００２３】
請求項６の発明は、心筋バイアビリティを残存心筋量と残存心筋の収縮能から評価できるようにした心筋断層画像処理装置を提供することを目的とする。
【００２４】
請求項７の発明は、請求項６の発明の目的に加え、拡張末期における最大カウントを統一基準として評価できるようにした心筋断層画像処理装置を提供することを目的とする。
【００２５】
請求項８の発明は、請求項７の発明の目的に加え、拡張末期における最大カウントを１００％にして評価できるようにした心筋断層画像処理装置を提供することを目的とする。
【００２６】
請求項９の発明は、請求項６の発明の目的に加え、心筋の各状態を視覚的に正確に表現できるようにした心筋断層画像処理装置を提供することを目的とする。
【００２７】
請求項１０の発明は、請求項９の発明の目的に加え、高レベルの分布により視覚性を向上させることができるようにした心筋断層画像処理装置を提供することを目的とする。
【００２８】
請求項１１の発明は、請求項６の発明の目的に加え、心筋の収縮または拡張に伴う回転のずれに対処して任意の心時相像の局所を一致させ、心筋バイアビリティを正確に表現できるようにした心筋断層画像処理装置を提供することを目的とする。
【００２９】
請求項１２の発明は、請求項６の発明の目的に加え、局所心筋の機能の遅れや進み具合を把握できるようにした心筋断層画像処理装置を提供することを目的とする。
【００３０】
請求項１３の発明は、請求項１２の発明の目的に加え、短軸断層における極座標または展開図に基づいて局所心筋の機能の遅れや進み具合を把握できるようにした心筋断層画像処理装置を提供することを目的とする。
【００３１】
請求項１４の発明は、請求項１３の発明の目的に加え、局所心筋の収縮の遅れを見落としせずに把握できるようにした心筋断層画像処理装置を提供することを目的とする。
【００３２】
請求項１５の発明は、患者ごとの拡張末期の最大カウントを統一基準として心筋バイアビリテイの診断や心機能を視覚的に理解し易くできるようにした心筋断層画像処理装置を提供することを目的とする。
【００３３】
請求項１６の発明は、請求項１５の発明の目的に加え、拡張末期の最大カウントを基準として心筋バイアビリテイの診断や心機能を視覚的に理解し易くできるようにした心筋断層画像処理装置を提供することを目的とする。
【００３４】
請求項１７の発明は、請求項１６の発明の目的に加え、分布を変更して診断を一層容易に行うことができるようにした心筋断層画像処理装置を提供することを目的とする。
【００３５】
【課題を解決するための手段】
請求項１の発明に係る心筋バイアビリティの評価方法は、カウント評価の統一基準として拡張末期の最大カウントを基準とするように設定し、この基準に基づいて任意の心時相の心筋断層データのカウントを評価するようにしたことを特徴とするものである。ここで統一基準は患者ごとに相違する相対的基準を含むが、いずれにしても断層ごとあるいは各時相ごとにばらばらではない統一した基準として用いられる。また、任意の心時相とは、拡張末期から収縮末期までの任意の収縮期と収縮末期から拡張末期までの任意の拡張期とに係わる全心周期をいい、この心時相における心筋断層データから得られるカウントが対象となるものである。ここで、拡張末期における全データ中から最大カウントが検索され、この最大カウントが任意の心時相に対する統一した基準として用いられる。本発明者は、このような統一基準を用いる評価方法を統一基準法と命名した。
【００３６】
削除
【００３７】
請求項２の発明に係る心筋バイアビリティの評価方法は、拡張末期の最大カウントを１００％に設定したことを特徴とするものである。
【００３８】
請求項３の発明に係る心筋バイアビリティの評価方法は、統一基準法による統一基準に基づいて任意の断層で各心筋断層像を順次並べ、これらを時相ごとに配列したことを特徴とするものである。任意の断層とは短軸断層や長軸面垂直断層や長軸面水平断層等をいい、これらの各心筋断層像を時系列的に配列する。本発明者は、このように各心筋断層像を順次並べて時系列的に配列する方法を断層像時相配列法と命名した。断層像時相配列法では心筋断層像として心筋断層データを含み、同一部位の局所心筋の配列も包含されるものである。
【００３９】
請求項４の発明に係る心筋バイアビリティの評価方法は、統一基準法により左室全体を網羅できるように各心筋断層像を順次並べて配列したことを特徴とするものである。断層像時相配列法では、心筋断層像として左室全体が網羅されるように順次配列することに主眼がある。
【００４０】
請求項５の発明に係る心筋バイアビリティの評価方法は、任意の心時相の短軸断層データを設定する場合、任意の一短軸断層面を基準として収縮または拡張に伴う他の短軸断層面における回転ずれを回転軸を中心として局所が一致するように補正するようにしたことを特徴とするものである。心筋は収縮に伴って回転（拡張ではその逆に回転）しているので、拡張末期と収縮末期では局所心筋が必ずしも一致しておらず、局所心筋の対応で一致させるために回転が補正される。本発明者はこのような回転の補正を回転補正法と命名した。
【００４１】
本発明の関連発明１に係る心筋バイアビリティの評価方法は、任意の心時相における任意の断層を評価する場合、複数の時相から時相間における局所の最高カウントを合成（以下「最高値合成」という）するようにしたことを特徴とするものである。収縮や拡張に伴う局所心筋の機能の遅れや進みは、複数の時相における最高値合成により把握することが可能となる。本発明者はこのような最高値合成を最高値合成法と命名した。
【００４２】
本発明の関連発明２に係る心筋バイアビリティの評価方法は、短軸断層から極座標または展開図が作成されるようにしたことを特徴とするものである。収縮や拡張に伴う局所心筋の機能の遅れや進みは、短軸断層からの複数時相間の最高値合成法により容易に把握することが可能となる。
【００４３】
本発明の関連発明３に係る心筋バイアビリティの評価方法は、収縮末期における極座標または展開図が作成されるようにしたことを特徴とするものである。局所心筋の収縮の遅れを漏れなく把握することができる。
【００４４】
本発明の関連発明４に係る心筋バイアビリティの評価方法は、患者ごとの局所的な心機能を
〔任意の局所心筋における｛（任意の心時相の心筋断層データのカウント
−他の心時相の心筋断層デ−タのカウント）
／他の心時相の心筋断層デ−タのカウント｝〕×１００％
で表されるカウントの変化率を用いて表示することを特徴とするものである。任意の心時相としては全心周期の心筋断層データを用いることができ、任意の心時相と他の心時相は異なる心時相である。
【００４５】
本発明の関連発明５に係る心筋バイアビリティの評価方法は、他の心時相の心筋断層デ−タのカウントが拡張末期の心筋断層デ−タのカウントであることを特徴とするものである。
【００４６】
本発明の関連発明６に係る心筋バイアビリティの評価方法は、カウントの変化率を任意の分布に基づいて変更表示可能にしたことを特徴とするものである。カウントの変化率は、僅かなカウントの変化や一定レベル以上の変化あるいは心室瘤のような変化率（負の増加率）のように目的に応じて表示レベルが変更される。
【００４７】
請求項６の発明に係る心筋断層画像処理装置は、予め記憶されたカウント評価の統一基準として患者ごとに設定される拡張末期の最大カウントに基づいて評価する統一基準データ設定手段と、このカウントに基づいて任意の心時相に対応する心筋断層データを設定する心時相データ設定手段と、上記各データに基づいて予め決められた分布に従って任意の各心時相の各心筋イメージを作成する心筋イメージ作成手段とを備えるものである。統一基準法による統一基準が任意の心時相に対して用いられ、この統一基準に基づいて任意の心時相における心筋断層像がイメージ表示される。
【００４８】
請求項７の発明に係る心筋断層画像処理装置は、統一基準データ設定手段が拡張末期の最大カウントを基準として評価する拡張末期データ設定手段であって、心筋イメージ作成手段が拡張末期と任意の心時相の各心筋イメージを作成することを特徴とするものである。拡張末期の最大カウントを基準として他の心時相のカウントを計測して画像処理される。
【００４９】
請求項８の発明に係る心筋断層画像処理装置は、拡張末期データ設定手段が、拡張末期における最大カウントを１００％に設定することを特徴とするものである。拡張末期の最大カウントを基準として他の心時相のカウントを計測して画像処理される。例えば健常部や軽度の梗塞部の収縮末期では１００％を越えて計測されることがある。
【００５０】
請求項９の発明に係る心筋断層画像処理装置は、心筋イメージ作成手段が、予め決められた分布の範囲内で任意の範囲を設定することを特徴とするものである。分布としては色分布の他に、濃度分布、輝度分布、図形パターン、カウントや領域毎のグラフによる分布等で表現することができる。なお、本発明においてイメージとはグラフによる表示を含む概念である。
【００５１】
請求項１０の発明に係る心筋断層画像処理装置は、任意の心時相の心筋イメージの分布の範囲に拡張末期の心筋イメージデータの分布の範囲よりもレベルの高い分布を追加することを特徴とするものである。任意の心時相では拡張末期にない色、濃度、輝度等の分布が追加表示され、拡張末期との視覚的な対比が表現される。
【００５２】
請求項１１の発明に係る心筋断層画像処理装置は、心時相データ設定手段が、心筋の心基部から心尖部までの任意の心時相の心筋断層データを設定する場合、任意の一短軸断層面を基準にして収縮または拡張に伴う他の短軸断層面の回転を補正するようにしたことを特徴とするものである。回転の補正は、短軸断層像の中心を基準として収縮や拡張に伴う回転によるずれを補正し、基準とする心時相像と局所を一致させる。
【００５３】
請求項１２の発明に係る心筋断層画像処理装置は、心時相データ設定手段が、任意の断層を設定する場合、複数の時相から時相間における局所の最高値合成を行うようにしたことを特徴とするものである。複数の時相の最高カウントを合成して心時相データが設定されるので、単一の時相との対比によって局所心筋の機能の遅れや進みがあっても表示することできる。
【００５４】
請求項１３の発明に係る心筋断層画像処理装置は、心時相データ設定手段が、任意の断層として短軸断層から極座標または展開図を設定することを特徴とするものである。複数の時相の最高カウントを合成して極座標表示または展開図表示される。
【００５５】
請求項１４の発明に係る心筋断層画像処理装置は、心時相データ設定手段が、収縮末期の極座標または展開図を設定する場合、複数の時相から時相間における局所の最高値合成をするようにしたことを特徴とするものである。局所心筋の収縮の遅れを表現した極座標または展開図で作成される。
【００５６】
請求項１５の発明に係る心筋断層画像処理装置は、予め記憶されたカウント評価の統一基準として患者ごとに設定される拡張末期の最大カウントに基づいて評価する統一基準データ設定手段と、このカウントに基づいて任意の心時相に対応する心筋断層データを設定する心時相データ設定手段と、任意の心時相の心筋断層データとその他の任意の心時相の心筋断層データとの間のカウントの変化率に基づいて局所的な心機能データを設定する心機能データ設定手段と、予め決められた分布に従って任意の心時相とその他の任意の心時相の各心筋イメージを作成するとともに心機能イメージを作成するイメージ作成手段とを備えるようにしたことを特徴とするものである。統一基準でカウント評価したうえでカウントの変化率を加え、全体と細部とがイメージ表示される。
【００５７】
請求項１６の発明に係る心筋断層画像処理装置は、統一基準データ設定手段が拡張末期の最大カウントを基準として評価する拡張末期データ設定手段であって、心機能データ設定手段が任意の心時相の心筋断層データとして拡張末期の心筋断層データを用い、イメージ作成手段が任意の心時相として拡張末期を用いることを特徴とするものである。統一基準として拡張末期における最大カウントに基づいたうえでカウントの変化率を加え、拡張末期を基準として全体と細部とがイメージ表示される。
【００５８】
請求項１７の発明に係る心筋断層画像処理装置は、イメージ作成手段が予め決められた分布の範囲内で任意の範囲に設定されることを特徴とするものである。診断の目的に応じて分布が変更される。
【００５９】
【作用】
請求項１の発明における心筋バイアビリティの評価方法は、カウント評価の統一基準として患者ごとの拡張末期の最大カウントに基づいて任意の心時相の心筋断層データのカウントを評価としたので、どの心時相の心筋断層も統一した基準に基づいて表現され、同じ条件で各心時相ごとの比較が可能となる。この統一基準法は心筋バイアビリティの基本原則として重用することができる。
【００６０】
削除
【００６１】
請求項２の発明における心筋バイアビリティの評価方法は、拡張末期における最大カウントを１００％に設定したので、任意の心時相は拡張末期の最大カウントを基準にして表現され、同じ条件で拡張末期との比較が明瞭となる。
【００６２】
請求項３の発明における心筋バイアビリティの評価方法は、統一基準法と断層像時相配列法とによって、統一基準で各心筋断層像を時相ごとに順次配列したので、立体的なあらゆる方向での心機能の把握が可能となる。
【００６３】
請求項４の発明における心筋バイアビリティの評価方法は、統一基準法と断層像時相配列法とによって、左室全体を網羅し、これらを時相ごとに順次表示したので、立体的なあらゆる方向での心機能の把握が可能となる。
【００６４】
請求項５の発明における心筋バイアビリティの評価方法は、任意の心時相の一短軸断層面を基準とし、回転補正法により心筋の収縮または拡張に伴う他の短軸断層の回転ずれを回転軸を中心として局所が一致するように補正するようにしたので、収縮または拡張に伴う回転を考慮して各心時相像の局所を一致させることが可能となる。
【００６５】
本発明の関連発明１における心筋バイアビリティの評価方法は、最高値合成法により複数の時相から局所の最高値合成を行って断層データを評価するようにしたので、局所心筋の機能の遅れや進みを把握可能となる。
【００６６】
本発明の関連発明２における心筋バイアビリティの評価方法は、最高値合成法により短軸断層から極座標または展開図を作成したので、局所心筋の機能の遅れや進み具合を短軸断層から把握可能となる。
【００６７】
本発明の関連発明３における心筋バイアビリティの評価方法は、最高値合成法により収縮末期における極座標または展開図が作成され、収縮の遅れが把握可能となる。
【００６８】
本発明の関連発明４における心筋バイアビリティの評価方法は、局所心筋の複数時相におけるカウントの変化率を表示するようにしたので、この変化率を患者ごとの各局所心筋の心機能として扱うことができる。
【００６９】
本発明の関連発明５における心筋バイアビリティの評価方法は、拡張末期と任意の心時相におけるカウントの変化率を表示するようにしたので、この変化率を拡張末期と他の心時相との局所心筋の心機能として扱うことができる。
【００７０】
本発明の関連発明６における心筋バイアビリティの評価方法は、各局所心筋における僅かなカウントの変化を明瞭に表現したり、一定のレベル以上の変化を観察することに主眼をおいたり、あるいは心室瘤のような負の増加率を示すものを把握したり、その診断の目的に応じて表示レベルが自由に設定される。
【００７１】
請求項５の発明における心筋断層画像処理装置は、統一基準データ設定手段により予め記憶された統一基準として患者ごとの拡張末期の最大カウントに基づいてカウントを評価し、心時相データ設定手段により上記カウントに基づいて任意の心時相の心筋断層に対応する断層データを設定し、設定されたデータに基づいて心筋イメージ作成手段により予め決められた分布に従って任意の各心時相の心筋イメージをそれぞれ画像処理する。
【００７２】
請求項６の発明における心筋断層画像処理装置は、統一基準データ設定手段として拡張末期データ設定手段が用いられ、この拡張末期データ設定手段により統一基準として拡張末期の最大カウントを用いて評価し、心筋イメージ作成手段により拡張末期と任意の心時相の各心筋イメージを作成する。
【００７３】
請求項７の発明における心筋断層画像処理装置は、拡張末期における最大カウントを１００％にし、この１００％を基準として他の心時相におけるカウントを評価する。
【００７４】
請求項８の発明における心筋断層画像処理装置は、心筋イメージ作成手段により、予め決められた分布の範囲内で任意の範囲を設定して視覚的に正確に表現される。
【００７５】
請求項９の発明における心筋断層画像処理装置は、任意の心時相の心筋イメージにレベルの高い分布が追加され、これによって視覚性が向上する。
【００７６】
請求項１０の発明における心筋断層画像処理装置は、心時相データ設定手段により心筋の心基部から心尖部までの任意の心時相の心筋短軸断層データを設定する場合、収縮または拡張に伴う短軸断面の回転を補正する。
【００７７】
請求項１１の発明における心筋断層画像処理装置は、心時相データ設定手段により複数の時相から局所の最高カウントを合成したので、収縮や拡張に伴う局所心筋の機能の遅れや進みが表現可能となる。
【００７８】
請求項１２の発明における心筋断層画像処理装置は、心時相設定手段により最高値合成を行ったので、極座標または展開図に基づいて収縮や拡張に伴う局所心筋の機能の遅れや進みが表現可能となる。
【００７９】
請求項１３の発明における心筋断層画像処理装置は、心時相データ設定手段により収縮末期における最高値合成を行ったので、収縮末期を特定したうえで極座標または展開図に基づいて局所心筋の収縮の遅れを表現可能となる。
【００８０】
請求項１４の発明に係る心筋断層画像処理装置は、統一基準データ設定手段により患者ごとの統一基準に基づいて評価し、心時相データ設定手段により上記カウントに基づいて任意の心時相に対応する心筋断層データを作成し、心機能データ設定手段により任意の心時相の心筋断層データとその他の任意の心時相の心筋断層データとの間のカウントの変化率に基づいて局所的な心機能データを作成し、イメージ作成手段により上記各データに基づいて任意の心時相とその他の任意の心時相の各心筋イメージを作成するとともに心機能イメージを画像処理する。
【００８１】
請求項１５の発明における心筋断層画像処理装置は、統一基準データ設定手段として拡張末期データ設定手段が用いられ、この拡張末期データ設定手段により統一基準として拡張末期の最大カウントを基準として評価し、心機能データ作成手段が任意の心時相の心筋断層データとして拡張末期の心筋断層データを用い、イメージ作成手段が任意の心時相として拡張末期の心筋イメージと心機能イメージとを用いて画像処理する。
【００８２】
請求項１６の発明における心筋断層画像処理装置は、心機能イメージ作成手段により予め決められた分布の範囲内で任意の分布に変更し、診断を容易に行う。
【００８３】
【実施例】
以下に添付図面を参照して本発明に係る好適な一実施例を詳細に説明する。
実施例１．
実施例１では、投与した放射性医薬品（「ＲＩ」という）の挙動を動態解析するための心筋バイアビリティの評価方法について説明する。
【００８４】
本実施例では、ＲＩとして一例であるが、 99mＴｃ（過テクネシウム酸塩，以下「Ｔｃ」という）−ＭＩＢＩまたは 99mＴｃ−ｓｅｓｔａｍｉｂｉ（以下単に「ミビ」という）が適用される。もちろん、ミビに限定されず、Ｔｃ−Ｔｅｔｒｏｆｏｓｍｉｎ（テトロフォスミン）やＴｃ−Ｔｅｂｏｒｏｘｉｍｅ（テボロキシム）等のＴｃ製剤や、またはＴｃ以外のＲＩを用いてもよく、心筋細胞に吸着あるいは取り込まれ若しくは間質にとどまるものであれば、動態解析用として適用が可能である。それ以外の核種や装置にも応用できる。例えばＰＥＴ（ポジトロンエミッショントモグラフィ）等も用いることが可能である。
【００８５】
このようなＲＩを用いた本実施例の心筋シンチグラフィでは、安静時にＴｃ製剤の適量の静脈投与がなされ、３０分以上経過してから心拍同期法により心拍同期断層像にて画像収集が行われる。なお、本発明では運動負荷や薬物負荷を施行させる心筋シンチグラフィも含んでいるが、以下の説明では安静時を例に挙げて説明する。
【００８６】
本実施例で用いた心拍同期法により、全心周期の心時相像が得られる。比較の基準を同一条件に設定するために、統一基準法によりカウント評価の統一基準を新たに設定した。この統一基準としては例えば拡張末期の最大カウントを用い、この最大カウントを１００％として、この基準を用いて任意の心時相における短軸断層像、長軸面垂直断層像、長軸面水平断層像等の任意の断層像を作成する。いずれの断層像も上記カウントに応じたレベルで表現される。なお、最大カウントを１００％以外の基準に変更して評価することを妨げるものではない。
【００８７】
そこで、レベルを視覚的に且つ客観的に理解できるように、色、濃度、輝度、図形パターン、数値、領域毎のグラク等の分布で表現する手法が使用される。いずれの手法を用いても、カウントの違いが把握しやすいという利点がある。色の分布で表現する方法としては、断層像時相配列法と、短軸断層像から得られる極座標表示方法や展開図表示方法等がある。ここでは、レベルを色の分布で表現した場合の一例を説明する。
【００８８】
図１は本発明に係る心筋バイアビリティの評価方法に用いる極座標表示の分布手法の一例を示す図である。カラー表示部２２の表示画面２２ａには、カウントの違いを色の分布で表現する場合に用いるテーブル（以下カラーバーという）が示されており、Ｃ１は拡張末期像用に使用されているカラーバーを示し、Ｃ２は収縮末期像等のその他の心時相像用に使用されるカラーバーを示している。
【００８９】
まず、拡張末期の全データ中から最大カウントを検索し、この最大カウントを基準とする。健常部の心筋では、拡張末期の最大カウントを基準にして他の心時相像を表現すると、拡張末期の最大カウントを越えてしまう。これに対応できるように、図１に示した如く拡張末期像のカラーバーＣ１よりも高いカウントを表示するようにレベル設定できる。図１に示した例では、カラーバー（Ｃ１，Ｃ２）をｎ％（ｎは正の数）ごとに段階的にレベルが設定されている。カラーバーＣ２は、カラーバーＣ１の最大カウントに対して例えば２段階レベル（１１０％、１２０％以上）が高く設定される。なお、各心時相像の最大カウントを表示できるように例えば１３０％、１４０％というように任意のレベルを追加することができる。カラーバーＣ１，Ｃ２は段階的表示に加えてリニアに表現することも妨げない。重要なことは、拡張末期の最大カウントを統一基準にして任意の心時相像のカウントを色を用いて視覚的に表現できることにある。
【００９０】
図１においては、極座標表示するために、回転補正法により予め収縮や拡張に伴う心筋の回転を補正して基準とする時相像の局所を一致させてある。図１において前壁中隔Ｐと下壁（後壁を含む）Ｑの梗塞例で、側壁Ｒは健常である。健常部Ｒでは拡張末期像３１での最大カウントを呈しており、このカウントを１００％として設定して収縮末期像３２を評価する。心機能イメージ（以下「ファンクショナルイメージ」という）３３の表示レベルは診断目的に応じて独立して自由に変更できる。健常部Ｒのカウントは収縮末期像３２では１００％を越えているが、ファンクショナルイメージ３３ではカウントの変化率はさほど高くなっていない。これは変化率を計算するときの分母が大きく、その結果小さい値をとるので、健常部のファンクショナルイメージは低いレベルを呈した色で表現されている。
【００９１】
また、前壁中隔Ｐと下壁Ｑでは、ともに拡張末期において低カウントを呈している（下壁Ｑは前壁中隔Ｐよりも低カウント）。収縮末期では前壁中隔Ｐのカウントはかなり上昇しているが、下壁Ｑのカウントは若干の上昇にとどまっている。ファンクショナルイメージ３３は、前壁中隔ＰにおいてスケールバーＣ２の最も高いスケール（カウントの変化率が極めて高いこと）を示し、十分な機能が維持されていることが理解される。さらに、下壁Ｑにおいてはカウントの変化率がかなり高いことが理解される。前壁中隔Ｐにバイアビリティがあることは収縮末期像３２とファンクショナルイメージ３３とから容易に診断できるが、下壁Ｑにおいては収縮末期像３２だけではバイアビリティの存在が分かり難い。この場合にファンクショナルイメージ３３を併用すれば、カウントの変化率が高いことが容易に理解でき、下壁Ｑにバイアビリティのあることが理解できる。このように拡張末期像３１で低カウントを呈した心筋局所のバイアビリテイの評価に非常に有用である。
【００９２】
ここで従来からＲＩを用いた心筋シンチグラフィは心拍同期法とは無関係に施行されていたが、心拍同期法を用いて従来法（Ｔｌ心筋シンチグラフィ）を駆使して評価したとしても、拡張末期像と収縮末期像はそれぞれの時相でのそれぞれの最大カウントによる相対的評価でしかなかった。本実施例によれば、拡張末期における最大カウントを統一基準として任意の心時相の評価を行えるものである。この統一基準は、局所心筋の状態や機能、さらには左室全体などの立体的なあらゆる方向への評価・診断に普及可能となり、臨床上必要とされる心筋情報の提供に応えられることとなる。
【００９３】
図２に示すように、展開図表示方法を用いる場合は、カラー表示部２２の表示画面２２ａには拡張末期と収縮末期の各心筋イメージ３１’，３２’及びファンクショナルイメージ３３’が展開されて表示される。この展開図表示では、各心筋イメージ３１’，３２’及びファンクショナルイメージ３３’について、注視したい領域を任意の位置で表示することが可能であるが、開く場所は健常部とするのが通例である。図２においては、健常部である側壁Ｒから開いた展開図表示がなされており、これにカラーバーＣ１，Ｃ２も表示することができ、診断の方法については極座標表示の場合と同様である。
【００９４】
図３に示すように、上記統一基準法に基づいて断層像時相配列法が用いられる。断層像時相配列法は、カラー表示部２２の表示画面２２ａに示すように、任意の断層で心拍同期法に基づいて左室全体を網羅し、これらを順次時相毎に配列することができる。各心筋断層像は、時相ごとに順次変化を示すので、局所の心筋血流のみならず、あらゆる方向での心機能（壁厚、回転、短縮、拡張、ひきつれ等）の把握が可能となる。これは、各心筋断層面での局所心筋のカウントの変化、カウントを表示する領域の面積や厚さや長さの変化を観察する等の方法により、局所心筋の状態や機能さらには左室全体の総合的な状態や機能までも把握可能となる。また、任意の時相における心筋断層像のスライス枚数から各心筋断層面に垂直方向での短縮または拡張も理解できることとなる。この断層像時相配列法を極座標表示法等と併用することによって、左室全体の総合的な状態や機能までも把握可能となる。断層像時相配列法としては短軸断層を前提として図示説明したが、短軸断層以外の任意の断層を時相配列することも容易である。なお、断層像時相配列法により注視したい局所心筋のみの同一部位を時相配列することも可能である。
【００９５】
図４に示すように、短軸断層面については収縮または拡張に伴う回転を考慮し、回転補正法により任意の一時相を基準に回転を補正する。任意の一時相として拡張末期Ａを基準とすることができる。拡張末期Ａに対して補正をする心時相Ｂの収縮または拡張に伴う回転ずれを検出し、この回転ずれに応じて心時相Ｂをその回転軸Ｏを中心として拡張末期Ａの回転度に補正し、補正された心時相Ｂ’が得られる。回転ずれは心時相Ｂの時相ごとに相違するが、拡張末期Ａに局所が一致するように補正される。このようにして得られた短軸断層データから極座標表示法や展開図表示法を行い、局所心筋の対比が明瞭となる。なお、回転補正を行わないときには、局所心筋の回転の様子を理解できるので、併用して評価を行うことが望ましい。
【００９６】
次に図５に示すように、最高値合成法により最高値合成を行うことができる。任意の心時相として収縮に係わる遅れを例示すると、健常部に比較して収縮が遅れる梗塞部がある。このために収縮末期像を単一の時相に限定すると、実際の収縮を見落としてしまうおそれが生じる。収縮末期と設定した時相よりも後の時相で最高カウントを呈する局所心筋が存在することがあるためである。図５に示すように健常部の収縮末期近傍から複数の心時相３２ａ，３２ｂを探しだす。先の心時相３２ａは後壁ｑｂよりも下壁ｑが低いカウントを示しているが、後の心時相３２ｂでは後壁ｑｂのカウントが減少しているにも拘らずに下壁ｑのカウントが増加しているので、下壁ｑでは収縮が遅れていることが把握される。そこで、複数の心時相３２ａ，３２ｂにおいて、下壁ｑの最高値を後の心時相３２ｂから用い、後壁ｑｂの最高値を先の心時相３２ａから用い、この両者の局所を合わせることで収縮末期３２の下壁（後壁を含む）Ｑの最高カウントが合成される。前壁中隔Ｐと側壁Ｒにおいても、同様の手法で最高値合成がなされる。このように収縮末期像の時相に幅をもたせ、複数の時相から各局所のカウントの最高値を合成して任意の断層としての短軸断層から極座標または展開図を作成することが望ましい。こうして得られた心筋断層像と単一の時相による心筋断層像とを対比すれば、収縮の遅れを呈する局所がどこであるかを理解できることとなる。なお、最高値合成法は、短軸断層以外の任意の断層または収縮末期以外の任意の心時相を対象とすることが可能である。
【００９７】
図６は図１，図２，図５と異なる例であるが、図６に示すように、各局所心筋の拡張末期像から任意の心時相像へのカウントの変化率をグラフで表現することができる。そこで変化率については、〔任意の局所心筋における｛（任意の心時相の心筋断層データのカウント−拡張末期の心筋断層データのカウント）／拡張末期の心筋断層データのカウント｝」×１００％で求められる。カウントの変化率を時系列的に表示することで、局所心筋の機能を推定できる。この場合に、僅かなカウントの変化を明瞭に表現したり、一定のレベル以上の変化を観察することに主眼をおいたり、心室瘤のような負の増加率を示すものを把握したり、診断の目的に応じて表示レベルを自由に設定できる。図１に示したカラーバーのように色分布を選択してレベルの範囲を設定することもでき、グラフ曲線表示以外の分布表示も可能である。
【００９８】
健常部８１においては全体的に高いカウントが表示され、急激な立ち上がりで高い振幅を示している。このようなカウントレベル、立ち上がりトレンド、立ち上がり時間、振幅の高さ、時間による微分解析等の各種要因から健常部８１をパターン化することができ、その評価が可能となる。次に、梗塞部８２，８３においてはカウントが高くならず、特に梗塞部８３では最高カウントをとる時相が遅れることが表示されている。さらに心室瘤８４においては、健常部の収縮による受動的延長という物理的ストレスに晒されるので、収縮期８４’のカウントは拡張期のカウントより減少する。また、負の増加率により心室瘤を呈する局所を特定することが容易になる。本発明者はこのようにカウントの変化率から心機能を把握することをカウント変化率法と命名した。なお、カウントの変化率は最高値合成後に行うことができる。
【００９９】
図７は、図６を時間で微分したグラフである。健常部８１は鋭い立ち上がりと高い頂点を示している。これに対し梗塞部８２では健常部８１よりは劣っている。また梗塞部８２より大きなダメージを受けた梗塞部８３は遅延収縮を呈している。梗塞部８３は、頂点の高さが梗塞部８２と大きな差がなく残存心筋量は比較的保たれているが、頂点に達する時間は大幅に遅れており、機能の低下が把握できる。心室瘤８４ではグラフの正と負とが逆になっており、その同定は容易である。このように頂点の高さとそれに達するまでの時間を健常部８１と比較することで残存心筋量と残存心機能を測定することが可能となる。本発明者はこのようなカウントの変化率を時間で微分して心機能を把握することをカウント変化率時間微分法と命名した。カウント変化率時間微分法によれば、心筋梗塞に限らず、心筋症や感染・免疫・代謝等による変性疾患やその他の疾患の局所心筋の状態や機能の把握が可能であり、また疾患の診断，治療の効果判定，予防等についての詳細な情報を提供できる。
【０１００】
以上説明したように実施例１によれば、統一基準法、断層像時相配列法、回転補正法、最高値合成法、カウント変化率法、カウント変化率時間微分法を単独または併用して駆使し、心筋バイアビリティおよび心機能を正確に評価することができる。
【０１０１】
実施例２．
次に上述した実施例１を具体的に実現する心筋断層画像処理装置について実施例２を用いて説明する。なお、以下の説明では、短軸断層像による表示を例に挙げて説明する。
【０１０２】
図８は本発明の一実施例による心筋断層画像処理装置の概略構成を示すブロック図である。心筋断層画像処理装置１は、ガンマカメラによって得られた断層データを入力する断層データ入力部１０と、画像処理にかかるキー操作を行うための操作部１２と、任意の断層の断層データを記憶する断層データ記憶部１４と、断層データから放射線（Ｔｃ）の分布を求める際にＩＣを用いて画像再生を行うとともに断層データの断層像時相配列処理，極座標処理，展開図処理等の画像処理を専用に行う画像処理部１６と、画像処理された表示用データを記憶する表示用記憶部１８と、画像処理過程をカラー表示するとともに任意の心時相の各心筋イメージや拡張末期像から任意の心時相像へのカウントの増加率を画像化した心機能イメージをカラー表示するカラー表示部２２と、一連の画像処理過程をカラー印刷するとともに前述の各心筋イメージや心機能イメージをカラー印刷するカラー印刷部２４と、本装置全体を制御する制御部２６と、各種のデータ、アドレス信号、制御信号を伝送するバス２８とを備えている。
【０１０３】
上記制御部２６は、プログラムによって制御動作を実行するＣＰＵ２６ａと、ＣＰＵ２６ａを動作させるための制御プログラムを格納したＲＯＭ２６ｂと、各種制御プログラムのワークエリアとして用いるＲＡＭ２６ｃとを備えている。上記ＲＯＭ２６ｂには、図９〜図１４に示す制御プログラム（最大カウント評価プログラム、心筋断層像構築プログラム、極座標や展開図等の表示処理プログラム、回転ずれ補正処理プログラム、最高値合成処理プログラム、心機能処理プログラム等を含む）及び図１〜図７に示す表示画面を形成するための表示プログラム等の制御プログラムが格納されている。なお、各図のフローに従う個々の動作は画像処理部１６等の各部で行われる。
【０１０４】
図９〜図１０に基づいて基本的動作を説明する。画像処理手順は、任意の心時相における任意の断層の断層データとして断層データ入力部１０から入力され、断層データ記憶部１４に格納される（ステップ１００）。任意の断層として短軸断層を選択し、任意の心時相として拡張末期と収縮末期を選択する場合を例示して説明する（ステップ１０１）。拡張末期においては、拡張末期の全データ中から最大カウントが検索されて登録され、この最大カウントを１００％として以後の処理過程において他のカウントが評価される（ステップ１０２）。この最大カウントの検索は本実施例の根幹をなすので、任意のステップに組み込むことが可能である。
【０１０５】
次に短軸断層の中で回転ずれの補正を必要とするか否かが選択される（ステップ１０３）。操作部１２の操作によって回転ずれの補正が指示されると、回転ずれ補正が行われ、画像処理部１６に登録される（ステップ１０４）。回転補正を必要とするか否かが選択されて処理された後に、短軸断層面を時相配列するか否かを選択する（ステップ１０５）。時相配列が指示されると、左室全体を網羅するように時相配列処理が行われ、画像処理部１６に登録される（１０６）。時相配列をするか否かが選択されて処理された後に、最高値合成を必要とするか否かが選択される（ステップ１０７）。最高値合成が指示されると、複数の時相間の最高値が合成されて画像処理部１６に登録される。任意の心時相として収縮末期を選択している場合には、収縮末期における最高値が合成される（ステップ１０８）。
【０１０６】
次に拡張末期と任意の心時相としての収縮末期のイメージ処理が行われ（ステップ１０９）、表示用記憶部１８に登録される。ファンクショナルイメージの処理を行う場合には、同様に表示用記憶部１８に登録される（ステップ１１０）。なお、ファンクショナルイメージの処理は後述する心機能表示を選択した後に行うステップに組みかえしても良い。
【０１０７】
表示方法としてイメージ表示する各種方法が選択されるが、ここでは極座標と展開図のみを図示して説明する（ステップ１１１）。極座標表示を選択した後、心機能も表示するか否かが選択される（ステップ１１２）。心機能も表示する場合には、表示用記憶部１８から各イメージが読み出され、図１に示すように拡張末期像３１と任意の心時相像としての収縮末期像３２とファンクショナルイメージ３３とが極座標によってカラー表示される。心機能を表示しない場合には、拡張末期像３１と任意の心時相像としての収縮末期像３２とがカラー表示される（ステップ１１４）。
【０１０８】
ステップ１１１にて展開図表示が選択された場合には上述した極座標表示とは表示方法が異なるだけである。表示用記憶部１８から各イメージが読み出され、図２に示すように拡張末期像３１’、収縮末期像３２’、ファンクショナルイメージ３３’の各展開図がカラー表示される（ステップ１１５〜１１７）。
【０１０９】
さらにステップ１０１にて短軸断層以外の他の断層データが選択された場合には、上述したステップ１０２と同様に最大カウント検索がなされる（ステップ１１８）。長軸面垂直断層や長軸面水平断層等の断層データは、ステップ１１９〜１２０に示すように断層像時相配列の処理がなされ、ステップ１０９〜１１０のイメージ処理がなされ、極座標や展開図以外の他の表示方法が選択されてカラー表示される。
【０１１０】
図１１及び図４に示す短軸断層面の回転ずれ補正処理について説明する。各短軸断層面による心時相像のずれを補正する場合、図４に示すように拡張末期の短軸断層面を基準にして他の全心周期の短軸断層面の回転を補正しようとするもので、極座標または展開図処理の前に行われる。勿論、回転ずれ補正を選択しない場合もある（ステップ１０３参照）。
【０１１１】
操作部１２の操作によって、心基部から心尖部までの短軸断層の中で回転ずれの補正を必要とする時相が選択される（ステップ２０１）。短軸断層の各時相ごとに補正に必要な回転度が設定されて、各時相の番号がＲＡＭ２６ｃに登録される（ステップ２０２〜２０４）。
【０１１２】
そして、操作部１２の操作によって回転ずれの補正が指示されると（ステップ２０５）、登録された各時相の番号に従い、各時相に対応した回転度に基づく回転ずれ補正が実行され、同時に選択されなかった時相にも登録された時相の回転度に応じてリニアに回転ずれ補正が実行される（ステップ２０６）。図４に示すように任意の心時相像Ｂは、拡張末期Ａの位置に補正されて心時相像Ｂ’として再記憶される（ステップ２０７）。なお、登録された時相のみ回転度を指定したが、必要なすべての時相に回転を指定するようにしても良い。
【０１１３】
図１２及び図５に示す最高値合成処理について説明する。前述のごとく梗塞部における収縮末期では心筋の収縮が遅れることがあるために、複数の時相から時相間における局所の最高値合成をすることが必要となる。収縮末期像については、時相を、例えば拡張末期像（Ｒ波からの最初の画像）から２００ｍｓｅｃ〜３００ｍｓｅｃとした場合、梗塞部心筋の一部では次の時相（３００ｍｓｅｃ〜４００ｍｓｅｃ）で最高カウントをとる場合がある（図６の梗塞部８３参照）。このようなときには、一部の心筋バイアビリティの評価を誤診してしまう危険性がある。この危険性を取り除くために、収縮末期像の時相を任意の期間（複数の時相間）に設定して、各局所の最高値合成を可能にしたものである。時間分解能の精度が向上するに伴って最高値合成は不可欠なものとなろう。
【０１１４】
図１２に基づいて最高値合成の設定手順を説明する。最高値合成は前述の回転補正を選択するか否かの処理後に行われる。操作部１２の操作によって、任意の複数の時相、例えば３２ａ，３２ｂが選択される。例えば収縮末期像の時相を、前述のように２００ｍｓｅｃ〜４００ｍｓｅｃに設定することができる（ステップ３０１）。もちろん、２００〜５００ｍｓｅｃや３００〜７００ｍｓｅｃ等のように任意に設定可能である。選択された時相に従って心基部から心尖部からまでの各短軸断層の局所ごとに最高カウントが検索される（ステップ３０２，３０３）。このようにして検出された短軸断層データは収縮末期像３２として登録される。
【０１１５】
図１３は心機能の画像処理を示すフローチャートである。まず断層データ記憶部１４より拡張末期と任意の心時相としての収縮末期の各短軸断層データから評価が必要となる局所心筋のカウントが検出される（ステップ４０１〜４０２）。そしてカウントの増加率（ファンクショナルスコア「Ｆ．Ｓ．」という）を下式により算出する（ステップ４０３）。すなわち、
Ｆ．Ｓ．＝｛（ＥＳ−ＥＤ）／ＥＤ｝×１００（％）
ＥＤ：拡張末期のカウント、
ＥＳ：収縮末期のカウント
である。
【０１１６】
以上のステップ４０１〜４０３の処理は短軸断層の心基部から心尖部までの断層面ごとに実施され（ステップ４０４）、全領域の処理が終了すると、カウントとカラーバーとによる極座標処理が選択される（ステップ４０５）。カラー表示部２２では、図１に示したように、その表示画面２２ａには拡張末期と収縮末期の各心筋イメージ３１，３２及びファンクショナルイメージ３３が表示されると同時に、カラーバーＣ１，Ｃ２による表示がなされる。
【０１１７】
次に図１４に示す色分布処理について説明する。拡張末期のデータ中のＴｃの最大カウントを検索し、この最大カウントを基準にして図１〜図３に示すごときカラーバーＣ１の色分布の範囲が設定される。カラーバーＣ１は上記最大カウントを１００％として均等あるいは任意に分布して表示され、レベルの高い分布を含むカラーバーＣ２を追加することができる。拡張末期ではカラーバーＣ１を用い、収縮末期ではカラーバーＣ２が使用される。なお、心機能については、拡張末期との比較を行うものではなく変化率を視覚化するので、異なるレベルに設定することができる（ステップ５０１）。上記最大カウントを基準として任意の心時相のカウントが評価され、拡張末期と任意の心時相の各心筋イメージの色分布が作成される（ステップ５０２〜５０３）。また、前述の計算式でファンクショナルイメージの色分布が作成される（ステップ５０３）。
【０１１８】
以上説明したように実施例２によれば、統一基準法、断層像時相配列法、回転補正法、最高値合成法、カウント変化率法、カウント変化率時間微分法を具現化して処理操作することができ、心筋バイアビリティおよび心機能を正確に評価する心筋断層画像処理装置を提供することができる。
【０１１９】
【発明の効果】
以上説明したように請求項１の発明に係る評価方法によれば、統一基準法により患者ごとに設定した拡張末期の最大カウントを基準としたので、どの心時相の心筋断層も拡張末期の最大カウントを基準にして表現され、拡張末期との比較が明瞭となり、どの心時相の心筋断層も同一の基準から表現され、心時相間における比較が明瞭となる。特に、心拍同期法を用いた心筋シンチグラフィを拡張末期の最大カウントから表現することができ、拡張末期を所望の心時相と対比することによって心筋バイアビリティを正確に評価することができる。
【０１２０】
削除
【０１２１】
請求項２の発明に係る評価方法によれば、拡張末期の最大カウントを１００％に設定したので、他の心時相との比較が明瞭となる。殊に、拡張末期と他の心時相を統一基準法によって評価する際、他の心時相で１００％以上にカウントを越える領域も表示することができ、視覚的表現を行った場合には一層明瞭に評価することができる。
【０１２２】
請求項３の発明に係る評価方法によれば、統一基準法と断層像時相配列法によって任意の断層の各心筋断層像を時相ごとに順次配列したので、立体的なあらゆる方向での心機能を評価することができる。
【０１２３】
請求項４の発明に係る評価方法によれば、断層像時相配列法によって左室全体を網羅して時相ごとに順次配列したので、局所の心筋血流のみならず、立体的なあらゆる方向での心機能を評価することができる。
【０１２４】
請求項５の発明に係る評価方法によれば、回転補正法により心筋の収縮または拡張に伴う短軸断面の回転を補正するようにしたので、任意の一心時相像を基準として他の心時相像の局所を一致させることができる。心基部から心尖部までの各局所の回転によるずれが補正されて位置を合わせることによって、各局所を直接対比して評価することができる。
【０１２５】
本発明の関連発明１に係る評価方法によれば、最高値合成法により時相間における局所の最高カウントを合成して断層データを得るので、収縮や拡張に伴う局所心筋の機能の遅れや進みを的確に評価することができる。
【０１２６】
本発明の関連発明２に係る評価方法によれば、最高値合成法に基づいて短軸断層データを得るので、極座標または展開図を作成して収縮や拡張に伴う局所心筋の機能の遅れや進みを的確に評価することができ、心筋バイアビリティの見落としを確実に防止することができる。
【０１２７】
本発明の関連発明３に係る評価方法によれば、最高値合成法により局所心筋に収縮の遅れがあっても、局所の最高値を合成することによって実際の収縮を見落さずに的確に把握することができるうえに、収縮末期を新たに定義することができる。これによって、梗塞部心筋の遅延収縮による心筋バイアビリティの見落としを確実に防止することができる。
【０１２８】
本発明の関連発明４に係る評価方法によれば、カウント変化率法により患者ごとの局所心筋の拡張末期から任意の心時相へのカウントの変化率を表示するようにしたので、この変化率を各局所心筋の心機能として扱うことができる。
【０１２９】
本発明の関連発明５は係る評価方法によれば、カウント変化率法により局所心筋の拡張末期から任意の心時相へのカウントの変化率を表示するようにしたので、この変化率を各局所心筋の心機能として扱うことができる。
【０１３０】
本発明の関連発明６に係る評価方法によれば、カウント変化率法によりカウントの増加を明瞭に表現したり、一定レベル以上の変化を観察することに主眼をおいたり、心室瘤のような負の増加を示すものも把握することができ、診断目的に応じて表示レベルを自由に設定することができる。
【０１３１】
請求項６の発明に係る心筋断層画像処理装置によれば、拡張末期の最大カウウントを患者ごとの統一基準として評価して各心時相を視覚的に表示したので、心筋バイアビリティの評価をするための資料を総合的且つ客観的見地から明確にすることができ、心筋バイアビリティを残存心筋量と残存心筋の収縮能から評価することが可能となる。
【０１３２】
請求項７の発明に係る心筋断層画像処理装置によれば、拡張末期と任意の心時相とを統一基準から視覚的に表示することができ、拡張末期を基準として他の任意の心時相との対比を行って心筋バイアビリティの評価をすることができる。
【０１３３】
請求項８の発明に係る心筋断層画像処理装置によれば、拡張末期における最大カウントを１００％にしてこれを基準に他の任意の時相との対比を容易に行うことができる。
【０１３４】
請求項９の発明に係る心筋断層画像処理装置によれば、心筋イメージが色分布やグラフ表示等で任意に設定されているので、心筋の各状態を視覚的に正確に表現することができる。
【０１３５】
請求項１０の発明に係る心筋断層画像処理装置によれば、任意の心時相の心筋イメージに高い分布を追加して表現することができるので、拡張末期との対比を視覚的に明瞭にすることができる。
【０１３６】
請求項１１の発明に係る心筋断層画像処理装置によれば、心筋の収縮や拡張に伴う回転に対処して任意の心時相像の局所を一致させることができ、極座標表示や展開図表示を行ったときには同一表示位置での局所対比を行うことができる。
【０１３７】
請求項１２の発明に係る心筋断層画像処理装置によれば、単一時相に固定せずに複数時相から任意の断層における局所ごとの最高値合成をしたので、局所心筋の機能の遅れや進みを把握することができ、心筋バイアビリティを正確に表現することができる。
【０１３８】
請求項１３の発明に係る心筋断層画像処理装置によれば、短軸断層における時相間の最高値合成をしたので、極座標や展開図によって心筋バイアビリティを正確に表現することができる。
【０１３９】
請求項１４の発明に係る心筋断層画像処理装置によれば、収縮末期における局所心筋の収縮の遅れを見落としすることなく収縮末期像として表現することができる。
【０１４０】
請求項１５の発明に係る心筋断層画像処理装置によれば、拡張末期の最大カウウントを患者ごとの統一基準として任意の心時相像を新たに表示したうえにこれらのカウントの変化率から局所の心機能イメージを併せて表示するので、全体的且つ局所的な心筋バイアビリティを表現することができる。
【０１４１】
請求項１６の発明に係る心筋断層画像処理装置によれば、拡張末期の最大カウントから任意の心時相像を新たに表示したうえにこれらのカウントの変化率から局所の心機能イメージを併せて表示することができる。
【０１４２】
請求項１７の発明に係る心筋断層画像処理装置によれば、診断目的に応じて心筋と心機能のイメージの分布を独立して変更することができる。
【図面の簡単な説明】
【図１】  本発明に係る心筋バイアビリティの評価方法を極座標表示した一例を示す説明図である。
【図２】  本発明に係る心筋バイアビリティの評価方法を展開図表示した一例を示す説明図である。
【図３】  本発明に係る心筋バイアビリティの評価方法を断層像時相配列法に基づいて表示した一例を示す説明図である。
【図４】  本発明に係る心筋バイアビリティの評価方法の回転ずれ補正の一例を示す説明図である。
【図５】  本発明の関連発明に係る心筋バイアビリティの評価方法の最高値合成の一例を示す説明図である。
【図６】  本発明の関連発明に係る心筋バイアビリティの評価方法のカウント変化率の一例を示すグラフである。
【図７】  本発明の関連発明に係る心筋バイアビリティの評価方法のカウント変化率時間微分法の一例を示すグラフである。
【図８】  本発明に係る心筋断層画像処理装置の概略構成を示すブロック図である。
【図９】  本発明に係る心筋断層画像処理装置による画像処理手順を説明するフローチャートである。
【図１０】  本発明に係る心筋断層画像処理装置による画像処理手順を図８に続いて説明するフローチャートである。
【図１１】  本発明に係る心筋断層画像処理装置による回転ずれ補正の手順を説明するフローチャートである。
【図１２】  本発明に係る心筋断層画像処理装置による最高カウントの合成手順を説明するフローチャートである。
【図１３】  本発明に係る心筋断層画像処理装置による心機能の画像処理手順を説明するフローチャートである。
【図１４】  本発明に係る心筋断層画像処理装置による色分布処理手順を説明するフローチャートである。
【符号の説明】
１０            断層データ入力部
１２            操作部
１４            断層データ記憶部
１６            画像処理部
１８            表示用記憶部
２２            カラー表示部
２４            カラー印刷部
２６            制御部
２６ａ          ＣＰＵ
２６ｂ          ＲＯＭ
２６ｃ          ＲＡＭ
２８            バス
３１，３１’    拡張末期像
３２，３２’    収縮末期像
３３，３３’    ファンクショナルイメージ（心機能像）
Ｃ１，Ｃ２      カラーバー
【図１】
0003724511.tif 000002
【図２】
0003724511.tif 000003
【図３】
0003724511.tif 000004
【図４】
0003724511.tif 000005
【図５】
0003724511.tif 000006
【図６】
0003724511.tif 000007
【図７】
0003724511.tif 000008
【図８】
0003724511.tif 000009
【図９】
0003724511.tif 000010
【図１０】
0003724511.tif 000011
【図１１】
0003724511.tif 000012
【図１２】
0003724511.tif 000013
【図１３】
0003724511.tif 000014
【図１４】
0003724511.tif 000015