JP6243580B2 - フォトンカウンティングctのための心臓再構成 - Google Patents

フォトンカウンティングctのための心臓再構成 Download PDF

Info

Publication number
JP6243580B2
JP6243580B2 JP2017519490A JP2017519490A JP6243580B2 JP 6243580 B2 JP6243580 B2 JP 6243580B2 JP 2017519490 A JP2017519490 A JP 2017519490A JP 2017519490 A JP2017519490 A JP 2017519490A JP 6243580 B2 JP6243580 B2 JP 6243580B2
Authority
JP
Japan
Prior art keywords
weighting
cardiac motion
projection data
flux
flux level
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
JP2017519490A
Other languages
English (en)
Other versions
JP2017529978A (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 JP2017529978A publication Critical patent/JP2017529978A/ja
Application granted granted Critical
Publication of JP6243580B2 publication Critical patent/JP6243580B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/405Source units specially adapted to modify characteristics of the beam during the data acquisition process
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5288Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

以下は、概してフォトンカウンティング検出器による心臓イメージングに関し、コンピュータ断層撮影(CT)への特定の適用で記述される。 しかしながら、以下は他イメージングアプリケーションにも適している。
コンピュータ断層撮影(CT)スキャナは、X線放射線を放射するX線管を含む。 放射線は視野内における被験体又は対象物を横切り、それによって減衰される。 検出器アレイは、視野を横切る放射線を検出して、それを示す信号を生成する。 再構成器は信号を再構成して、ボリュメトリック画像データを生成する。
スペクトルイメージングのために構成されるコンピュータ断層撮影スキャナは、フォトンカウンティング検出器ピクセルを含んでいる。 このようなピクセルにより、生成信号は電流又は電圧を含み、各々は、検出されるフォトンのエネルギーを示すピーク振幅又はピーク高さを有する。 信号は、エネルギー分解されて、その後処理される。 再構成器は、スペクトル及び非スペクトルボリュメトリック画像データを再構成する。
心臓CTは、スキャンの間、心臓運動を仮想的に止めるようにすることが必要とされるので、高速ガントリ回転の間、十分な信号対雑音比を得るため、高いフラックスを必要とする。高いフラックス心臓CTスキャンの場合、検出器のいくつかは、高いフラックスを受ける。 たとえば、患者が小さな減衰のみを有する、患者の周辺及び/又は患者の肺を横切る放射線を受ける検出器は高いフラックスを受ける。
残念なことに、フォトンカウンティング検出器ピクセルは、限られたカウントレートを有する。 結果として、フォトンカウンティング検出器ピクセルで構成されるCTスキャナは、高いフラックス心臓CT調査のために適切でない。
一態様において、本方法は、フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、第一のフラックスレベルは第二のフラックスレベルより低くなり、第一の心臓運動は第二の心臓運動よりも大きくなるように、心周期信号と協調して第一のフラックスレベルと第二の異なるフラックスレベルとの間で放出放射線のフラックスを変調するステップを含む。本方法は、変調される放出放射線を検出するステップと、検出される変調放出放射線を示す投影データを生成するステップを更に含む。本方法は、第一の再構成ウィンドウで投影データを再構成するステップであって、第一のボリュメトリック画像データを生成するために、第一の心臓運動フェーズ及びより低い第一のフラックスレベルに対応する投影データの第一のサブセットに第一の重み付けを適用し、第二の心臓運動フェーズ及びより高い第二のフラックスレベルに対応する投影データの第二のサブセットに第二の異なる重み付けを適用し、第一の重み付けは第二の重み付けより大きくなる、ステップを更に含む。
別の態様において、イメージングシステムは、検査領域の周りを回転し、検査領域を横切る放射線を放出するように構成される放射線源を含む。放射線源コントローラは、フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、第一のフラックスレベルは第二のフラックスレベルより低くなり、第一の心臓運動は第二の心臓運動よりも大きくなるように、心周期信号と協調して第一のフラックスレベルと第二の異なるフラックスレベルとの間で放出放射線のフラックスを変調するために放射線源を制御するように構成される。画素の放射線感受性のアレイは、検査領域を横切る放射線を検出し、検出される放射線を示す信号を生成するように構成される。再構成器は、第一の再構成ウィンドウを用いて投影データを再構成し、第一のボリュメトリック画像データを生成するために、第一の心臓運動フェーズ及びより低い第一のフラックスレベルに対応する投影データの第一のサブセットに第一の重み付けを適用し、第二の心臓運動フェーズ及びより高い第二のフラックスレベルに対応する投影データの第二のサブセットに第二の異なる重み付けを適用し、第一の重み付けは第二の重み付けより大きくなるように構成される。
更に他の態様において、コンピュータ読取り可能な記憶媒体がコンピュータ読取り可能な命令でエンコードされ、プロセッサで実行されるとき、プロセッサに、フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、第一のフラックスレベルは第二のフラックスレベルより低くなり、第一の心臓運動は第二の心臓運動よりも大きくなるように、心周期信号と協調して第一のフラックスレベルと第二の異なるフラックスレベルとの間で放出放射線のフラックスを変調させ、変調される放出放射線を検出する検出器によって生成される投影データを受けさせ、第一の再構成ウィンドウを用いて投影データを再構成させ、それは、第一のボリュメトリック画像データを生成するために、第一の心臓運動フェーズ及びより低い第一のフラックスレベルに対応する投影データの第一のサブセットに第一の重み付けを適用し、第二の心臓運動フェーズ及びより高い第二のフラックスレベルに対応する投影データの第二のサブセットに第二の異なる重み付けを適用し、第一の重み付けは第二の重み付けより大きくなる。
本発明は、様々な構成要素及び構成要素の構成、並びに様々なステップ及びステップの構成の形態を取ることができる。図面は、好ましい実施形態を例示する目的のものにすぎず、本発明を限定するものと解釈されるべきではない。
フォトンカウンティング検出器ピクセルのアレイを備えるイメージングシステムの例を図示する。 ECG信号の例を図示する。 ECG信号例に関連するフラックス変調パターンの例を図示する。 ECG信号例に関連する心臓再構成ウィンドウの例を図示する。 ECG信号例に関連する肺再構成ウィンドウの例を図示する。 ECG信号例に関連する肺再構成ウィンドウの他の例を図示する。 ECG信号例に関連する肺再構成ウィンドウの他の例を図示する。 実施例による方法の例を図示する。
図1は、コンピュータ断層撮影(CT)スキャナのようなイメージングシステム100の例を図示する。 イメージングシステム100は、回転ガントリ102及び静止ガントリ104を含む。 回転ガントリ102は、静止ガントリ104によって回転可能に支持される。 回転ガントリ102は、縦又はz軸について検査領域106のまわりを回転するように構成される。 イメージングシステム100は、スキャンの前に、その間、及び/又はその後に、検査領域106の被験体又は対象物を支持する被験体支持部107を更に含む。 被験体支持部107は、検査領域106に、又はそれから被験体又は対象物をロード及び/又はアンロードするように用いられることもできる。
図の例において、心周期モニタ108と特徴識別子110は、心臓スキャンのためにイメージングシステム100例で使用される。 1つの例において、心周期モニタ108及び/又は特徴識別子110は、イメージングシステム100の部分である。 他の例において、心周期モニタ108及び/又は心周期機能識別子110は、イメージングシステム100から分離される。 心周期モニタ108は被験体上及び/又はその近くに配置されるセンサーから電気信号を受け取って、心電図(ECG又はEKG)のような心臓の電気的活動を示す心周期信号を生成する。
図2は、2つの心周期202及び204を備えるECG信号200の副部分の例を示す。 図2において、第一又はy軸206はミリボルト(mV)の単位の電圧を表し、第二又はx軸208は時間を表す。 心周期202及び204の各々は、心室収縮期の期間210及び心室拡張期の期間212を含む。 通常、心室収縮期の期間210は、心周期202及び204のQRS群216のRピーク214の直後に始まって、T波218の終端の近くで終わる。 収縮期の期間210の間に、心室は収縮して(QRS群216)、それから再分極化する(T波218)。
心室拡張期期間は、U波220(乳頭筋再分極化)及びP波222(心房収縮)を含む。 心室拡張期期間212の間、心臓は弛緩し、心室は循環する血液で充満する。 R-R間隔224は、心周期202及び204のRピーク214に対して、心周期の期間を表す。 領域226は、心臓の運動が心周期の他の領域に対して最小で、「静的フェーズ」と称されている領域に対応する。通常、このような領域は、R-R間隔224におよそ40%、70%などに位置される。
図 1及び2に関して、特徴識別子110は、心周期モニタ108によって生成される心周期信号の所定の特徴を特定するように構成される。 たとえば、特徴識別子110は、特定波(例えば、QRS群216、T波218、U波220、P波222及び/又は他の波)のピークを特定するように構成される。 1つの非限定的な例において、特徴識別子110がRピーク214を特定するように構成される。これは、電圧レベル閾値化及び/又は他のアプローチを通じて実現されることができる。特定されると、特徴識別子110は心周期の所定の特徴が特定されたことを示すトリガー信号を生成する。
図1に関して、イメージングシステム100は回転ガントリ102によって回転可能に支持されるX線管のような放射線源112を更に含む。 放射線源112は検査領域106のまわりで回転ガントリ102により回転して、検査領域106を横切る放射線を発生させて、放射するように構成される。 イメージングシステム100は、放射線源コントローラ114を更に含む。 放射線源コントローラ114は、生成される放射線のフラックスを変調するように構成される。 たとえば、放射線コントローラ114は、フラックスを変調するために、放射線源112のカソード熱電流を選択的に変化させ、放射線源112の電子の流れを妨げるようにチャージし、放射される放射線をフィルタリングすることなどができる。 図の例において、放射線源コントローラ114は、所定の変調パターンに基づいてフラックスを変調する。
図3は、ECG信号200に関連する変調パターン300を示す。図3において、第一又はy軸302は、ミリアンペア(mA)の単位の電流を表す。 図示される変調パターン300は、x-(時間)軸208をECG信号200と共有する。 この例において、変調パターン300は、より高いフラックスレベル304及びより低い(非ゼロの)フラックスレベル306の間の生成される放射線のフラックスをスイッチする(切り替える)。 更に、変調パターン300は、静的フェーズ226に対応する時間ウィンドウ308のための生成される放射線のフラックスをより高いフラックスレベル304及び静的フェーズ226の外部のより低いフラックスレベル306にスイッチする。他の実施例において、変調パターン300は、異なる及び/又はさらなるフェーズのために、個別及び/又は連続態様等で、二つよりも多くの異なるフラックスレベルの間でスイッチすることができる。バリエーションにおいて、変調パターンは、管電流及びそれゆえにフラックスはゼロである一つ又はそれより多くの領域も含む。 これは、ビーム経路において、又はそこから位置されることができる管及び/又は物理的なフィルタをスイッチングするグリッドで実現されることができる。
図1に戻って、イメージングシステム100は放射線感受性検出器ピクセル116の1又は2次元アレイ115を更に含む。 ピクセル116は、検査領域106の間に放射線源112に対向して位置され、検査領域106を横切る放射線を検出して、それを示す電気信号(投影データ)を生成する。 図の例において、ピクセル116は直接変換フォトンカウンティング検出器ピクセルを含む。 このようなピクセルにより、生成信号は、検出されるフォトンのエネルギーを示すピーク振幅又はピーク高さを有する電流又は電圧を含む。 直接変換フォトンカウンティング検出器ピクセルは、CdTe、CdZnTe、Si、Ge、GaAs又は他の直接変換物質のような何れかの適切な直接変換物質も含む。
イメージングシステム100は、検出器ピクセル116によって出力される電気信号を処理して、検出されるフォトンのエネルギーを示す電圧のようなパルス又は他のパルスを生成するパルスシェーパー118を更に含む。 バリエーションにおいて、検出器ピクセル116によって出力される電気信号は前処理(例えば、増幅等)され、パルスシェーパー118は前処理される信号を処理する。 イメージングシステム100は、エネルギーがパルスを識別するエネルギー識別器120を更に含む。 図の例において、エネルギー識別器120は、パルスの振幅を関心エネルギーに対応する少なくとも1つのエネルギー閾値と比較する、少なくとも1つのコンパレーター122を含む。コンパレーター122は、検出されるフォトンのエネルギーが閾値の上か、又は下かを示す出力信号を生成する。
イメージングシステム100は、各々の閾値のためのカウント値をインクリメント(又はデクリメント)するカウンター124を更に含む。 たとえば、特定の閾値のためのコンパレーター122の出力は、パルスの振幅が対応する閾値を上回ることを示すとき、その閾値のためのカウント値はインクリメントされる。 イメージングシステム100は、異なるエネルギー範囲に対応する、カウントされるパルスをエネルギービンに割り当てるビナー126を更に含む。たとえば、ビンは2つの閾値の間のエネルギー範囲のために規定される。 この例の場合、ビナー126は、より高い閾値のためでなく低い閾値のためのカウントをもたらすフォトンを、2つの閾値の間のエネルギー範囲のために規定されるビンに割り当てる。
イメージングシステム100は、スペクトル及び/又は非スペクトルアルゴリズムを使用してビニングされるデータを再構成して、スペクトル及び/又は非スペクトルボリュメトリック画像データを生成する再構成器128を更に含む。 再構成器128は、フィルタ補正逆投影法、反復法及び/又は他の再構成アルゴリズムを使用することができる。 これは、重み付け及び非重み付け再構成アルゴリズムを含む。 以下に詳細に説明するように、重み付け再構成アルゴリズム例は、フラックスがより高いフラックスレベル304にあるときに収集される投影データのために、より大きく投影データを重み付けする。 別の例において、重み付け再構成アルゴリズムは、フラックスがより低いフラックスレベル306にあるときに収集される投影データのために、より大きく投影データを重み付けする。
このようなアルゴリズムは、図3に示される変調パターン300及び/又はそれが心周期の静的フェーズの間、高いレベルにあるようにフラックスを制御する他のフラックス変調パターンを使用する心臓調査によく適している。たとえば、フラックスが静的フェーズ226のためにより高いフラックスレベル304にあるときに収集される投影データのために、より大きく投影データを重み付けする重み付け再構成アルゴリズムは心臓のボリュメトリック画像データを生成するように使用されることができ、フラックスが他のフェーズのためにより低いフラックスレベル306にあるときに収集される投影データのために、より大きく投影データを重み付けする重み付け再構成アルゴリズムは、肺(及び随意に被験体の周辺)のボリュメトリック画像データを生成するのに用いられることができる。
これは、以下に記載されているように、実現されることができる。 スキャンの前に、管電流変調パターンは決定される。 たとえば、パターンは、最後のRピーク後の300ミリ秒(300ミリ秒)電流を強化し、更に300ミリ秒後に、再び電流を減らす。このパターンは、特定の長さのほぼ安定した鼓動を仮定する。 再構成のために、画像は心周期の70%で再構成される。 より高いフラックスレベル304の期間は、再構成のために必要なゲーティングウィンドウ幅より長くなりうる。 その場合、フラックスがより高いフラックスレベル304にあるとき、重み付けは必ずしも高くなる必要はない。 しかしながら、フラックスがより高いフラックスレベル304にあるとき、重み付けは常に高くなることができる。
このように、心臓のボリュメトリック画像データは、診断目的のために適切な信号対雑音比を有する。 更に、肺のボリュメトリック画像データは、ピクセル116の放射線感度の高いアレイのフォトンカウンティング検出器ピクセルの限られたカウントレートによってもたらされるアーチファクトを含まない。更に、被験体ドーズは、フラックス304のより高いレベルで完全にスキャンされることに対して減らされる。
心臓のボリュメトリック画像データ及び肺のボリュメトリック画像データは、ディスプレイスクリーンの、異なる、オーバラップする、又は同じディスプレイウィンドウにおいて並列に視覚的に表示されることができる。更に、心臓のボリュメトリック画像データ及び肺のボリュメトリック画像データは、代わりに、視覚的に表示されることができる。 画像プロセッサ130は、単一のボリュメトリック画像データセットを形成するために、心臓のボリュメトリック画像データ及び肺のボリュメトリック画像データを結合するように用いられることができる。このために、1つの例において、心臓又は肺は、セグメント化され、肺のボリュメトリック画像データ又は心臓のボリュメトリック画像データと結合されることができる。
イメージングシステム100は、オペレーターコンソール132としての役割を果たし、ディスプレイのような出力装置及びキーボード、マウス等のような入力装置を含む計算システムを更に含む。コンソール132に常駐するソフトウェアは、ユーザーがシステム100の操作を制御することを可能にする。 たとえば、コンソール132は、フラックスを変調し、特定の再構成アルゴリズムを使用するために、特定されるイメージングプロトコルに基づいて、イメージングプロトコル(例えば、フラックス変調及び重み付け再構成による心臓イメージングプロトコル)を特定し、スキャンを開始し、放射線源コントローラ114及び再構成器128を制御するように用いられることができる。
心臓モニタ108、特徴識別子110、パルスシェーパー118、識別器120、コンパレーター122、カウンター124、ビナー126、再構成器128、画像プロセッサ130及びコンソール132が電気回路及び/又はソフトウェア/ファームウェアによって実施されることができることは評価されるべきである。 たとえば、コンソール132は一つ又はそれより多くのプロセッサ(例えば、中央演算処理ユニット、マイクロプロセッサなど)及び一つ又はそれより多くのプロセッサで実行されるとき、一つ又はそれより多くのプロセッサに、ここに記述される機能の一つ又はそれより多くを実行させる、コンピュータ実行可能な命令でエンコードされる(一時的なメモリを除外し、物理メモリを含む)コンピュータ読取り可能な媒体を含む。
上で簡潔に述べられるように、再構成器128は重み付け再構成アルゴリズムを使用することができる。 以下は、重み付け再構成アルゴリズムの例を記述する。 概して、各々の再構成ウィンドウの幅は、全てのボクセルが少なくとも180度(180°)プラスファン角度に対して(又は平行コーンビームリビニングへのファンビームの後の180°に対して)源によって放射されるように、再構成の完全性条件によって決定される。
図4は、ECG 200及び変調パターン300に関連して再構成ウィンドウ400を用いる再構成アルゴリズム例を示す。 この例において、再構成ウィンドウ400は、心臓のための心臓再構成ウィンドウである。 この例において、心臓再構成ウィンドウ400は、静的フェーズ226に対応し、放射線源電流がより高いフラックスレベル304で変調される期間内にある。他の実施例において、心臓再構成ウィンドウ400は、より狭くなるか、又はより広くなることができ、放射線源電流がより低いフラックスレベル306で変調される時間をカバーすることを含む。
図 5、6及び7はそれぞれ、ECG 200及び変調パターン300に関連して再構成ウィンドウ500、600及び700を用いる再構成アルゴリズム例を示す。これらの図において、再構成ウィンドウ500、600及び700は、肺のための呼吸再構成ウィンドウである。
図5において、呼吸再構成ウィンドウ500は、放射線源電流がより低いフラックスレベル306で変調される全時間をカバーする。 図6において、呼吸再構成ウィンドウ600は、放射線源電流がより低いフラックスレベル306で変調される時間の部分のみをカバーする。 図7において、呼吸再構成ウィンドウ700は、心臓再構成ウィンドウ400(図4)のすぐ近くにあるが、放射線源電流がより低いフラックスレベル306で変調される場合だけである。
結果として、図7の呼吸再構成ウィンドウ700の場合、心臓及び肺画像の動き状態は最小の差を有する。 呼吸再構成ウィンドウ700は、(図7で示されるように)心臓再構成ウィンドウ400の前又は後の何れかに位置される。バリエーションにおいて、呼吸再構成ウィンドウ700は、心臓再構成ウィンドウ400の前の部分及び後の部分に分割され得る。
非反復再構成の非限定的な例は、「コーン-ビームCTのためのアパーチャ重み付け心臓再構成」(Koken, 物理医療生物学51 (2006) 3433-3448)に記述される。 反復再構成の非限定的な例は、「ARTを使う心臓コーン-ビームCTボリューム再構成」(ニールセン、医療物理、32(4) (2005) 851-860)に記述される。 Kokenから、適切な再構成が、等式 1で示されるように公式化されることができる。
Figure 0006243580
ここで
Figure 0006243580
はオブジェクトポイントを表し、
Figure 0006243580
はオブジェクト関数を表し、l及び
Figure 0006243580
はインデックスを表し、
Figure 0006243580
は重み付け関数を表し、P4(φ、u、h)は前処理される投影データを表し、
Figure 0006243580
はxy-平面に投影される平行線の角度を表し、
Figure 0006243580
及び
Figure 0006243580
はオブジェクトポイントの投影される検出器座標を表す。
心臓再構成に対して、
Figure 0006243580
となり、ここでωcは心臓重み付けを表し、ωapはアパーチャ重み付けを表し、
Figure 0006243580
となる。肺再構成に対して、
Figure 0006243580
となり、ここで全ての
Figure 0006243580
に対して
Figure 0006243580
になり、さもなければ0になる。
他の実施例において、検出器へのX線フラックスが検出器のカウントレート限界を示す所定の閾値を上回らない、すべての投影に対して
Figure 0006243580
になり、検出器へのX線フラックスが検出器のカウントレート限界を示す所定の閾値を上回る、すべての投影に対して
Figure 0006243580
になる。この実施例は、過度のフラックスレベルがハードウェア、ソフトウェア、及び/又はその他を通じて検出器自体の範囲内で検出されることができる構成によく適している。
他の実施例において、
Figure 0006243580
は境界上のcos2スムージング、心臓フェーズポイントまでの最小時間距離又はデータの180度プラスファン角度(それより大きくない)を保証する時間幅を含む。
Figure 0006243580
となり、ここで
Figure 0006243580
は全ての
Figure 0006243580
に対して1になり、さもなければ0になる。この例において、Isatは検出器が飽和する強度であり、ここに記憶されるように、
Figure 0006243580
は重み付け関数である。
図8は、ここに記述される実施例による方法を例示する。
これらの方法の操作の順序は限定的なものではないことを理解されたい。従って、本明細書では、他の順序も企図される。更に、1つ若しくは複数の操作が省かれてもよく、及び/又は1つ若しくは複数の追加の操作が含まれてもよい。
802において、患者はフォトンカウンティング検出器ピクセルを備えるイメージングシステムの被験体支持部にロードされ、スキャンのためにイメージング領域に関して位置される。
804において、放射線放出フラックスは少なくとも第一のフラックスレベル及び第二のフラックスレベルの間で変調され、第一のフラックスレベルは第二のフラックスレベルより高く、静的フェーズの間の放射線放出フラックスがより高いフラックスなるように、変調は患者の心周期の静的フェーズと協調する。
806において、スキャンが実行され、投影データが生成される。 スキャンは環状、ヘリカルステップのためのゲーティングされる再構成、(例えば、軸の)シュート、及び/又は、他はスキャンになり得ることは評価される。
808において、投影データは、静的フェーズの外側の一つ又はそれより多くのフェーズに対応する投影データの第二の部分より高い、静的フェーズに対応する投影データの第一の部分を重み付けして心臓ボリュメトリック画像データを生成する心臓重み付け再構成アルゴリズムで再構成される。
810において、投影データは、より高い静的フェーズに対応する投影データの第一の部分より高い静的フェーズの外側の一つ又はそれより多くのフェーズに対応する投影データの第二の部分を重み付けして呼吸ボリュメトリック画像データを生成する肺重み付け再構成アルゴリズムで再構成される。
812において、ここに、及び/又は他において議論されるように、心臓ボリュメトリック画像データ及び呼吸ボリュメトリック画像データが表示される。
ここに議論される少なくとも部分の方法は、コンピュータ可読記憶媒体(即ち物理メモリ及び他の非一時的媒体)にエンコードされる又は埋め込まれたコンピュータ可読命令によって実施されてよく、それらのコンピュータ可読命令は、コンピュータプロセッサによって実行されるときに、プロセッサに上述した操作を実施させる。追加又は代替として、コンピュータ可読命令の少なくとも1つは、信号、搬送波、及び他の一時的な媒体によって搬送される。
好ましい実施形態を参照して本発明を述べてきた。上記の詳細な説明を読んで理解すれば、修正形態及び変形形態が想到されよう。本発明は、添付の特許請求の範囲又はその均等物の範囲内にある限り、全てのこのような修正形態及び変形形態を含むものと解釈されることが意図される。

Claims (21)

  1. フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、前記第一のフラックスレベルは前記第二のフラックスレベルより低くなり、前記第一の心臓運動は前記第二の心臓運動よりも大きくなるように、心周期信号と協調して前記第一のフラックスレベルと前記第二の異なるフラックスレベルとの間で放出放射線のフラックスを変調するステップと、
    前記変調放出放射線を検出するステップと、
    前記検出される変調放出放射線を示す投影データを生成するステップと、
    第一の再構成ウィンドウを用いて前記投影データを再構成するステップであって、第一のボリュメトリック画像データを生成するために、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの第一のサブセットに第一の重み付けを適用し、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの第二のサブセットに第二の異なる重み付けを適用し、前記第一の重み付けは前記第二の重み付けより大きくなる、ステップと
    を有する、方法。
  2. 前記第一の重み付けは非ゼロ値を有し、前記第二の重み付けはゼロの値を有する、請求項1に記載の方法。
  3. 前記第一の心臓運動フェーズに対応する前記投影データの全セットに前記第一の重み付けを適用するステップ
    を更に有する、請求項1又は2に記載の方法。
  4. 前記第一の心臓運動フェーズに対応する前記投影データの第一のサブセットに前記第一の重み付けを適用するステップ
    を更に有する、請求項1又は2に記載の方法。
  5. 第二の再構成ウィンドウを用いて前記投影データを再構成するステップであって、第二のボリュメトリック画像データを生成するために、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの前記第二のサブセットに第三の重み付けを適用し、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの前記第一のサブセットに第四の重み付けを適用し、前記第三の重み付けは前記第四の重み付けより大きくなる、ステップ
    を更に有する、請求項3又は4に記載の方法。
  6. 1に等しい前記第一の重み付けを、前記検出器のカウントレート限界を示す所定の閾値を超えない、前記検出器に照射される前記X線フラックスに応答して前記投影データに適用するステップと、
    ゼロに等しい前記第二の重み付けを、前記検出器のカウントレート限界を示す前記所定の閾値を超える、前記検出器に照射される前記X線フラックスに応答して前記投影データに適用するステップと
    を更に有する、請求項1に記載の方法。
  7. 第二の再構成ウィンドウで前記投影データを再構成するステップであって、前記第二のボリュメトリック画像データを生成するために、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの前記第二のサブセットに第三の重み付けを適用し、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの前記第一のサブセットに第四の重み付けを適用し、前記第三の重み付けは前記第四の重み付けより大きくなる、ステップ
    を更に有する、請求項2に記載の方法。
  8. 前記第二の心臓運動フェーズに隣接する一つ又はそれより多くの心臓運動フェーズに対応する前記投影データに前記第一の重み付けを適用するステップ
    を更に有する、請求項7に記載の方法。
  9. 前記一つ又はそれより多くの心臓運動フェーズは前記第二の心臓運動フェーズの前に発生する、請求項8に記載の方法。
  10. 前記一つ又はそれより多くの心臓運動フェーズは前記第二の心臓運動フェーズの後に発生する、請求項8に記載の方法。
  11. 前記一つ又はそれより多くの心臓運動フェーズの前記第一のサブセットは前記第二の心臓運動フェーズの前に発生し、前記一つ又はそれより多くの心臓運動フェーズの前記第二のサブセットは前記第二の心臓運動フェーズの後に発生する、請求項8に記載の方法。
  12. 検査領域の周りを回転し、前記検査領域を横切る放射線を放出するように構成される放射線源と、
    フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、前記第一のフラックスレベルは前記第二のフラックスレベルより低くなり、前記第一の心臓運動は前記第二の心臓運動よりも大きくなるように、心周期信号と協調して前記第一のフラックスレベルと前記第二の異なるフラックスレベルとの間で放出放射線の前記フラックスを変調するために前記放射線源を制御するように構成される放射線源コントローラと、
    前記検査領域を横切る放射線を検出し、前記検出される放射線を示す信号を生成するように構成される、画素の放射線感受性のアレイと、
    第一の再構成ウィンドウを用いて前記投影データを再構成するように構成される再構成器であって、前記第一のボリュメトリック画像データを生成するために、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの第一のサブセットに第一の重み付けを適用し、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの第二のサブセットに第二の異なる重み付けを適用し、前記第一の重み付けは前記第二の重み付けより大きくなる、再構成器と
    を有する、イメージングシステム。
  13. 前記第一の重み付けは非ゼロ値を有し、前記第二の重み付けはゼロである、請求項12に記載のイメージングシステム。
  14. 前記第一の重み付けは前記第一の心臓運動フェーズに対応する前記投影データの全セットに適用される、請求項12又は13に記載のイメージングシステム。
  15. 前記第一の重み付けは前記第一の心臓運動フェーズに対応する前記投影データの第一のサブセットに適用される、請求項12又は13に記載のイメージングシステム。
  16. 前記再構成器は、第二の再構成ウィンドウを用いて前記投影データを再構成するように更に構成され、第二のボリュメトリック画像データを生成するために、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの前記第二のサブセットに第三の重み付けを適用し、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの前記第一のサブセットに第四の重み付けを適用し、前記第三の重み付けは前記第四の重み付けより大きくなる、請求項13又は14に記載のイメージングシステム。
  17. 前記再構成器は、第二の再構成ウィンドウを用いて前記投影データを再構成するように更に構成され、第二のボリュメトリック画像データを生成するために、前記第二の心臓運動フェーズに対応する前記投影データの前記第二のサブセットに第三の重み付けを適用し、前記第一の心臓運動フェーズに対応する前記投影データの前記第一のサブセットに第四の重み付けを適用し、前記第三の重み付けは前記第四の重み付けより大きくなる、請求項13に記載のイメージングシステム。
  18. 前記第一の重み付けは前記第二の心臓運動フェーズに隣接する一つ又はそれより多くの心臓運動フェーズに対応する前記投影データに適用される、請求項17に記載のイメージングシステム。
  19. 前記一つ又はそれより多くの心臓運動フェーズは前記第二の心臓運動フェーズの前又は後の少なくとも一つに発生する、請求項18に記載のイメージングシステム。
  20. 前記一つ又はそれより多くの心臓運動フェーズの前記第一のサブセットは前記第一の心臓運動フェーズの前に発生し、前記一つ又はそれより多くの心臓運動フェーズの前記第二のサブセットは前記第一の心臓運動フェーズの後に発生する、請求項18に記載のイメージングシステム。
  21. コンピュータ読取り可能な命令でエンコードされるコンピュータ読取り可能な記憶媒体であって、プロセッサによって実行されるとき、前記プロセッサに、フラックスが、第一の心臓運動を有する第一の心臓運動フェーズの間、第一のフラックスレベルにあり、第二の心臓運動を有する第二の心臓運動フェーズの間、第二のフラックスレベルにあり、前記第一のフラックスレベルは前記第二のフラックスレベルより低くなり、前記第一の心臓運動は前記第二の心臓運動よりも大きくなるように、心周期信号と協調して前記第一のフラックスレベルと前記第二の異なるフラックスレベルとの間で放出放射線のフラックスを変調させ、
    前記変調される放出放射線を検出する検出器によって生成される投影データを受けさせ、
    第一の再構成ウィンドウを用いて投影データを再構成させ、第一のボリュメトリック画像データを生成するために、前記第一の心臓運動フェーズ及び前記より低い第一のフラックスレベルに対応する前記投影データの第一のサブセットに第一の重み付けを適用し、前記第二の心臓運動フェーズ及び前記より高い第二のフラックスレベルに対応する前記投影データの第二のサブセットに第二の異なる重み付けを適用し、前記第一の重み付けは前記第二の重み付けより大きくなる、コンピュータ読取り可能な記憶媒体。
JP2017519490A 2014-10-20 2015-10-12 フォトンカウンティングctのための心臓再構成 Active JP6243580B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201462065939P 2014-10-20 2014-10-20
US62/065,939 2014-10-20
PCT/IB2015/057791 WO2016063171A1 (en) 2014-10-20 2015-10-12 Cardiac reconstruction for photon counting ct for heart and lung images

Publications (2)

Publication Number Publication Date
JP2017529978A JP2017529978A (ja) 2017-10-12
JP6243580B2 true JP6243580B2 (ja) 2017-12-06

Family

ID=54396923

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017519490A Active JP6243580B2 (ja) 2014-10-20 2015-10-12 フォトンカウンティングctのための心臓再構成

Country Status (5)

Country Link
US (1) US9775580B2 (ja)
EP (1) EP3209214B1 (ja)
JP (1) JP6243580B2 (ja)
CN (1) CN107072627B (ja)
WO (1) WO2016063171A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11200709B2 (en) * 2016-12-27 2021-12-14 Canon Medical Systems Corporation Radiation image diagnostic apparatus and medical image processing apparatus
JP7066359B2 (ja) * 2017-09-19 2022-05-13 キヤノンメディカルシステムズ株式会社 X線ct装置
CN108065952B (zh) * 2018-01-03 2021-08-13 东软医疗系统股份有限公司 一种剂量调制方法及装置
EP3745959B1 (en) * 2018-01-31 2023-12-20 Koninklijke Philips N.V. Non-spectral computed tomography (ct) scanner configured to generate spectral volumetric image data

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6931094B2 (en) * 2003-10-10 2005-08-16 Ge Medical Systems Global Technology Company, Llc Methods and systems for smoothing
JP2006101926A (ja) * 2004-09-30 2006-04-20 M & C:Kk 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法
WO2006067671A2 (en) * 2004-12-22 2006-06-29 Philips Intellectual Property & Standards Gmbh Method and apparatus for cardiac computed tomography
EP2238906B1 (en) * 2005-09-07 2015-04-08 Kabushiki Kaisha Toshiba X-ray computed tomography apparatus
JP4945203B2 (ja) * 2005-09-07 2012-06-06 株式会社東芝 X線コンピュータ断層撮影装置
US7532702B2 (en) 2005-11-23 2009-05-12 General Electric Company Method and system for performing CT image reconstruction with motion artifact correction
US20070147579A1 (en) * 2005-12-23 2007-06-28 De Man Bruno K B Method and system for radiographic imaging with organ-based radiation profile prescription
WO2008021661A2 (en) * 2006-08-01 2008-02-21 Koninklijke Philips Electronics, N.V. Stereo tube computed tomography
WO2010070583A1 (en) * 2008-12-17 2010-06-24 Koninklijke Philips Electronics N.V. X-ray examination device and method
DE102008063311A1 (de) * 2008-12-30 2010-07-08 Siemens Aktiengesellschaft Verfahren zur Rauschreduktion von CT-Bilddaten und Bildbearbeitungssystem
JP2011056167A (ja) * 2009-09-14 2011-03-24 Toshiba Corp 放射線画像撮影装置及び放射線画像撮影方法
JP2011067333A (ja) * 2009-09-25 2011-04-07 Fujifilm Corp 放射線画像撮影装置及び撮影制御装置
EP2372703A1 (en) * 2010-03-11 2011-10-05 Fraunhofer-Gesellschaft zur Förderung der Angewandten Forschung e.V. Signal processor, window provider, encoded media signal, method for processing a signal and method for providing a window
JP5869001B2 (ja) * 2011-01-31 2016-02-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 検出値処理装置
WO2013076613A1 (en) * 2011-11-23 2013-05-30 Koninklijke Philips Electronics N.V. Image domain de-noising
CN104220899A (zh) 2012-03-27 2014-12-17 皇家飞利浦有限公司 利用具有光子计数探测器的成像系统的常规成像
JP6261915B2 (ja) * 2012-08-30 2018-01-17 東芝メディカルシステムズ株式会社 X線ct装置、画像処理装置、及び画像処理方法
US9111337B2 (en) * 2012-10-12 2015-08-18 Mayo Foundation For Medical Education And Research Low dose cardiac CT imaging with time-adaptive filtration
JP6329490B2 (ja) * 2013-02-05 2018-05-23 株式会社日立製作所 X線ct装置及び画像再構成方法
US20140236488A1 (en) * 2013-02-19 2014-08-21 Mindways Software, Inc. Computed Tomography Calibration Systems and Methods
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
KR102205353B1 (ko) * 2013-09-26 2021-01-20 삼성전자주식회사 엑스선 촬영 장치 및 그 제어 방법
US9974495B2 (en) * 2014-01-20 2018-05-22 Hitachi, Ltd. X-ray CT apparatus, image processing device, and image reconstruction method

Also Published As

Publication number Publication date
US9775580B2 (en) 2017-10-03
WO2016063171A1 (en) 2016-04-28
US20170231590A1 (en) 2017-08-17
EP3209214A1 (en) 2017-08-30
CN107072627B (zh) 2020-10-16
EP3209214B1 (en) 2019-06-05
CN107072627A (zh) 2017-08-18
JP2017529978A (ja) 2017-10-12

Similar Documents

Publication Publication Date Title
US9848843B2 (en) Dynamic perfusion imaging
US6370217B1 (en) Volumetric computed tomography system for cardiac imaging
EP1522045B1 (en) Motion artifact correction of tomographical images
JP5551436B2 (ja) ステレオ管コンピュータ断層撮像
US9235907B2 (en) System and method for partial scan artifact reduction in myocardial CT perfusion
JP6243580B2 (ja) フォトンカウンティングctのための心臓再構成
JP2009545395A5 (ja)
US10420519B2 (en) Computed tomography system
KR102205353B1 (ko) 엑스선 촬영 장치 및 그 제어 방법
US11179122B2 (en) Bolus imaging
KR20190079371A (ko) 조영제를 주입하여 컴퓨터 단층 촬영하는 방법 및 장치
JP7159359B2 (ja) 可動容積についての放射性画像を改善するシステム
US20220071578A1 (en) Improved method of acquiring a radiographic scan of a region-of-interest in a metal containing object
US10984564B2 (en) Image noise estimation using alternating negation
WO2018087049A1 (en) Dose reduction in dynamic radiography
US20220180576A1 (en) Cone beam artifact correction for gated imaging
KR20180063753A (ko) 의료 영상 장치 및 동작 방법
Sawall et al. Low-dose phase-correlated cone-beam micro-CT of small animals
Sawall et al. Functional phase-correlated micro-CT imaging of small rodents with low dose

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170411

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20170411

TRDD Decision of grant or rejection written
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20171016

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20171102

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171109

R150 Certificate of patent or registration of utility model

Ref document number: 6243580

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250