特開 2023-148281 心筋血流の循環指標算出方法、プログラムおよび情報処理装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023148281

(43)【公開日】2023-10-13

(54)【発明の名称】心筋血流の循環指標算出方法、プログラムおよび情報処理装置

(51)【国際特許分類】

   G01T 1/161 20060101AFI20231005BHJP

【ＦＩ】

G01T1/161 B

G01T1/161 D

【審査請求】未請求

【請求項の数】8

【出願形態】ＯＬ

(21)【出願番号】P 2022056215

(22)【出願日】2022-03-30

(71)【出願人】

【識別番号】000230250

【氏名又は名称】日本メジフィジックス株式会社

(74)【代理人】

【識別番号】100113549

【弁理士】

【氏名又は名称】鈴木 守

(74)【代理人】

【識別番号】100115808

【弁理士】

【氏名又は名称】加藤 真司

(72)【発明者】

【氏名】為重 喜行

(72)【発明者】

【氏名】宮村（布村） 梢

(72)【発明者】

【氏名】白井 求

(72)【発明者】

【氏名】久保田 雅博

【テーマコード（参考）】

4C188

【Ｆターム（参考）】

4C188EE02

4C188FF04

4C188KK24

4C188KK33

4C188MM04

(57)【要約】      （修正有）

【課題】簡便に心筋血流の循環指標を求める技術を提供する。
【解決手段】負荷状態にある被験者にＴｌ製剤を投与してダイナミック収集によるＳＰＥＣＴ撮像を行うステップと、負荷状態にある被験者をＳＰＥＣＴ撮像して負荷後像を取得するステップと、Ｔｌ製剤の再分布時に安静状態にある被験者をＳＰＥＣＴ撮像して再分布像を取得するステップと、ダイナミック収集によって得られたダイナミック画像に基づいて時間放射能曲線を求め、ｓＭＢＦを求めるステップと、負荷後像および再分布像に基づいてＴｌのＷＯＲを求めるステップと、次の関係式（１）に基づき、ｓＭＢＦとＷＯＲからｒＭＢＦを求めるステップと、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めたｒＭＢＦとｓＭＢＦとに基づいてＣＦＲを求めるステップとを備える。
【選択図】図５

【特許請求の範囲】

【請求項1】

被験者の心筋血流の循環指標を算出する方法であって、
負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するステップと、
前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるステップと、
前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるステップと、
次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるステップと、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求めるステップと、
を備える、循環指標算出方法。

【請求項2】

負荷時心筋血流量（ｓＭＢＦ）を求める前記ステップにおいて、ダイナミック収集によって得られた複数の時間点における前記負荷時の心筋ＳＰＥＣＴ画像データに基づいて時間放射能曲線を求め、前記時間放射能曲線に基づいて負荷時心筋血流量（ｓＭＢＦ）を求める、請求項１に記載の循環指標算出方法。

【請求項3】

負荷時心筋血流量（ｓＭＢＦ）を求める前記ステップにおいて、前記負荷時の心筋ＳＰＥＣＴ画像データと前記Ｔｌ製剤の投与量と被験者の体重とに基づいてＳＵＶ（Standard Uptake Value）を求め、予め取得されたＳＵＶと負荷時心筋血流量（ｓＭＢＦ）との相関関係に基づいて、負荷時心筋血流量（ｓＭＢＦ）を求める、請求項１に記載の循環指標算出方法。

【請求項4】

前記負荷時心筋血流量（ｓＭＢＦ）と、前記冠動脈血流予備能（ＣＦＲ）とに基づいて、ＣＦＣ（Coronary Flow Capacity）を算出するステップをさらに備える、請求項１乃至３のいずれか１項に記載の循環指標算出方法。

【請求項5】

前記循環指標の少なくとも１つをポーラーマップで表示するステップをさらに備える請求項１乃至４のいずれか１項に記載の循環指標算出方法。

【請求項6】

前記ポーラーマップで表示するステップにおいて、負荷時心筋血流量（ｓＭＢＦ）および洗い出し率（ＷＯＲ）の少なくとも１つまたは両方をポーラーマップで表示し、ポーラーマップで表示された前記循環指標と並べて表示することを含む、請求項５に記載の循環指標算出方法。

【請求項7】

被験者の心筋血流の循環指標を算出するためのプログラムであって、コンピュータに、
負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するステップと、
前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるステップと、
前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるステップと、
次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるステップと、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求めるステップと、
を実行させるプログラム。

【請求項8】

被験者の心筋血流の循環指標を算出するための情報処理装置であって、
負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するデータ取得部と、
前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるｓＭＢＦ算出部と、
前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるＷＯＲ算出部と、
次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるｒＭＢＦ算出部と、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求める循環指標算出部と、
求めた循環指標を出力する出力部と、
を備える情報処理装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、心筋血流の循環指標を算出する技術に関する。

【背景技術】

【0002】

心筋血流検査に用いられるＴｌ製剤（塩化タリウム（^２０１Ｔｌ））は、虚血の検出や心筋生存能（心筋Viability）の評価に優れている。また、Ｔｌ製剤は、再分布現象があるため、１回の投与で負荷時と再分布時の撮像が可能である。運動負荷Ｔｌ検査で利用される心筋からの洗い出し率（ＷＯＲ：Washout rate）は、次の式で計算される。
｛（負荷後像カウント）－（再分布像カウント）｝／（負荷後像カウント）×１００％

【0003】

なお、再分布像は負荷後像より３～４時間後に撮像されるので、その時間分だけカウントが減衰する。この減衰分を補正するために、上式において、再分布像カウントに時間補正係数を乗じることもある。

【0004】

ＷＯＲは、３枝病変の診断に有用である。正常部位ではＷＯＲは４０～５０％であるが、心筋虚血部位ではＷＯＲは低下し、４０％以下を異常値とすることが多い。

【0005】

また、心筋血流検査に用いられる指標として、冠動脈血流予備能（ＣＦＲ：Coronary Flow Reserve）が知られている。ＣＦＲは、安静時の心筋血流量（ｒＭＢＦ：rest Myocardial Blood Flow）と負荷時の心筋血流量（ｓＭＢＦ：stress Myocardial Blood Flow）の比である。ＣＦＲの低下は、冠微小循環障害や機能的虚血を反映するだけでなく、冠動脈の独立した予後予測因子である。さらに、心筋血流検査に用いられる指標として、Coronary Flow Capacity（ＣＦＣ）が注目されている。ＣＦＣは、ｓＭＢＦおよびＣＦＲとの組み合わせにより算出される指標であり、ｒＭＢＦ、ｓＭＢＦ又はＣＦＲ単独で心筋血流の循環指標を評価するよりも優れた評価が可能となる。

【先行技術文献】

【非特許文献】

【0006】

【非特許文献1】Shinya Shiraishi他「Prediction of Left Main or 3-Vessel Disease Using Myocardial Perfusion Reserve on Dynamic Thallium-201 Single-Photon Emission Computed Tomography With a Semiconductor Gamma Camera」 Circulation Journal 2015; 79:623-631

【非特許文献2】Robert M Bober他「The impact of revascularization on myocardial blood flow as assessed by positron emission tomography」 Eur J Nucl Med Mol Imaging 2019;46(6):1226-1239

【発明の概要】

【発明が解決しようとする課題】

【0007】

上述したとおり、ＣＦＲは、ｓＭＢＦとｒＭＢＦの比であるので、ＣＦＲを算出するためには、負荷時と安静時の両方において、Ｔｌ，Ｔｃ等の心筋ＳＰＥＣＴ製剤を被験者に投与する必要があり、患者の負担が大きいと共に算出手順も煩雑であった（非特許文献１）。また、ＣＦＣは心筋血流の循環指標の包括的な評価指標であり、ＣＦＲとｓＭＢＦの組み合わせから算出されるため、ＣＦＣの算出手順も煩雑であった（非特許文献２）。

【0008】

そこで、本発明は、上記背景に鑑み、簡便にＣＦＲを始めとする心筋血流の循環指標を定量する技術を提供することを目的とする。

【課題を解決するための手段】

【0009】

本発明の心筋血流の循環指標算出方法は、被験者の心筋血流の循環指標を算出する方法であって、負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するステップと、前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるステップと、前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるステップと、次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるステップと、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求めるステップとを備える。

【0010】

本発明のプログラムは、被験者の心筋血流の循環指標を算出するためのプログラムであって、コンピュータに、負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するステップと、前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるステップと、前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるステップと、次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるステップと、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求めるステップと、を実行させる。

【0011】

本発明の情報処理装置は、被験者の心筋血流の循環指標を算出するための情報処理装置であって、負荷状態の被験者に放射性タリウム製剤（Ｔｌ製剤）を投与した後早期に収集された負荷時の心筋ＳＰＥＣＴ画像データ、および、前記Ｔｌ製剤の再分布時に収集された安静時の心筋ＳＰＥＣＴ画像データを取得するデータ取得部と、前記負荷時の心筋ＳＰＥＣＴ画像データから負荷時心筋血流量（ｓＭＢＦ）を求めるｓＭＢＦ算出部と、前記負荷時の心筋ＳＰＥＣＴ画像データおよび前記安静時の心筋ＳＰＥＣＴ画像データに基づいて放射能の洗い出し率（ＷＯＲ）を求めるＷＯＲ算出部と、次の関係式（１）に基づき、負荷時心筋血流量（ｓＭＢＦ）と洗い出し率（ＷＯＲ）から安静時心筋血流量（ｒＭＢＦ）を求めるｒＭＢＦ算出部と、
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
求めた安静時心筋血流量（ｒＭＢＦ）と負荷時心筋血流量（ｓＭＢＦ）とに基づいて、冠動脈血流予備能（ＣＦＲ）を含む前記循環指標を求める循環指標算出部と、求めた循環指標を出力する出力部とを備える。

【発明の効果】

【0012】

本発明によれば、Ｔｌ製剤を１回投与して得られた心筋ＳＰＥＣＴ画像に基づいて、簡便にＣＦＲを始めとする心筋血流の循環指標を求めることができる。

【図面の簡単な説明】

【0013】

【図1】実施の形態の情報処理装置の構成を示す図である。

【図2】検査プロトコルおよび情報処理装置に入力されるＳＰＥＣＴデータを説明するための図である。

【図3】ＣＦＣのステージを求める方法の一例を示す図である。

【図4】（ａ）ｓＭＢＦとＣＦＲのデータをポーラーマップ表示した例を示す図である。（ｂ）ｓＭＢＦとＣＦＲのデータを３Ｄ表示した例を示す図である。

【図5】実施の形態の循環指標算出方法の例を示すフローチャートである。

【図6】実施の形態の循環指標算出方法の別の例を示すフローチャートである。

【図7】Ｔｌ製剤を２回投与する従来の方式について説明するための図である。

【発明を実施するための形態】

【0014】

以下、本実施の形態に係る心筋血流の循環指標の算出方法および情報処理装置について図面を参照して説明する。なお、以下の説明はあくまでも好ましい態様の一例を示したものであり、特許請求の範囲に記載された発明を限定する意図ではない。

【0015】

図１は、実施の形態の情報処理装置１の構成を示す図である。情報処理装置１は、心筋ＳＰＥＣＴ画像データを取得するデータ取得部１０と、心筋ＳＰＥＣＴ画像データに基づいてＣＦＲやＣＦＣといった心筋血流の循環指標を求める演算部１１と、求めた心筋血流の循環指標を出力する出力部１２とを有している。

【0016】

図２は、情報処理装置１で処理するデータを取得するための検査プロトコルおよび取得する心筋ＳＰＥＣＴ画像データについて説明する図である。まず、被験者に対して負荷薬剤を投与し、心臓に負荷をかける。ここで用いられる負荷薬剤は、アデノシン（adenosine）やジピリダモール（Dipyridamole）等である。なお、本実施の形態では、負荷薬剤により心臓に負荷をかけているが、被験者に運動をさせることで心臓に負荷をかけてもよい。

【0017】

続いて、被験者にＴｌ製剤（塩化タリウム（^２０１Ｔｌ））を投与し、プレダイナミック収集により心臓の位置を確認した後、ダイナミック収集によるＳＰＥＣＴ撮像を行い、ダイナミック画像を取得する。続いて、心電図同期ＳＰＥＣＴにより負荷後像（負荷時の心筋ＳＰＥＣＴ画像データ）を取得する。ダイナミック収集の時間は一例として、Ｔｌ製剤の投与から約６分間であり、心電図同期ＳＰＥＣＴの時間は、一例として７分間である。

【0018】

次に、Ｔｌ製剤の投与から３～４時間後に、心電図同期ＳＰＥＣＴにより再分布像を取得する。ここで、心電図同期ＳＰＥＣＴの時間は一例として７分間である。図２に示す手順で取得したダイナミック画像と負荷後像と再分布像のデータを情報処理装置１に入力する。

【0019】

図１に戻って、情報処理装置１の構成について説明する。演算部１１は、ｓＭＢＦ算出部２０と、ＷＯＲ算出部２１と、ｒＭＢＦ算出部２２と、ＣＦＲ算出部２３と、ＣＦＣ算出部２４とを有している。ｓＭＢＦ算出部２０は、ダイナミック画像から、時間放射能曲線（ＴＡＣ：Time activity curve）を求め、時間放射能曲線に基づいてｓＭＢＦを求める機能を有する。

【0020】

なお、本実施の形態では、ｓＭＢＦを求めるためにダイナミック画像を用いる例を取り上げて説明しているが、ダイナミック画像を用いない方法でｓＭＢＦを求めてもよい。一例として、ＳＵＶ（standardized uptake value）とｓＭＢＦとの相関関係を予め求め、その相関関係のデータを記憶部に記憶しておく。被験者へのＴｌ製剤の投与量、Ｔｌ製剤の心筋カウントおよび被験者の体重からＳＵＶを算出し、Ｓｔａｇｅ分類図によりｓＭＢＦを求めることとしてもよい。この場合、図２で説明した検査プロトコルにおいて、ダイナミック収集は不要である。

【0021】

ＷＯＲ算出部２１は、データ取得部１０にて取得した負荷後像と再分布像に基づいて、ＴｌのＷＯＲを算出する機能を有する。ＷＯＲは上述した公知の算出式に基づいて、算出することができる。

【0022】

ｒＭＢＦ算出部２２は、ｓＭＢＦとＷＯＲからｒＭＢＦを求める機能を有する。本実施の形態では、次の関係式が成り立つと仮定する。
ＷＯＲ＝（ｓＭＢＦ－ｒＭＢＦ）／ｒＭＢＦ・・・（１）
ｓＭＢＦ算出部２０およびＷＯＲ算出部２１にて、ｓＭＢＦとＷＯＲが求められているので、ｒＭＢＦ算出部２２は、上記の関係式（１）に基づいて、ｒＭＢＦの値を求める。

【0023】

ＣＦＲ算出部２３は、求めたｒＭＢＦとｓＭＢＦとに基づいてＣＦＲを求める機能を有する。ＣＦＣ算出部２４は、ＣＦＣ（Coronary Flow Capacity）を求める機能を有する。図３は、ＣＦＣのステージを求める方法の一例を説明するための図である。図３は、心筋血流ＰＥＴ検査で予め求めたｓＭＢＦとＣＦＲの相関関係からＳｔａｇｅ分類図としてそれぞれのカットオフ値を決められた図表である。情報処理装置１は、図３に示すテーブルを記憶部に記憶している。情報処理装置１は、図３に示すテーブルを参照し、心臓の各部位について求めたｓＭＢＦとＣＦＲの値に対応するＣＦＣのステージを求める。

【0024】

出力部１２は、演算部１１にて求めたＣＦＲおよびＣＦＣのデータを出力する機能を有する。出力部１２は、ＣＦＲおよびＣＦＣだけでなく、ｓＭＢＦのデータを出力してもよく、本実施の形態では、ｓＭＢＦとＣＦＲとＣＦＣを並べて表示する。

【0025】

図４（ａ）および図４（ｂ）は、出力画面の例を示す図である。図４（ａ）は、算出されたｓＭＢＦとＣＦＲとＣＦＣのデータをポーラーマップ表示した例を示し、図４（ｂ）は３Ｄ表示した例を示す。図４は白黒で記載しているが、実際の画面では、血流量に応じた色が付与されている。このようにポーラーマップあるいは３Ｄ表示によりＣＦＲ、ＣＦＣを表示することにより、診断に役立てることができる。

【0026】

本実施の形態の情報処理装置１のハードウェアの例は、ＣＰＵ、ＲＡＭ、ＲＯＭ、ハードディスク、ディスプレイ、キーボード、マウス、通信インターフェース等を備えたコンピュータである。上記した各機能を実現するモジュールを有するプログラムをＲＡＭまたはＲＯＭに格納しておき、ＣＰＵによって当該プログラムを実行することによって、上記した情報処理装置１が実現される。このようなプログラムも本発明の範囲に含まれる。

【0027】

図５は、本実施の形態の循環指標算出方法を示すフローチャートである。循環指標算出方法は、ＳＰＥＣＴ撮像を行う段階（Ｓ１０～Ｓ１３）と、ＳＰＥＣＴデータに基づく演算処理を行う段階（Ｓ１４～Ｓ２０）とに分けられる。演算処理は、上述した情報処理装置１によって行われる。

【0028】

ＳＰＥＣＴ撮像を行う段階では、まず、被験者に対して、アデノシンやジピリダモール等の負荷薬剤を投与し、心臓に負荷をかける（Ｓ１０）。被験者にＴｌ製剤を静注し、ダイナミック収集によるＳＰＥＣＴ撮像を行い（Ｓ１１）、続いて、負荷後像のＳＰＥＣＴ撮像を行う（Ｓ１２）。その後、しばらく時間をおいて（例えば３時間程度）、再分布像のＳＰＥＣＴ撮像を行う（Ｓ１３）。以上のようにして、１回のＴｌ製剤の投与により、ダイナミック画像と、負荷後像および再分布像を取得する。

【0029】

情報処理装置１は、上述の手順により取得されたＳＰＥＣＴデータを取得する（Ｓ１４）。情報処理装置１は、ダイナミック画像に基づいてｓＭＢＦを求め（Ｓ１５）、負荷後像と再分布像とに基づいてＷＯＲを算出する（Ｓ１６）。

【0030】

続いて、情報処理装置１は、上述した関係式（１）にｓＭＢＦとＷＯＲを代入して、ｒＭＢＦを算出し（Ｓ１７）、算出されたｒＭＢＦとｓＭＢＦとに基づいてＣＦＲを算出する（Ｓ１８）。次に、情報処理装置１は、ｓＭＢＦとＣＦＲとに基づいて、ＣＦＣを算出する（Ｓ１９）。情報処理装置１は、求めたｓＭＢＦとＣＦＲとＣＦＣのデータを出力する（Ｓ２０）。なお、Ｓ１７とＳ１８の工程を省略し、ｓＭＢＦとＷＯＲから直接ＣＦＣを算出してもよい。

【0031】

図６は、ｓＭＢＦをＳＵＶから求めてＣＦＲおよびＣＦＣを算出する方法を示すフローチャートである。循環指標算出方法は、ＳＰＥＣＴ撮像を行う段階（Ｓ３０～Ｓ３２）と、ＳＰＥＣＴデータに基づく演算処理を行う段階（Ｓ３３～Ｓ４０）とに分けられる。演算処理は、上述した情報処理装置１によって行われる。

【0032】

ＳＰＥＣＴ撮像を行う段階では、まず、被験者に対して、アデノシンやジピリダモール等の負荷薬剤を投与し、心臓に負荷をかける（Ｓ３０）。被験者にＴｌ製剤を静注し、負荷後像のＳＰＥＣＴ撮像を行う（Ｓ３１）。その後、しばらく時間をおいて（例えば３～４時間程度）、再分布像のＳＰＥＣＴ撮像を行う（Ｓ３２）。以上のようにして、１回のＴｌ製剤の投与により、負荷後像および再分布像を取得する。

【0033】

情報処理装置１は、まず、被験者に投与したＴｌ製剤の投与量と被験者の体重のデータの入力を受け付ける（Ｓ３３）。続いて、情報処理装置１は、上述の手順により取得されたＳＰＥＣＴデータを取得する（Ｓ３４）。情報処理装置１は、負荷後像とＴｌ製剤の投与量と被験者の体重とに基づいてＳＵＶを求め、ＳＵＶとｓＭＢＦとの相関関係からｓＭＢＦを求める（Ｓ３５）。また、情報処理装置１は、負荷後像と再分布像とに基づいてＷＯＲを算出する（Ｓ３６）。

【0034】

続いて、情報処理装置１は、上述した関係式（１）にｓＭＢＦとＷＯＲを代入して、ｒＭＢＦを算出し（Ｓ３７）、算出されたｒＭＢＦとｓＭＢＦとに基づいてＣＦＲを算出する（Ｓ３８）。次に、情報処理装置１は、ｓＭＢＦとＣＦＲとに基づいて、ＣＦＣを算出する（Ｓ３９）。情報処理装置１は、求めたｓＭＢＦとＣＦＲとＣＦＣのデータを出力する（Ｓ４０）。なお、Ｓ３７とＳ３８の工程を省略し、ｓＭＢＦとＷＯＲから直接ＣＦＣを算出してもよい。

【0035】

本実施の形態の循環指標算出方法によれば、１回のＴｌ製剤の投与によって取得したＳＰＥＣＴデータに基づいてＣＦＲを求めることができる。図７は、従来のＳＰＥＣＴデータの取得方法を示す図である。従来は、安静時のダイナミック収集を行うために、２回目のＴｌ投与を行う必要があったが、本実施の形態によれば、２回目のＴｌ投与が要らなくなり、簡便にＣＦＲを算出することができる。

【0036】

上記した実施の形態では、循環指標としてＣＦＲ，ＣＦＣを求める例を挙げたが、本発明の循環指標算出方法は、実施の形態で示した指標以外の指標を用いることができる。例えば、循環指標として、ＳＳＳ（Summed Stress Score）とＳＲＳ（Summed Rest Score）との差分であるＳＤＳ（Summed Difference Score）を求めてもよい。

【0037】

上記した実施の形態では、算出されたｓＭＢＦとＣＦＲとＣＦＣのデータをポーラーマップ表示する例を挙げたが、虚血心筋量と負荷心筋量と左室駆出率（Ejection Fraction、ＥＦ）の指標を同時に表示してもよい。

【0038】

３つの指標を同時に表示し、予後の良さを表すStageで分類してもよい。Stageは、虚血心筋量と負荷心筋量とＥＦのそれぞれのカットオフ値によって決まる。負荷心筋量（Stress flow）と虚血心筋量（ＣＦＲ）にＥＦを加えることにより予後層別化でき、予後評価を適切に行うことができる。

【実施例0039】

１．ＳＰＥＣＴデータの取得
複数の被験者に対し、アデノシン等を投与して薬剤負荷をかける。
被験者に塩化タリウム（²⁰¹Ｔｌ）注ＮＭＰ（商品名、日本メジフィジックス株式会社製）を静注して、ＳＰＥＣＴ撮像により、被験者の心筋のダイナミック画像を取得し、続いて、心電図同期ＳＰＥＣＴ画像（負荷後像）を取得する。薬剤負荷の影響がなくなるまで待ち、再び心電図同期ＳＰＥＣＴ画像（再分布像）を取得する。

【0040】

さらに、安静状態にある被験者に対して２回目の塩化タリウム（²⁰¹Ｔｌ）注ＮＭＰを静注して、ＳＰＥＣＴ撮像により、被験者の心筋のダイナミック画像を取得する。

【0041】

２．従来法によるＣＦＲの算出
負荷時のダイナミック画像からｓＭＢＦを算出すると共に、安静時のダイナミック画像からｒＭＢＦを算出し、求めたｓＭＢＦとｒＭＢＦの比によりＣＦＲを算出する。これは従来法によるＣＦＲの算出方法であり、この方法で求めたＣＦＲの値を比較例とする。また、ｓＭＢＦとＣＦＲとに基づいてＣＦＣの値を求め、比較例とする。

【0042】

３．本実施例によるＣＦＲの算出
上記１．で求めたＳＰＥＣＴデータのうち、２回目のダイナミック画像を用いないで、１回目のダイナミック画像と、負荷時画像と、再分布画像とを用いて本実施例の方法によりＣＦＲを算出する。すなわち、ダイナミック画像からｓＭＢＦを算出し、負荷時画像と再分布画像に基づいてＷＯＲを算出する。続いて、上述した関係式（１）にｓＭＢＦとＷＯＲを代入して、ｒＭＢＦを算出し、算出されたｒＭＢＦとｓＭＢＦとに基づいてＣＦＲを算出する。また、ｓＭＢＦとＣＦＲとに基づいてＣＦＣの値を求め、比較例とする。

【0043】

４．本実施例により算出されるＣＦＲおよびＣＦＣと比較例との比較
本実施例により算出されるＣＦＲと比較例とを比較すると、本実施例により算出されるＣＦＲは比較例のＣＦＲと精度良く一致する。また、ＣＦＣについても、比較例と本実施例の結果は精度良く一致する。

【符号の説明】

【0044】

１ 情報処理装置
１０ データ取得部
１１ 演算部
１２ 出力部
２０ ｓＭＢＦ算出部
２１ ＷＯＲ算出部
２２ ｒＭＢＦ算出部
２３ ＣＦＲ算出部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】