JP6684909B2 - 造影剤濃度マップを生成する方法 - Google Patents

造影剤濃度マップを生成する方法 Download PDF

Info

Publication number
JP6684909B2
JP6684909B2 JP2018530859A JP2018530859A JP6684909B2 JP 6684909 B2 JP6684909 B2 JP 6684909B2 JP 2018530859 A JP2018530859 A JP 2018530859A JP 2018530859 A JP2018530859 A JP 2018530859A JP 6684909 B2 JP6684909 B2 JP 6684909B2
Authority
JP
Japan
Prior art keywords
contrast
computed tomography
agent concentration
contrast agent
map
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2018530859A
Other languages
English (en)
Other versions
JP2018537226A (ja
JP2018537226A5 (ja
Inventor
レイズ カルミ
レイズ カルミ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2018537226A publication Critical patent/JP2018537226A/ja
Publication of JP2018537226A5 publication Critical patent/JP2018537226A5/ja
Application granted granted Critical
Publication of JP6684909B2 publication Critical patent/JP6684909B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/484Diagnostic techniques involving phase contrast X-ray imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5223Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/37Determination of transform parameters for the alignment of images, i.e. image registration using transform domain methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Description

本発明は、造影剤濃度マップを生成する方法に関する。本発明は、医療用撮像装置と通信するように構成された装置、コンピュータ断層撮影スキャナ及びコンピュータ読取可能な記憶媒体にも関する。
コンピュータ断層撮影応用分野において、ヨード造影剤に基づく“ヨードマップ”又はガドリニウム、バリウム、ビスマス若しくは金等の余り一般的でない薬剤に基づく他のマップ等の造影剤濃度マップは、多くの疾病及び臨床状態の診断精度を向上させる際に非常に有用であり得る。このようなヨードマップを生成するための2つの基本的アプローチが知られており、各アプローチは賛否両論を持つ。第1のアプローチは、造影剤を投与し又は投与しない2つの通常のコンピュータ断層撮影スキャンを利用する。これらスキャンは、造影前(pre-contrast)及び造影後(post-contrast)コンピュータ断層撮影スキャンとも呼ばれる。画像減算が後続する体積的空間位置合わせを適用することにより、造影剤マップが生成される(即ち、造影後画像−造影前画像である)。第2のアプローチは、単一のスペクトルコンピュータ断層撮影スキャンから造影剤マップを直接生成するものである。デュアルエネルギ又は光子計数ベース等のスペクトルコンピュータ断層撮影は、造影剤を生体物質から定量的に区別する能力を有している。
スペクトルコンピュータ断層撮影技術に関する位置合わせ減算法の幾つかの利点が存在する。先ず第1に、同じ放射線量に対するコントラスト対雑音比(CNR)は、関連するデュアルエネルギヨードマップにおけるよりも約3倍高い。この一層高いCNRは、例えば、血液又は軟組織HUに対するヨード強調の間のハウンズフィールド単位(HU)差分が、デュアルレイヤコンピュータ断層撮影(120kVpで動作される)における高エネルギ第2レイヤ画像に対する低エネルギ第1レイヤ画像におけるヨード強調の間の差分よりも約3倍高いという事実から理解することができる。この比較は、2つの技術の間の全体の線量の等化が、従来のコンピュータ断層撮影画像(2スキャン減算プロトコルにおける)及びシングルレイヤ画像(デュアルエネルギシングルスキャンプロトコルにおける)の各々のオリジナルの信号対雑音をおおよそ同じレベルにするという考慮を含んでいる。
更に、2スキャンプロトコルにおいては、臨床診断のために真の非造影画像が利用可能である。この画像は、特に低線量コンピュータ断層撮影スキャンにおいては、デュアルエネルギ仮想非造影画像より通常は大幅に良好である。
加えて、位置合わせ減算技術においては、骨及びカルシウム識別及び削除を、時には、シングルスキャンスペクトルコンピュータ断層撮影におけるよりも大幅に正確に実施することができる。このことは、造影剤濃度が、骨及びカルシウムが同じに留まるのに反して、造影前及び造影後スキャンの間で大きく変化することによるものである。
結局のところ、デュアルエネルギヨードマップ及び仮想非造影(VNC)画像は、特に低線量スペクトルコンピュータ断層撮影スキャンにおいては著しい不正確さ及びアーチファクトを被る傾向がある。
上記2スキャン技術の欠点は、非常に正確な体積的空間位置合わせが必要とされることである。相当に優秀な位置合わせアルゴリズムが市場及び研究において既に利用可能であるが、位置合わせ結果は依然として多くの場合において十分に正確ではなく、減算結果に誤位置合わせアーチファクトが存在する。
加えて、2つのスキャンを実行することは、スペクトルコンピュータ断層撮影で単一の造影剤スキャンのみを実行することより、幾らか一層複雑で時間の掛かる臨床作業及び患者計画を必要とする。
両技術の限界は、ヨードマップの外見を改善するために、両技術において、取得されたヨードマップからノイズ及びアーチファクトを低減するために相対的に強いフィルタを適用することが普通であり、これらの技術は、時には、当該マップの空間解像度及び相対的に低いヨード濃度の小さな特徴構造及び構造の検出可能性を大幅に悪化させる。
本発明の目的は、既存の各方法を、これらの各欠点を克服するように組み合わせることである。
本発明は、非造影(非コントラスト強調)コンピュータ断層撮影スキャン、造影(コントラスト強調)コンピュータ断層撮影スキャン及び対応するスペクトルコンピュータ断層撮影データから造影剤濃度マップを生成する方法に関するもので、該方法は、
a.前記非造影コンピュータ断層撮影スキャン、前記造影コンピュータ断層撮影スキャン及び前記対応するスペクトルコンピュータ断層撮影データから少なくとも2つの異なる一次造影剤濃度マップを生成するステップと、
b.各一次造影剤濃度マップの局部的品質分析を実行するステップと、
c.各一次造影剤濃度マップに対する局部的な体積的重みを前記局部的品質分析に基づいて決定するステップと、
d.二次造影剤濃度マップを前記2つの一次造影剤濃度マップ及びこれらマップの対応する局部的な体積的重みに基づいて生成するステップと、
を有する。
好ましい実施態様において、少なくとも前記造影スキャンは、通常の造影コンピュータ断層撮影画像及びスペクトルコンピュータ断層撮影データの両方を発生させるためにスペクトルコンピュータ断層撮影モードで実行される。
“局部的分析”とは、単一ボクセルの寸法を含む如何なる寸法をも有することが可能な幾つかのボクセル群に対して実行される分析と理解されるべきである。分析されたボクセル群は、全体の造影剤濃度マップの区画を形成するか、又は形成しないものとすることができる。即ち、局部的分析は全体のマップ又は該マップの一部のみに関するものであり得る。当該分析は、“局部的な体積的重み”を対応するボクセル群に割り当てることができるように実行される。該対応するボクセル群は、当該“局部的分析”を受けたものと同一のボクセル群又は別の組のボクセル群とすることができる。最後のステップは、典型的に、両一次マップの加重平均である二次マップを生成し、該加重平均演算が算出された重みを考慮に入れることにある。
好ましい実施態様において、前記一次造影剤濃度マップのうちの少なくとも1つは、前記非造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンに対して体積的空間位置合わせステップ及び画像減算ステップを実行することにより得られる。当該位置合わせアルゴリズムは、弾性、アファイン又は剛性位置合わせに関する技術を利用することができる。この場合、当該方法は、前記体積的空間位置合わせステップを、前記造影コンピュータ断層撮影スキャン及び/又は前記スペクトルコンピュータ断層撮影データから仮想非造影画像ボリュームを計算し、前記非造影コンピュータ断層撮影スキャン画像ボリュームを前記仮想非造影画像ボリュームに変形することを可能にする変形関数を計算し、該変形関数を前記造影コンピュータ断層撮影スキャン画像ボリュームに適用することにより改善するステップを更に有することができる。確かなことに、空間位置合わせ精度は造影剤濃度マップを生成する位置合わせ−減算方法の重大な制約であり、利用可能なスペクトルデータから賢い方法で利益を得ることは興味深いことである。前記位置合わせステップを改善する他の方法は、
a.体積的空間位置合わせステップ及び画像減算ステップを実行することにより得られる一次造影剤濃度マップを、他の一次造影剤濃度マップからの情報を用いて改善するステップと、
b.該造影剤濃度マップを前記造影コンピュータ断層撮影スキャン画像ボリュームから減算して、改変された非造影コンピュータ断層撮影画像ボリュームを得るステップと、
c.前記改変された非造影コンピュータ断層撮影画像ボリューム及び前記造影コンピュータ断層撮影スキャン画像ボリュームに対して体積的空間位置合わせステップ及び画像減算ステップを実行することにより新たな造影剤濃度マップを生成するステップと、
からなるステップを更に追加することである。
好ましい実施態様において、前記一次造影剤濃度マップのうちの少なくとも1つは、前記スペクトルコンピュータ断層撮影データのみから得られる。
この場合、前記スペクトルコンピュータ断層撮影データのみから得られた一次造影剤濃度マップの前記局部的品質分析は、局部的放射線量レベルを考慮するステップを有することができる。例えば、該局部的放射線量レベルは、スキャンパラメータ及びツールから、好ましくは計算された線量マップから得ることができる。
典型的に、前記非造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンの少なくとも一方は低線量スキャンとすることができる。実際のところ、一次スキャンは、通常、粗い分析のためにのみ必要とされ、従って一層低い線量ほど良い。
前記局部的品質分析及び前記対応する局部的な体積的重み付けは、前記一次造影剤濃度マップの各ボクセルに対して実行することができる。これは、実際に、分析されるボクセル群が単一のボクセルの寸法を持つ前述したケースに相当する。
各一次造影剤濃度マップの尺度又は単位は、前記二次造影剤濃度マップの生成の前に適応化及び/又は正規化することができる。
前記一次造影剤濃度マップのうちの少なくとも1つの前記局部的品質分析は、局部的画像ノイズ推定、一群の隣接するピクセルに対する標準偏差測定、微細構造若しくは形状解析、局部的画像アーチファクト分析、空間解像度分析、局部的空間周波数若しくは波長解析、又は予め定められた限界外のマップ値若しくは値勾配の分析等の技術を有することができる。
前記一次造影剤濃度マップの各々における前記局部的な体積的重みの間の滑らかな移行を形成することができる。このことは、最終的マップの可読性を向上させ得るが、幾つかの血管の壁等の幾つかの不連続性をぼかし得る。かくして、このような機能を、関連する平滑化の量を加重することができる分割アルゴリズムと結合することが興味深い。
本発明は、医療用撮像装置と通信するように構成された装置にも関するものであり、該装置は、
a.非造影コンピュータ断層撮影スキャン、造影コンピュータ断層撮影スキャン及び対応するスペクトルコンピュータ断層撮影データから少なくとも2つの異なる一次造影剤濃度マップを生成し;各一次造影剤濃度マップの局部的品質分析を実行し;各一次造影剤濃度マップに対する局部的な体積的重みを前記局部的品質分析に基づいて決定し;二次造影剤濃度マップを前記2つの一次造影剤濃度マップ及びこれらマップの対応する局部的な体積的重みに基づいて生成するように構成されたプロセッサと、
b.前記二次造影剤濃度マップを表示する手段と、
を有する。
本発明は、本発明による方法を実施するように構成された装置を有するコンピュータ断層撮影スキャナにも関するものである。
本発明は、プロセッサにより実行された場合に該プロセッサに本発明による方法を実施させるコンピュータ読取可能な命令によりコード化されたコンピュータ読取可能な記憶媒体にも関するものである。
本発明の更なる態様及び実施態様は、当業者であれば、後述する詳細な説明を精読及び理解することにより認識されるであろう。数々の付加的利点及び利益は、当業者によれば、好ましい実施態様の後述する詳細な説明を精読すれば明らかとなるであろう。
本発明は、本発明の実施態様の後述する詳細な説明を精読することにより、且つ、添付図面を精査することにより一層良好に理解されるであろう。
本発明は、種々の構成要素及び構成要素の配置で、並びに種々の処理動作及び処理動作の構成で具現化し得る。図面は、好ましい実施態様を図示する目的のためだけのものであり、本発明を限定するものと見なされるべきではない。より良く描くために、特定のフィーチャは削除され得るか又は寸法は実寸通りではない場合がある。
図1は、典型的又はスペクトル型コンピュータ断層撮影スキャナである。 図2は、本発明の好ましい実施態様による方法の概略フローチャートを示す。 図3は、図2の一次造影剤濃度マップの画像強度を示す。 図4は、図2の方法により取得された一次造影剤濃度マップの臨床例を示す。 図5は、図2の方法の一次造影剤濃度の1つを得るために空間位置合わせ精度を改善する1つの方法を示すフローチャートである。 図6は、図2の方法の一次造影剤濃度の1つを得るために空間位置合わせ精度を改善する他の方法を示すフローチャートである。 図7は、図2による方法の組み合わせステップを更に詳細に示すフローチャートである。 図8は、図2による方法の組み合わせステップを更に詳細に示すフローチャートである。
図1は、コンピュータ断層撮影(CT)スキャナ等の例示的撮像システム100を概略的に示している。撮像システム100は、回転ガントリ102及び静止ガントリ104を含んでいる。回転ガントリ102は静止ガントリ104により回転可能に支持されている。回転ガントリ102は、長軸又はz軸を中心として検査領域106の周りに回転するよう構成される。撮像システム100は更に被検者支持部107を有し、該支持部107は走査の前、中及び/又は後において被検者又は被写体を検査領域106内で支持する。被検者支持部107は当該被検者又は被写体を検査領域106に装填し及び/又は検査領域106からとりだすために使用することもできる。撮像システム100はX線管等の放射源112を更に含み、該放射源は回転ガントリ102により回転可能に支持される。放射源112は、回転ガントリ102と共に検査領域106の周りを回転すると共に、検査領域106を横切る放射を発生及び放出するよう構成される。一例において、放射源112は積分期間及び/又はそれ以外内において2以上の放出電圧(例えば、80及び140kVp、80、100及び120kVp等)の間で放出電圧を切り換えるように構成される。変形例において、撮像システム100は異なる放出電圧において放射を放出する複数の放射源112を含む。他の変形例において、放射源112は単一の広域スペクトルX線管を含む。撮像システム100は、更に、放射源コントローラ114を含むことができる。この場合、該放射源コントローラ114は、発生される放射束を変調するように構成される。例えば、放射源コントローラ114は、当該放射束を変調するために、選択的に放射源112のカソード加熱電流を変化させ、放射源112の電子流を阻止するために電荷を供給し、放出された放射をフィルタし、等々を行うことができる。図示された例において、放射源コントローラ114は当該放射束を予め定められた変調パターンに基づいて変調する。
撮像システム100は、放射感知検出器ピクセル116の一次元又は二次元アレイ115を更に含んでいる。ピクセル116は、検査領域106を跨いで放射源112と対向して配置され、検査領域106を横切った放射を検出し、該放射を示す電気信号(投影データ)を発生する。一実施態様において、ピクセル116は直接変換光子計数検出器ピクセルを含む。このようなピクセルによれば、発生される信号は、検出された光子のエネルギを示すピーク振幅又はピーク高を有する電流又は電圧を含む。上記直接変換光子計数検出器ピクセルは、CdTe,CdZnTe,Si,Ge,GaAs又は他の直接変換物質等の如何なる好適な直接変換物質を含むこともできる。
他の実施態様において、検出器アレイ115は、放射源112に対して検査領域106の反対側で或る角度の弧にわたって延びる。該検出器アレイ115は、検査領域106を横切る放射を検出し、該放射を示す信号(投影データ)を発生する。当該放射源電圧が少なくとも2つの放出電圧の間で切り換えられ、及び/又は2以上のX線管が2つの異なる放出電圧において放射を放出する場合、検出器アレイ115は、これら放射源電圧の各々に関する信号を発生する。単一の広域スペクトルX線管の場合、検出器アレイ115は、斯かる信号を生成するエネルギ分解検出器(例えば、多層シンチレータ/フォトダイオード、直接変換光子計数、等)を含む。
再構成プロセッサは、上記信号を、再構成アルゴリズムメモリ又は何れかに記憶された1以上のスペクトルベースの分解アルゴリズムにより再構成する。上記アルゴリズムのうちの少なくとも1つを使用する該再構成プロセッサは、異なる画像ベースに対応する2組以上の体積的画像データを含むスペクトル体積的画像データを生成する。例えば、デュアルエネルギの場合、これらのベースは、光電/コンプトン散乱対、水/ヨウ素対(又は他の物質ベースの対)、2つの異なる有効keVX線エネルギ対、等であり得る。他の例において、光子計数CTによれば、当該再構成プロセッサは、例えば3以上のエネルギ窓が存在する場合、kエッジ画像ベースを含む2以上の画像ベースを発生することができる。該再構成プロセッサは非スペクトル体積的画像データを発生することもできる。
図2は、本発明による好まし実施態様の概略フローチャートである。該実施態様は、特に完全に診断的な非造影画像に加えて高品質ヨードマップが必要とされる臨床状況においてスペクトルCTスキャナを最適な態様で使用するための方法である。該方法は、スペクトルCT結果が品質の劣化を被り得る低線量スキャンに特に関係するものであり得る。キーとなるアイデアは、当該2つの方法の何れかから別々に達成され得るものより優れたヨードマップを提供するために、位置合わせ減算技術からの及びスペクトル解析技術からの相補的情報を組み合わせることである。入力される非造影コンピュータ断層撮影スキャン及び造影コンピュータ断層撮影スキャン並びに対応する利用可能なスペクトルデータから、2つの造影剤(例えば、ヨード)濃度マップが得られる。第1マップは従来の減算−位置合わせ方法を介してのものである一方、第2マップは上記スペクトルデータから導出される。両方法は、導入部分で説明されたものである。最終的に、これら2つのマップは二次造影剤濃度マップを生成するために組み合わされる。このステップは、図7及び図8において更に詳細に説明される。
図3は、上記2つの技術から導出された造影剤マップの間の根本的差異を示す(明瞭化のために、1D画像プロファイルのみが示されている)。曲線301は、位置合わせ−減算マップの結果であり、局所化された領域に大きなエラーを生じさせる残存した誤位置合わせアーチファクト330を伴っている。位置合わせが良好である他の領域において、該造影剤マップは良好な品質を示している。点線曲線311及び312は、曲線301を計算するために使用された造影及び非造影画像に各々対応している。曲線312は誤って位置合わせされているので、空間的誤位置合わせずれが重大なアーチファクト330を生じさせている。
曲線302は、造影スキャンのみから導出されたスペクトル解析の結果である。局部的な鋭いアーチファクトは存在しないが、全体的ノイズは、実際の構造として解釈され得るような幾らかの低周波数ノイズを含み、著しく高い。
提案される方法における狙いは、この図示に関しては、上記2つの不十分な品質の曲線から最適化され且つ改善された曲線を発生させることである。
造影剤濃度マップの対応する臨床例が図4に示されている。図4のa、b及びcは、従来のコンピュータ断層撮影スキャンから位置合わせ−減算技術を介して取得されたヨード濃度マップの例を示す。図4のd及びeは、スペクトルデュアルエネルギコンピュータ断層撮影スキャンから取得されたヨード濃度マップの例を示す。
本発明の一態様においては、位置合わせ−減算法を介して取得された造影剤濃度マップの空間位置合わせ精度を改善することができる。造影前及び造影後スキャンの間の空間位置合わせ問題は、時には、造影剤投与後の局部的ハウンズフィールド単位(HU)の変化による画像不一致により影響される。スペクトルコンピュータ断層撮影は、実質的に非コントラスト強調を示す画像を用いることにより一層正確な位置合わせを実行することを支援することができ、その場合、上記画像は実際に造影されたスキャンから、又は、そうでなければ、反復方式で位置合わせアーチファクトを補正することにより導出される。これら2つのオプションは、図5及び図6のフローチャートに示される。
図5においては、位置合わせ減算マップの処理のための位置合わせ精度を改善するために、スペクトル仮想非造影画像(VNC)ボリューム(最初に計算される)が使用される。このオプションは、使用されているスペクトルコンピュータ断層撮影システムが高品質のVNC結果を提供する場合、又は高造影剤濃度により造影前スキャン及び造影後スキャンの間に非常に大きなHU変化が存在する臨床プロトコルに特に適し得る。
図6において、初期位置合わせ−減算マップはスペクトルマップからの情報を用いて補正される。補正された位置合わせ−減算マップが造影画像から減算されて、一種の“従来の”仮想非造影マップを形成し、該マップは造影前画像ボリュームと位置合わせされる。造影剤強調の仮想的削除は、より良好な位置合わせ結果を助けることができる。改善された変形フィールドが、改善された位置合わせ−減算マップを生成するために使用され、これが最終マップとして使用される。この方式においては空間位置合わせが2つの反復処理として実行されること指摘することが興味深い。
図7及び図8は、2つのマップ型を最適に組み合わせる方法ステップを示す。図7においては、位置合わせ−減算マップ及びスペクトルマップが先ず独立に計算される。組み合わせ造影剤濃度マップは、上記位置合わせ−減算マップからの及びスペクトルマップからの相補的情報を用いて直接生成される。太枠内のブロックは、図8において更に詳述される。
図8のフローチャートは、位置合わせ−減算マップからの及びスペクトルマップからの値を重み付けすることにより組み合わせ造影剤マップを如何にして生成するかのオプションを示す。キーポイントは、上記2つのマップの各々に関して局部的マップ品質の分析がなされ、それに従って体積的重みが導出されることである。例えば、1つのマップにおいて局部的マップ品質が低く、第2のマップにおいて局部的品質が高い場合、組み合わせ重みは高品質のマップに有利にする。“局部的”なる文言は、各画像ピクセル又は一群のピクセル若しくはボクセル、又はそれ以外で決定された関心領域を記述する。
オプションとして、当該2つの方法からのマップの尺度又は単位は、組み合わせの前に適応化し又は正規化することができる。
CTにおける空間的位置合わせの局部的信頼レベルは、文献“Supervised quality assessment of medical image registration: Application to intra-patient CT lung registration”, S.E. Muenzing etal., Medical Image Analysis, Dec. 2012に公開された方法等の幾つかの既知の技術により評価することができる。
スペクトル分析の局部的信頼レベルは、放射線レベル及び造影剤濃度の範囲に関する既知で事前測定されたシステム限界を考慮することにより評価することができる。局部的線量レベルは、計算された線量マップ等のスキャンパラメータ及びツールから得ることができる。
オプションとして、2つのマップの重みの間の滑らかな移行は、2つの重み分布の特別な平滑化により発生させることができる。
以上、本発明を図面及び上記記載において詳細に図示及び説明したが、このような図示及び説明は解説的又は例示的であって、限定するものでないと見なされるべきである。即ち、本発明は論じられた実施態様に限定されるものではない。
開示された実施態様に対する他の変形例は、当業者によれば、請求項に記載された本発明を実施する際に図面、本開示及び添付請求項の精査から理解し、実施することができる。請求項において、“有する”なる文言は他の要素又はステップを排除するものではなく、単数形は複数を排除するものではない。単一のプロセッサ又は他のユニットは、請求項に記載された幾つかの項目の機能を満たすことができる。特定の手段が互いに異なる従属請求項に記載されているという単なる事実は、これら手段の組み合わせを有利に使用することができないということを示すものではない。請求項における如何なる符号も、当該範囲を限定するものと見なしてはならない。

Claims (15)

  1. 非造影コンピュータ断層撮影スキャン、造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンと同時に取得されるスペクトルコンピュータ断層撮影データから造影剤濃度マップを生成する方法であって、
    a.前記非造影コンピュータ断層撮影スキャン、前記造影コンピュータ断層撮影スキャン及び前記スペクトルコンピュータ断層撮影データから少なくとも2つの異なる一次造影剤濃度マップを生成するステップと、
    b.各一次造影剤濃度マップの局部的品質分析を実行するステップと、
    c.各一次造影剤濃度マップに対する局部的な体積的重みを前記局部的品質分析に基づいて決定するステップと、
    d.新たな造影剤濃度マップを前記少なくとも2つの一次造影剤濃度マップ及びこれらマップの対応する局部的な体積的重みに基づいて生成するステップと、
    を有する、方法。
  2. 前記一次造影剤濃度マップのうちの少なくとも1つが、前記非造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンの画像データに対して、体積的空間位置合わせステップ及び画像減算ステップを実行することにより得られる、請求項1に記載の方法。
  3. 前記造影コンピュータ断層撮影スキャンで得た画像ボリュームび前記スペクトルコンピュータ断層撮影データから仮想非造影画像ボリュームを計算し、前記非造影コンピュータ断層撮影スキャンで得た画像ボリュームを前記仮想非造影画像ボリュームに変形することを可能にする変形関数を計算し、該変形関数を前記造影コンピュータ断層撮影スキャンで得た画像ボリュームに適用することにより、前記体積的空間位置合わせステップを改善するステップを有する、請求項2に記載の方法。
  4. 前記新たな造影剤濃度マップを生成するステップは、
    a.前記体積的空間位置合わせステップ及び前記画像減算ステップを実行することにより得られる一次造影剤濃度マップを他の一次造影剤濃度マップからの情報を用いて改善するステップと、
    b.前記改善を行った一次造影剤濃度マップを前記造影コンピュータ断層撮影スキャンの画像データから減算して、改変された非造影コンピュータ断層撮影スキャンを得るステップと、
    c.前記改変された非造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンの画像データに対して、体積的空間位置合わせステップ及び画像減算ステップを実行するステップと、
    を有する、請求項2に記載の方法。
  5. 前記一次造影剤濃度マップのうちの少なくとも1つが、前記スペクトルコンピュータ断層撮影データのみから得られる、請求項1に記載の方法。
  6. 前記スペクトルコンピュータ断層撮影データから得られた一次造影剤濃度マップの前記局部的品質分析が、局部的放射線量レベルを考慮するステップを有する、請求項1に記載の方法。
  7. 前記局部的放射線量レベルが、スキャンパラメータから計算される線量マップから得られる、請求項6に記載の方法。
  8. 前記非造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンの少なくとも一方が低線量スキャンである、請求項1に記載の方法。
  9. 前記局部的品質分析及び前記対応する局部的な体積的重み付けが前記一次造影剤濃度マップの各ボクセルに対して実行される、請求項1に記載の方法。
  10. 各一次造影剤濃度マップのマップ値が、前記新たな造影剤濃度マップの生成の前に適応化及び/又は正規化される、請求項1に記載の方法。
  11. 前記一次造影剤濃度マップのうちの少なくとも1つの前記局部的品質分析が、局部的画像ノイズ推定、一群の隣接するピクセルに対する標準偏差測定、微細構造若しくは形状解析、局部的画像アーチファクト分析、空間解像度分析、局部的空間周波数若しくは波長解析、又は予め定められた範囲から外れるマップ値又はマップ値勾配の分析から選択される技術を有する、請求項1に記載の方法。
  12. 前記一次造影剤濃度マップの各々における前記局部的な体積的重みの間の滑らかな移行が形成される、請求項1に記載の方法。
  13. 医療用撮像装置と通信する装置であって、
    a.非造影コンピュータ断層撮影スキャン、造影コンピュータ断層撮影スキャン及び前記造影コンピュータ断層撮影スキャンと同時に取得されるスペクトルコンピュータ断層撮影データから少なくとも2つの異なる一次造影剤濃度マップを生成し;各一次造影剤濃度マップの局部的品質分析を実行し;各一次造影剤濃度マップに対する局部的な体積的重みを前記局部的品質分析に基づいて決定し;新たな造影剤濃度マップを前記少なくとも2つの一次造影剤濃度マップ及びこれらマップの対応する局部的な体積的重みに基づいて生成するプロセッサと、
    b.前記新たな造影剤濃度マップを表示する手段と、
    を有する、装置。
  14. 請求項1に記載の方法を実施する装置を有する、コンピュータ断層撮影スキャナ。
  15. プロセッサにより実行された場合に該プロセッサに請求項1に記載の方法を実施させるコンピュータ読取可能な命令によりコード化された、非一時的コンピュータ読取可能記憶媒体。
JP2018530859A 2015-12-17 2016-12-08 造影剤濃度マップを生成する方法 Active JP6684909B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP15200863 2015-12-17
EP15200863.7 2015-12-17
PCT/EP2016/080198 WO2017102529A1 (en) 2015-12-17 2016-12-08 Method for generating contrast agent concentration map

Publications (3)

Publication Number Publication Date
JP2018537226A JP2018537226A (ja) 2018-12-20
JP2018537226A5 JP2018537226A5 (ja) 2019-03-28
JP6684909B2 true JP6684909B2 (ja) 2020-04-22

Family

ID=55023925

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018530859A Active JP6684909B2 (ja) 2015-12-17 2016-12-08 造影剤濃度マップを生成する方法

Country Status (5)

Country Link
US (1) US10217247B2 (ja)
EP (1) EP3250127B1 (ja)
JP (1) JP6684909B2 (ja)
CN (1) CN107530037B (ja)
WO (1) WO2017102529A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110072457A (zh) 2016-12-15 2019-07-30 皇家飞利浦有限公司 对脉管结构进行可视化
US10580132B2 (en) * 2017-04-13 2020-03-03 Canon Kabushiki Kaisha Medical image processing apparatus, control method therefor, and non-transitory storage medium storing program
US11813101B2 (en) * 2018-01-31 2023-11-14 Koninklijke Philips N.V. Image quality improved virtual non-contrast images generated by a spectral computed tomography (CT) scanner
CN109934888B (zh) * 2019-04-24 2020-11-17 清华大学 非对比剂增强的磁共振动态血管成像方法及系统
DE102020206729A1 (de) * 2020-05-28 2021-12-02 Siemens Healthcare Gmbh Bewegungskorrekturverfahren

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002032292A2 (en) * 2000-10-18 2002-04-25 The Board Of Trustess Of The Leland Stanford Junior University Methods for development and use of diagnostic and therapeutic agents
US20040019275A1 (en) * 2002-07-23 2004-01-29 Maria Iatrou Method and apparatus for positioning a CT reconstruction window
US7787927B2 (en) * 2003-06-20 2010-08-31 Merge Cad Inc. System and method for adaptive medical image registration
US7894880B2 (en) * 2003-10-23 2011-02-22 The Board Of Trustees Of The Leland Stanford Junior University Measurement of renal extraction fraction using contrast enhanced computed tomography
JP4679068B2 (ja) * 2004-04-26 2011-04-27 株式会社東芝 X線コンピュータ断層撮影装置
DE102004021965B4 (de) * 2004-05-04 2009-05-07 Siemens Ag Verfahren zur Erstellung von tomographischen Aufnahmen eines schlagenden Herzens
US20080108894A1 (en) * 2004-11-15 2008-05-08 Elgavish Gabriel A Methods and Systems of Analyzing Clinical Parameters and Methods of Producing Visual Images
DE102006009222B4 (de) * 2006-02-28 2008-02-28 Siemens Ag Verfahren und Vorrichtung zur Bestimmung der Konzentration einer Substanz in einem Körpermaterial mittels Mehr-Energie-Computertomographie
US20070238968A1 (en) * 2006-04-05 2007-10-11 Vitaliy Rappoport Method and apparatus for providing low dose dual energy for PET attenuation correction
DE102007024158B4 (de) * 2007-05-24 2017-09-14 Bayer Intellectual Property Gmbh Auswahlverfahren für zwei Kontrastmittel zur Verwendung in einer Dual-Energy-CT-Untersuchung, Kontrastmittelkombination und Erzeugung von CT-Aufnahmen mit einer Kontrastmittelkombination mit unterschiedlichen Energiespektren
JP2009022450A (ja) 2007-07-18 2009-02-05 Ge Medical Systems Global Technology Co Llc X線ct装置および画像作成方法
EP2274727A1 (en) * 2008-04-08 2011-01-19 Vereniging voor Christelijk Hoger Onderwijs, Wetenschappelijk Onderzoek en Patiëntenzorg A pet data processing system, an arrangement, a method and a computer program product for determining a distribution of a tracer uptake
JP5539342B2 (ja) * 2008-07-18 2014-07-02 コーニンクレッカ フィリップス エヌ ヴェ スペクトルイメージング
US9342882B2 (en) * 2012-09-06 2016-05-17 Siemens Medical Solutions Usa, Inc. Methods and systems for determining the concentration of a contrast agent
CN104871207B (zh) * 2012-12-20 2018-07-17 皇家飞利浦有限公司 图像数据处理方法及系统
CN104700390B (zh) 2013-12-09 2019-07-23 通用电气公司 双能ct造影剂增强扫描图像中的钙化部分识别方法
JP6667999B2 (ja) * 2014-05-16 2020-03-18 キヤノンメディカルシステムズ株式会社 画像処理装置、画像処理方法、及びプログラム

Also Published As

Publication number Publication date
WO2017102529A1 (en) 2017-06-22
CN107530037B (zh) 2018-11-27
US10217247B2 (en) 2019-02-26
CN107530037A (zh) 2018-01-02
JP2018537226A (ja) 2018-12-20
EP3250127A1 (en) 2017-12-06
US20180276853A1 (en) 2018-09-27
EP3250127B1 (en) 2018-10-31

Similar Documents

Publication Publication Date Title
JP6684909B2 (ja) 造影剤濃度マップを生成する方法
Yu et al. Noise performance of low-dose CT: comparison between an energy integrating detector and a photon counting detector using a whole-body research photon counting CT scanner
Clark et al. Micro-CT of rodents: state-of-the-art and future perspectives
US10013779B2 (en) Metal artifact reduction for 3D-digtial subtraction angiography
US9025815B2 (en) System and method for multi-material correction of image data
US11328391B2 (en) System and method for controlling noise in multi-energy computed tomography images based on spatio-spectral information
US10111638B2 (en) Apparatus and method for registration and reprojection-based material decomposition for spectrally resolved computed tomography
US9547889B2 (en) Image processing for spectral CT
US20170023496A1 (en) Spectral x-ray imaging
US9367903B2 (en) Method and apparatus for correcting artifacts during generation of X-ray images, in particular computed tomography, or radiography by means of temporal modulation of primary radiation
US20110110573A1 (en) Image generation device with optimized dose control
US9374565B2 (en) System and method for multi-material correction of image data
JP7346429B2 (ja) スペクトルコンピュータ断層撮影(ct)スキャナによって生成される画質が向上したバーチャル非造影画像
Feng et al. An experimental method to correct low-frequency concentric artifacts in photon counting CT
JP2004237076A (ja) マルチ・モダリティ・イメージング方法及び装置
Stein et al. Photon-counting computed tomography–basic principles, potenzial benefits, and initial clinical experience
Wiegert et al. Projection Extension for Region Of Interest Imaging in Cone-Beam CT1
Dinten et al. Dual-energy x-ray absorptiometry using 2D digital radiography detector: application to bone densitometry
US9262825B2 (en) Image reconstruction in interleaved multi-energy imaging
JP7442455B2 (ja) 分光イメージングデータからの分光炎症マップ
US20210304458A1 (en) Method for providing automatic adaptive energy setting for ct virtual momochromatic imaging
CN116172594B (zh) 用于生成患者的结果图像数据集的方法及装置
Haase et al. Estimation of statistical weights for model-based iterative CT reconstruction
Puett et al. Contrast enhanced imaging with a stationary digital breast tomosynthesis system
Feng Overcoming Spectral Inconsistency-induced Image Artifacts in Photon Counting CT Systems

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190212

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190212

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20190212

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20190411

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190423

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190722

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191007

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191219

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200303

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200330

R150 Certificate of patent or registration of utility model

Ref document number: 6684909

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250