(57)【特許請求の範囲】
【請求項１】
管状構造内の部分ボリューム効果の存在下で前記管状構造のマルチスペクトルＸ線画像を区分化するための方法であって、
管状構造を表すボクセルを含むコントラスト強調スペクトル画像データを取得するステップであって、前記管状構造が血管である、ステップと、
取得された前記コントラスト強調スペクトル画像データに基づいてコントラストマップを生成するステップと、
前記コントラストマップ内の前記管状構造を局所化するステップと、
局所化された前記管状構造の中心線を抽出するステップであって、前記中心線が、前記局所化された前記管状構造の部分ボリュームを構成して前記中心線に沿う領域を含む、ステップと、
前記コントラストマップと、前記局所化された管状構造と、前記局所化された前記管状構造の前記抽出された中心線に沿う造影剤スペクトル応答と、対応する複数の減衰値を有する物質に隣接する複数の関心のある測定スペクトルポイントを用いてスペクトルモデルに基づいて取得されたボクセルごとの最もあり得る組織型とに基づいて、これらを用いて各ボクセル内の相対的な造影剤の比を計算し、更新済みコントラストマップを生成するステップと、
前記更新済みコントラストマップに基づいて前記管状構造を区分化するステップと、を含み、
前記中心線は、基準スペクトルポイントに基づいて生成され、前記基準スペクトルポイントは、前記管状構造の部分ボリュームがない若しくは殆どない管腔の領域に対応する異なるスペクトル画像のボクセルの１組の値を示すものであり、
前記スペクトルモデルが、複数の所定の異なる種類の物質に関するスペクトル曲線を有するマルチエネルギグラフである、方法。
【請求項２】
前記更新済みコントラストマップが、部分ボリューム領域内に部分的な造影剤を含む、請求項１に記載の方法。
【請求項３】
前記取得されたコントラスト強調スペクトル画像データに基づいてヨウ素／カルシウム分離情報を生成するステップと、
前記コントラストマップ、前記局所化された管状構造、前記抽出済みの中心線、及び前記スペクトルモデルに基づいて更新済みヨウ素／カルシウム分離情報を生成するステップと
を更に含む、請求項１又は２に記載の方法。
【請求項４】
前記取得されたコントラスト強調スペクトル画像データに基づいて仮想非造影画像を生成するステップと、
前記コントラストマップ、前記局所化された管状構造、前記抽出済みの中心線、及び前記スペクトルモデルに基づいて更新済み仮想非造影画像を生成するステップと
を更に含む、請求項３に記載の方法。
【請求項５】
停止基準が満たされているかどうかを判定し、前記停止基準が満たされていないと判定することに応答して、
前記更新済みコントラストマップ内の前記管状構造を局所化するステップと、
前記更新済みコントラストマップ内で局所化された前記管状構造から第２の中心線を抽出するステップと、
前記更新済みコントラストマップ内で局所化される前記管状構造、抽出済みの前記第２の中心線、及び前記スペクトルモデルに基づいて第２の更新済みコントラストマップを生成するステップと
を更に含む、請求項１乃至４の何れか一項に記載の方法。
【請求項６】
停止基準が満たされているかどうかを判定するステップと、
前記停止基準が満たされていると判定することに応答して、前記更新済みコントラストマップを洗練されたコントラストマップとして出力するステップと
を更に含む、請求項１乃至５の何れか一項に記載の方法。
【請求項７】
前記管状構造の管腔を区分化するために能動輪郭汎関数を適用するステップ
を更に含む、請求項１乃至６の何れか一項に記載の方法。
【請求項８】
前記更新済みコントラストマップ及び物質解析に基づいて前記管状構造の管腔の境界を推定するステップ
を更に含む、請求項１乃至６の何れか一項に記載の方法。
【請求項９】
部分ボリューム補正モデルを求めるステップと、
前記部分ボリューム補正モデルを適用し、前記更新済みコントラストマップ内の部分ボリューム効果を補正するステップと
を更に含む、請求項１乃至８の何れか一項に記載の方法。
【請求項１０】
管状構造を表すボクセルを含むコントラスト強調スペクトル画像データを取得するスペクトルアナライザと、
管状構造内の部分ボリューム効果の存在下で前記管状構造のマルチスペクトルＸ線画像を区分化するためのスペクトル解析データプロセッサと、を含み、
前記スペクトル解析データプロセッサは、
複数の所定の異なる種類の物質に関するスペクトル曲線を有するマルチエネルギグラフであるスペクトルモデルに基づいて更新済みコントラストマップを生成し、
血管である前記管状構造を前記更新済みコントラストマップに基づいて区分化し、
前記コントラストマップ内の前記管状構造を局所化し、
局所化された前記管状構造の中心線を抽出し、前記中心線が、前記局所化された前記管状構造の部分ボリュームを構成して前記中心線に沿う領域を含み、
前記コントラストマップ、前記局所化された管状構造、前記抽出された中心線に沿う造影剤スペクトル応答に基づいて前記更新済みコントラストマップを生成するように構成されており、
前記中心線は、基準スペクトルポイントに基づいて生成され、
前記基準スペクトルポイントは、撮像されたボリューム内の単一のボクセルに対応する異なるスペクトル画像の１組の値、対応する複数の減衰値を有する物質に隣接する複数の関心のある測定スペクトルポイントを用いてスペクトルモデルに基づいて取得されたボクセルごとの最もあり得る組織型、及び前記スペクトルモデルを示すものであり、これらを用いて各ボクセル内の相対的な造影剤の比が計算される、計算システム。
【請求項１１】
前記スペクトル解析データプロセッサが、
前記スペクトル解析データ内の管状構造を局所化するローカライザと、
局所化された前記管状構造の中心線を抽出するエクストラクタと、
前記スペクトル解析データ、前記局所化された管状構造、抽出済みの前記中心線、及びモデルに基づいて更新済みスペクトル解析データを生成するスペクトル解析データ更新器と、
前記更新済みスペクトル解析データに基づいて前記管状構造を区分化するセグメンタと、
所定の停止基準に基づいて前記更新済みスペクトル解析データの更なる更新を引き起こす決定ロジックと
を含む、請求項１０に記載の計算システム。
【請求項１２】
前記スペクトル解析データ更新器が、
所定の最小直径閾値を上回る直径を有する領域に対応する前記管状構造に沿って２つ以上の領域を選択する領域セレクタと、
前記被選択領域のそれぞれについてスペクトルポイントを識別するスペクトルポイントアイデンティファイヤと、
前記識別されたスペクトルポイントに基づいて基準スペクトルポイントを決定する基準スペクトルポイントデターミナと、
前記スペクトル画像データ内の関心のある測定スペクトルポイントを識別する測定スペクトルポイントアイデンティファイヤと、
エネルギスペクトルグラフ内に前記基準スペクトルポイント及び前記測定スペクトルポイントをプロットするスペクトルポイントプロッタと、
ピクセルの最もあり得る組織型に対応するスペクトル線をピクセルごとに選択するスペクトル線セレクタと、
前記基準スペクトルポイントから前記測定ポイントのそれぞれを通り、対応する被選択スペクトル線に至る前記スペクトルポイントをピクセルごとに投影するスペクトルポイントプロジェクタと、
投影される線が交差する前記被選択スペクトル線上のポイントをピクセルごとに割り当てるスペクトルポイントアサイナと、
前記基準スペクトルポイント、前記測定スペクトルポイント、及び前記割り当てられたポイントに基づいてピクセル内の造影剤の相対比を算出する相対比デターミナと
を含み、
前記スペクトル解析データ更新器が前記相対比に基づいて前記更新済みスペクトル解析データを生成する、
請求項１１に記載の計算システム。
【請求項１３】
コンピュータ可読命令で符号化されるコンピュータ可読記憶媒体であって、プロセッサによって実行されるとき、請求項１乃至９の何れか一項に記載の方法の各ステップを前記プロセッサに実行させる、コンピュータ可読記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
以下の内容は、概してスペクトル画像データ処理に関し、コンピュータ断層撮影法（ＣＴ）への特定の応用例と共に説明されている。但し、以下の内容は他の撮像様式にも適している。
【背景技術】
【０００２】
コンピュータ断層撮影法（ＣＴ）とは、非侵襲的な撮像手続きである。ＣＴスキャナは、検出器アレイに対向し、検査領域にわたる回転可能なガントリ上に設けられるｘ線管を概して含む。回転可能なガントリ、従ってｘ線管は検査領域の周りを回転する。ｘ線管は、検査領域を横断し、検出器アレイによって検出される放射線を放つように構成される。それに応答して、検出器アレイは、検出される放射線を示す投影データを生成して出力する。ボリューム画像データを生成するために投影データが再構築される。ボリューム画像データは、相対的な放射線濃度に対応するグレースケール強度値として表わされるボクセルを含む。
【０００３】
グレースケール値は、走査対象の減衰特性を概して反映し、走査対象内の解剖学的構造等の構造を示す。物質による光子の吸収率は物質を横断する光子のエネルギに依存するので、検出される放射線は、対象及び／又は目的物の走査物質の元素組成又は物質組成（例えば原子番号）を示す追加情報を提供するスペクトル情報も含む。不都合なことに非スペクトルＣＴスキャナでは、検出器アレイによって出力される投影データがエネルギスペクトルにわたって積分されるエネルギフルエンスに比例するので、スペクトル特性を反映しない。
【０００４】
スペクトル（マルチエネルギ）撮像向けに構成されるＣＴスキャナは上記のスペクトル特性を捕捉する。概して、そのように構成されるＣＴスキャナは、少なくとも２つの異なる光子エネルギ範囲にそれぞれ対応する少なくとも２つの減衰値を生成する。かかるスキャナは、異なる平均スペクトルを有する放射線を放つように構成される複数のｘ線管、走査中に少なくとも２つの異なる放出電圧（例えば８０ｋＶｐ及び１４０ｋＶｐ）を制御可能に切り替えられるように構成される単一のｘ線管、単一の広域スペクトルｘ線管、並びに（例えば光子計数検出器、スペクトル感度又は応答が異なる少なくとも２組のシンチレータ等を有する）エネルギ分解検出器及び弁別エレクトロニクスを有するエネルギ分解検出器アレイを含み得る。
【０００５】
冠動脈疾患（ＣＡＤ）は、心臓の動脈の内壁上にプラークが蓄積することに応じて引き起こされ得る。ＣＡＤに由来する合併症は、心臓への血流を減らす冠動脈の血管狭窄（狭窄）又は閉塞や、不安定プラークの破裂等を含み得る。冠動脈血管再生術を行う判定は、誘発性心筋虚血の客観的証拠を視覚的に示すことによって最も良く知らされる。胸痛に苦しむＣＡＤ患者における血管再生の現在の主要な決定要因は、冠動脈血管の血管狭窄率である。不都合なことに、冠動脈血管再生術を施すかどうかを判定するのに十分明確な情報を与えるには、狭窄率の測定は不十分であり得る。
【０００６】
血流予備量比（ＦＦＲ：fractional flow reserve）とは、通常の動脈の仮想症例における血流に対する病的動脈を通る最大、充血性、心筋血流の比率である。ＦＦＲは、心臓カテーテル検査法における狭窄後の圧力と狭窄前の圧力との比率によって測定される。文献は、多枝ＣＡＤを有する１，００５名の患者のＦＡＭＥ（Fractional Flow Reserve Versus Angiography for Multivessel Evaluation）研究において、とりわけＦＦＲ＜０．８の虚血性病変の血管再生を確保するＦＦＲガイド下の血管再生が、血管造影ガイド下方式を受ける患者に比べて著しく低い有害事象発生率及び少ない冠状動脈ステントの設置に関連したことを示している。
【０００７】
近年、マルチスライスＣＴ画像データから計算流体力学を使用してＦＦＲを推定する方法が提案された。ＣＴを検討することに基づいてＦＦＲを推定することは、医師にとって貴重な機能診断情報を非侵襲的に導出することを可能にする。心臓のＣＴデータからＦＦＲを推定するのに必要な条件は、とりわけ狭窄領域内の管腔の形状寸法を正確に区分化することである。不都合なことに、心臓のＣＴ画像内の限られた分解能、及び石灰化や他のプラーク成分等の他の生物学的物質に対する造影剤の部分ボリューム効果が管腔の区分化の質を低下させる。
【０００８】
急性冠症候群（ＡＣＳ）とは、冠動脈の閉塞に起因する一群の症状を指す。心臓専門医にとっての臨床的課題は、急性心筋梗塞等のＡＣＳを検出し予防することである。組織病理学的研究は、急性冠症候群による急死の原因は破裂性プラークがあることだと示している。従って、破裂し易いアテローム性動脈硬化症のプロファイルを明らかにする努力がなされている。剖検調査に基づき、責任病変の大半は概して非狭窄性であり、覆い被さる薄い線維性被膜を有する大きい壊死性コアを含み、血管区分の広範囲の再モデリングを引き起こす。
【０００９】
マルチスライスコンピュータ断層撮影（ＭＳＣＴ）は、動脈硬化性プラークを検出すること、並びにその組成及び機械的性質を評価することを可能にする。従って心臓のＣＴは、ＡＣＳを検出し予防する有効な手段を提供し得る。しかし、ＭＳＣＴを用いてプラークを正確に特徴付けることは依然として課題である。１つの制限は、石灰化プラーク等の超高密度物質から生じるブルーミングアーティファクトである。カルシウムブルーミングは、主にＣＴシステムの限られた空間分解能及びビーム硬化効果から生じる。このぼやけ効果は、石灰化領域の過度な表現、及びそれに対応して他のプラーク成分の過小評価を招く。
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本明細書に記載の態様は上記の問題及び他の問題に対処する。
【課題を解決するための手段】
【００１１】
以下、コントラスト強調スペクトルＣＴ画像データを使用して管状構造（例えば冠動脈内腔）の区分化を改善する手法を説明する。この手法は、構造の局所化と区分化との間の相互関係、及びスペクトル解析による定量造影剤マッピングを活用する。これらのそれぞれは、他のものの結果を過去の情報として使用することで改善された結果をもたらすことができ、この手法は１回又は複数回の反復を行ったり来たりし、区分化を最適化し得る。その結果生じる区分化は、少なくともＣＴ血流予備量比解析（ＣＴ−ＦＦＲ）及び急性冠症候群のリスク評価を改善するために使用され得る。
【００１２】
一態様では、方法は、管状構造を表すボクセルを含むコントラスト強調スペクトル画像データを取得するステップを含む。この方法は、取得されたコントラスト強調スペクトル画像データに基づいて少なくともコントラストマップを生成するステップを更に含む。この方法は、スペクトルモデルに基づいて更新済みコントラストマップを生成するステップを更に含む。この方法は、更新済みコントラストマップに基づいて管状構造を区分化するステップを更に含む。
【００１３】
別の態様では、計算システムが、コントラスト強調スペクトル画像データを受け取り、それに基づいてスペクトル解析データを生成するスペクトルアナライザを含み、スペクトル解析データはコントラストマップを含む。この計算システムは、スペクトル解析データを洗練させ、洗練されたスペクトル解析データを生成するスペクトル解析データプロセッサを更に含む。
【００１４】
別の態様では、コンピュータ可読記憶媒体がコンピュータ可読命令で符号化される。コンピュータ可読命令は、プロセッサによって実行されるとき、管腔の局所化と区分化との間の相互関係及びスペクトル解析による定量造影剤マッピングを利用する反復的更新に基づいてヨウ素マップを洗練させ、それにより洗練されたヨウ素マップを生成すること、及び関心のある構造を洗練されたヨウ素マップから区分化することをプロセッサに実行させる。
【００１５】
本発明は、様々なコンポーネント及びコンポーネントの構成、並びに様々なステップ及びステップの構成の形を取り得る。図面は好ましい実施形態を例示するためのものに過ぎず、本発明を限定するものとして解釈すべきではない。
【図面の簡単な説明】
【００１６】
【図１】画像データプロセッサに関連して撮像システムの一例を概略的に示す。
【図２】スペクトル解析データプロセッサを含む画像データプロセッサの一例を概略的に示す。
【図３】スペクトル解析データアップデータを含むスペクトル解析データプロセッサの一例を概略的に示す。
【図４】スペクトル解析データアップデータの一例を概略的に示す。
【図５】狭窄を有する血管の下位部分を示し、この中では血液のスペクトル応答曲線内の基準造影剤を決定するための基準スペクトルポイントを決定する領域が識別されている。
【図６】図５の血管の下位部分を示し、この中では血液のスペクトル応答曲線内の基準造影剤を決定するための基準スペクトルポイントが識別されている。
【図７】測定スペクトルポイントが識別されている、図５の血管の下位部分を示す。
【図８】１組のスペクトル応答曲線の一例を有するスペクトル応答グラフの一例を示す。
【図９】図７の基準及び測定スペクトル応答ポイント、並びに１組のスペクトル応答曲線を含むスペクトル応答グラフの一例を示す。
【図１０】スペクトル応答ポイントから測定スペクトル応答ポイントを通ってスペクトル応答曲線に至る線が投影されている図９の例を示す。
【図１１】本明細書に記載の実施形態による方法の一例を示す。
【発明を実施するための形態】
【００１７】
まず図１を参照し、コンピュータ断層撮影（ＣＴ）スキャナ等の撮像システム１００が示されている。図示の撮像システム１００は、スペクトル撮像向けに構成される。本明細書で使用されるとき、スペクトル撮像とは、撮像システム１００が少なくとも２つの異なる所定のエネルギ範囲にそれぞれ対応する少なくとも２つの減衰値を生成することを指す。撮像システム１００は、概して固定されたガントリ１０２及び回転するガントリ１０４を含み、回転するガントリ１０４は固定されたガントリ１０２によって回転自在に支持され、ｚ軸を中心に検査領域１０６の周りを回転する。
【００１８】
ｘ線管等の放射線源１０８が回転するガントリ１０４によって回転自在に支持され、回転するガントリ１０４と共に回転し、検査領域１０６を横断する多エネルギ放射線を放つ。一例では、放射線源１０８が単一の広域スペクトル管を含む。改変形態では、放射線源１０８が、走査中に少なくとも２つの異なる放出電圧（例えば８０ｋＶｐ、１４０ｋＶｐ等）を制御可能に切り替えられるように構成される。切り替えは、積分期間（integration period）内とすることができ、積分期間を行ったり来たりすることができ、走査間（例えば連続した走査ごとに異なるｋＶが使用されるデュアルスピンスキャン）等とすることができる。更に別の改変形態では、放射線源１０８が、異なる平均エネルギスペクトルで放射線を放つように構成される少なくとも２つの管を含む。別の改変形態は上記のものの組合せを含む。
【００１９】
放射線感知検出器アレイ１１０が、検査領域１０６にわたり放射線源１０８の反対側に角度をなした弧を張る。放射線感知検出器アレイ１１０は、ｚ軸に沿って互いに平行に配置される１つ又は複数の検出器の並びを含む。放射線感知検出器アレイ１１０は、検査領域１０６を横断する放射線を検出し、それを示すスペクトル投影データを生成する。この場合もやはり、かかる投影データは、少なくとも２つの異なる所定のエネルギ範囲及び／又は重みにそれぞれ対応する少なくとも２つの減衰値を含む。一例では、放射線感知検出器アレイ１１０が、ｘ線エネルギ感度が異なる少なくとも２つのシンチレータと、対応する光感度を有する少なくとも２つの対応する光センサとを有するマルチエネルギ積分検出器、直接変換材料（例えばＣｄＴｅ、ＣｄＺｎＴｅ等）を有する光子計数検出器、他のエネルギ分解検出器等のエネルギ分解検出器アレイを含む。改変形態では、放射線感知検出器アレイ１１０が、１０８にある任意選択的な２つの異なるｘ線管に対応する２つの異なる検出アレイを含む。
【００２０】
リコンストラクタ１１２は、放射線感知検出器アレイ１１０によって出力されるスペクトル投影データを再構築する。この再構築は、例えばエネルギ依存性の成分を再構築し、異なるエネルギ値又は範囲にそれぞれ対応する１つ又は複数のスペクトル画像を生成することを含む。リコンストラクタ１１２は、エネルギ依存性の成分を組み合わせて非スペクトル画像データを生成することもできる。診察台等の対象支持部１１４は、検査領域１０６内で対象又は目的物を支持する。オペレータコンソール１１６は、モニタ等の人間が読むことができる出力装置、及びキーボードやマウス等の入力装置を含む。コンソール１１６上にあるソフトウェアは、スペクトル撮像のためにグラフィカルユーザインタフェース（ＧＵＩ）又は他の態様によってオペレータが撮像システム１００と対話し且つ／又は撮像システム１００を操作することを可能にする。
【００２１】
計算システム１２０は、少なくとも１個のマイクロプロセッサ１２２及びコンピュータ可読記憶媒体（「メモリ」）１２４を含む。メモリ１２４は一時的記憶媒体を除き、物理メモリ及び／又は他の非一時的記憶媒体を含む。マイクロプロセッサ１２２は、メモリ１２４内に記憶される画像データプロセッサ１２６の命令を少なくとも実行する。マイクロプロセッサ１２２は、搬送波、信号、又は他の一時的媒体によって運ばれるコンピュータ可読命令も実行することができる。計算システム１２０はコンソール１１６の一部とすることができ、且つ／又は図１に示されているようにコンソール１１６とは別であり得る。
【００２２】
画像データプロセッサ１２６の命令は、少なくともコントラスト強調スペクトル画像データを処理するための命令を含む。以下でより詳細に説明されるように、一例では画像データ処理命令が、コントラスト強調スペクトル画像データ内に示される関心のある組織を区分化することに関連して、コントラスト強調スペクトル画像データを処理する命令を含む。管状構造（例えば冠動脈等の血管や、脳、頸動脈、肝臓、肺の血管等の他の血管等の管腔）の区分化に関し、画像データプロセッサ１２６は、くびれ（例えば狭窄）領域を含む形状寸法を正確に区分化できるようにする。以下で更に説明されるように、命令は、区分化を利用する他の処理のための命令も含み得る。
【００２３】
コントラスト強調スペクトル画像データは、撮像システム１００及び／又は他の撮像システムによって及び／又はそれらから生成され得る。スペクトル画像データは、画像保管通信システム（ＰＡＣＳ）、電子医療記録（ＥＭＲ）、放射線医学情報システム（ＲＩＳ）、病院情報システム（ＨＩＳ）等のデータリポジトリ１２８、及び／又は他のデータリポジトリから代わりに得られても良い。区分化及び／又は他の処理の結果が、コンソール１１６、データリポジトリ１２８、及び／又は他の記憶装置内に記憶され、別の機器１３０等に伝達され得る。
【００２４】
図２は、画像データプロセッサ１２６の一例を概略的に示す。
【００２５】
画像データプロセッサ１２６は、コントラスト強調スペクトル画像データを入力として受け取る。この例では、受け取られるコントラスト強調スペクトル画像データが、心臓のスペクトルコントラスト強調ＣＴ画像データを含む。
【００２６】
画像データプロセッサ１２６の命令は、スペクトルアナライザ２０２のための命令を含む。スペクトルアナライザ２０２は、受け取られるスペクトルコントラスト強調ＣＴスキャン画像データを処理し、１組又は複数組の初期スペクトル解析データをもたらす。かかるスペクトル解析は、全走査ボリューム又はその下位ボリュームに対して行われ得る。スペクトルアナライザ２０２は、受け取られるスペクトルコントラスト強調ＣＴスキャン画像データを処理するために知られている及び／又は他のスペクトル解析技法を使用することができる。
【００２７】
かかる技法は、定量造影剤（例えばヨウ素やガドリニウム等）マップに対応するボリューム画像データセットを生成すること、仮想非造影（ＶＮＣ：virtual non-contrasted）データ、仮想非強調データ、ヨウ素／カルシウム分離、並びに／又は他の物質解析及び／若しくは識別に対応するボリューム画像データセットを生成することを含む。技法の一例が、Thomas et al., “Automatic lumen segmentation in calcified plaques: dual-energy CT versus standard reconstructions in comparison with digital subtraction angiography,” AJR (2010)の中で論じられている。
【００２８】
改変形態では、スペクトルアナライザ２０２が省略され、様々な機器によって生成されるスペクトル解析データを画像データプロセッサ１２６が受け取る。本明細書で論じられるスペクトル解析データは、管状構造のくびれ領域内を含め、例えば心臓のＣＴ画像内の限られた分解能、石灰化や他のプラーク等の他の生物学的物質に対する造影剤の部分ボリューム効果、石灰化プラーク等の超高密度物質から生じるブルーミングアーティファクト等による画質の劣化を起こしやすい場合がある。
【００２９】
画像データプロセッサ１２６の命令は、スペクトル解析データプロセッサ２０４のための命令を含む。以下でより詳細に説明されるように、スペクトル解析データプロセッサ２０４は、スペクトル解析データを反復プロセスによって更新し洗練させる。一例ではこのプロセスは、反復ごとに、スペクトル解析データを少なくとも造影剤のスペクトル応答曲線に基づいて、及びスペクトル解析データのピクセル又はボクセルごとの最もあり得る組織型に基づいて更新すること、及びピクセルごとの更新値を与える各ピクセル内の相対的造影剤（relative contrast agent）を計算することを含む。最後の反復の更新済みスペクトル解析データは、洗練されたスペクトル解析データを与える。
【００３０】
画像データプロセッサ１２６の命令は、他の画像処理プロセッサ２０４のための命令を含む。他の画像データ処理プロセッサは、他の画像データ処理に関連して洗練されたスペクトル解析データを利用する。例えば、洗練されたスペクトル解析データが冠動脈の管腔の区分化を含む場合、他の画像データ処理アルゴリズムの例は、ＣＴ血流予備量比解析（ＣＴ−ＦＦＲ）アルゴリズム、急性冠症候群（ＡＣＳ）リスク評価アルゴリズム、及び／又は他のアルゴリズムを含み得る。
【００３１】
洗練されたスペクトル解析データは、ＣＴ血流予備量比解析（ＣＴ−ＦＦＲ）及び／又は急性冠症候群（ＡＣＳ）リスク評価を改善することができる。例えば洗練された区分化は、管腔の区分化の精度を改善することができ、そのことはＣＴ−ＦＦＲの精度を改善し得る。別の例では、洗練された区分化が、従来の撮像に関連するぼやけ効果を最小限にし、ＡＣＳリスクを推定するためのプラーク成分のより正確な定量化及び特徴付けを可能にすることができる。
【００３２】
図３は、スペクトル解析データプロセッサ２０４の一例を概略的に示す。
【００３３】
説明目的で、以下の内容は冠動脈の区分化に関連して説明される。しかし、この例は限定的ではなく、他の構造も本明細書で考えられることが理解されるべきである。例えば以下の内容は、解剖学的（例えば血管等）及び／又は非解剖学的な任意の管状構造と共に利用され得る。
【００３４】
スペクトル解析データプロセッサ２０４は、コントラスト強調スペクトル画像データ、及びスペクトルアナライザ２０２によって生成される初期スペクトル解析データを入力として受け取る。この例では、初期スペクトル解析データは、造影剤がヨウ素である場合はヨウ素マップ等の造影剤マップを少なくとも含む。改変形態では、造影剤マップがガドリニウム、ビスマス、金等の高減衰物質を含む他の造影剤に基づき得る。
【００３５】
スペクトル解析データプロセッサ２０４は、ローカライザ３０２を含む。ローカライザ３０２は、（最初の反復の）初期スペクトル解析データ及び（その後の反復ごとの）更新済みスペクトル解析データを入力として受け取る。ローカライザ３０２は、（初期又は更新済み）スペクトル解析データ内の冠動脈及び／又は他の構造を局所化する。知られている手法、並びに／又は他の（ユーザ対話なしの）自動及び／若しくは（ユーザ対話ありの）手動の手法が使用され得る。
【００３６】
スペクトル解析データプロセッサ２０４は、エクストラクタ３０４を更に含む。エクストラクタ３０４は、局所化された冠動脈を有するスペクトル解析データを入力として受け取る。エクストラクタ３０４は、スペクトル解析データ内の冠動脈及び／又は他の構造ごとに中心線を抽出する。知られている手法、並びに／又は他の（ユーザ対話なしの）自動及び／若しくは（ユーザ対話ありの）手動の手法が使用され得る。
【００３７】
スペクトル解析データプロセッサ２０４は、スペクトル解析データアップデータ３０６を更に含む。スペクトル解析データアップデータ３０６は、スペクトル解析データ、局所化された冠動脈、及び抽出済みの中心線を入力として受け取る。スペクトル解析データアップデータ３０６は、スペクトル解析データを更新する。上記で簡潔に論じられ、以下でより詳細に論じられるように、かかる更新は、中心線に沿った造影剤のスペクトル応答に基づいて、及びボクセルごとの最もあり得る組織型に基づいてスペクトル解析データを更新することを含み、各ボクセル内の相対的造影剤を計算するために使用される。
【００３８】
スペクトル解析データプロセッサ２０４は、セグメンタ３０８を更に含む。セグメンタ３０８は、更新済みスペクトル解析データを入力として受け取る。セグメンタ３０８は、更新済みスペクトル解析データから関心のある構造を区分化する。知られている及び／又は他の区分化手法が使用され得る。或る手法では、スペクトル解析データ上に１組の輪郭が適用され、推定される中心線を取り囲む面として初期設定される。輪郭は、Chan-Vese最小分散及び測地面モデルに基づき得る。
【００３９】
上記の手法は、Holtzman-Gazit et al., “Segmentation of thin structures in volumetric medical images,” Image Processing, IEEE Transactions on 15.2 (2006): 354-363の中でより詳細に論じられている。別の手法が、Lesage et al., “A review of 3D vessel lumen segmentation techniques: Models, features and extraction schemes,”Medical image analysis 13.6 (2009): 819-845の中で説明されている。他の手法も本明細書で考えられる。これらの手法から、冠動脈の管腔の正確な区分化が実現され得る。
【００４０】
冠動脈の区分化では、セグメンタ３０８が、冠動脈の管腔を区分化するために上記の及び／又は他の手法を利用することができる。セグメンタ３０８は、例えば動脈壁のプラーク及び／又は他の組織形状を識別し、区分化することもできる。これは、造影剤の強調がコントラスト強調画像データから除去される仮想非造影（ＶＮＣ）画像データを使用して実現され得る。
【００４１】
一実施形態では、セグメンタ３０８が管腔の境界の推定を改善することもできる。例えば、スペクトル解析画像データ及び区分化は、動脈内の各ボクセル内の造影剤の相対量のマッピングを提供する。相対量は、部分ボリューム効果及び限られた空間分解能により、動脈壁の近くで１００％又は０％と異なり得る。より鮮明な壁及び一層改善された管腔の区分化を推定するために幾つかのモデルが使用され得る。
【００４２】
或る手法は、知られている管の不鮮明な形状についての過去の知識を（通常は較正によって）使用する。この手法は、肺気道の管腔及び壁の測定に使用されている。この手法の一例が、Saba et al., “Maximizing quantitative accuracy of lung airway lumen and wall measure obtained from X-ray CT imaging,” J. Appl. Physiol. (2003)の中で論じられている。別の技法は、撮像プロトコル及びパラメータに対応する既知の３Ｄ点広がり関数を用いた３Ｄデコンボリューションを適用する。
【００４３】
スペクトル解析データプロセッサ２０４は、決定ロジック３１０を更に含む。決定ロジック３１０は、停止基準に基づいて別の反復が行われるべきかどうかを判定する。例えば、決定ロジック３１０は、更新済みスペクトル解析データが停止基準を満たしていないことに応答して別の反復を引き起こす。さもなければそれ以上の反復は行われず、更新済みスペクトル解析データが洗練されたスペクトル解析データとして出力される。
【００４４】
停止基準は、経験的な反復回数を選択し、連続した反復間の変化を十分小さくなるまで確認し、全体的な性能指数を既定の閾値と比較することによって、及び／又は他の手法によって決定され得る。他の臨床応用では、この手法が、脳、頸動脈、肝臓、肺等の他の臓器内の狭い血管の区分化を改善するために拡張され得る。
【００４５】
図４は、スペクトル解析データアップデータ３０６の一例を概略的に示す。
【００４６】
説明目的で、以下の内容は図３で論じられた冠動脈の区分化に関連して説明される。同様に、この例も限定的ではなく、他の組織も本明細書で考えられることが理解されるべきである。
【００４７】
スペクトル解析データアップデータ３０６は、スペクトル解析データ、局所化された冠動脈、及び抽出済みの中心線を入力として受け取る。この例では、図２と同様にスペクトル解析データがヨウ素マップ等の造影剤マップを含む。
【００４８】
スペクトル解析データアップデータ３０６は、領域セレクタ４０２を含む。領域セレクタ４０２は、自動及び／又は手動の手法を使用し、スペクトル解析データ内の冠動脈に沿って２つ以上の領域を選択する。概して、被選択領域は、比較的広幅の断面を有する管状構造の領域に対応する。この領域は、所定の最小直径閾値に基づいて及び／又は他の態様で決定され得る。
【００４９】
この例では、それらの領域が、管腔の中心に部分ボリュームがない若しくは殆どない、又は空間分解能の壁との重なりがない若しくは殆どない動脈部分を与える。これらの領域内では、入力されるスペクトル解析画像データを用いて造影剤の定量化が正確に求められ得る。図５は、２つの領域５０２及び５０４の選択を示す。２つの領域５０２及び５０４は、狭窄５０６の前の領域及び狭窄５０６の後の領域にそれぞれ対応する。別の実施形態では３つ以上の領域が選択される。更に領域は、領域５０２及び５０４を含むことができ、且つ／又はそれらの領域と異なることができる。
【００５０】
図４を参照し、スペクトル解析データアップデータ３０６は、スペクトルポイントアイデンティファイヤ４０４を更に含む。スペクトルポイントという用語は、撮像されたボリューム内の単一のピクセル又はボクセルに対応する様々なスペクトル画像の１組の値を指す。例えばスペクトルポイントは、以下で論じられる図８〜図１０内のグラフ内の点として表され得る。スペクトルポイントアイデンティファイヤ４０４は、被選択領域５０２及び５０４について、被選択領域５０２及び５０４それぞれのスペクトルポイントを識別する。図６は、図５の２つの領域５０２及び５０４に関する２つの点（ＡＬ）６０２及び（ＡＵ）６０４をそれぞれ示す。
【００５１】
図４に戻り、スペクトル解析データアップデータ３０６は、基準スペクトルポイントデターミナ４０６を更に含む。基準スペクトルポイントデターミナ４０６は、図５のスペクトル解析データ内で識別されるスペクトルポイント６０２及び６０４から基準スペクトルポイントを決定する。一例では基準スペクトルポイントデターミナ４０６が、補間スペクトル応答線（例えば（ＡＵ＋ＡＬ）／２）を中心線に沿って生成する。これは部分ボリューム及び壁との重なりを有し得る狭領域を含み、かかる領域では、コントラストの定量化を局所データだけでは正確に求めることができない可能性がある。
【００５２】
スペクトル解析データアップデータ３０６は、測定スペクトルポイントアイデンティファイヤ４０８を更に含む。測定スペクトルポイントアイデンティファイヤ４０８は、スペクトル画像データ内の関心のある測定スペクトルポイントを識別する。図７は、関心のある測定スペクトルポイントの複数の例（Ｂ１）７０２、（Ｂ２）７０４、及び（Ｂ３）７０６を示す。図７では、関心のある測定スペクトルポイント（Ｂ１）７０２が、高減衰値を有する物質（図７ではハウンスフィールド単位に関して示されている）に隣接する。関心のある測定スペクトルポイント（Ｂ２）７０４は、中程度の減衰値を有する物質に隣接する。関心のある測定スペクトルポイント（Ｂ３）７０６は、低減衰値を有する物質に隣接する。例えば高減衰値は、カルシウムによって、又は血管腔内の造影剤とは異なるプラークを標的とする造影剤によって引き起こされる場合があり、中程度の減衰値は柔組織によって引き起こされる場合があり、低減衰値は脂質によって引き起こされる場合がある。
【００５３】
スペクトル解析データアップデータ３０６は、スペクトルポイントプロッタ４１０を更に含む。スペクトルポイントプロッタ４１０は、エネルギスペクトルグラフ上に基準スペクトルポイント及び測定スペクトルポイントをプロットする。デュアルエネルギ低エネルギ−高エネルギグラフの一例が図８に示されている。図８では、グラフ８００が、スペクトル画像データの第１の高エネルギ画像内のピクセルの減衰を表す第１の軸８０２と、スペクトル画像データの第２の低エネルギ画像内の同一ピクセルの減衰を表す第２の軸８０４とを含む。この例では、低エネルギ軸及び高エネルギ軸がハウンスフィールド単位（ＨＵ）を単位とする。別の例では、低エネルギ軸及び高エネルギ軸が他の単位に関することができ、光電効果とコンプトン効果との分解、仮想単色画像、物質ベースの表現等の値を表し得る。
【００５４】
図８では、第１のスペクトル曲線８０６が血液−水曲線を表し、第２のスペクトル曲線８０８がカルシウム−水曲線を表し、第３のスペクトル曲線８１０がカルシウム／線維性−水曲線を表し、第４のスペクトル曲線８１２が線維性組織−水曲線を表し、第５の曲線８１４が脂肪、脂質−水スペクトル曲線を表し、第６の曲線８１６がカルシウム／脂質−水曲線を表す。このグラフの特定の曲線は異なることができ、関心のある組織を表し得る。グラフ８００は２つのエネルギに関するが、３つ以上のエネルギが検討されても良いことが理解されるべきである。例えば、Ｎのエネルギ（Ｎは正の整数）が検討され得る。かかる例では、グラフ８００はＮ次元（Ｎ−Ｄ）グラフを含み得る。
【００５５】
図９は、図７のスペクトルポイントに関連したデュアル（低及び高）エネルギグラフ９００を示す。図９では、第１の軸９０２が高エネルギ画像内のピクセルの減衰を表し、第２の軸９０４が低エネルギ画像内の同一ピクセルの減衰を表す。図９では、第１のスペクトル曲線９０６がヨウ素−水曲線を表し、第２のスペクトル曲線９０８がカルシウム−水曲線を表し、第３のスペクトル曲線９１０が血液−水曲線を表し、第４のスペクトル曲線９１２が柔組織−水曲線を表し、第５の曲線９１４が脂肪−水スペクトル曲線を表す。ヨウ素−血液スペクトルポイント（Ａ）９１６は、基準スペクトルポイントを表す。スペクトルポイントＢ１、Ｂ２、及びＢ３は、図７の関心のある測定スペクトルポイント（Ｂ１）７０２、（Ｂ２）７０４、及び（Ｂ３）７０６に対応する。
【００５６】
図４に戻り、スペクトル解析データアップデータ３０６は、スペクトル線セレクタ４１２を更に含む。スペクトル線セレクタ４１２は、ピクセルの最もあり得る組織型に対応するスペクトル線をピクセルごとに選択する。これは、一例では仮想非造影画像及びモデルに基づいて行われ得る。仮想非造影画像はデータを提供し、かかるデータ内では、特定の領域内でカルシウムが関与しているかどうかを自動で判定するのがより容易である。
【００５７】
スペクトル解析データアップデータ３０６は、スペクトルポイントプロジェクタ４１４を更に含む。スペクトルポイントプロジェクタ４１４は、スペクトルポイントをピクセルごとに投影する。一例では、そのためにスペクトルポイントプロジェクタ４１４が、基準スペクトルポイント（Ａ）９１６から測定ポイント（Ｂ１）７０２、（Ｂ２）７０４、及び（Ｂ３）７０６のそれぞれを通り、対応する被選択スペクトル曲線に至る線を投影する。これが図１０に示されており、図１０は、基準スペクトルポイント（Ａ）９１６から測定ポイント（Ｂ１）７０２、（Ｂ２）７０４、及び（Ｂ３）７０６を通り被選択曲線９０８、９１２、及び９１４まで投影される線１００２、１００４、１００６を示す。
【００５８】
図４に戻り、スペクトル解析データアップデータ３０６は、スペクトルポイントアサイナ４１６を更に含む。スペクトルポイントアサイナ４１６は、投影される線が交差する被選択曲線上のポイントをピクセルごとに割り当てる。これも図１０に示されており、図１０ではスペクトルポイントアサイナ４１６が、基準スペクトルポイント（Ａ）９１６から測定ポイント（Ｂ１）７０２、（Ｂ２）７０４、及び（Ｂ３）７０６を通り被選択曲線９０８、９１２、及び９１４に至る投影１００２、１００４、１００６に基づき、ポイント（Ｃ１）１００８、（Ｃ２）１０１０、及び（Ｃ３）１０１２を割り当てている。この例では、スペクトルポイント（Ｂ１）７０２がカルシウム（骨）−水の線９０８上に投影されている。別の例では、スペクトルポイント（Ｂ１）７０２が血液−水の線９１０、柔組織−水の線９１２、及び／又は他のスペクトル応答線上に投影され得る。
【００５９】
図４に戻り、スペクトル解析データアップデータ３０６は、相対比デターミナ４１８を更に含む。相対比デターミナ４１８は、（例えば仮定上の血中造影剤値に対する）ボクセル内の造影剤の相対比Ｒを計算する。一例では、この比率が絶対距離の比率であるＲ＝｜Ｂ−Ｃ｜／｜Ａ−Ｃ｜によって算出され、Ａ、Ｂ、及びＣはベクトルであり、それぞれＮの要素を含み、Ｎは解析データ内の異なるスペクトル画像の数である。スペクトル解析データアップデータ３０６は、更新済みスペクトル解析画像データとして相対的な血中造影剤値を出力する。
【００６０】
図３に関して論じられたように、決定ロジック３１０は、別の反復が行われるか、又は洗練されたスペクトル解析データとして更新済みスペクトル解析データが出力されるかを決定する。図２に関して論じられたように、洗練されたスペクトル解析データは、ＣＴ−ＦＦＲやＡＣＳリスク評価等の他の処理に関連して利用され得る。
【００６１】
図１１は、本明細書に記載の実施形態による方法の一例を示す。
【００６２】
本明細書に記載の方法の中の行為の順序は限定的ではないことを理解すべきである。そのため、他の順序も本明細書で予期される。加えて、１つ若しくは複数の行為が省略され且つ／又は１つ若しくは複数の追加の行為が含められても良い。
【００６３】
１１０２で、ボクセルが管状を表す画像データを含む、コントラスト強調スペクトル画像データが取得される。
【００６４】
１１０４で、コントラスト強調スペクトル画像データから少なくともコントラストマップが生成される。このステップはモデルを使用して実現され得る。ヨウ素／カルシウム分離やＶＮＣ画像等の他の情報も生成され得る。
【００６５】
１１０６で、管状構造が局所化され、局所化された管状構造の中心線が抽出される。最初の反復では、このステップは初期ヨウ素マップに基づく。その後の反復では、このステップは直近の更新済みヨウ素マップに基づく。
【００６６】
１１０８で、局所化された管状構造、抽出済みの中心線、及びマルチエネルギグラフモデルに基づいてコントラストマップが更新される。更新済みコントラストマップは部分ボリューム領域内の部分的なヨウ素（partial iodine）を識別するために、物質蓄積及び管状構造壁の解析を行うために、及びＶＮＣ画像を洗練させるために使用され得る。
【００６７】
１１１０で、更新済みコントラストマップ及び管状構造区分化モデルに基づいて管状構造が区分化される。このステップは、洗練されたコントラストマップ上に適用される能動輪郭汎関数を使用した管腔の区分化、洗練されたコントラストマップ及び物質解析に基づく管腔境界推定の改善、部分ボリューム効果及び限られた分解能効果のための補正モデル、区分化のためのプラーク領域の局所化等を含み得る。
【００６８】
１１１２で、停止基準が満たされているかどうかが判定される。停止基準の例は、経験的な反復回数、連続した反復の結果の差異、全体的な性能指数との比較、及び／又は他の手法を含む。
【００６９】
停止基準が満たされていない場合、行為１１０６から１１１０が繰り返される。しかし、停止基準が満たされている場合、１１１４で、洗練されたコントラストマップが出力される。本明細書で論じられたように、洗練されたコントラストマップは他の画像処理に関連して使用され得る。
【００７０】
上記の内容は、コンピュータプロセッサによって実行されるときに記載された行為をそのプロセッサに実行させる、コンピュータ可読記憶媒体上に符号化され又は埋め込まれるコンピュータ可読命令によって実施され得る。加えて、又は或いは、コンピュータ可読命令の少なくとも１つが信号、搬送波、又は他の一時的媒体によって運ばれる。
【００７１】
本発明が好ましい実施形態に関して説明されてきた。上記の詳細な説明を読んで理解するとき、他者は修正形態及び改変形態に気付くことがある。添付の特許請求の範囲又はその均等物の範囲内に含まれる限り、本発明はそのような全ての修正形態及び改変形態を含むものとして構成されることを意図する。
【図１】
0006559678.tif 000002
【図２】
0006559678.tif 000003
【図３】
0006559678.tif 000004
【図４】
0006559678.tif 000005
【図５】
0006559678.tif 000006
【図６】
0006559678.tif 000007
【図７】
0006559678.tif 000008
【図８】
0006559678.tif 000009
【図９】
0006559678.tif 000010
【図１０】
0006559678.tif 000011
【図１１】
0006559678.tif 000012