(57)【要約】【課題】被験者の個人差の如何に関わらず、冠動脈走行を反映した心臓の正しい診断をより容易に支援することのできる心疾患診断支援プログラム、心疾患診断支援方法及び心疾患診断支援装置を提供する。【解決手段】三次元冠動脈画像、左心室の三次元形態画像、三次元心臓ダイナミック核医学画像、各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得部１０２と、三次元形態画像を複数の領域に分割する第１の分割処理部１０３と、その分割した複数の領域を仕切る境界線を抽出し、三次元心臓ダイナミック核医学画像を分割する第２の分割処理部１０６と、その分割した領域毎に被験者のＣＦＲを算出するＣＦＲ算出部１０７と、ＣＦＲと、領域に対応する指標値との比較に基づき、治療方針を決定するための指標を求める指標算出部１０８と、指標算出部１０８で算出された指標を表示部１３０に表示する表示処理部１０９とを備える。【選択図】図１
2019191079.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータに、
被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する処理と、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する処理と、
前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する処理と、
前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出する処理と、
前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める処理と、
前記指標を表示部に表示する処理を実行させるための心疾患診断支援プログラム。
【請求項２】
請求項１に記載の心疾患診断支援プログラムにおいて、
前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする心疾患診断支援プログラム。
【請求項３】
被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得ステップと、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割ステップと、
前記第１の分割ステップで前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割ステップと、
前記第２の分割ステップで前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出ステップと、
前記ＣＦＲ算出ステップで算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出ステップと、
前記指標算出ステップで算出された前記指標を表示部に表示する表示ステップとを含む心疾患診断支援方法。
【請求項４】
請求項３に記載の心疾患診断支援方法において、
前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする心疾患診断支援方法。
【請求項５】
被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得部と、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割処理部と、
前記第１の分割処理部で前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割処理部と、
前記第２の分割処理部で前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出部と、
前記ＣＦＲ算出部で算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出部と、
前記指標算出部で算出された前記指標を表示部に表示する表示処理部とを備えることを特徴とする心疾患診断支援装置。
【請求項６】
請求項５に記載の心疾患診断支援装置において、
前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする心疾患診断支援装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心疾患診断支援プログラム、心疾患診断支援方法及び心疾患診断支援装置に関する。
【背景技術】
【０００２】
従来、心筋シンチグラフィは、被験者に放射性同位元素で標識された薬剤（放射性医薬品）を注射し、注射した薬剤から放出される放射線を撮影して、放射線量をコンピュータ処理して画像にし、心筋の血流やエネルギー代謝などをイメージングする検査である。この検査は、侵襲性が低く、適切な放射性医薬品を選択することで心臓に関する様々な生理学的、生化学的情報が得られるといった優れた特徴を有しているため、広く臨床において利用されている。例えば、薬物又は運動負荷時の心筋血流量との比によって求められるいわゆる冠動脈血流予備能（Coronary Flow Reserve：CFR）は、心疾患によって冠動脈がどの程度ダメージを受けているかを表す指標として用いられている。
【０００３】
上述した様に、CFRは、疾患により冠動脈の受けたダメージを表す指標として用いられているが、部分冠血流予備能比（fractional flow reserve：FFR）との乖離がみられる場合があることが指摘されている（非特許文献１、２参照）。そして、この乖離の態様を、治療方針の決定に利用することが提案されている。
【０００４】
一方、心筋シンチグラムの極座標表示の領域分割基準として、AHA（American Heart Association）の推奨する17segmentや、冠動脈の大まかな支配領域を表す3segmentなどが従来利用されてきた。
【先行技術文献】
【特許文献】
【０００５】
【非特許文献１】Nils P. Johnson et al.著、Is Discordance of Coronary Flow Reserve and Fractional Flow Reserve Due to Methodology or Clinically Relevant Coronary Pathophysiology?、2012年、p.193-202
【非特許文献２】Tim P. van de Hoef et al.著、 Physiological Basis and Long-Term Clinical Outcome of Discordance Between Fractional Flow Reserve and Coronary Flow Velocity Reserve in Coronary Stenoses of Intermediate Severity.、2014年、 p.301-11
【発明の概要】
【発明が解決しようとする課題】
【０００６】
一般的に用いられている心筋シンチグラムの診断では、上述した画一的に分類されたセグメントでの評価が行われている。しかし、これらの分割法は一律の基準で設定されるため、個人の冠動脈走行とは一致しないことがある。
【０００７】
すなわち、心臓の冠動脈走行には個人差があるので、17segmentや3segmentといった分割を行っても、各分割領域は、個人の冠動脈の支配領域を反映していない場合がある。例えば3Segmentで右冠動脈（Right Coronary Artery：以下、RCA）領域として表された場所が、実際のところは左前下行枝（Left Anterior Descending Coronary Artery：以下、LAD）領域であったというような乖離が生じることがある。よって、この様な方法では、心臓の冠動脈走行の個人差を考慮しながら冠動脈の疾患部位を特定することは困難である。
【０００８】
また、何らかの疾患のために冠動脈の1つに障害が起きた場合、従来の分割方法では、障害が起きた冠動脈が複数の分割領域にまたがっている場合がある。従って、従来の領域分割基準のみでは障害が起きている部位の同定が困難になるといった問題もある。
【０００９】
本発明は、上記課題に鑑みてなされたもので、被験者の個人差の如何に関わらず、冠動脈走行を反映した心臓の正しい診断をより容易に支援することのできる心疾患診断支援プログラム、心疾患診断支援方法及び心疾患診断支援装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明の一側面に係る心疾患診断支援プログラムは、コンピュータに、被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する処理と、前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する処理と、前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する処理と、前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出する処理と、前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める処理と、前記指標を表示部に表示する処理を実行させる。
【００１１】
本発明の別の一側面に係る心疾患診断支援方法は、被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得ステップと、前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割ステップと、前記第１の分割ステップで前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割ステップと、前記第２の分割ステップで前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出ステップと、前記ＣＦＲ算出ステップで算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出ステップと、前記指標算出ステップで算出された前記指標を表示部に表示する表示ステップとを含む。
【００１２】
本発明のまた別の一側面に係る心疾患診断支援装置は、被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得部と、前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割処理部と、前記第１の分割処理部で前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割処理部と、前記第２の分割処理部で前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出部と、前記ＣＦＲ算出部で算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出部と、前記指標算出部で算出された前記指標を表示部に表示する表示処理部とを備えることを特徴とする。
【発明の効果】
【００１３】
このように構成された本発明に係る心疾患診断支援プログラム、心疾患診断支援方法及び心疾患診断支援装置によれば、心臓疾患の診断をより正しく支援し得る。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態に係る心筋画像処理装置及びその周辺装置を示す機能構成図である。
【図２】左心室の三次元形態画像をボロノイ分割する方法を説明する説明図である（その１）。
【図３】左心室の三次元形態画像をボロノイ分割する方法を説明する説明図である（その２）。
【図４】左心室の三次元形態画像をボロノイ分割する方法を説明する説明図である（その３）。
【図５】左心室の三次元形態画像をボロノイ分割する方法を説明する説明図である（その４）。
【図６】心筋画像処理装置及びその周辺装置で行われる全体処理のフローチャートである。
【図７】第１の分割処理部によるボロノイ分割処理の結果を示す説明図である。
【図８】境界線極座標展開部による境界線の極座標展開処理の結果を示す説明図である。
【図９】冠動脈と左室心筋とを抽出する操作に供される画面の構成を示す画面構成図である。
【図１０】抽出済みの冠動脈から領域分割の起点となる部位を選択する操作に供される画面の構成を示す画面構成図である。
【図１１】支配領域の計算操作に供される画面の構成を示す画面構成図である。
【図１２】従来の3segmentで支配領域を決める方法と本法との比較結果を示す説明図である。
【図１３】従来の17segmentで支配領域を決める方法と本法との比較結果を示す説明図である。
【図１４】形態画像の一例である三次元ＣＴ画像から得られた、三次元冠動脈画像と三次元左心室画像を重ね合わせた画像の一例である。
【発明を実施するための形態】
【００１５】
以下、本発明の実施形態について、図面を用いて説明する。なお、すべての図面において、同様の構成要素には同一の符号を付し、適宜に説明を省略する。なお、以下の説明はあくまでも好ましい態様の一例を示したものであり、本発明の範囲を限定する意図ではない。
【００１６】
はじめに、本実施形態の概要について説明する。図１は本発明の実施形態に係る心疾患診断支援装置１００及びその周辺装置を示す機能構成図である。ここでは周辺装置として、撮像装置１１０、操作受付部１２０、表示装置１３０（表示部）及びデータベース１４０を例示するが、これに限られない。
【００１７】
撮像装置１１０は、被験者の核医学画像を撮像し、これを心疾患診断支援装置１００に提供する。このような撮像装置１１０としては、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置又はＰＥＴ(Positron Emission Tomography)装置等の核医学撮像装置がある。これらの装置により撮像された核医学画像は、特定の放射性同位元素でラベルされた薬剤（以下、放射性医薬品）を投与し、放射性医薬品より直接的または間接的に放出されたγ線を専用のカメラによって検出し、再構築することによって取得される。
【００１８】
なお、放射性医薬品としては、ＳＰＥＣＴ用製剤として塩化タリウム（２０１TlCl）注射液、テトロホスミンテクネチウム（９９ｍTc Tetrofosmin）注射液、ヘキサキス（２−メトキシイソブチルイソニトリル）テクネチウム（９９ｍTc MIBI）注射液、15-(4-ヨードフェニル)- 3(R,S)-メチルペンタデカン酸（123I）注射液や、ＰＥＴ用製剤として１３Ｎ−アンモニアや塩化ルビジウム（82RbCl）注射液などがある。また、単一の検査において、異なる２核種からなる放射性医薬品を投与してもよい。
【００１９】
心疾患診断支援装置１００は、情報取得部１０２、第１の分割処理部１０３、境界線極座標展開部１０４、心臓ダイナミック核医学画像正規化部１０５、第２の分割処理部１０６、ＣＦＲ算出部１０７、指標算出部１０８及び表示処理部１０９を備える。
【００２０】
情報取得部１０２は、被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する。ここでの取得には、予め生成されてデータベース１４０に保持されている画像や指標値をデータベース１４０から読み出すことが含まれる。三次元冠動脈画像とは、好ましくは、ＣＴや造影又は非造影のＭＲＩといった、撮像装置１１０とは別のモダリティ（図示せず）を使用した検査法で冠動脈を撮像して得られた冠動脈画像に対して三次元画像処理を施した画像のことをいうが、これに限定されない。つまり、被験者の冠動脈が描出されている画像であれば良い。なお、三次元冠動脈画像は、後述する三次元形態画像と別に撮像したものであっても良いが、冠動脈が描出されている三次元形態画像上で冠動脈を抽出して三次元画像処理を施したものであっても良い。後者の場合には、三次元冠動脈画像は、三次元形態画像と合わせて一つの画像として生成される。三次元画像処理としては、ＭＰＲ（Multi Planar Reconstruction）処理、三次元冠動脈抽出処理、三次元左心室抽出処理等を挙げることができる。
【００２１】
「少なくとも左心室」には、左心室の他に、左心室及び右心室の組み合わせが含まれる。「少なくとも左心室の三次元形態画像」には、左心室のみを撮像して得られた形態画像に対して三次元画像処理を施した画像と、左心室及び右心室の双方を撮像して得られた形態画像に対して三次元画像処理を施した画像とが含まれる。なお、形態画像は、心臓の形態を認識できる画像であれば、特に限定されない。形態画像としては、例えば、ＣＴや造影又は非造影のＭＲＩといった、撮像装置１１０とは別のモダリティ（図示せず）を使用した検査法で撮像された画像を三次元画像処理した画像の他、撮像装置１１０で得られた核医学画像に対して三次元画像処理を施した画像を、用いることができる。図１４に、形態画像の一例である三次元ＣＴ画像から得られた画像を示す。この画像では、ＣＴ画像から抽出された左心室の三次元形態画像上に、同じくＣＴ画像から抽出された左心室の三次元冠動脈画像が重ね合わせて表示されている。
【００２２】
心臓ダイナミック核医学画像は、心臓を検査対象として遂行される核医学測定により得られる画像であり、核医学測定の方法としてＳＰＥＣＴ又はＰＥＴを用いて得られた画像である。三次元心臓ダイナミック核医学画像とは、心臓ダイナミック核医学画像データに対して三次元画像処理を施した画像のことをいう。三次元冠動脈画像、三次元形態画像及び三次元心臓ダイナミック核医学画像は、データベース１４０に予め格納しておき、必要に応じてデータベース１４０から読み出せるようにしておいても良いし、心疾患診断支援装置１００の外部（例えば撮像装置１１０）からネットワーク等を介して読み出せるようにしておいても良い。
【００２３】
被験者における各冠動脈の局所的な狭窄状態を表す指標値は、例えば、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴから選択された指標値とする事ができる。ＦＦＲは、血流障害全体に対する冠動脈狭窄の寄与の割合を示す指標である。ＦＦＲは、センサーと大動脈圧の間の血管抵抗がない理想的な血管であった場合の何％の最大冠血流供給能力を有しているかを示している。ＦＦＲはその多くのエビデンスから、いずれの指標も心筋虚血の指標として有用な指標であると示されている。ＦＦＲはカテーテル検査室内では評価できる最も簡便かつ信頼性のある指標と言える。
【００２４】
第１の分割処理部１０３は、三次元冠動脈画像に基づいて、少なくとも左心室の三次元形態画像を複数の領域に分割する。好ましい態様において、この分割はボロノイ分割の手法により行う事ができる（詳細についてはhttp://gihyo.jp/dev/serial/01/geometry/0011を参照）。ボロノイ分割とは、与えられた点の近傍をその点の勢力圏と定義する分割方法であり、物理学、生態学、地域問題など様々な分野で用いられている主要な概念である。例えば、空間内にｎ個の点が与えられたとき、各点の勢力圏を定義し、これをボロノイ多角形で表すといった手法が用いられる。
【００２５】
分割の対象となる点は母点と呼ばれる。ボロノイ多角形の頂点はボロノイ点、辺はボロノイ辺と呼ばれる。ボロノイ辺はその両側にある母点から等距離にある点の軌跡、すなわち、その母点を結ぶ線分の垂直二等分線の一部である。また、ボロノイ点は三つの母点を頂点とする三角形の外心である。図２〜図５に示した例において、ボロノイ分割の対象となる直線Ａ，Ｂは冠動脈を模式的に示している（図２）。
【００２６】
直線Ａでは、所定の間隔で４つの母点Ａ１，Ａ２，Ａ３，Ａ４が特定される（図３〜図５）。同様に、直線Ｂでは、所定の間隔で４つの母点Ｂ１，Ｂ２，Ｂ３，Ｂ４が特定される（図３〜図５）。以下、本実施形態によって行われる画像処理における第１の分割処理の手順の一例を図３〜図５を参照しつつ列記する。
【００２７】
まず、図３に示すように、母点Ａ１，Ｂ１，Ｂ２を頂点とする三角形を形成する。そして、辺Ａ１Ｂ１の中点を点Ｍ１と定義し、辺Ａ１Ｂ２の中点を点Ｍ２と定義する。
次に、母点Ａ２，Ｂ１，Ｂ２を頂点とする三角形を形成する。そして、辺Ａ２Ｂ１の中点を点Ｍ３と定義し、辺Ａ２Ｂ２の中点を点Ｍ４と定義する。
次に、母点Ａ２，Ｂ２，Ｂ３を頂点とする三角形を形成する。そして、辺Ａ２Ｂ２の中点を点Ｍ４と定義し、辺Ａ２Ｂ３の中点を点Ｍ５と定義する。
次に、母点Ａ３，Ｂ２，Ｂ３を頂点とする三角形を形成する。そして、辺Ａ３Ｂ２の中点を点Ｍ６と定義し、辺Ａ３Ｂ３の中点を点Ｍ７と定義する。
次に、母点Ａ３，Ｂ３，Ｂ４を頂点とする三角形を形成する。そして、辺Ａ３Ｂ３の中点を点Ｍ７と定義し、辺Ａ３Ｂ４の中点を点Ｍ８と定義する。
次に、母点Ａ４，Ｂ３，Ｂ４を頂点とする三角形を形成する。そして、辺Ａ４Ｂ３の中点を点Ｍ９と定義し、辺Ａ４Ｂ４の中点を点Ｍ１０と定義する。
【００２８】
続いて、図４に示すように、母点Ｂ１，Ａ１，Ｂ２を頂点とする三角形を形成する。そして、辺Ｂ１Ａ１の中点を点Ｎ１と定義し、辺Ｂ２Ａ１の中点を点Ｎ２と定義する。
次に、母点Ｂ１，Ａ２，Ｂ２を頂点とする三角形を形成する。そして、辺Ｂ１Ａ２の中点を点Ｎ３と定義し、辺Ｂ２Ａ２の中点を点Ｎ４と定義する。
次に、母点Ｂ２，Ａ２，Ｂ３を頂点とする三角形を形成する。そして、辺Ｂ２Ａ２の中点を点Ｎ４と定義し、辺Ｂ３Ａ２の中点を点Ｎ５と定義する。
次に、母点Ｂ２，Ａ３，Ｂ３を頂点とする三角形を形成する。そして、辺Ｂ２Ａ３の中点を点Ｎ６と定義し、辺Ｂ３Ａ３の中点を点Ｎ７と定義する。
次に、母点Ｂ３，Ａ３，Ｂ４を頂点とする三角形を形成する。そして、辺Ｂ３Ａ３の中点を点Ｎ７と定義し、辺Ｂ４Ａ３の中点を点Ｎ８と定義する。
次に、母点Ｂ３，Ａ４，Ｂ４を頂点とする三角形を形成する。そして、辺Ｂ３Ａ４の中点を点Ｎ９と定義し、辺Ｂ４Ａ４の中点を点Ｎ１０と定義する。
図５は、中点群Ｍ１〜Ｍ１０，Ｎ１〜Ｎ１０を結んだ折れ線（ボロノイ辺）を基に得られた直線であって、直線Ａ，Ｂ間に仮想的に引かれた近似直線Ｌを示す。近似直線Ｌは、例えば、公知の最小二乗法により導き出される。冠動脈の様に蛇行している線の間に境界線を引く場合は、その境界線も蛇行した形となる。この場合、好ましい態様として、各中点（本例におけるＭ１〜Ｍ１０、Ｎ１〜Ｎ１０に相当）を結んだ線を作成し、必要に応じてスムージング処理を行う事によって得る事ができる。スムージングは公知の方法を用いることができ、例えば３点スムージングや５点スムージングの手法を用いることができる。
また、境界線の位置は、冠動脈の太さを考慮した補正を行っても良い。例えば、左右の冠動脈の直径が６：４の関係となっている場合は、上記の手法により得られる分割点を、母点間の中点ではなく、母点を結んだ線分を６：４に分割する点として定義することができる。
【００２９】
境界線極座標展開部１０４は、第１の分割処理部１０３で左心室の三次元形態画像を分割した複数の領域を仕切る境界線を抽出し、当該境界線を極座標展開する。心臓ダイナミック核医学画像正規化部１０５は、三次元心臓ダイナミック核医学画像を極座標展開する。なお、このような極座標展開する方法は、例えば特開２０１７−６２２１３号公報等で開示されているなど、既に周知であるため、詳細な説明は省略するが、例えば、撮像装置１１０により心臓を撮像して得られた断層画像に基づき、その断層画像から得られた画素の最大カウントを、心尖部を中心として心基部に向かって同心円上に放射状に展開する方法がある。
【００３０】
第２の分割処理部１０６は、極座標展開された三次元心臓ダイナミック核医学画像に同じく極座標展開された境界線を重畳した合成画像を生成する。第２の分割処理部１０６は、その生成した合成画像に基づいて極座標展開された三次元心臓ダイナミック核医学画像を第１の分割処理部１０３と同様にボロノイ分割で複数の領域に分割する。ここで分割された各領域を、その領域内に走行している冠動脈の支配領域とする。
【００３１】
ＣＦＲ算出部１０７は、第２の分割処理部１０６で三次元心臓ダイナミック核医学画像を分割した領域毎に被験者のＣＦＲを算出する。ＣＦＲは下記の式（１）で表される。
［数１］
ＣＦＲ＝薬物又は運動負荷時の心筋血流量／安静時の心筋血流量・・・（１）
【００３２】
なお、本例において、心筋血流量は、三次元心臓ダイナミック核医学画像を用いて算出される。心筋血流量の定量法は、公知の種々の方法を用いることができる。公知の方法としては、文献（Masao Miyagawa et al.著、Estimation of myocardial flow reserve utilizing an ultrafast cardiac SPECT: Comparison with coronary angiography, fractional flow reserve, and the SYNTAX score、2017年、p.347-353）記載の2コンパートメントモデルに基づく方法の他、例えば特許5695003号又は特許5700712号に開示された方法を用いることができる。
【００３３】
指標算出部１０８は、ＣＦＲ算出部１０７で算出された領域におけるＣＦＲと、情報取得部１０２が取得した領域に対応する、被験者における各冠動脈の局所的な狭窄状態を示す指標値との比較に基づき、治療方針を決定するための指標を求める。
この指標は、例えば、ＣＦＲとＦＦＲとの不一致の態様から、被験者の状態を大きく４つに分ける事ができる（指標Ａ〜Ｄ）。各指標Ａ〜Ｄにおける治療方針は、例えば下記の様に決定される。
Ａ：ＣＦＲの値もＦＦＲの値も保たれている状態（例えば、ＣＦＲ＞２．０、ＦＦＲ＞０．８）は、正常領域と判断される。
Ｂ：ＣＦＲの値は保たれているが、ＦＦＲの値が低下している状態（例えば、ＣＦＲ＞２．０、ＦＦＲ＜０．８）は、血管性病変が疑われるので、ＰＣＩ（Percutaneous Coronary Intervention）による治療で効果があると判断される。
Ｃ：ＣＦＲの値もＦＦＲの値も低下している状態（例えば、ＣＦＲ＜２．０、ＦＦＲ＜０．８）は、血管性病変と末梢性病変が共に疑われるので、バイパス術等の選択が検討される。
Ｄ：ＣＦＲの値が低下しているが、ＦＦＲの値は保たれている状態（例えば、ＣＦＲ＜２．０、ＦＦＲ＞０．８）は、末梢性病変が疑われるので、ＰＣＩによる治療の対象とはされず、薬物療法が選択される。
【００３４】
表示処理部１０９は、指標算出部１０８で算出された指標を表示装置１３０に表示する。例えば、表示処理部１０９は、指標算出部１０８で算出された指標を、第２の分割処理部１０６で生成された合成画像に対応付けて表示装置１３０に表示しても良い。
【００３５】
（心疾患診断支援装置１００及びその周辺装置で行われる全体処理の流れ）
図６は、心筋画像処理装置及びその周辺装置で行われる全体処理のフローチャートである。
まず、情報取得部１０２（図１）は、被験者における三次元冠動脈画像、左心室の三次元形態画像、三次元心臓ダイナミック核医学画像及びＦＦＲを取得する（ステップＳ１０１）。三次元冠動脈画像、左心室の三次元形態画像及び三次元心臓ダイナミック核医学画像は、それぞれ対応する機器を用いて生成された画像である。それぞれの画像やＦＦＲは、予めデータベース１４０に格納したものを読みだしても良いし、ネットワークを通じて各機器から直接取り込まれるものであっても良い。またこれらの情報は、ハードディスクや半導体メモリといった記憶媒体に記憶されたものを、周辺機器インターフェース（図示せず）等を介して取り込むものであっても良い。
【００３６】
次に、第１の分割処理部１０３（図１）は、三次元冠動脈画像に基づいて、左心室の三次元形態画像を複数の領域（図７）にボロノイ分割する（ステップＳ１０２）。ボロノイ分割処理で行われる一連の流れを図９，図１０，図１１に示す。図９は、冠動脈と左室心筋とを抽出する操作に供される画面構成を示す画面構成図である。図１０は、冠動脈から領域分割の起点となる部位を選択する操作に供される画面構成を示す画面構成図である。図１１は、三次元形態画像上で支配領域の境界線の計算操作に供される画面構成を示す画面構成図である。境界線の作成方法は、上述した通りである。
【００３７】
次に、境界線極座標展開部１０４（図１）は、第１の分割処理部１０３で左心室の三次元形態画像を分割した複数の支配領域を仕切る境界線を抽出し、当該境界線を公知の方法を用いて極座標展開する（ステップＳ１０３）。
【００３８】
次に、心臓ダイナミック核医学画像正規化部１０５（図１）は、三次元心臓ダイナミック核医学画像を公知の方法を用いて極座標展開する（ステップＳ１０４）。
【００３９】
次に、第２の分割処理部１０６は、心臓ダイナミック核医学画像正規化部１０５で極座標展開された三次元心臓ダイナミック核医学画像に、境界線極座標展開部１０４で極座標展開された境界線を重畳した合成画像（図８）を生成する（ステップＳ１０５）。ここで、図８に示されているSep,Ant,Latは、それぞれ、中隔、前壁、側壁を示している。第２の分割処理部１０６は、その生成した合成画像に基づいて、心臓ダイナミック核医学画像正規化部１０５で極座標展開された三次元心臓ダイナミック核医学画像を複数の支配領域（図８）に分割する。
【００４０】
図８の上段には、心尖部を中心として心基部に向かって同心円上に放射状に極座標展開した画像を示している。図８の下段には、複数の支配領域に仕切られた左心室の各部位に対応させるように、三次元心臓ダイナミック核医学画像に境界線を重畳した画像を示している。
【００４１】
次に、ＣＦＲ算出部１０７は、第２の分割処理部１０６で三次元心臓ダイナミック核医学画像を分割した領域毎に被験者のＣＦＲを算出する（ステップＳ１０６）。ここで算出されるＣＦＲの具体例については、上述した通りである。
【００４２】
次に、指標算出部１０８は、ＣＦＲ算出部１０７で算出された領域におけるＣＦＲと、領域に対応するＦＦＲとの比較に基づき、治療方針を決定するための指標を求める（ステップＳ１０７）。ここで算出される指標の具体例については、上述した通りである。
【００４３】
表示処理部１０９は、指標算出部１０８で算出された指標を表示装置１３０に表示する（ステップＳ１０８）。
【実施例１】
【００４４】
次に、図１２及び図１３を用いて、従来の3segmentや17segmentで冠動脈の支配領域を決める方法との比較により、本実施形態に係る心筋画像表示方法（以下、単に本法と略記する）の有用性を検討した実験結果について説明する。なお、本発明は、本実施例に限定されるものではない。
【００４５】
まず、図１２を参照しつつ、3segmentで支配領域を決める方法と本法との比較を行う。線ＢＬ１は、3segmentで使用される境界線である。線ＢＬ２は、本法によるボロノイ分割で得られた境界線である。図１２から明らかなように、3segmentで定めた支配領域は、本法で定めた境界線（BL2）によって画された支配領域と異なり複数の冠動脈の支配を受けていると言える。よって、本法で分割した複数の支配領域は、3segmentで支配領域を決める方法と比べてより患者の冠動脈走行を反映した結果を得る事ができると言える。
【００４６】
次に、図１３を参照しつつ、17segmentで支配領域を決める方法と本法との比較を行う。線ＢＬ３は、17segmentで使用される境界線である。線ＢＬ４は、本法によるボロノイ分割で得られた境界線である。図１３から明らかなように、17segmentで定めた支配領域は、本法で定めた境界線（BL4）によって画された支配領域と異なり複数の冠動脈の支配を受けていると言える。よって、本法で分割した複数の支配領域は、17segmentで支配領域を決める方法と比べてより患者の冠動脈走行を反映した結果を得る事ができると言える。
【００４７】
したがって、本実施形態によれば、CFRとFFRとの不一致の態様に基づく治療方針の決定を、冠動脈走行に基づく支配領域ごとに行うことができる。
【００４８】
以上、本実施形態によれば、被験者の個人差の如何に関わらず、心臓の正しい診断をより容易に支援し得ることのできる技術を提供することができる。
【００４９】
ここまで実施形態を示して本発明を説明したが、これらは一例である。また、本発明の各種の構成要素は、個々に独立した存在である必要はなく、複数の構成要素が単一の構成要素として構成されていること、一つの構成要素が複数の構成要素に分割されて形成されていること、ある構成要素が他の構成要素の一部であること、ある構成要素の一部と他の構成要素の一部とが重複していること、等を許容する。
【００５０】
また、上述した各種の構成要素は、必ずしも必須の構成要素ではなく、本発明の効果を阻害しない程度に省いても構わないし、同等に機能又は作用する他の構成要素に代えてもよい。
【００５１】
また、上記実施形態において示すフローも本発明の実施の一例であって、ここに示した手順に限られない。従って、上記フローチャートに図示した一のステップが複数に分離されて実行される、複数のステップが一のステップとして実行される、複数のステップが並行して実行される等の態様を許容する。
【００５２】
なお、上記実施形態では、第１の分割処理部１０３は、複数の三角形（図５）の各辺の中点を検出しその中点を結んでボロノイ辺を形成する例を示した。本発明はこの態様に限らず、例えば、第１の分割処理部１０３は、直線Ａ（図２）上の任意の点から直線Ｂ（図２）に向けて引いた直線のうち、最短となる直線の中点を結んでボロノイ辺を形成しても良い。また、そのように求めた中点につき、上述した方法により血管の太さの違いを反映した補正を行っても良い。
【００５３】
なお、上記実施形態では、近似直線Ｌ（図５）を公知の最小二乗法で導き出す例を示した。本発明はこの態様に限らず、近似直線Ｌは、一次相関などの手法により統計解析して導き出しても良い。要するに、近似直線Ｌは、昨今、一般的に行われている何れの方法を用いて導出されても良い。
【００５４】
なお、上記実施形態では、第２の分割処理部１０６が三次元心臓ダイナミック核医学画像を複数の領域に分割する前に、境界線極座標展開部１０４や心臓ダイナミック核医学画像正規化部１０５が極座標展開を行う例を示した。本発明はこの態様に限らず、第２の分割処理部１０６が三次元心臓ダイナミック核医学画像を複数の領域に分割する前に、必ずしも極座標展開を行う必要はない。
【００５５】
本実施形態は以下の技術思想を包含する。
（１）コンピュータに、
被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する処理と、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する処理と、
前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する処理と、
前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出する処理と、
前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める処理と、
前記指標を表示部に表示する処理を実行させるための心疾患診断支援プログラム。
（２）前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする（１）に記載の心疾患診断支援プログラム。
（３）被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得ステップと、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割ステップと、
前記第１の分割ステップで前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割ステップと、
前記第２の分割ステップで前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出ステップと、
前記ＣＦＲ算出ステップで算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出ステップと、
前記指標算出ステップで算出された前記指標を表示部に表示する表示ステップとを含む心疾患診断支援方法。
（４）前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする（３）に記載の心疾患診断支援方法。
（５）被験者における三次元冠動脈画像、少なくとも左心室の三次元形態画像、安静時及び負荷時の三次元心臓ダイナミック核医学画像、前記被験者における各冠動脈の局所的な狭窄状態を表す指標値を取得する情報取得部と、
前記三次元冠動脈画像に基づいて、前記左心室の前記三次元形態画像を複数の領域に分割する第１の分割処理部と、
前記第１の分割処理部で前記左心室の前記三次元形態画像を分割した複数の前記領域を仕切る境界線を抽出し、当該境界線を用いて前記三次元心臓ダイナミック核医学画像を分割する第２の分割処理部と、
前記第２の分割処理部で前記三次元心臓ダイナミック核医学画像を分割した前記領域毎に前記被験者のＣＦＲを算出するＣＦＲ算出部と、
前記ＣＦＲ算出部で算出された前記領域におけるＣＦＲと、前記領域に対応する前記指標値との比較に基づき、治療方針を決定するための指標を求める指標算出部と、
前記指標算出部で算出された前記指標を表示部に表示する表示処理部とを備えることを特徴とする心疾患診断支援装置。
（６）前記指標値は、ＦＦＲ（Fractional Flow Reserve）、ＦＦＲＣＴのうちいずれか一つであることを特徴とする（５）に記載の心疾患診断支援装置。
【符号の説明】
【００５６】
１００  心疾患診断支援装置
１０２  情報取得部
１０３  第１の分割処理部
１０６  第２の分割処理部
１０７  ＣＦＲ算出部
１０８  指標算出部
１０９  表示処理部
１３０  表示装置（表示部）
【図１】
2019191079.tif 000003
【図２】
2019191079.tif 000004
【図３】
2019191079.tif 000005
【図４】
2019191079.tif 000006
【図５】
2019191079.tif 000007
【図６】
2019191079.tif 000008
【図７】
2019191079.tif 000009
【図８】
2019191079.tif 000010
【図９】
2019191079.tif 000011
【図１０】
2019191079.tif 000012
【図１１】
2019191079.tif 000013
【図１２】
2019191079.tif 000014
【図１３】
2019191079.tif 000015
【図１４】
2019191079.tif 000016