JP2016531630A - 心筋灌流の病理学的特性評価のためのシステム及び方法 - Google Patents

Abstract

医用イメージングモダリティ１０によって連続的に取得された、関心被験者２０の心臓の少なくとも一部の複数の医用画像を解析することによって、心筋灌流病変の特性評価をする方法であって、複数の医用画像から、選択された心筋像位置の強度をサンプリングし、心筋像位置の対応するサンプリングされた強度に、取得順序を表す指標を割当て、強度曲線６０を得るステップと、得られた強度曲線６０に基づいて、複数の心筋セグメントのうちの少なくとも心筋セグメントのサブセットのうち、時空間的灌流不均一性又は灌流脱位相を示す指数６４、６６を計算するステップとを含む方法と、上記方法を実行することによって、複数の医用画像を解析することによって、心筋灌流病変の特性評価をするシステム５２と、上記システム５２を含む医用イメージングモダリティ１０と、上記方法を実行するためのソフトウェアモジュール４８とに関する。

Description

本発明は、心臓の医用イメージングの分野、特に、心臓の医用画像解析の分野に関する。

心臓磁気共鳴（ＣＭＲ）イメージング、より最近では、心臓コンピュータ断層撮影（ＣＣＴ）イメージングによる心臓の初回通過強調イメージングは、心筋灌流の定量化を可能とする。この定量化は、時間強度曲線の半定量分析又は真の定量分析を含む。半定量分析は、例えばピーク強度、最大アップスロープ、平均通過時間等である時間強度曲線の幾つかの特徴の定量化を含む。真の定量分析では、心筋において得られる動脈入力関数（ＡＩＦ）及び時間強度曲線の数理解析から実際の心筋血流が計算される。Jerosch-Heroldによる「Quantification of myocardial perfusion by cardiovascular magnetic resonance」（Journal of Cardiovascular Magnetic Resonance ２０１０、１２：５７）に半定量分析方法及び真の定量分析方法の両方についての詳細な考察が記載されている。

国際特許公開公報ＷＯ２００５／００４０６６Ａ１は、心灌流の定量的評価方法について説明している。この方法は、一連の心臓画像上に表現される心筋を、少なくとも１つの画像画素を含む画像セグメントに分割することと、各画像セグメントについて、心灌流パラメータを決定することとを含む。更に、正常灌流パラメータ値を有する少なくとも１つの画像セグメントが選択され、次に、残りの画像セグメントの心灌流パラメータは、当該正常灌流パラメータ値に基づいている。一実施形態では、灌流パラメータは、心筋における造影剤の分布についての時間強度プロファイルの最大アップスロープである。正常最大アップスロープが、少なくとも１つの画像セグメントについて導出され、相対最大アップスロープが、正常最大アップスロープに関して、各セグメントについて計算される。これらの値に基づいて、負荷時に導出される心筋灌流パラメータと、安静時に導出される心筋灌流パラメータとの比率が、各セグメントについて計算される。一例として、安静時及び負荷時に導出された相対最大アップスロープの比率として規定される心筋灌流予備能指数（Myocardial Perfusion Reserve Index：ＭＰＲＩ）が各セグメントについて計算されることが説明されている。

現在、心筋灌流は、特性評価及び分析のために、侵襲的方法を適用しなければならないという欠点がある。非侵襲的イメージング方法は、冠状動脈微小血管機能障害（ＭＶＤ）と多岐冠動脈疾患（左主冠動脈疾患ＣＡＤを含む）とを確実に区別することができない。というのは、これらの状態は、共に、広範性の心筋虚血をもたらし、視覚分析及び定量分析では、同様の外観を有するからである。したがって、患者に対して侵襲的血管造影が行われ、心筋虚血を有する患者における正常な冠状動脈が示された後、冠状動脈微小血管機能障害が診断される。

初回通過強調画像の解析は、心筋血流の真の定量化をもたらす。灌流欠損は、心筋血流の真の定量化によって解析可能であるが、非侵襲的方法は、虚血の分布及び重大度のみ評価可能であり、心筋の様々なセグメントにおける灌流の時空間的均質性を測定することができない。例えばＭＶＤ及びＣＡＤは、共に、通常、心筋の心内膜層への造影剤の延着に関連付けられる重大で空間的に広がった虚血によって特徴付けられる。

したがって、本発明は、関心の被験者の心臓の少なくとも一部の複数の医用画像を解析することによって心筋灌流病変の特性評価をする改良方法を提供することを目的とする。当該複数の医用画像は、医用イメージングモダリティによって連続的に取得される。

この目的は、以下のステップ：
複数の医用画像における関心被験者の心臓の選択された部分の輪郭を描写し、選択された部分を複数のセグメントにセグメント化するステップと、
複数の医用画像から、選択された心筋像位置の強度をサンプリングし、心筋像位置の対応するサンプリングされた強度に、医用画像の各医用画像の取得の順序を表す指標を割当て、選択された心筋像位置のそれぞれに対し、強度曲線を得るステップと、
得られた強度曲線に基づいて、複数の心筋セグメントのうちの少なくとも心筋セグメントのサブセットのうち、時空間的灌流不均一性又は灌流脱位相を示す指数を計算するステップとを含む方法によって達成される。

「心臓の選択された部分」との表現は、本願において使用される場合は、具体的に、心臓の任意の血管腔と理解されるものとし、また、左心室及び上行大動脈を含むと理解されるものとする。

「医用イメージングモダリティ」との表現は、本願において使用される場合は、具体的に、心臓磁気共鳴（ＣＭＲ）イメージングデバイス、心臓コンピュータ断層撮影（ＣＣＴ）イメージングデバイス、冠動脈造影（ＣＡ）デバイス、ＣＣＴ血管造影（ＣＣＴＡ）デバイス、血管内超音波（ＩＶＵＳ）デバイス、単光子放出コンピュータ断層撮影（ＳＰＥＣＴ）デバイス、ポジトロン放出断層撮影（ＰＥＴ）デバイス及び心エコー検査デバイスを含むものとする。

「時空間的灌流不均一性」との表現は、本願では、「時空間的脱位相」とも表現され、具体的に、病理学的異常がある場合に、不均一な心筋血流の時間的及び空間的分布と理解されるものとする。

心臓の選択された部分の輪郭を描写するステップは、手動で、半自動的に又は全自動的に行われてよい。適切なセグメント化技術は、当技術分野において知られており、市販されているので、本明細書では、より詳細には説明されない。

心筋セグメントのサブセットは、灌流領域といった複数の心筋セグメントの厳密なサブセットを含んでも、心筋全体を含んでもよい。

各医用画像の取得の順序を表す指標は、時間的尺度を参照してもよいし、時間そのものであってもよい。

方法の１つの利点は、特定の心筋病変を、高速及び単純に特徴付けるために使用することができる、心臓の選択された部分の様々な領域への灌流の分布に関する追加の時間的情報が提供される点である。具体的には、計算された指数を、冠動脈疾患と微小血管疾患とを特性評価において区別するために使用することができる。

方法のもう１つの利点は、非侵襲的方法によって、追加情報が提供される点であり、これにより、多くの場合、侵襲的な血管造影評価を回避することができる。

心臓の選択された部分の輪郭を描写するステップを適用する前に、複数の医用画像は、関心被験者の呼吸動作に対し補正するために、画像位置合わせ技術が施されていてもよい。適切な位置合わせ技術は、当技術分野においてよく知られているため、本明細書では、詳細には説明されない。

更に、得られた強度曲線は、より優れた結果のためにフィルタリングが施されていてもよい。当業者に適している任意のフィルタリング技術が採用されてよい。

好適な実施形態では、方法は更に、複数のセグメントの各セグメントにおける心筋血流の定量化を行うステップを含む。心筋血流の定量化を行うステップは、例えばデコンボリューション方法といった背景技術のセクションにおいて提示された検討のための文献に説明されている従来技術の１つによって、又は、当業者が適切であると見なす任意の他の技術によって行われてよい。各セグメントにおける心筋血流の真の定量化は、心筋灌流病変の特性評価をするための補助的情報を提供する。

方法の別の実施形態では、複数の医用画像は、関心被験者に造影剤を投与した後に、医用イメージングモダリティによって取得される。「造影剤」との表現は、本願において使用される場合は、具体的に、医用イメージングモダリティによって取得される場合に、造影剤を囲む関心被験者の組織よりも、ベースラインに比べてより大きい信号を生成する任意の造影剤と理解されるものとする。医用イメージングモダリティが、ガンマ線を使用する動作原理に基づいている場合、「造影剤」との表現は、放射性トレーサー物質も含むものとする。このようにすると、信号対雑音比が向上し、心筋灌流病変のより正確な特性評価が達成できる。

一実施形態では、複数の医用画像は、関心被験者に造影剤を投与した後に、心臓を通る造影剤の初回通過の間に取得される。

別の実施形態では、複数の医用画像は、関心被験者に造影剤を投与した後に、造影剤の濃度が、関心被験者内で定常状態に到達する時間の間に取得される。これは、負荷エコー灌流イメージングによって医用画像を取得する場合に、特に重要である。

一実施形態では、複数の医用画像は、例えば動脈スピンラベリングといった内因性のコントラストを使用して、心臓を通る造影剤の初回通過の間に取得される。

別の好適な実施形態では、方法は更に、心臓の選択された部分における基準場所を特定するステップを含み、指数を計算するステップにおいて、強度曲線は、特定された基準場所によって決定される基準時間を参照して評価される。これにより、指数を計算するための正確な基準が提供される。好適な基準場所は、左心室である。これは、左心室が、造影剤を、初回通過時に受け取り、選択された心筋像位置のそれぞれについて得られた強度曲線におけるアップスロープに先行する場所だからである。「アップスロープ」との表現は、本願において使用される場合は、具体的に、得られた強度曲線の強度が、造影剤の投与によって、所定の閾値を、造影剤がない強度に比べて、例えば１０％超えた瞬間と理解されるものとする。

別の好適な実施形態では、方法は更に、選択された心筋位置のそれぞれについて、心筋像位置のそれぞれのサンプリングされた強度の特徴の発生まで、特定された基準場所によって決定される基準時間に対する個々の期間を自動的に決定するステップを含む。個々の期間は、指数を計算するステップにおいて使用される。特徴は、強度上昇の開始の時点、ピーク強度の時点、又は、当業者に適切と思われる任意の他の特徴であってよい。このようにすると、時空間的灌流不均一性を示す指数が容易に自動的に計算できる。

好適には、指数を計算するステップは、特定された基準場所における特徴の発生の時間に対する個々の心筋位置のそれぞれにおける特徴の発生までの時間の変動を示す統計的尺度の計算を含む。これにより、心筋セグメントのうち、時空間的灌流不均一性を記述する指数が、非常に顕著な方法で提供される。統計的尺度は、一組の数がどれだけ遠くまで広がっているかを示す尺度として統計学において習慣的に使用される分散の形を有してよい。この意味で、分散は、一組の数の標準偏差の二乗であってよい。一般に、統計的尺度は、個々の心筋位置のそれぞれにおける特徴の発生までの時間の変動を示すのに適していると当業者が考える任意の他の形式であってもよい。

好適には、関心被験者の心臓の少なくとも一部の複数の医用画像を取得することは、関心被験者の心臓の周期動作と少なくとも部分的に同期される。例えば医用画像は、ＱＲＳ群のＲピークといった心電図における基準イベントの前又は後の固定時間量において取得される。方法のこの実施形態の利点は、複数の医用画像のすべての医用画像が、同様の心臓状態において取られるので、医用画像間に心筋心室の動きがほとんどなく、心筋が比較的静止した状態でレンダリングされる点である。

方法の更に別の実施形態では、心筋像位置の強度をサンプリングするステップにおいて、心筋像位置は、心筋と平行な方向及び心筋を横断する方向において選択される。このようにすると、計算された指数は、心筋セグメントのうちの時空間的灌流不均一性を特に適切に表す。

方法は更に、灌流図を生成し、灌流図をユーザに表示するステップを含む。「灌流図（perfusogram）」との用語は、本願において使用される場合は、具体的に、位置及び時間の関数として、心筋における像位置の強度のカラー表現と理解されるものとする。灌流図の水平軸は、時間を表し、垂直軸は、心筋における位置（即ちセグメント）を表す。「灌流図」は、例えばMarcel Breeuwerによる「Comprehensive visualization of first-pass myocardial perfusion: The uptake movie and the perfusogram」（International Society for Magnetic Resonance in Medicine：ＩＳＭＲＭ２００２）である先の出版物において、心筋灌流を可視化するための有用なツールであると説明されている。例えば灌流図は、複数の医用画像を取得するために用いた医用イメージングモダリティに通常あるモニタユニット上に表示される。

更に、方法は、少なくとも１つの特徴位置及び／又は少なくとも１つの特徴時点を示す少なくとも１つのマーカーを灌流図内に組み込むステップを含むことが好適である。ユーザに素早く灌流図内に含まれる情報を伝達できるこのようなマーカーの例は、冠動脈灌流領域の位置、（例えばＡＨＡ１７セグメントモデルによる）左心室のセグメントの境界、並びに、開始時間及びピーク時間といった造影剤の通過時間における特徴的な瞬間である。マーカーは、直線、曲線、又は、心筋の領域をマーク付けする閉ループとして形成されてよい。

別の好適な実施形態では、方法は更に、複数のコンピュータリンクを組み込むステップを含み、複数のコンピュータリンクの各コンピュータリンクは、灌流図内のある場所に割当てられ、複数のコンピュータリンクの各コンピュータリンクは、複数の医用画像の１つの医用画像を表すデータセットにリンク付けされる。このようにすると、心筋血流に関する詳細情報が、ユーザ、通常は、医療スタッフメンバーに容易に提供される。

本発明は、関心被験者の心臓の少なくとも一部の複数の医用画像を解析することによって、心筋灌流病変の特性評価をするシステムを提供することをもう１つの目的とする。複数の医用画像は、医用イメージングモダリティによって連続的に取得される。当該システムは、
複数の医用画像における関心被験者の心臓の選択された部分の輪郭を描写し、選択された部分を複数のセグメントにセグメント化する描写ユニットと、
強度サンプリング及び解析ユニットとを含む。

強度サンプリング及び解析ユニットは、
複数の医用画像から、心筋像位置の強度をサンプリングし、心筋像位置の対応するサンプリングされた強度に、複数の医用画像の各医用画像の取得の順序を表す指標を割当て、
複数の心筋セグメントのうちの少なくとも心筋セグメントのサブセットのうち、時空間的灌流不均一性又は灌流脱位相を示す指数を計算する。システムは、方法について上記された利点と同じ利点を提供できる。

本発明の別の態様では、心筋灌流病変の特性評価をするシステムは、複数の医用画像を取得するために用いられた医用イメージングモダリティの一体部分である。好適には、医用イメージングモダリティは、磁気共鳴イメージング装置としてデザインされる。磁気共鳴イメージング装置は、有利には、医用画像の取得を、関心被験者の心臓の周期動作に同期させる同期化手段を含む。

本発明の更に別の態様では、上記された心筋灌流病変の特性評価する方法の何れかの一実施形態又はその組み合わせを実行するためのソフトウェアモジュールが提供される。行われるべき方法のステップは、ソフトウェアモジュールのプログラムコードに変換され、プログラムコードは、医用イメージングモダリティの制御ユニットのメモリユニットに実装可能であり、医用イメージングモダリティの制御ユニットのプロセッサユニットによって実行可能である。

制御ユニットは、医用イメージングモダリティの制御機能に通常使用される制御ユニットであってよい。或いは、制御ユニットは、方法ステップを実行するように、特別に割当てられた追加の制御ユニットであってもよい。

ソフトウェアモジュールは、方法のロバストで信頼できる実行を可能にし、また、方法ステップの高速変更及び／又は画像位置合わせアルゴリズムの適応を可能にする。

本発明のこれらの及び他の態様は、以下に説明される実施形態を参照して明らかとなろう。

図１は、磁気共鳴イメージングシステムとしてデザインされる本発明による医用イメージングモダリティの一実施形態の一部の概略図である。 図２は、対照個体の心臓を通る（上）、微小血管機能障害を有する個体の心臓を通る（中央）、及び、３枝冠動脈疾患を有する個体の心臓を通る（下）造影剤の初回通過時の心筋灌流の、本発明によって取得される強度曲線の概略図である。 図３は、図２の強度曲線から計算される時空間的灌流不均一性を示す計算された指数を含む図を示す。 図４は、灌流図（perfusogram）の一例を示す。

図１は、人間である関心被験者２０の心臓の少なくとも一部の医用画像（磁気共鳴画像で表される）を取得する、磁気共鳴イメージングシステムとしてデザインされる本発明による医用イメージングモダリティ１０の一実施形態の一部の概略図を示す。磁気共鳴イメージングシステムは、静磁場を生成するために設けられる主磁石１４を有する磁気共鳴スキャナ１２を含む。主磁石４は、関心被験者２０がその中で位置付けられるべき検査空間１６を中心軸１８の周りに提供する中心ボアを有する。図１では、明瞭にする理由から、関心被験者２０を支えるための従来の台は省略されている。実質的に静的な磁場は、中心軸１８に並列に位置合わせされる検査空間１６の軸方向を規定する。更に、磁気共鳴イメージングシステム１０は、静磁場に重ねられる勾配磁場を生成するために設けられる磁場勾配コイルシステム２２を含む。磁場勾配コイルシステム２２は、当技術分野において知られているように、主磁石１４のボア内に同心円状に配置される。

原理上、本発明は、静磁場内に検査領域を提供する任意の他のタイプの磁気共鳴イメージングシステムにも適用可能である。更に、当然ながら、本発明は、関心被験者の心臓の少なくとも一部の医用画像を取得する任意の他の医用イメージングモダリティと共に使用されてもよい。本発明を適用可能である医用イメージングモダリティの例としては、心臓コンピュータ断層撮影（ＣＣＴ）イメージングデバイス、冠動脈造影（ＣＡ）デバイス、ＣＣＴ血管造影（ＣＣＴＡ）デバイス、血管内超音波（ＩＶＵＳ）デバイス、単光子放出コンピュータ断層撮影（ＳＰＥＣＴ）デバイス、ポジトロン放出断層撮影（ＰＥＴ）デバイス及び心エコー検査デバイスが挙げられる。

更に、磁気共鳴イメージングシステムは、関心被験者２０の又は関心被験者２０内の原子核を励起させるために、無線周波数送信位相時に、検査空間１６に無線周波数電磁場を印加するように設けられる全身コイル２４としてデザインされる無線周波数アンテナを含む。全身コイル２４は更に、無線周波数受信位相時に、関心被験者２０の又は関心被験者２０内の励起された原子核からの磁気共鳴信号を受信するためにも設けられている。磁気共鳴イメージングシステムの動作状態において、無線周波数送信位相及び無線周波数受信位相は、連続的に起きる。全身コイル２４は、中心軸を有し、動作状態では、全身コイル２４の中心軸と磁気共鳴イメージングシステムの中心軸１８とが一致するように、主磁石１４のボア内に同心円状に配置される。当技術分野においてよく知られているように、円筒状の金属製無線周波数シールド２６が、磁場勾配コイルシステム２２と全身コイル２４との間に同心円状に配置される。

磁気共鳴イメージングシステムは更に、少なくとも、磁気共鳴スキャナ１２及び磁場勾配コイルシステム２２の機能を制御するために設けられる制御ユニット２８を含む。制御ユニット２８は、通常のモニタユニット３６を含む。

更に、磁気共鳴イメージングシステムは、制御ユニット２８に接続され、制御ユニット２８によって制御される無線周波数送信器ユニット３０を含む。無線周波数送信器ユニット３０は、無線周波数送信位相時に、無線周波数スイッチングユニット３２を介して、磁気共鳴無線周波数の無線周波数電力を全身コイル２４に供給するように設けられる。無線周波数受信位相時、無線周波数スイッチングユニット３２は、全身コイル２４からの磁気共鳴信号を、制御ユニット２８内にある画像処理ユニット３４に送る。画像処理ユニット３４は、取得された磁気共鳴信号から被験者２０の一部の磁気共鳴画像を生成するように、取得した磁気共鳴信号を処理する。この技術の多くの様々な変形例が当業者にはよく知られているので、本明細書では、更に詳細に説明する必要はない。

関心被験者２０の心臓の磁気共鳴画像を取得するために、磁気共鳴イメージングシステムは、医用画像の取得を、関心被験者の心臓の周期動作に同期させる同期化手段を含む。同期化手段は、心電図デバイス３８及び同期化ユニット４２として形成される。

心電図デバイス３８は、関心被験者２０の心臓の心電図データを測定するために設けられる。このために、心電図デバイス３８の複数の電極４０が、関心被験者２０上に配置される。更に、心電図デバイス３８は、勾配磁場によって生成されるアーチファクトを減らすために、心電図データをフィルタリングする手段を含む。当業者には、適切なフィルタリング手段が知られているので、本明細書では、より詳細に説明しない。

心電図デバイス３８は、同期化ユニット４２に結合される。同期化ユニット４２は、心臓活動のＱＲＳ群のＲピークの検出から磁気共鳴信号を取得する取得期間をトリガするトリガ信号５０を生成する。同期化ユニット４２は、次に、制御ユニット２８に結合される。制御ユニット２８は、勾配磁場を生成する磁場勾配コイルシステム２２の制御信号を生成するために、同期化ユニット４２から提供されるトリガ信号５０によって同期される。このために、制御ユニット２８は、トリガ信号５０を受信すると、複数のシーケンスを生成する。各シーケンスは、無線周波数磁場及び勾配磁場を含む。

医用イメージングモダリティ１０は、関心被験者２０の心臓の少なくとも一部の複数の医用画像を解析することによって心筋灌流を病理学的に特性評価する心筋灌流の病理学的特性評価システム５２を含む。心筋灌流の病理学的特性評価システム５２は、制御ユニット２８のハウジング内にあり、描写ユニット５４、血流アナライザ５６、及び、強度サンプリング及び解析ユニット５８を含む。これらのデバイスの機能及び相互作用については、以下に詳細に説明する。

関心被験者２０に造影剤を投与した後、関心被験者２０の心臓を造影剤が初回通過する際に、磁気共鳴イメージングシステムによって、関心被験者２０の心臓の磁気共鳴画像の形である複数の医用画像が、連続的に取得される。複数の医用画像が、３人の異なる関心被験者、即ち、対照個体、３枝冠動脈疾患を有する個体及び微小血管機能障害を有する個体について、取得される。関心被験者２０の心臓の複数の医用画像の取得は、上記されたように、関心被験者２０の各自の心臓の周期動作に同期されている。

各関心被験者２０の心臓の複数の医用画像を取得した後、医用画像は、本発明による心筋灌流病変の特性評価をする方法の一実施形態によって、解析される。この方法は、一個体の複数の医用画像に適用されて、以下に詳細に説明される。当然ながら、この方法は、残りの２つの個体の複数の医用画像にも同様に適用される。

この方法を実行可能とするために、心筋灌流の病理学的特性評価システム５２は、ソフトウェアモジュール４８（図１）を含む。行われるべき方法ステップは、ソフトウェアモジュール４８のプログラムコードに変換される。プログラムコードは、心筋灌流の病理学的特性評価システム５２のメモリユニット４４に実装可能であり、心筋灌流の病理学的特性評価システム５２のプロセッサユニット４６によって実行可能である。

準備ステップにおいて、複数の医用画像の各医用画像は、画像登録アルゴリズムに提出される。

方法の第１のステップにおいて、複数の医用画像内に、関心被験者２０の心臓の左心室の輪郭が描写され、左心室は、心筋灌流の病理学的特性評価システム５２の描写ユニット５４によって、複数のセグメントにセグメント化される。

方法の次のステップにおいて、心筋灌流の病理学的特性評価システム５２の血流アナライザ５６が、複数のセグメントの各セグメントにおける心筋血流の真の定量化を行う。

次に、心筋灌流の病理学的特性評価システム５２の強度サンプリング及び解析ユニット５８によって、心筋像位置が選択される。心筋像位置は、心筋に平行する方向と心筋を横断する方向とにおいて選択される。次に、強度サンプリング及び解析ユニット５８は、複数の医用画像から、選択された心筋像位置の強度をサンプリングし、心筋像位置の対応するサンプリングされた強度に、各医用画像の取得時間によって形成される取得順序を表す指標を割当て、選択された心筋像位置のぞれぞれに対し、経時的な強度曲線６０が得られる。

図２に、３人の異なる関心被験者について得られた経時的な強度曲線６０が概略的に示される（上：対照個体、中央：微小血管機能障害（ＭＶＤ）を有する個体、下：３枝冠動脈疾患（ＣＡＤ）を有する個体）。

図２における経時的な強度曲線６０について、適切な時間的尺度を提供するために、強度サンプリング及び解析ユニット５８によって、左心室における基準場所が特定される。この基準場所は、造影剤の初回通過時に造影剤を最初に受け取る場所であり、選択された心筋像位置の強度の上昇の開始の前に来る。基準場所における強度の上昇の開始の時点が、基準時間の役割を果たす。経時的な強度曲線６０は、特定された基準場所及び基準時間を参照して評価される。

方法の次のステップにおいて、選択された心筋位置のそれぞれについて、特定された基準場所及び対応する基準時間に対する個々の期間が、各心筋像位置のサンプリングされた強度曲線６０において特徴が生じるまで、自動的に決定される。この実施形態では、特徴とは、経時的な強度曲線のピーク強度である。ピーク強度までの個々の期間ＴＴＰＩは、指数６４、６６を計算する次のステップにおいて使用される。

方法の後続のステップにおいて、強度サンプリング及び解析ユニット５８は、得られた経時的な強度曲線６０に基づいて、複数の心筋セグメントのうちの心筋セグメントのサブセットのうち、時空間的灌流不均一性又は脱位相を示す第１の指数６４を計算する。

左心室灌流の時間的脱位相を示す第１の指数６４は、次のステップにおいて、基準時間と、経時的な強度曲線６０のピーク強度の発生までの時間との間のピーク強度までの期間ＴＴＰＩの分散として計算される。ここでは、分散とは、期間ＴＴＰＩの標準偏差の二乗であると理解される。

左心室灌流の時間的脱位相を示す第２の指数６６は、基準時間と、経時的な強度曲線６０のピーク強度の発生までの時間との間のピーク強度までの期間ＴＴＰＩの変動係数として計算される。

図３に示される、計算された指数６４、６６についての結果は、対照個体（左列ＣＴＲＬ）と、３枝冠動脈疾患を有する個体（中央列ＣＡＤ）と、微小血管機能障害を有する個体（右列ＭＶＤ）との間には、著しい差があることを明らかに示す。

図２における経時的な強度曲線６０から明らかとなるように、ＭＶＤ及び３枝ＣＡＤは、共に、通常、心筋の心内膜層への造影剤の延着に関連付けられる重大で空間的に広がった虚血によって特徴付けられる。しかし、時空間的脱位相解析を追加することによって、左心室の様々な部位への灌流の時間的分布に関する追加情報が提供される。図３は、計算された指数６４、６６が、元になっている病態生理学の違いによるＣＡＤとＭＶＤとの主要な差を構成し、特定の特徴付けを可能にすることを示す。

ＣＡＤを有する個体では、心筋における心筋血流の時空間的分布は、大部分が、不均一である。ＭＶＤを有する個体では、病理学的変化は、心筋収縮を有する心収縮における微視的循環及びその相互作用が伴う。この場合、心外膜層への灌流は、妨げられることなく、非常に均一である。しかし、最初通過波の経壁伝播に遅延がある。ＣＡＤと同様に、これは、強度の上昇開始が遅延した広がった虚血を引き起こす。しかし、この特徴は、心筋全体に亘って時間的領域において均一であるため、ＣＡＤとＭＶＤとの非侵襲的な区別が可能となる。

方法の別のステップにおいて、灌流図６２が生成され、モニタユニット３６によって、ユーザに対して表示される。図４は、グレイスケール又はカラーコーディングを有する灌流図６２の一例を示す。灌流図６２には、特徴位置及び／又は特徴時点を示す直線のようなマーカー（図示せず）が組み込まれる。閉ループとして形成され、心筋があまりよく灌流されていない時間及び空間のレジームを示す他のマーカー（図示せず）が灌流図６２に挿入されてもよい。

本発明の更に別のステップにおいて、灌流図６２内に、複数のコンピュータリンクが組み込まれる。複数のコンピュータリンクの各コンピュータリンクは、灌流図６２内のある場所に割当てられ、複数のコンピュータリンクの各コンピュータリンクは、複数の医用画像のうちの１つの医用画像を表すデータセットにリンク付けされる。図４に示されるコンピュータマウスカーソルは、心筋の特定の時点及び特定の場所を表す灌流図６２の場所を指し示している。マウスボタンをクリックすることによって、心筋の当該時点及び場所に対応する、医用画像を表すデータセットへのコンピュータリンクがアクティブにされ、当該医用画像が、磁気共鳴イメージングシステムのモニタユニット３６上に表示される。

本発明は、図面及び上記説明において詳細に例示され、説明されたが、当該例示及び説明は、例示的に見なされるべきであり、限定的に見なされるべきではない。本発明は、開示される実施形態に限定されない。開示された実施形態の他の変形態様は、図面、開示内容及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解され、実施される。請求項において、「含む」との用語は、他の要素又はステップを排除するものではなく、また、「ａ」又は「ａｎ」との不定冠詞も、複数形を排除するものではない。特定の手段が相互に異なる従属請求項に記載されることだけで、これらの手段の組み合わせを有利に使用することができないことを示すものではない。請求項における任意の参照符号は、範囲を限定するものと解釈されるべきではない。

１０ 医用イメージングモダリティ
１２ 磁気共鳴スキャナ
１４ 主磁石
１６ 検査空間
１８ 中心軸
２０ 関心被験者
２２ 磁場勾配コイルシステム
２４ 全身コイル
２６ 無線周波数シールド
２８ 制御ユニット
３０ 無線周波数送信器ユニット
３２ 無線周波数スイッチングユニット
３４ 画像処理ユニット
３６ モニタユニット
３８ 心電図デバイス
４０ 電極
４２ 同期化ユニット
４４ メモリユニット
４６ プロセッサユニット
４８ ソフトウェアモジュール
５０ トリガ信号
５２ 心筋灌流の病理学的特性評価システム
５４ 描写ユニット
５６ 血流アナライザ
５８ 強度サンプリング及び解析ユニット
６０ 経時的な強度曲線
６２ 灌流図
６４ 第１の指数
６６ 第２の指数
ＡＩＦ 動脈入力関数
ＬＶ 左心室
ＴＴＰＩ ピーク強度までの時間

Claims (15)

1. 関心被験者の心臓の少なくとも一部の複数の医用画像を解析することによって、心筋灌流病変の特性評価をする方法であって、前記複数の医用画像は、医用イメージングモダリティによって連続的に取得され、当該方法は、
   前記複数の医用画像における前記関心被験者の前記心臓の選択された部分の輪郭を描写し、選択された前記部分を複数のセグメントにセグメント化するステップと、
   前記複数の医用画像から、選択された心筋像位置の強度をサンプリングし、前記心筋像位置の対応するサンプリングされた強度に、前記複数の医用画像の各医用画像の取得の順序を表す指標を割当て、選択された前記心筋像位置のそれぞれに対し、強度曲線を得るステップと、
   得られた前記強度曲線に基づいて、前記複数の心筋セグメントのうちの少なくとも心筋セグメントのサブセットのうち、時空間的灌流不均一性又は灌流脱位相を示す指数を計算するステップと、
   を含む、方法。
2. 前記複数のセグメントの各セグメントにおける心筋血流の定量化を行うステップを更に含む、請求項１に記載の方法。
3. 前記複数の医用画像は、前記関心被験者に造影剤を投与した後に、前記医用イメージングモダリティによって取得される、請求項１又は２に記載の方法。
4. 前記心臓の選択された前記部分における基準場所を特定するステップを更に含み、前記指数を計算する前記ステップにおいて、前記強度曲線は、特定された前記基準場所によって決定される基準時間を参照して評価される、請求項１乃至３の何れか一項に記載の方法。
5. 選択された前記心筋位置のそれぞれについて、前記心筋像位置のそれぞれのサンプリングされた強度の特徴の発生まで、特定された前記基準場所によって決定される基準時間に対する個々の期間を自動的に決定するステップを更に含み、
   前記個々の期間は、前記指数を計算する前記ステップにおいて使用される、請求項４に記載の方法。
6. 前記指数を計算する前記ステップは、特定された前記基準場所における前記特徴の発生の時間に対する個々の前記心筋位置のそれぞれにおける特徴の発生までの時間の変動を示す統計的尺度の計算を含む、請求項４又は５に記載の方法。
7. 前記関心被験者の前記心臓の少なくとも一部の前記複数の医用画像を取得することは、前記関心被験者の前記心臓の周期動作と少なくとも部分的に同期される、請求項１乃至６の何れか一項に記載の方法。
8. 心筋像位置の強度をサンプリングする前記ステップにおいて、前記心筋像位置は、心筋と平行な方向及び前記心筋を横断する方向において選択される、請求項１乃至７の何れか一項に記載の方法。
9. 灌流図を生成し、前記灌流図をユーザに表示するステップを更に含む、請求項１乃至８の何れか一項に記載の方法。
10. 少なくとも１つの特徴位置及び／又は少なくとも１つの特徴時点を示す少なくとも１つのマーカーを前記灌流図内に組み込むステップを更に含む、請求項１乃至９の何れか一項に記載の方法。
11. 複数のコンピュータリンクを組み込むステップを更に含み、前記複数のコンピュータリンクの各コンピュータリンクは、前記灌流図内のある場所に割当てられ、前記複数のコンピュータリンクの各コンピュータリンクは、前記複数の医用画像の１つの医用画像を表すデータセットにリンク付けされる、請求項１乃至１０の何れか一項に記載の方法。
12. 関心被験者の心臓の少なくとも一部の複数の医用画像を解析することによって、心筋灌流病変の特性評価をするシステムであって、前記複数の医用画像は、医用イメージングモダリティによって連続的に取得され、当該システムは、
    前記複数の医用画像における前記関心被験者の前記心臓の選択された部分の輪郭を描写し、選択された前記部分を複数のセグメントにセグメント化する描写ユニットと、
    前記複数の医用画像から、心筋像位置の強度をサンプリングし、前記心筋像位置の対応するサンプリングされた強度に、前記複数の医用画像の各医用画像の取得の順序を表す指標を割当て、前記複数の心筋セグメントのうちの少なくとも心筋セグメントのサブセットのうち、時空間的灌流不均一性又は灌流脱位相を示す指数を計算する強度サンプリング及び解析ユニットと、
    を含む、システム。
13. 請求項１２に記載のシステムを含む、医用イメージングモダリティ。
14. 磁気共鳴イメージング装置としてデザインされる、請求項１３に記載の医用イメージングモダリティ。
15. 請求項１乃至１１の何れか一項に記載される心筋灌流病変の特性評価する方法を実行するためのソフトウェアモジュールであって、行われるべき前記方法のステップは、前記ソフトウェアモジュールのプログラムコードに変換され、前記プログラムコードは、前記医用イメージングモダリティの制御ユニットのメモリユニットに実装可能であり、前記医用イメージングモダリティの前記制御ユニットのプロセッサユニットによって実行可能である、ソフトウェアモジュール。