JP2006517825A - ボリュメトリック画像再構成のための方法及び装置 - Google Patents

ボリュメトリック画像再構成のための方法及び装置 Download PDF

Info

Publication number
JP2006517825A
JP2006517825A JP2006502445A JP2006502445A JP2006517825A JP 2006517825 A JP2006517825 A JP 2006517825A JP 2006502445 A JP2006502445 A JP 2006502445A JP 2006502445 A JP2006502445 A JP 2006502445A JP 2006517825 A JP2006517825 A JP 2006517825A
Authority
JP
Japan
Prior art keywords
projections
rebinned
marker
labeled
detector
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.)
Granted
Application number
JP2006502445A
Other languages
English (en)
Other versions
JP4431136B2 (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 Electronics 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 Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of JP2006517825A publication Critical patent/JP2006517825A/ja
Application granted granted Critical
Publication of JP4431136B2 publication Critical patent/JP4431136B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

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/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by 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/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • 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
    • 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/027Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/412Dynamic
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00X-ray or gamma ray systems or devices
    • Y10S378/901Computer tomography program or processor

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pulmonology (AREA)
  • Physiology (AREA)
  • General Physics & Mathematics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Generation (AREA)

Abstract

コンピュータトモグラフィデータの再構成は、多くの場合、フィルタリングされた逆投影に先立つ投影データの平行リビニングステップに基づく。心臓再構成の場合、ゲーティング技法が、時間依存の画像再構成用の適当な投影を選択するために適用される。通常、ゲーティングウィンドウは、ファンビーム/コーンビーム投影だけの関数である。本発明によれば、個々の投影が検出された検出器表面上の位置に対応する関数が、個々の投影が本当に選択されたゲーティングウィンドウ内にあるか否かを更に検証するために使用される。

Description

本発明は、画像再構成の分野に関する。本発明は、関心対象に対して放射線のコーンビーム又はファンビームを移動させるコンピュータトモグラフィ(CT)スキャナからのX線伝送データを再構成することに関して特定のアプリケーションを見いだし、特にそれに関して記述される。しかしながら、本発明は、CTスキャナ、核カメラ及び非平行な軌跡を表すデータを処理する他の診断スキャナからのデータの再構成に関してもアプリケーションを見いだすことが認められるべきである。特に、本発明は、ボリュメトリックな画像の再構成方法と、ボリュメトリックな画像再構成のための画像処理装置と、ボリュメトリックな画像再構成のために使用できるコンピュータ読み取り可能な媒体を含むコンピュータプログラム製品とに関する。
従来、スパイラル又はヘリカルCTスキャナは、透過放射線の薄いスライス又はビームを放出するX線源を有する。X線源は、回転軸に沿って平行移動される対象(object)の周りを回転運動するように取り付けられる。放射線検出器の円弧又はリングは、対象を通り抜けた放射線を受け取る。放射線検出器からのデータは、対象を通じた単一のらせん状スライスを表す。検出器からのデータは、3次元画像表現に再構成される。
CTデータの再構成は、多くの場合、フィルタリングされた逆投影に先立つ投影データの平行リビニングステップに基づく。これは、シングルライン及びコーンビーム投影データにそれぞれ当てはまる。運動する対象、例えば人間の心臓から取得されるCTデータの再構成の場合、状況は、心臓がほぼ周期的な運動を行うという事実によって複雑である。心臓のシングルフェーズのより高い解像度の再構成を実施するために、ゲーティング方法が適用される。ゲーティングウィンドウは、一般に、心電図(ECG)から決定され、利用可能な投影データの量を、関心のある心臓フェーズに属するものに制限する。
これまで知られているすべてのケースにおいて、ゲーティングウィンドウは、再構成のために使用される平行リビニングされた投影をラベリングするために使用される。
平行リビニング自体は、さまざまな異なる放射線源位置において取得される複数のファンビーム/コーンビーム投影から、平行な投影を生成する。従って、複数のファンビーム/コーンビーム投影が、さまざまな異なる時点に取得される。平行リビニングされた投影は、扇角の角度レンジ内において測定された一組のファンビーム/コーンビーム投影からの情報を含む。同じ投影角において取得されるファンビーム投影と同じ時点に測定されるのは、平行リビニングの平行投影の中心線だけである。従って、ECGに基づくゲーティングウィンドウから生じる値によって、平行リビニングされた投影を直接的にラベリングすることは、ゲーティングウィンドウに属さなかった投影データを再構成に使用することになり、従って、画像の不所望の動きアーチファクトをもたらす。
本発明の目的は、画像内の不所望の動きアーチファクトを低減することである。
本発明の例示的な実施形態によれば、上述の目的は、対象の複数の投影からボリュメトリックな画像を再構成する方法であって、複数の投影が、例えばヘリカル走査パスのような走査パスに沿って対象に対して移動される放射線源及び対応する検出器によって決定され、方法が、複数のリビニングされた投影を決定するために複数の投影についてリビニングステップを実施するステップと、複数のリビニングされた投影の個々のものが、ゲーティングウィンドウ内にあるかどうかを、複数のリビニングされた投影の個々のものが検出された検出器上の位置に関する位置情報に基づいて、決定するステップと、を含む方法によって解決されることができる。有利には、本発明のこの例示的な実施形態によれば、ボリュメトリック画像再構成ののち、画像内の不所望の動きアーチファクトは低減される。このように、例えば、高時間分解能の心臓CT再構成が実施されるとき、分解能が改善されることができる。
請求項2に記載の本発明の例示的な実施形態によれば、複数のリビニングされた投影の個々のものが検出された検出器上の位置に依存する関数が決定され、関数は、複数の投影のラベリングを、複数のリビニングされた投影のリビニングされた投影空間に移す。この関数は、リビニングされた投影の複数の個々のものに対応するボクセルが、本当にゲーティングウィンドウ内にあるか否かを検証するために使用される。有利には、本発明のこの例示的な実施形態は、ボクセルがゲーティングウィンドウ内にあり、最終的な画像を再構成するために使用されるべきか否かの検証される信頼できる決定を可能にする。
請求項3に記載の本発明の別の例示的な実施形態によれば、重み付け、フィルタリング及び重み付けの正規化が、ボクセルがリビニングされた投影空間の許容される領域上に投影されるかどうかを決定するために、実行される。
請求項4に記載の本発明の別の例示的な実施形態によれば、複数の投影が、CTスキャナによって決定され、方法が、高時間分解能のCT再構成に関し、ゲーティングウィンドウが心電図から決定される。代替例として、アコースティック又は光学センサが、ECGの代わりに使用されてもよい。
請求項5に記載の本発明の別の例示的な実施形態によれば、リビニングステップが、平行リビニング及び放射線源の移動を含み、検出器が、ヘリカルパスに沿う。
請求項6に記載の本発明の別の例示的な実施形態によれば、計算ユニットを有する画像処理装置であって、計算ユニットが、複数のリビニングされた投影を決定するために複数の投影についてリビニングを実施するように構成され、計算ユニットが、更に、複数のリビニングされた投影の個々のものがゲーティングウィンドウ内にあるかどうかを、複数のリビニングされた投影の個々のものが検出された検出器上の位置に関する位置情報に基づいて決定するように構成される、画像処理装置が提供される。
請求項6の画像処理装置の他の例示的な実施形態が、請求項7乃至9に与えられている。
請求項10に記載の本発明の別の例示的な実施形態によれば、リビニングステップと、複数のリビニングされた投影の個々のものがゲーティングウィンドウ内にあるかどうかについて位置情報に基づく決定と、を、コンピュータプログラムがコンピュータ上で実行されるとき、コンピュータに実行させるためのコンピュータプログラム手段を有するコンピュータ読み取り可能な媒体を具えるコンピュータプログラム製品が提供される。有利には、この例示的な実施形態によるコンピュータプログラムは、計算の最小限の量が実施されることを必要とし、従って、必要とされる計算パワーの量を低減する。更に、コンピュータプログラムを実行するのに必要なストレージが最小化されることができる。
特定の投影が検出された検出器上の位置に対応する位置情報が、この投影が本当にゲーティングウィンドウ内にあるかどうかを決定するために使用されることが、本発明の例示的な実施形態の要旨とみなされることができる。
本発明のこれら及び他の見地は、以下に記述される実施例から明らかになり、それらを参照して説明される。
以下、本発明の例示的な実施形態について、添付の図面を参照して説明する。
以下の記述において、本発明は、平行リビニングの例示的な実施形態に沿って説明される。しかしながら、本発明は、新たな人為的な投影が、他の時点に他のジオメトリにおいて取得された複数の他の投影から再構成される各タイプのリビニングに適用されることができる。言い換えると、本発明は、他の投影ジオメトリへのリビニングが実施されるいかなる適切な走査軌跡にも適用されることができる。
図1は、本発明の例示的な実施形態による、入力/出力ユニット4及びCTスキャナ6に接続される画像処理装置2の簡略化された概略図を示す。画像処理装置2は、計算ユニット8を有する。コンピュータトモグラフィCTスキャナ6は、検査領域12を規定する静止ガントリ10を有する。参照符号14は、静止ガントリ10に取り付けられる回転ガントリを示す。例えばX線管のような透過放射線源16が、それとともに回転する回転ガントリ14上に配置される。放射線源16は、放射線のコーン形状のビーム18を生成するようにコリメートされ、該ビームは、回転ガントリ14の回転中、検査領域12を通過する。
更に、放射線ビーム18をシェーピングし、ビーム18のオン及びオフを選択的にゲート制御するためのコリメータ及びシャッターシステム20が設けられる。更に、放射ビーム18は、放射線源16を電子的に切り替えることによって、オン及びオフをゲート制御されることができる。
参照符号22は、対象の関心領域が少なくとも部分的にCTスキャナ6の検査領域12内にあるように、走査されるべき対象を支持する対象支持体、例えば患者用寝台を示す。回転ガントリ14には、放射線源16の反対側の位置に、放射線検出器24が設けられる。放射線検出器24は、それが検査領域12の対象を通過した後の放射線源16によって発されたコーンビーム18の実質的な量を検出するように構成される。好適には、放射線検出器は、2次元の検出器アレイである。更に、放射線検出器24は、放射線検出器24が検査領域12の周囲を放射線源16と共に回転するように、それとともに回転する回転ガントリ14に取り付けられる。
回転ガントリ14が検査領域12の周囲を回転する間、関心対象をその上に保持する対象支持体22は、検査領域12の中央の水平軸に沿って平行移動する。これにより、放射線源16及び放射線検出器24は、関心対象に対してヘリカルな経路をたどる。任意には、代替実施例において、支持体22が静止したままであり、その一方で、静止ガントリ10が関心対象に対して平行移動され又は他の方法で移動され、それによって、放射線源16及び放射線検出器24が、走査中、関心対象に対するヘリカルなパスに沿って移動されるように、静止ガントリ10が構成される。
走査中、放射線検出器24は、予め決められた時間間隔でサンプリングされる。放射線検出器24から読み取られるサンプリング結果は、電気信号、すなわち電気的なデータであり、以下において投影と呼ばれる。従って、関心対象の走査全体の全データセットは、複数の投影からなり、投影の数は、放射線検出器24がサンプリングされる時間間隔に対応する。複数の投影は、画像処理装置2に伝送され、画像処理装置2において、計算ユニット8によって読み取られる。複数の投影に基づいて、計算ユニット8は、検査領域12の中央の水平軸に対して垂直な平面に沿って、関心対象のスライス画像を再構成する。スライス画像を生成する代わりに、計算ユニット8は、3次元画像又は動画化された表現を決定するように構成されてもよい。後者の場合、例えば、個別の画像は、人間の心臓の周期的な動きを示す画像の動画化されたシーケンスが生成されるように、一緒に組み立てられる。画像は、入力/出力ユニット4を介して出力される。更に、入力/出力ユニット4は、操作者が、CTスキャナ6及び画像処理装置2における走査のパラメータを手動で制御することを可能にする。別の実施例において、絞りのようなアパーチャシステムが、コーンビーム18に変わってファンビームが生成されるように、放射線源16と検査領域12との間に設けられてもよい。
図1の参照符号26は、検査されるべき患者の胸郭に取り付けられるセンサを示す。センサ26は、心電計(ECG)28に接続され、心電計28は、画像処理装置2に接続される。代替例として、ECGの代わりに、ゲーティング関数は、投影から直接的に決定されてもよく(キモグラム)、又は他のセンサ、例えばアコースティックセンサ若しくは光学センサによって決定されてもよい。
以下、本発明について、心臓のシングルフェーズの高解像度再構成が行われる高時間分解能心臓CT再構成の例に関して詳しく記述する。しかしながら、本発明は、高時間分解能心臓CT再構成に限定されず、運動、例えば周期的な運動をする対象が検査されるCTスキャナ、核カメラ及び他の診断スキャナからのデータの再構成に関連して一般に適用されることができることに留意しなければならない。
以下、図1に示されるシステムの動作の例示的な実施形態について、図2を参照して説明する。図2は、図1のシステムの動作のための本発明の方法による例示的な実施形態の簡略化されたフローチャートを示す。
ステップS1における開始後、例えば、操作者が、入力/出力ユニット4への入力を介して図1のシステムの動作を開始したのち、ECGデータが、ステップS2において、ECG28から読み込まれる。ECGデータは、センサ26によって測定される患者の心拍を表す。
方法は、ステップS3へ続く。ステップS3において、画像処理装置2、すなわち画像処理装置2の計算ユニット8に読み込まれる測定されたECGデータに基づいて、計算ユニット8は、ECGのR−R間隔の範囲内にゲーティングウィンドウ幅及び位置を決定する。これは、例えば心臓医のような操作者の対話によって行われることができる。R−R間隔は、人間の心臓の拍動のアノテーション間の間隔である。このウィンドウ幅及び位置は、総取得間隔を通じて一定であり、又は変化しうる。総取得間隔は、CTスキャナにおける走査中に決定される投影のすべて又はいくつかを含む、最初の投影から最後の投影にわたるものとして規定される。
ステップS3においてゲーティングウィンドウが決定されたのち、方法は、ステップS4へ続く。ステップS4において、画像処理装置2は、CTスキャンが患者の完全な投影データセットを取得するように患者に実施されるように、CTスキャナ6を制御する。
ECGデータ及びCT投影を続けて読み込む代わりに、有利には、ECGデータ及びCT投影が、並行して決定される。回顧的(retrospective)なゲーティング及び先見的(prospective)なゲーティングの2つのゲーティング方法が、適用されうる。
回顧的なゲーティングにおいて、投影及びECGデータは連続的に決定され、そののち、すなわちあとから振り返ってみて、ゲーティングウィンドウの位置及び投影のどれが画像の再構成のために使用されるべきかに関する決定が、ECGデータからなされる。
先見的なゲーティングにおいて、ECGデータ及びCT投影は、ここでも並行して決定される。しかしながら、心臓の各サイクルの初め(Rピーク)に、投影が各サイクルのどの部分に決定されるべきかについての決定がなされる。取得間隔は、真空管(放射線源)がサイクルのこの短い間隔においてのみオンにされるとともに、投影がこの間隔においてのみ決定されるように、制御されることができる。
ステップS4における走査後、方法は、ステップS5へ続く。ステップS5において、画像処理装置2の計算ユニット8は、ステップS3において規定されたゲーティングウィンドウ内に測定された各投影を、第1のマーカによってラベリングする。言い換えると、ステップS5において、取得された投影(コーンビーム又は扇形ビーム)の各々は、それが、走査が行われる間に心臓サイクルの1つのゲーティングウィンドウ内に測定された場合、第1のマーカによってラベリングされる。
ラベリングステップS5の後、計算ユニット8は、リビニングステップを実施する。ステップS6において、計算ユニット8は、投影データセット全体について平行リビニングを実施する。これについては、図3a乃至図4を参照してより詳しく説明する。
既に図1を参照して示したように、放射線源16及び放射線検出器24は、共に、患者の周囲を回転する。図3a乃至図3fにおいて、参照符号30は、ヘリカルパスを示し、放射線源16及び放射線検出器24は、このヘリカルパスに沿って回転される。これは、図3a乃至図3fに示されており、これらの図には、連続する投影が、放射線源16及び放射線検出器24のヘリカルパスに沿って示されている。図3a乃至図3fにおいて、放射線検出器24は、複数のライン32乃至36及び複数の列37乃至46を有する2次元検出器アレイとして描かれている。
平行リビニングステップS6の間、平行な投影、すなわち複数の投影の平行な放射線又は列が、最終の画像を再構成するために使用される。これらの平行な投影は、参照符号50によって示されている。図3a乃至図3fに示される検出器構造において、リビニングステップの間に選択される平行投影50の各々は、2次元検出器アレイ24の列37乃至44に対応する。図3aから図3fから分かるように、平行な投影50は、互いに平行な平面を有する。
しかしながら、図3a乃至図3fから分かるように、平行な投影50の各々は、2次元検出器アレイ24の異なる位置に当たる。言い換えると、図3a乃至図3fに示される検出器構造24において、リビニングの間に選択される平行な投影50の各々は、2次元検出器アレイ24の異なる列37乃至44によって検出される。
以下、リビニングステップの他の部分について、図4を参照して説明する。図4に示されるように、複数の平行な投影50は、複数の平行な投影50からなるデータセットに一緒にまとめられる。このデータセットは、以下において、平行リビニングされた投影と呼ばれる。
次のステップS7において、第1のマーカによってラベリングされた投影のうちの1つからの情報を含む平行リビニングされた投影の各々は、第2のマーカによってラベリングされる。次のステップS8において、列に依存する関数が、第2のマーカによってラベリングされたすべての平行リビニングされた投影について、決定される。言い換えると、図3a乃至図3fを参照して、関数が、検出器アレイ24の列37乃至44のうちどれが個々の平行投影50を検出したかに依存して、決定される。こうして、本発明によれば、個々の平行な放射線50が検出された放射線検出器24の位置に関する位置情報に基づいて、列に依存する関数が、第1のマーカによってラベリングされた投影からの情報を含むすべての平行リビニングされた投影について、決定される。
次のステップS9において、再構成プロセスの間、第2のマーカによってラベリングされた平行リビニングされた投影及び列に依存する関数に基づいて、後続の画像表現のためのボクセルが、再構成される。ステップS9においてボクセルが再構成されたのち、方法は、ステップS10へ続き、ステップS10において終了する。言い換えると、本発明により、元々ファンビーム/コーンビーム投影のみの関数であるゲーティングウィンドウが、投影角度及び平行リビニング後の放射線検出器24の列座標の関数になる。有利には、逆投影の間、放射線検出器24の列座標に関するこのような位置情報を取り入れることによって、選ばれたゲーティングウィンドウに本当に属するデータのみが、平行な投影50からの画像を再構成するために使用されるので、再構成方法の時間分解能が、大幅に増加されることができる。
以下、図2のステップS9の例示的な実施形態について、図5を参照して説明する。図5は、図2のステップS9の例示的な実施形態のフローチャートを示す。ステップS20における開始後、CTスキャナ6の放射線ジオメトリに依存する余弦タイプの重み付け関数が、ステップS21において決定される。放射線ジオメトリは、ヘリカルパス、コーンビーム又はファンビームのディメンジョン、放射線源16に対する放射線検出器24の位置、並びに放射線検出器24のライン及び列に関する位置情報からなるグループから選択される少なくとも1つの要素に関する情報を含む。ステップS22の後、余弦タイプの重み付け関数が、ステップS22において第2のマーカによってラベリングされた平行リビニングされた投影に適用される。ステップS22のラベリング後、方法はステップS23へ続く。ステップS23において、重み付けされた、第2のマーカによってラベリングされた、平行リビニングされた投影が、適切なハイパスフィルタによってフィルタリングされる。ステップS23におけるフィルタリングののち、フィルタリングされた、重み付けされた、第2のマーカをもつ、平行リビニングされた投影の各々についてのボクセルの位置が、ステップS24において計算される。方法は、ステップS25へ続く。
ステップS25において、第2のマーカによってマークされた投影の領域上に投影されるボクセルが、選択される。言い換えると、ステップS27において、決定されたボクセルが、本当にゲーティングウィンドウ内にあるかどうかのテストがなされる。このために、ボクセルが、検出器24上に「投影され」、ボクセルが、検出器24上のゲーティングウィンドウ領域に投影されるかどうかの決定がなされる。上述したように、ゲーティングウィンドウ領域は、放射線検出器24の表面上のゲーティングウィンドウの「投影領域」に対応する。
ステップS25の後、方法は、ステップS26へ続き、ステップS26において、正規化ステップが実行される。ステップS26において、特定の角度からのすべての投影の合計の重みが1になることを確実にするために、各投影及びステップS25において選択されたボクセルの各々の重み付け関数が、正規化される。ステップS26における正規化の後、方法は、ステップS27へ続き、そこで終了する。
図6は、0.5秒の回転時間をもつCTシステムの平行リビニングされた検出器について、回転ガントリ14の回転中心からの距離の関数としての時間値を示す。実際の検出器システムの場合、投影のレンジにわたって68ミリ秒のバリエーションが見られる。回転中心を中心とする288mmの円の場合、20ミリ秒のバリエーションがある。
上述したように、平行リビニング自体は、ヘリカルパス30上の異なる放射線源位置、従って異なる時点において取得された複数のファンビーム又はコーンビーム投影から、平行な投影を生成する。平行リビニングされた投影は、扇角の角度レンジ内で測定された一組のファンビーム又はコーンビーム投影からの情報を含む。しかしながら、同じ投影角において適用されたファンビーム投影と同じ時点において測定されるのは、平行な投影の中心線だけであることに留意すべきである。従って、ECGに基づくゲーティングウィンドウに起因する値の範囲内の平行リビニングされた投影を直接的にラベリングすることは、再構成において、ゲーティングウィンドウに本当は属さなかった投影データを使用することになり、それによって、画像に不所望の動きアーチファクト又は欠陥を生じさせる。本発明によれば、このような不所望の動きアーチファクト又は欠陥が、大幅に低減されることができる。本発明によれば、最終的な画像の再構成のために最終的に使用される投影データが、ゲーティングウィンドウ内にあるか否かのテストが、位置情報を使用することによって行われる。画像を最終的に再構成するために使用される投影データがゲーティングウィンドウ内にあるか否かを更に調べるために、個々の投影データが検出器24のどの位置で検出されたかに関する情報が、使用される。
既に上述したように、高時間分解能心臓CT再構成において使用されることに加えて、本発明は、運動する対象が走査されうる核カメラ及び他の診断スキャナに関連するアプリケーションにも適用されることができる。
本発明による、CTスキャナに接続される画像処理装置の例示的な一実施形態を示す図。 図1の装置を動作させるための方法の例示的な実施形態を示すフローチャート。 a乃至fは、図2の方法をより詳しく説明するための簡略図。 図2の方法をより詳しく説明するためのリビニングプロセスの簡略化されたグラフィック表現を示す図。 本発明による図2の方法のステップS9の例示的な実施形態を示すフローチャート。 0.5秒の回転時間をもつCTシステムについて、放射線源及び平行リビニングされた検出器に関する検出器からの回転中心からの距離の関数としての時間値を示す図。

Claims (10)

  1. 対象の複数の投影からボリュメトリックな画像を再構成する方法であって、
    前記複数の投影が、前記対象に対して移動される放射線源及び対応する検出器によって決定され、
    前記方法が、
    複数のリビニングされた投影を決定するために、前記複数の投影についてリビニングステップを実施するステップと、
    前記複数のリビニングされた投影の個々のものがゲーティングウィンドウ内にあるかどうかを、前記複数のリビニングされた投影の個々のものが検出された前記検出器上の位置に関する位置情報に基づいて、決定する決定ステップと、
    を含む方法。
  2. 前記ゲーティングウィンドウを規定するステップと、
    前記ゲーティングウィンドウ内に検出された前記複数の投影の各々を、第1のマーカによってラベリングするステップと、
    前記第1のマーカによってラベリングされた前記複数の投影のうちの1つの情報を含む前記複数のリビニングされた投影の各々を、第2のマーカによってラベリングするステップと、
    を更に含み、前記決定ステップが、
    前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものが検出された前記検出器上の位置に関する前記位置情報を決定するステップと、
    前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものが検出された前記検出器上の前記位置に依存する関数を決定するステップであって、前記関数が、前記複数の投影のラベリングを、前記複数のリビニングされた投影のリビニングされた投影空間に移す、ステップと、
    前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものに対応するボクセルを決定するステップと、
    前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものに対応する前記ボクセルが、前記ゲーティングウィンドウ内にあるかどうかを、前記関数に基づいて検証するステップと、
    を更に含む、請求項1に記載の方法。
  3. 前記放射線源によって発せられる放射線の放射線ジオメトリに依存する重み関数を決定するステップと、
    前記重み関数を、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影に適用するステップと、
    前記重み付けされた、前記第2のマーカによってラベリングされた、複数のリビニングされた投影をフィルタリングするステップと、
    前記フィルタリングされた、重み付けされた、前記第2のマーカによってラベリングされた、複数のリビニングされた投影の個々のものについて、前記リビニングされた投影空間における前記ボクセルの位置を計算するステップと、
    前記ボクセルが、前記リビニングされた投影空間の許容される領域に投影されるかどうかを、前記関数に基づいて決定するステップと、
    特定の角度からのすべての投影の重みが合計で1になることを確実にするために、前記ボクセルに関する前記重み関数を正規化するステップと、
    前記重み関数を乗じられた許容される検出器値を、前記ボクセルに逆投影するステップと、
    を含む、請求項2に記載の方法。
  4. 前記方法が、高時間分解能CT再構成に関し、前記複数の投影が、CTスキャナによって決定され、前記ゲーティングウィンドウが、心電図から決定される、請求項1に記載の方法。
  5. 前記リビニングステップが、平行リビニングを含み、前記放射線源及び前記検出器の移動が、ヘリカルパスに沿う、請求項1に記載の方法。
  6. 対象の複数の投影からボリュメトリックな画像を再構成するための画像処理装置であって、
    前記複数の投影が、前記対象に対して移動される放射線源及び対応する検出器によって決定され、
    前記画像処理装置が、計算ユニットを有し、前記計算ユニットが、複数のリビニングされた投影を決定するために、前記複数の投影についてリビニングを実施するように構成され、前記計算ユニットが更に、前記複数のリビニングされた投影の個々のものがゲーティングウィンドウ内にあるかどうかを、前記複数のリビニングされた投影の個々のものが検出された前記検出器上の位置に関する位置情報に基づいて、決定するように構成される、画像処理装置。
  7. 前記計算ユニットが更に、前記ゲーティングウィンドウを決定するように構成され、
    前記計算ユニットが更に、前記ゲーティングウィンドウ内に検出された前記複数の投影の各々を、第1のマーカによってラベリングするように構成され、
    前記計算ユニットが更に、前記第1のマーカによってラベリングされた前記複数の投影のうちの1つの情報を含む前記複数のリビニングされた投影の各々を、第2のマーカによってラベリングするように構成され、
    前記計算ユニットが更に、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものが検出された前記検出器上の前記位置に関する前記位置情報を決定するように構成され、
    前記計算ユニットが更に、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものが検出された前記検出器上の前記位置に依存する関数であって、前記複数の投影のラベリングを、前記複数のリビニングされた投影のリビニングされた投影空間に移す関数を決定するように構成され、
    前記計算ユニットが更に、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の個々のものに対応するボクセルを決定するように構成され、
    前記計算ユニットが更に、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影の複数の個々のものに対応する前記ボクセルが、前記ゲーティングウィンドウ内にあるかどうかを、前記関数に基づいて検証するように構成される、請求項6に記載の画像処理装置。
  8. 前記計算ユニットが更に、前記放射線源によって発せられる放射線の放射線ジオメトリに依存する重み付け関数を決定するように構成され、
    前記計算ユニットが更に、前記重み関数を、前記第2のマーカによってラベリングされた前記複数のリビニングされた投影に適用するように構成され、
    前記計算ユニットが更に、前記重み付けされた、前記第2のマーカによってラベリングされた、複数のリビニングされた投影をフィルタリングするように構成され、
    前記計算ユニットが更に、前記フィルタリングされた、重み付けされた、前記第2のマーカによってラベリングされた、複数のリビニングされた投影の個々のものについて、前記リビニングされた投影空間における前記ボクセルの位置を計算するように構成され、
    前記計算ユニットが更に、前記ボクセルが、前記リビニングされた投影空間の許容される領域に投影されるかどうかを、前記関数に基づいて決定するように構成され、
    前記計算ユニットが更に、特定の角度からのすべての投影の重みが合計で1になることを確実にするように、前記ボクセルに関する前記重み関数を正規化するように構成される、請求項7に記載の画像処理装置。
  9. 前記画像処理装置が、コンピュータ断層撮影機の一部であり、前記ゲーティングウィンドウが、心電図から決定される、請求項6に記載の画像処理装置。
  10. コンピュータ読み取り可能な媒体を含むコンピュータプログラムであって、
    複数のリビニングされた投影を決定するために、複数の投影についてリビニングステップを実施するステップと、
    前記複数のリビニングされた投影の個々のものが、ゲーティングウィンドウ内にあるかどうかを、前記複数のリビニングされた投影の個々のものが検出された検出器上の位置に関する位置情報に基づいて、決定するステップと、
    を、前記コンピュータプログラムがコンピュータ上で実行されるとき、前記コンピュータに実行させるコンピュータプログラム手段を有する、コンピュータプログラム。
JP2006502445A 2003-02-19 2004-02-10 ボリュメトリック画像再構成のための方法及び装置 Expired - Fee Related JP4431136B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP03100386 2003-02-19
PCT/IB2004/000414 WO2004075115A1 (en) 2003-02-19 2004-02-10 Method and device for volumetric image reconstruction

Publications (2)

Publication Number Publication Date
JP2006517825A true JP2006517825A (ja) 2006-08-03
JP4431136B2 JP4431136B2 (ja) 2010-03-10

Family

ID=32892953

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006502445A Expired - Fee Related JP4431136B2 (ja) 2003-02-19 2004-02-10 ボリュメトリック画像再構成のための方法及び装置

Country Status (6)

Country Link
US (1) US7340028B2 (ja)
EP (1) EP1599838B1 (ja)
JP (1) JP4431136B2 (ja)
AT (1) ATE341046T1 (ja)
DE (1) DE602004002574T2 (ja)
WO (1) WO2004075115A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012527287A (ja) * 2009-05-18 2012-11-08 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 補間不要な、扇形平行ビーム・リビニング

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10304662A1 (de) * 2003-02-05 2004-08-19 Siemens Ag Verfahren zur Erzeugung von Bildern in der Computertomographie mit einem 3D-Bildrekonstruktionsverfahren
JP2008519636A (ja) * 2004-11-15 2008-06-12 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 周期運動物体の検査のためのコンピュータ断層撮影方法
EP2238906B1 (en) * 2005-09-07 2015-04-08 Kabushiki Kaisha Toshiba X-ray computed tomography apparatus
EP1948022A2 (en) * 2005-11-10 2008-07-30 Philips Intellectual Property & Standards GmbH Signal processing unit for producing images
EP2073710B1 (en) * 2006-09-29 2016-05-04 Koninklijke Philips N.V. Fly-by scanning for computed tomography
US20090310740A1 (en) * 2008-06-16 2009-12-17 General Electric Company Computed tomography method and system
CN105093342B (zh) * 2014-05-14 2017-11-17 同方威视技术股份有限公司 螺旋ct系统及重建方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6104775A (en) * 1998-10-29 2000-08-15 Picker International, Inc. 3D image reconstruction for helical partial cone beam scanners using wedge beam transform
US6343108B1 (en) 1999-06-18 2002-01-29 Philips Medical Systems (Cleveland), Inc. Cone beam scanner using oblique surface reconstructions
US6522712B1 (en) * 1999-11-19 2003-02-18 General Electric Company Reconstruction of computed tomographic images using interpolation between projection views

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012527287A (ja) * 2009-05-18 2012-11-08 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 補間不要な、扇形平行ビーム・リビニング

Also Published As

Publication number Publication date
DE602004002574D1 (de) 2006-11-09
US20060062346A1 (en) 2006-03-23
JP4431136B2 (ja) 2010-03-10
WO2004075115A1 (en) 2004-09-02
ATE341046T1 (de) 2006-10-15
EP1599838A1 (en) 2005-11-30
DE602004002574T2 (de) 2007-06-14
US7340028B2 (en) 2008-03-04
EP1599838B1 (en) 2006-09-27

Similar Documents

Publication Publication Date Title
US8184883B2 (en) Motion compensated CT reconstruction of high contrast objects
JP4487095B2 (ja) 撮像システムの回顧的心臓ゲーティングを用いた冠状動脈の石灰化の検出
EP1723607B1 (en) Motion compensation
US7221728B2 (en) Method and apparatus for correcting motion in image reconstruction
JP4347807B2 (ja) 投影データから運動情報を導く方法及び装置
US7054475B2 (en) Apparatus and method for volumetric reconstruction of a cyclically moving object
US20080267455A1 (en) Method for Movement Compensation of Image Data
JP4602146B2 (ja) 運動する対象の構造データの取得装置
JP2004160222A (ja) コンピュータ断層撮影方法および装置
US7702063B2 (en) CT method for the examination of cyclically moving object
JP2009534079A (ja) 複数の部分的な走査軌道を有するコーンビーム・コンピュータ断層撮影
EP1663004A2 (en) Computer tomography method using a cone-shaped bundle of rays
EP1940294A1 (en) Short scan cardiac ct on a quasi axial trajectory
JP5539719B2 (ja) 画像形成システム
US20070183640A1 (en) Gated cone-beam computed tomography reconstruction with motion compensation
JP4431136B2 (ja) ボリュメトリック画像再構成のための方法及び装置
CN1839412A (zh) 产生时间内插断层摄影图像的方法和装置
JP3878788B2 (ja) トモシンセシスを行う方法及び装置
JP4691505B2 (ja) 周期的に動作する対象のためのコンピュータ断層撮影方法
EP2076885B1 (en) Imaging system, imaging method and computer program for imaging a region of interest
EP1696796A1 (en) Computer tomography method for objects moving periodically
EP4002279A1 (en) Method for use in x-ray ct image reconstruction

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070209

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

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20091218

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121225

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees