JP7340549B2 - kVpスイッチングマルチエネルギー取得から非スペクトルボリュメトリック画像データを生成するイメージングシステム - Google Patents

kVpスイッチングマルチエネルギー取得から非スペクトルボリュメトリック画像データを生成するイメージングシステム Download PDF

Info

Publication number
JP7340549B2
JP7340549B2 JP2020571549A JP2020571549A JP7340549B2 JP 7340549 B2 JP7340549 B2 JP 7340549B2 JP 2020571549 A JP2020571549 A JP 2020571549A JP 2020571549 A JP2020571549 A JP 2020571549A JP 7340549 B2 JP7340549 B2 JP 7340549B2
Authority
JP
Japan
Prior art keywords
spectral
kvp
projection data
image data
sets
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
JP2020571549A
Other languages
English (en)
Other versions
JP2021529027A (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 JP2021529027A publication Critical patent/JP2021529027A/ja
Application granted granted Critical
Publication of JP7340549B2 publication Critical patent/JP7340549B2/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
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • 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
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/408Dual energy

Description

本発明は、概して、イメージング(撮像)に関し、より詳細には、kVpスイッチングマルチエネルギー取得からスペクトルボリュメトリックデータ及び非スペクトルボリュメトリックデータを生成するように構成されたイメージングシステムに関し、特に、コンピュータトモグラフィ(CT)への応用について説明される。
非スペクトルコンピュータトモグラフィ(CT)スキャナは一般に、検出器の1つ又は複数の行の反対側の回転可能なガントリ上に取り付けられたX線管を有する。X線管は、X線管と検出器の1又は複数の行との間に位置する検査領域の周りを回転し、検査領域を横切る広帯域放射線を放出する。例えば、120keVのピーク(最大)管電圧(すなわち、120kVp)では、放出される放射線のエネルギースペクトル(低エネルギー光子のフィルタリングを伴う)は、40keV~120keVであり得る。検出器の1又は複数の行は検査領域を横断する放射線を検出し、それを示す投影データ(線積分)を生成する。投影データは、ボリュメトリック画像データを生成するように再構成される。
再構成されたボリュメトリック画像データのボクセルは、相対放射線密度に対応するグレースケール値を使用して表示される。グレースケール値は、スキャンされた被検体の減衰特性を反映し、一般に、スキャンされた被検体内の解剖学的構造などの構造を示す。材料による光子の減衰は、材料を横切る光子のエネルギーに依存するので、検出された放射線は、被検体のスキャンされた材料の元素又は材料組成(例えば、原子番号)を示す追加の情報を提供するスペクトル情報も含む。しかしながら、投影データの値は、エネルギースペクトル(例えば、40keV~120keV)にわたって積分されるエネルギーフルエンスに比例し、ボリュメトリックデータは、エネルギー依存の情報を反映しない。
スペクトル(マルチエネルギー)CTスキャナは、異なるエネルギーバンドに対する投影データを生成するように構成される。一例では、これはkVpスイッチングによって達成される。例えば、デュアルエネルギー構成では、第1の電圧(例えば、より低いkVp)が第1の積分期間にわたってX線管電圧に印加され、第2の電圧(例えば、より高いkVp)が第2の積分期間にわたってX線管電圧に印加され、第1の電圧は第3の積分期間にわたってX線管電圧に印加され、第2の電圧は第4の積分期間にわたってX線管電圧に印加される。より低い及びより高いkV投影データは、例えば光電効果及びコンプトン散乱成分にデコンポーズ(分解)され、これらは、個々に再構成され、次いで、組み合わされて、仮想的な単エネルギーのボリュメトリックデータを生成する。
このように構成され、kVpスイッチングモードで動作するCTスキャナは、仮想単エネルギーボリュメトリックデータを再構成するが、一般に、非スペクトルボリュメトリックデータを再構成しない。残念ながら、仮想単エネルギーボリュメトリックデータは、非スペクトルボリュメトリックデータに反映された全ての情報を含まない。例えば、仮想単エネルギーボリュメトリックデータは、非スペクトルボリュメトリックデータと同じ軟組織コントラストを示さない。しかしながら、放射線科医は、非スペクトルボリュメトリック画像データを使用して、微妙な軟組織コントラスト変化を病理として識別するように訓練されてきた。その結果、仮想単エネルギーボリュメトリックデータでは、病理による軟組織のコントラスト変化を見つけることがより困難になり、病理が見逃される可能性がある。
本明細書で説明される態様は、上記で参照された問題及び他の問題に対処する。
1つの態様において、再構成システムは、少なくとも2つの放射線源電圧間のkVpスイッチングを通じて生成された少なくとも2組の投影データをデコンポーズするように構成されるデコンポーザを有する。各組は、少なくとも2つの放射線源電圧のそれぞれ別個の放射線源電圧に対応する。システムは、少なくとも2組の投影データを処理しスペクトル画像データを生成するように構成されるスペクトルチャネルを更に有する。システムは、少なくとも2組の投影データを処理し予め決められた基準kVpの非スペクトル画像データを生成するように構成される非スペクトルチャネルを更に有する。基準kVpは、ユーザによって、又はシステムから、例えば、スカウトスキャンの分析によって、予め決定されることができる。
別の態様において、コンピュータ可読媒体が、コンピュータ実行可能命令で符号化され、前記コンピュータ実行可能命令は、コンピュータのプロセッサによって実行されるとき、プロセッサに、少なくとも2つの放射線源電圧間のkVpスイッチングを通じて生成された少なくとも2組の投影データをデコンポーズするステップであって、各組の投影データは少なくとも2つの放射線源電圧のうちのそれぞれ別個の放射線源電圧に対応する、ステップと、少なくとも2組の投影データを処理し、予め決められた基準kVpについて非スペクトル画像データを生成するステップと、を実行させる。
別の態様において、方法は、少なくとも2つの放射線源電圧間のkVpスイッチングを通じて生成される少なくとも2組の投影データをデコンポーズするステップを有する。各組の投影データは、少なくとも2つの放射線源電圧のうちのそれぞれ別個の放射線源電圧に対応する。本方法は、少なくとも2組の投影データに基づいて、予め決められた基準kVpの非スペクトル画像データを生成するステップを更に有する。
本発明は、さまざま構成要素及び構成要素の取り合わせ、並びにさまざまなステップ及びステップの取り合わせの形をとることができる。図面は、好ましい実施形態を例示する目的のためだけのものであり、本発明を限定するものとして解釈されるべきではない。
kVpスイッチングのために構成され、kVpスイッチングマルチエネルギー取得からスペクトルボリュメトリックデータ及び非スペクトルボリュメトリックデータを生成するよう構成された再構成器を有する例示的なイメージングシステムを概略的に示す図。 図1の再構成器の一例を概略的に示す図。 図1の再構成器の別の例を概略的に示す図。 本明細書で説明される実施形態による例示的な方法を示す図。 本明細書で説明される実施形態による別の例示的な方法を示す図。 例示的な仮想単エネルギー画像を示す図。 本明細書に記載されるように生成された例示的なkVp画像を示す図。 図6の例示的な仮想単エネルギー画像と図7のkVp画像との間の差を示す図。
以下では、kVpスイッチングのために構成され、kVpスイッチングマルチエネルギー取得からスペクトル及び非スペクトルボリュメトリック画像データを生成するように構成されるイメージングシステムの例について説明する。
最初に図1を参照して、コンピュータトモグラフィ(CT)スキャナのようなイメージングシステム100が概略的に示されている。イメージングシステム100は、ほぼ静止したガントリ102と、静止ガントリ102によって回転可能に支持され、z軸の周りの検査領域106の周りを回転する回転ガントリ104とを有する。
X線管のような放射線源108が、回転ガントリ104によって回転可能に支持され、回転ガントリ104と共に回転し、検査領域106を横切る広帯域放射線を放出する。放射線源電圧コントローラ(kVp CTRL)110は、放射線源108のピーク放出電圧を制御する。kVp CTRL110は、X線管電圧を少なくとも2つの電圧(例えば、80kVp及び140kVp等)の間で切り換えるように構成される。その結果、放射線源108は、第1のエネルギースペクトルを有する少なくとも第1の放射線ビームと、第2の異なるエネルギースペクトルを有する第2の放射線ビームとを交互に生成する。
検出器112は、検出器素子114の1次元又は2次元アレイを有し、各行は、xy平面内に延在し、複数の行が、z方向に沿って配置される。検出器112は、検査領域106をはさんで放射線源108に対向する角弧に沿って、回転ガントリ104によって回転可能に支持される。検出器112は放射線源108と協調して回転し、検査領域106を横断する放射線を検出し、異なるエネルギースペクトルの各々について異なる組の投影データを生成する。
一例では、kVp CTRL110が、積分期間(IP)と積分期間(IP)の間にX線管電圧を交替させる。一般に、積分期間IPは、検出器112が測定のために所定の角度インクリメントにわたって回転する間に放射線を検出する期間である。各積分期間IPごとに、各検出器素子114は線積分を生成する。積分期間IP/角度インクリメントにわたる線積分の組が、ビューである。投影データは、少なくとも180°にわたって取得された1組のビューに、それぞれ異なるエネルギースペクトルの各々についてファン角度を加えたものを含む。図示の実施形態では、異なるエネルギースペクトルのビューがインタリーブされる。
再構成器116は、それぞれ異なるエネルギースペクトルの投影データを再構成し、ボリュメトリックデータを生成する。これは、以下に詳述するように、投影データからスペクトルボリュメトリックデータ及び非スペクトルボリュメトリックデータを生成することを含む。スペクトルボリュメトリック画像データの例には、低エネルギー及び高エネルギー、単エネルギー/モノクロ、仮想非コントラスト、有効Z(原子番号)、ヨウ素のみ、及び/又は他のスペクトルボリュメトリック画像データなどが含まれる。非スペクトルボリュメトリック画像データは、予め決められたX線管kVpに対応する。
一例では、再構成器116は、物理メモリ及び/又は他の非一時メモリのようなコンピュータ可読記憶媒体(一時媒体を除く)上に記憶され、埋め込まれ、符号化されるなどしたコンピュータ実行可能命令を実行するよう構成されるプロセッサ(例えば、中央処理ユニット、マイクロプロセッサなど)を用いて実現される。いくつかの実施形態では、再構成器116は、グラフィックス処理ユニット(GPU)、特定用途向け集積回路(ASIC)などの特化したハードウェアを有する。再構成器116は、図示されるように装置100の一部であり、及び/又は装置から離れた場所にありうる。
寝台のような被検体支持体118は、検査領域106内の物体又は被検体を支持する。被検体支持体118は、被検体又は物体を、ロードし、スキャンし、及び/又はアンロードするために、検査領域106に対して被検体又は物体をガイドするように、撮像プロシージャを実行することと協調して移動可能である。
オペレータコンソール120は、オペレータがkVpスイッチングイメージングプロトコル、再構成アルゴリズムを選択するなど、システム100の動作を制御することを可能にする。オペレータコンソール120は、マウス、キーボードなどの入力装置と、ディスプレイモニタなどの出力装置とを有する。
図2は、再構成器116の一例を概略的に示す。説明の目的で、この例は、放射線源108の電圧が、ビュー間に2つの電圧(より低い、及びより高い)間で切り換えられ、検出器112は、低kV投影データ及び高kV投影データを生成するデュアルエネルギースキャンに関連して記載される。
角度リビナ202は、低kVスペクトル投影データの組及び高kVスペクトル投影データの組を受け取り、それらを角度に関してリビニングし、低kVの平行線積分の組及び高kVの平行線積分の組を生成する。1つの非限定的な例では、これは既知の又は他の角度に関するリビニングアプローチによって達成される。適切なアプローチの一例は、2002年10月25日出願の「Dynamic Detector Interlacing for Computed Tomography」という名称の米国特許第6,963,631B2号に記載されており、その全体が参照により本明細書に組み込まれる。
デコンポーザ204は、低kV及び高kV平行線積分の組を、光電効果及びコンプトン散乱又は他の基準のようなそれぞれ異なる基底又は寄与(成分)にデコンポーズする。非限定的なアプローチは、エネルギーの関数として各基底(物質)の寄与を有するルックアップテーブル(LUT)を作成することと、LUTを記憶することと、線積分を2つの基底の寄与にデコンポーズするためにLUT内の値を使用することとを含む。デコンポーズのためにLUTを作成し使用するための適切なアプローチの例は、2013年12月4日出願の「マルチエネルギーイメージング」という名称の米国特許第9,324,142B2号に記載されており、その全体が参照により本明細書に組み込まれる。
別の手法では、Ih,l=∫Sh,l(E)D(E)exp(-∫μ(E)dl)dEとして低kV及び高いkVの平行線積分をモデル化し、ここで、添え字h,lは、高(high)及び低(low)を指し、Sh,l(E)は、高低の管スペクトルであり、exp(-∫μ(E)dl)dEは、は線積分を表す。デュアルエネルギー処理は、∫μ(E)dlを基底成分に分離し、従って例えばIh,l=∫Sh,l(E)D(E)exp(-LB1B1(E)-LB2B2(E))dEであり、ここでLB1及びLB2は、2つの基底対の寄与であり、D(E)は、検出器のスペクトル応答を表し、B1(E)及びB2(E)は、それらの対応する減衰エネルギー依存である。デコンポーズは画素ごとに、2つの未知数LB1及びLB2を測定値Ih、Ilにマッピングする非線形系を解く。ここで、2つの基底対の寄与は、光電効果及びコンプトン散乱寄与、Ih,l=∫Sh,l(E)D(E)exp(-LpP(E)-LsS(E))dEであり、ここで、Lp及びLsは、光電効果及びコンプトン散乱寄与であり、P(E)及びS(E)は、それらの対応する減衰のエネルギー依存であり、デコンポーズは、各画素について、2つの未知数Lp及びLsを、測定値Ih、Ilにマッピングする非線形システムを解く。本明細書において、他の基底対も企図される。別の例は、Alverez, et al., "Energy-selective Reconstructions in X-ray Computerized Tomography," Phys. Med. Biol. 1976, Vol. 21, No. 5, 733-744に記載されている。他のアプローチもまた、本明細書において意図される。
再構成器116は、両方の寄与を処理する2つの処理チェーン206及び208を有する。処理(スペクトル)チェーン206は寄与を処理し、スペクトルボリュメトリックデータを生成する。処理(非スペクトル)チェーン208は寄与を処理し、非スペクトルボリュメトリック画像データを生成する。処理チェーン206は、2つの寄与を受け取り、等間隔の平行線積分を生成するために、平行線積分を(例えば補間を介して)半径方向にリビニングする半径方向リビナ210を有する。1つの非限定的な例において、これは既知の又は他の角度に関するリビニングアプローチによって達成される。適切なアプローチの例は、米国特許第6,963,631号B2に記載されている。処理チェーン206は、デコンポーズされ半径方向にリビニングされた投影データを再構成し、第1及び第2の寄与について第1及び第2のボリュメトリック画像データを生成する画像プロセッサ212を更に有する。画像プロセッサ212は、フィルタード逆投影、反復、及び/又は他の再構成アプローチを使用する。
処理チェーン206は、スペクトルボリュメトリック画像データを生成するために第1及び第2のボリュメトリック画像データを組み合わせる(例えば、重み付けされた合計)スペクトルプロセッサ214を更に有する。例えば、画像プロセッサ212は、単エネルギーボリュメトリックデータ画像を生成することができる。本明細書に記載されるように、他のスペクトルボリュメトリックデータは、仮想的な非コントラスト、有効Z、ヨウ素のみ等のスペクトルボリュメトリックデータを含む。
処理チェーン208は、予め決められた基準kVp218について、デコンポーズされた投影データを再コンポーズする再コンポーザ216を有する。予め決められた基準kVp218は、デフォルトであるか、ユーザ指定されるか、プログラム可能な基準電圧である。再コンポーズの例は、Ireference_kVp=∫D(E)Sreference_kVp(E)exp(-LpP(E)-LsS(E))dEである。例えば、低kVp及び高kVpは、80kVp及び140kVpであり、基準kVp218が120kVpである場合、再コンポーズは、I120_kVp=∫D(E)120_kVp(E)exp(-LpP(E)-LsS(E))dEである。一例では、これはエネルギースペクトルの全てのエネルギーについて、特定のエネルギーにおける材料の減衰に、基準kVp218についてのその特定のエネルギーにおけるエネルギースペクトルの部分を乗じ、次いで、積を合計することによって実現される。
処理チェーン208は更に、例えば、半径方向リビナ210に関連して説明したように、再コンポーズされた投影データを受け取り、平行線積分を半径方向にリビニングして、等間隔の平行線積分を生成する半径方向リビナ220を有する。1つの例において、半径方向リビナ210及び220は図示のように、別個の半径方向リビナである。変形例において、半径方向リビナ210及び220は、同じ半径方向リビナの一部であるか、又は同じ半径方向リビナである。
処理チェーン208は更に、例えば画像プロセッサ212に関連して説明したように、再コンポーズされ半径方向にリビニングされた投影データを受け取り、データを再構成する画像プロセッサ222を有する。一例において、画像プロセッサ212及び222は、図示されるように、それぞれ別個の画像プロセッサである。変形例において、画像プロセッサ212及び222は、同一の画像プロセッサの一部であるか、又は同一の画像プロセッサである。
一般に、処理チェーン208は、異なるソーススペクトルで取得されたインタリーブされた取得を処理し、特定のkVp取得からの投影データで再構成された非スペクトルボリュメトリック画像データのように見える非スペクトルボリュメトリック画像データを生成する。例として、非限定的な一例において、処理チェーン208が、80kVp及び140kVpでインタリーブされた取得を処理し、単一の120kVp取得からの投影データで再構成された非スペクトルボリュメトリック画像データのように見えるボリュメトリック画像データを生成する。
図3は、再構成器116の別の例を概略的に示す。この例は、角度リビナ202、デコンポーザ204、第1のチェーン206(半径方向リビナ210、画像プロセッサ212、及びスペクトルプロセッサ214)を有する。これらのコンポーネントは上述されており、ここでの説明は繰り返さない。この例では、スペクトルプロセッサ214が、例えば40keVから120keVの放射線ビームのエネルギースペクトルにわたって、単エネルギー画像を生成する。
この例では、チャネル208は、(上述したように)基準kVp218及び非スペクトルプロセッサ302を有する。非スペクトルプロセッサ302は、放射線ビームエネルギースペクトル(例えば、40keVから120keVまで)にわたって単エネルギー画像の加重平均を計算することによって、非スペクトルボリュメトリックデータを生成する。この場合、相対加重は、被検体又は平均的な被検体によってフィルタリングされた基準kVp(例えば、120kVp)における実効X線管スペクトルから導出される。この加重は、非スペクトルボリュメトリック画像データにおいて示されるコントラストを模倣することになる。
図4は、本明細書で説明される実施形態による例示的な方法を示す。以下の動作の順序は限定的なものではなく、他のシリアル処理及び/又はパラレル処理など、他の順序も本明細書で企図されることを理解されたい。
402において、放射線源108の電圧は、CTスキャンのために検査領域106の周りで放射線源108を回転させながら、kVp CTRL110を通じて、少なくとも2つの異なる電圧の間で切り替えられる。
404において、検出器112は、少なくとも2つの異なるX線管電圧のそれぞれについて別々に放射線を検出し、2つの電圧のうち第1の電圧についての第1の投影データと、2つの電圧のうち第2の電圧についての第2の投影データとを生成する。
406において、処理チェーン206は、第1及び第2の投影データを処理して、本明細書で説明されるように及び/又は他のやり方で、スペクトルボリュメトリック画像データを生成する。
408において、処理チェーン208は、第1及び第2の投影データを処理して、本明細書で説明されるように及び/又は他のやり方で、非スペクトルボリュメトリック画像データを生成する。
410において、スペクトルボリュメトリック画像データ及び/又は非スペクトルボリュメトリック画像データが出力され、例えば、表示され、保存され、別の装置に伝達される。
図5は、本明細書で説明される実施形態による別の例示的な方法を示す。
502において、放射線源108の電圧は、CTスキャンのために検査領域106の周りで放射線源108を回転させながら、kVp CTRL110を通じて、少なくとも2つの異なる電圧の間で切り替えられる。
504において、検出器112は、少なくとも2つの異なるX線管電圧のそれぞれについて別々に放射線を検出し、2つの電圧のうち第1の電圧についての第1の投影データと、2つの電圧のうち第2の電圧についての第2の投影データとを生成する。
506において、処理チェーン206は、第1及び第2の投影データを処理し、本明細書で説明されるように及び/又は他のやり方で、スペクトルボリュメトリック画像データを生成する。
508において、スペクトルボリュメトリック画像データは、本明細書で説明されるように及び/又は他のやり方、非スペクトルボリュメトリック画像データを生成するよう処理される。
510において、スペクトルボリュメトリック画像データ及び/又は非スペクトルボリュメトリック画像データが出力され、例えば、表示され、保存され、他の装置に伝達される。
上記はコンピュータ可読記憶媒体(一時媒体を除く)上に符号化され又は埋め込まれたコンピュータ可読命令であって、コンピュータプロセッサ(例えば、中央処理装置(CPU)、マイクロプロセッサなど)によって実行されるとき、プロセッサに本明細書に記載の動作を実行させるコンピュータ可読媒体によって実現されることができる。加えて又は代替として、コンピュータ可読命令のうちの少なくとも1つは、コンピュータ可読記憶媒体ではない信号、搬送波、又は他の一時媒体によって担持される。
図6-図8は、単エネルギー画像が、非スペクトル画像とは異なるコントラストを反映することを示している。図6は、単エネルギー画像を示し、図7は、本明細書に記載されたアプローチを用いて生成されたkVp非スペクトル画像を示し、図8は、それらの間の差分画像を示す。図8の差分画像は、図6の単エネルギー画像と図7の非スペクトル画像との間のコントラスト差を表す。
本発明が図面及び前述の説明において詳細に図示及び説明されてきたが、そのような図示及び説明は説明的又は例示的であり、限定的ではないと考えられるべきである。本発明は開示された実施形態に限定されるものではない。開示された実施形態に対する他の変形は、図面、開示、及び添付の請求項の検討から、請求項に記載された発明を実施する際に当業者によって理解され、及び実施されることができる。
請求項において、単語「有する、含む(comprising)」は、他の構成要素又はステップを排除せず、不定冠詞「a」又は「an」は複数性を排除しない。単一のプロセッサ又は他のユニットが、請求項に列挙されるいくつかの項目の機能を果たすことができる。特定の手段が相互に異なる従属請求項に記載されているという単なる事実は、これらの手段の組み合わせが有利に使用されることができないことを示すものではない。
コンピュータプログラムは、他のハードウェアと一緒に、又はその一部として供給される光記憶媒体又はソリッドステート媒体などの適切な媒体上に記憶され/配布されることができるが、インターネット又は他の有線又は無線電気通信システムを通じてなど、他の形態で配布されることもできる。請求項におけるいかなる参照符号も、その範囲を限定するものとして解釈されるべきではない。

Claims (17)

  1. 再構成システムであって、
    少なくとも2つの放射線源電圧の間のkVpスイッチングを通じて生成される少なくとも2組の投影データをデコンポーズするデコンポーザであって、各組の投影データが前記少なくとも2つの放射線源電圧のそれぞれ別個の放射線源電圧に対応する、デコンポーザと、
    前記少なくとも2組の投影データを処理してスペクトル画像データを生成するスペクトルチャネルと、
    前記少なくとも2組の投影データを処理して予め決められた基準kVpの非スペクトル画像データを生成する非スペクトルチャネルと、
    を有
    前記非スペクトルチャネルは、前記予め決められた基準kVpについて前記少なくとも2組の投影データを再コンポーズする再コンポーザを有し、
    前記再コンポーザは、前記エネルギースペクトルの全てのエネルギーについて、特定のエネルギーにおける減衰に前記基準kVpにおけるエネルギースペクトルの一部を乗じ、その積を合計することによって、前記少なくとも2組の投影データを再コンポーズするよう構成される、
    再構成システム。
  2. 前記非スペクトルチャネルは更に、
    前記再コンポーズされた投影データをリビニングする半径方向リビナと、
    前記リビニングされた投影データを再構成して非スペクトル画像データを生成する画像プロセッサと、
    を有する、請求項1に記載の再構成システム。
  3. 前記再コンポーザは、Ireference_kVp=∫D(E)Sreference_kVp(E)exp(-LpP(E)-LsS(E))dEを使用して前記少なくとも2組の投影データを再コンポーズするよう構成され、ここで、Ireference_kVpは基準kVpにおける強度測定値を表し、Sreference_kVpは基準kVpにおける放射線源のエネルギースペクトルを表し、Lpは光電効果寄与を表し、Lsはコンプトン散乱寄与を表し、D(E)は検出器のスペクトル応答を表し、P(E)は光電効果寄与の減衰のエネルギー依存を表し、S(E)は、コンプトン散乱の減衰のエネルギー依存を表す、請求項1又は2に記載の再構成システム。
  4. 前記非スペクトルチャネルは更に、所定の重みの組により、放射線源のエネルギースペクトルにわたる単エネルギースペクトル画像データの加重平均を計算することによって、非スペクトルボリュメトリックデータを生成する非スペクトルプロセッサを有する、請求項1に記載の再構成システム。
  5. 前記所定の重みの組は、被検体又は平均的な被検体によってフィルタリングされた前記基準kVpにおける有効管スペクトルから導出される、請求項4に記載の再構成システム。
  6. 前記スペクトルチャネルは更に、
    前記デコンポーズされた投影データをリビニングする半径方向リビナと、
    リビニングされた投影データを再構成してスペクトル基底画像データを生成する画像プロセッサと、
    スペクトル基底画像データを処理してスペクトルボリュメトリック画像データを生成するスペクトルプロセッサと、
    を有する、請求項1に記載の再構成システム。
  7. 前記非スペクトルチャネルは、前記スペクトルボリュメトリック画像データを処理して非スペクトルボリュメトリック画像を生成する非スペクトルプロセッサを有する、請求項6に記載の再構成システム。
  8. 前記予め決められた基準kVpは、デフォルト値である、請求項1乃至7のいずれか1項に記載の再構成システム。
  9. 前記予め決められた基準kVpは、ユーザにより定められた値である、請求項1乃至7のいずれか1項に記載の再構成システム。
  10. 請求項1に記載の再構成システムを有するコンピュータトモグラフィイメージングシステム。
  11. コンピュータ実行可能命令により符号化されたコンピュータ可読媒体であって、前記コンピュータ実行可能命令は、プロセッサによって実行されるとき、前記プロセッサに、
    少なくとも2つの放射線源電圧間のkVpスイッチングを通じて生成される少なくとも2組の投影データをデコンポーズするステップであって、各組の投影データは前記少なくとも2つの放射線源電圧のうちのそれぞれ別個の放射線源電圧に対応する、ステップと、
    前記少なくとも2組の投影データを処理して、スペクトル画像データを生成するステップと、
    前記少なくとも2組の投影データを処理して、予め決められた基準kVpの非スペクトル画像データを生成するステップと、
    前記予め決められた基準kVpについて前記少なくとも2組の投影データを再コンポーズするステップと、
    前記エネルギースペクトルの全てのエネルギーについて、特定のエネルギーにおける減衰に前記基準kVpにおけるエネルギースペクトルの一部を乗じ、その積を合計することによって、前記少なくとも2組の投影データを再コンポーズするステップと
    を実行させる、コンピュータ可読媒体。
  12. 前記コンピュータ実行可能命令は更に、前記プロセッサによって実行されるとき、前記プロセッサに、
    前記スペクトル画像データを処理し、非スペクトル画像データを生成するステップ、
    を実行させる、請求項11に記載のコンピュータ可読媒体。
  13. 前記コンピュータ実行可能命令は更に、前記プロセッサによって実行されるとき、前記プロセッサに、
    再コンポーズされた少なくとも2組の投影データに基づいて非スペクトル画像データを生成するステップと、
    を実行させる、請求項11に記載のコンピュータ可読媒体。
  14. 少なくとも2つの放射線源電圧の間のkVpスイッチングを通じて生成される少なくとも2組の投影データをデコンポーズするステップであって、各組の投影データは前記少なくとも2つの放射線源電圧のうちのそれぞれ別個の放射線源電圧に対応する、ステップと、
    前記少なくとも2組の投影データを処理して、スペクトル画像データを生成するステップと、
    前記少なくとも2組の投影データに基づいて、予め決められた基準kVpの非スペクトル画像データを生成するステップと、
    前記予め決められた基準kVpについて前記少なくとも2組の投影データを再コンポーズするステップと、
    前記エネルギースペクトルの全てのエネルギーについて、特定のエネルギーにおける減衰に前記基準kVpにおけるエネルギースペクトルの一部を乗じ、その積を合計することによって、前記少なくとも2組の投影データを再コンポーズするステップと
    を有する方法。
  15. 再コンポーズされた前記少なくとも2組の投影データに基づいて前記非スペクトル画像データを生成するステップと、
    を更に有する、請求項14に記載の方法。
  16. 前記少なくとも2組の投影データを処理してスペクトル画像データを生成するステップを更に有する、請求項15に記載の方法。
  17. 前記スペクトル画像データを処理して前記非スペクトル画像データを生成するステップを更に有する、請求項16に記載の方法。
JP2020571549A 2018-06-29 2019-06-24 kVpスイッチングマルチエネルギー取得から非スペクトルボリュメトリック画像データを生成するイメージングシステム Active JP7340549B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201862691999P 2018-06-29 2018-06-29
US62/691,999 2018-06-29
PCT/EP2019/066570 WO2020002174A1 (en) 2018-06-29 2019-06-24 Imaging system configured to generate non-spectral volumetric image data from a kvp switching multi-energy acquisition

Publications (2)

Publication Number Publication Date
JP2021529027A JP2021529027A (ja) 2021-10-28
JP7340549B2 true JP7340549B2 (ja) 2023-09-07

Family

ID=67003512

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020571549A Active JP7340549B2 (ja) 2018-06-29 2019-06-24 kVpスイッチングマルチエネルギー取得から非スペクトルボリュメトリック画像データを生成するイメージングシステム

Country Status (5)

Country Link
US (1) US20210272338A1 (ja)
EP (1) EP3815046A1 (ja)
JP (1) JP7340549B2 (ja)
CN (1) CN112334950A (ja)
WO (1) WO2020002174A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3967233A1 (en) * 2020-09-09 2022-03-16 Koninklijke Philips N.V. Controller for a ct scanner

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006503631A (ja) 2002-10-25 2006-02-02 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コンピュータ断層撮影用のダイナミック検出器インターレーシング
US20120045109A1 (en) 2009-05-12 2012-02-23 Koninklijke Philips Electronics N.V. Perfusion imaging
JP2014036836A (ja) 2012-07-17 2014-02-27 Toshiba Corp X線ct装置およびx線ct装置の制御方法
US20140056503A1 (en) 2011-04-28 2014-02-27 Koninklijke Philips N.V. Multi-energy imaging
WO2016147844A1 (ja) 2015-03-19 2016-09-22 株式会社日立製作所 X線ct装置及びマルチエネルギー像作成方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005009206A2 (en) * 2003-06-25 2005-02-03 Besson Guy M Dynamic multi-spectral imaging system
US8199875B2 (en) * 2009-12-11 2012-06-12 General Electric Company System and method of acquiring multi-energy CT imaging data
JP5869001B2 (ja) * 2011-01-31 2016-02-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 検出値処理装置
CN105408940B (zh) * 2013-07-23 2019-09-03 皇家飞利浦有限公司 混合式(谱/非谱)成像探测器阵列和对应的处理电子设备
JP6479967B2 (ja) * 2014-09-19 2019-03-06 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. スペクトル投影拡張
JP6402282B1 (ja) * 2015-09-16 2018-10-10 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. オブジェクト用のx線画像装置
JP6894928B2 (ja) * 2016-06-07 2021-06-30 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. スペクトル及び非スペクトル物質分解を組み合わせることによる定量的撮像の精度及び分解能の向上

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006503631A (ja) 2002-10-25 2006-02-02 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コンピュータ断層撮影用のダイナミック検出器インターレーシング
US20120045109A1 (en) 2009-05-12 2012-02-23 Koninklijke Philips Electronics N.V. Perfusion imaging
US20140056503A1 (en) 2011-04-28 2014-02-27 Koninklijke Philips N.V. Multi-energy imaging
JP2014036836A (ja) 2012-07-17 2014-02-27 Toshiba Corp X線ct装置およびx線ct装置の制御方法
WO2016147844A1 (ja) 2015-03-19 2016-09-22 株式会社日立製作所 X線ct装置及びマルチエネルギー像作成方法

Also Published As

Publication number Publication date
CN112334950A (zh) 2021-02-05
WO2020002174A9 (en) 2021-01-28
EP3815046A1 (en) 2021-05-05
WO2020002174A1 (en) 2020-01-02
US20210272338A1 (en) 2021-09-02
JP2021529027A (ja) 2021-10-28

Similar Documents

Publication Publication Date Title
JP4401751B2 (ja) アーティファクト低減を容易にする方法及び装置
CN106462985B (zh) 利用锐度驱动的正则化参数的迭代图像重建
US9074986B2 (en) System and method for reducing high density artifacts in computed tomography imaging
US8705822B2 (en) Method for creating images indicating material decomposition in dual energy, dual source helical computed tomography
US7995702B2 (en) System and method of data interpolation in fast kVp switching dual energy CT
JP2007521906A (ja) 動きアーチファクトの補償
US10420519B2 (en) Computed tomography system
US10290122B2 (en) Multi-spectral CT imaging
JP7340549B2 (ja) kVpスイッチングマルチエネルギー取得から非スペクトルボリュメトリック画像データを生成するイメージングシステム
JP2021511875A (ja) スペクトルボリューム画像データを生成するように構成された非スペクトルコンピュータ断層撮影(ct)スキャナ
JP7123927B2 (ja) 交互打ち消しを使用する画像ノイズ推定
US20200342639A1 (en) Single ct backprojector with one geometry calculation per voxel for multiple different types of projection data
CN117425918A (zh) 用于ct重建的方法
EP2864964B1 (en) Image reconstruction in interleaved multi-energy imaging
CN111201452B (zh) 体积图像数据的材料选择性自适应混合
JP2022533583A (ja) 3dプレスキャンボリュメトリック画像データに基づくプロトコル依存2dプレスキャン投影画像
EP3992912A1 (en) Methods and systems for generating a spectral computed tomography image
JP6718448B2 (ja) コンピュータ断層撮影のための画像再構成
US20230320687A1 (en) Controller for a ct scanner
JP2023067768A (ja) 計算機式断層写真法イメージング・システムにおいて二重エネルギ撮像を利用するためのシステム及び方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220425

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230208

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230314

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230529

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: 20230817

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230828

R150 Certificate of patent or registration of utility model

Ref document number: 7340549

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150