JP5199081B2 - 心臓ct撮影のバンドアーチファクトの抑制 - Google Patents

心臓ct撮影のバンドアーチファクトの抑制 Download PDF

Info

Publication number
JP5199081B2
JP5199081B2 JP2008517643A JP2008517643A JP5199081B2 JP 5199081 B2 JP5199081 B2 JP 5199081B2 JP 2008517643 A JP2008517643 A JP 2008517643A JP 2008517643 A JP2008517643 A JP 2008517643A JP 5199081 B2 JP5199081 B2 JP 5199081B2
Authority
JP
Japan
Prior art keywords
region
window
time window
radiation
processor
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.)
Expired - Fee Related
Application number
JP2008517643A
Other languages
English (en)
Other versions
JP2008546464A (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
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 NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2008546464A publication Critical patent/JP2008546464A/ja
Application granted granted Critical
Publication of JP5199081B2 publication Critical patent/JP5199081B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4085Cone-beams
    • 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/027Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral

Description

本発明は診断用の撮影技術に関する。本発明は、特に被検体の心臓のコンピュータ断層撮影に適用されるものであり、特にそれを参照して説明される。しかしながら、本発明はまた、その他の種類のコンピュータ断層撮影、単一光子放出型コンピュータ断層撮影(SPECT)、陽電子放出型断層撮影(PET)、及び3次元x線撮影などにも適用可能なものである。
一般に、コンピュータ断層撮影システムは、検査されるべき対象物の周りを回転するx線源及びx線検出器を有している。対象物は複数の方向から、x線源からのx線ビームで照射される。x線検出器はそれぞれの方向で対象物を通過したx線を受け取り、問題としている方向に関する減衰プロファイルを形成する。減衰プロファイルは、問題としている方向で対象物を貫通するx線形路に沿ってx線が吸収又は散乱されることに起因する、入射x線の減衰を表す。
ヘリカル心臓円錐ビーム画像は、位相選択アルゴリズムを用いて再構成される。一般的に、心臓画像の生成には特定の心位相が選択される。選択された位相、すなわち、相異なる心サイクルではあるが同一の心位相に対応する時点、に時間的に近接して収集されたデータのみが、マルチスライス再構成処理で同時に使用される。スキャンパラメータ、患者の心拍とその変動性、及び心臓ゲーティングのウィンドウ幅と位置に応じて、各々のボクセルの再構成に使用されるサイクル数は異なる。ボクセルは、一般的に、所与のボクセルすなわち照射ウィンドウを通過した全ての心サイクルにわたる全ての利用可能な光線から再構成される。
CT画像の品質は、しばしば、互いに空間的に近接するが異なる軸方向位置すなわちZ座標を有するボクセルが、異なる心サイクルの位相点付近からの投影を用いて再構成されることにより生じるバンドアーチファクトによって低下されてしまう。相異なるZ座標にあるボクセル群は、相異なるグループの心サイクル群又は照射ウィンドウ群で照射され得る。第1のボクセルは或る1つの照射ウィンドウで収集されたデータを用いて再構成され、その一方で、第2のボクセルは異なる1つの照射ウィンドウで収集されたデータを用いて再構成される。心サイクル数は照射ウィンドウが異なると変わってくる。このことは、これら2つのボクセルに関して計算されるCT数の間に相違を生じさせ得るものである。
例えばサジタル(矢状)断像又はコロナル(冠状)断像などのZ方向に細長い断像は、例えば非周期的な心臓運動、各心サイクル内での位相点の決定の一貫性の欠如、様々なボクセルに寄与する心サイクル数の相違、及びその他などの幾つかの原因に帰せられ得るストリーク(streak)やアーチファクトを含んでいることがある。
バンドアーチファクトの発生源を是正することに依存しない、心臓の円錐ビーム撮影におけるバンドアーチファクトを抑制する技術が望まれる。
本発明は、心臓の円錐ビーム撮影におけるバンドアーチファクトを抑制する方法及び装置を提供することを目的とする。
本発明の一態様に従って、被検体が検査領域内で動くときに1つ又は複数の連続する心サイクル内に現れる予め選択された位相点で撮影を行う診断用撮影システムが開示される。照射ウィンドウが、予め選択された位相点を取り囲む位相点ウィンドウ内で被検体の再構成されるボクセルの各々を照射する。時間ウィンドウ決定プロセッサ又はアルゴリズムが、Nを整数として、位相点ウィンドウが照射ウィンドウ内に完全に含まれるN個の連続する心サイクルを含む時間ウィンドウを決定し、且つ決定された時間ウィンドウに従って前記照射ウィンドウを短く切り詰める。Z領域分割プロセッサ又はアルゴリズムが、軸方向に沿って第1及び第2の領域を決定する。第1及び第2の領域は互いに交互にされる。重み付けプロセッサが、時間ウィンドウの第1及び第2の領域にある再構成されるボクセルの測定値に、正規化された加重を適用する。逆投影プロセッサ又はアルゴリズムが、重み付けされた測定値を立体的な画像表現へと3次元的に逆投影する。
本発明の他の一態様に従って撮影方法が開示される。予め選択された位相点を取り囲む位相点ウィンドウ内で被検体の再構成されるボクセルの各々が照射される。この位相点の各々は被検体が検査領域内で動くときに1つ又は複数の連続する心サイクル内に現れるものである。Nを整数として、位相点ウィンドウが照射ウィンドウ内に完全に含まれるN個の連続する心サイクルを含む時間ウィンドウが決定される。決定された時間ウィンドウに従って照射ウィンドウが短く切り詰められる。軸方向に沿って時間ウィンドウ内で第1及び第2の領域が決定される。この第1及び第2の領域は互いに交互にされる。時間ウィンドウの第1及び第2の領域にある再構成されるボクセルの測定値に重み付けが適用される。重み付けを適用された測定値は立体的な画像表現へと3次元的に逆投影される。
本発明の1つの利点は、バンドアーチファクトが抑制されることにある。他の1つの利点は、バンドアーチファクトの発生源を決定することを回避する技術にある。他の1つの利点は、心臓及びその他のゲーティング撮影技術による画像が改善されることにある。
以下の好適な実施形態についての詳細な説明を読むことにより、数多くの更なる効果が当業者に明らかになるであろう。
本発明は様々な構成要素とその配置、及び様々な処理段階及びその編成の形態を取り得る。図面は、好適な実施形態を例示するためだけのものであり、本発明を限定するものとして解釈されるべきではない。
図1を参照するに、撮像システム10は、放射線源14を有するコンピュータ断層撮影スキャナ12を含んでおり、放射線源14は検査領域16に向けて、好ましくは円錐ビーム又はウェッジビームである放射線ビームを生成する。放射線ビームは、検査領域16に配置された撮影対象の関心領域を横切るとき相互作用して部分的に吸収され、放射線ビームが検査領域を通過すると、空間的に変化する放射線吸収が生成される。放射線ビームが検査領域16を通過した後、好ましくは2次元検出器である放射線検出器18が、吸収により減衰された放射線を検出する。放射線源14と検出器18の各放射線検出素子との間の経路は光線と呼ばれる。
好ましくは、放射線源14は円錐ビーム状のx線を生成する。放射線源14及び検出器18は、検出器18が放射線源14からのx線を継続して受光するように、回転式ガントリー20上の反対側で対向し合うようにして取り付けられる。放射線源14及び検出器18が回転式ガントリー20上で検査領域16の周りを回転するとき、複数回の回転にまたがって視界(view)群が収集される。各視界すなわちデータの2次元配列は、検出器18の検出素子の同時サンプリングによって収集された、放射線源14の位置に頂点を有する円錐状の光線群を表す。ヘリカル円錐ビームによるコンピュータ断層撮影においては、被検体支持台又は寝台26はモータ駆動28によって軸方向すなわちZ方向に直線的に移動される。
選択的に、円錐ビームによるコンピュータ断層撮影の投影データは、(i)被検体支持台26を、各軸方向スキャンの最中は静止させ且つ軸方向スキャン同士の間に線形に一段進めるか、(ii)螺旋軌道を規定するように被検体支持台26を連続的に移動させるか、の何れかで複数の回転にまたがって収集される。放射線検出器18の検出素子の出力は、電気的な収集・積分された減衰投影値μdに変換されて、データメモリ30に記録される。各投影データμdは、放射線源14から検出器18の検出素子群のうちの対応する1つへの直線に沿っての減衰の線積分に相当する。
典型的な円錐ビームの幾何学構成では、一般的に、この線積分の添字はその測定値を測定するのに使用された検出素子に対応する。しかしながら、この線積分の添字が検出素子の番号との直接的な対応関係を欠いている場合も意図される。このような直接的な対応関係の欠如は、例えば、リビニングされた(re-binned)投影間の補間によって生じ得る。
マルチスライス式スキャナにおける放射線源を焦点とする収集の幾何学構成では、データメモリ30に記録された減衰の線積分の測定値群、すなわち、投影データセットの投影値群は、P(α,β,n)とパラメータ化されることができる。ただし、αは回転式ガントリー20の位置によって決定される放射線源14の光源角、βは扇内の角度(Φを扇角としてβ∈[−Φ/2,Φ/2])、nは検出器の行番号である。
心臓モニタ32は患者の心サイクルを監視し、典型的に各サイクルのR波、すなわち、各R−R間隔におけるR波に対する位相点34を検出する。位相点34の位置は、心臓運動の特徴及び要求される診断情報に従って臨床医によって選択される。減衰データは、ソート手段38によって、選択された心位相の各々の間に収集されたデータセット、すなわち、心位相特定データセットに仕分けされる。円錐から平行ビーム配置への心位相特定データは、リビニングプロセッサ40によって一組の平行な視界へとリビニングされる。この平行な視界は、回転軸に垂直な平面である軸平面に投影される。各視界は等間隔のπライン群を含んでいる。πラインとは、スキャン視野(FOV)と交差する軸平面内に含まれ、且つ正準座標θπ,lによって特徴付けられる線積分として定義されるものである。ただし、θπは伝播角∈[0,π)、lは治療中心からの距離である。特に、短い時間ウィンドウによって定められる心位相では、1つの心位相でのデータは、複数の回転及び心サイクルの各々における短い円弧部分で収集されたデータに相当する。この円弧部分のデータは単独では少なすぎ、完全なデータセットとはなり得ない。完全なデータセットを生成するため、データは幾つかの心サイクルにわたって収集され、必要であれば補間される。心位相特定データセットは対応する位相メモリ42に記録される。再構成プロセッサ44は、各々の選択された心位相に関するデータから、画像メモリ46に記録される対応する3次元画像へと視界データを処理する。或る1つの再構成技術においては、並列な投影の各々内の測定値がフィルタ48を用いてフィルタリングされる。コンボルバ(convolver)50は、例えば1次元ヤコビアンカーネル(Jacobian kernel)等のランプカーネルを用いた1次元コンボリューションを実行する。このコンボリューションは平行な測定値セットに沿って実行される。データは角度によって畳み込まれ、角度範囲θ∈[0,π)をカバーする2次元データセットが完成される。完全なデータセットのために補間が必要とされることがあり得る。逆投影プロセッサ52は、詳細に後述されるように、畳み込まれたデータの3次元画像表現への正規化重み付け逆投影を実行する。ビデオプロセッサ54は、画像メモリ46の内容物の一部又は全てを処理し、例えば3次元レンダリング、選択画像スライス、最大強度投影、CINEアニメーション、又はこれらに類するもの等の、人間が視認可能な画像表現を作り出す。この人間が視認可能な画像表現はユーザインターフェース58のディスプレー56上に表示される。ユーザインターフェース58は好ましくは、パーソナルコンピュータ、ワークステーション、ノート型コンピュータ、又はこれらに類するものである。必要に応じて、画像メモリ46に含まれる選択された内容物は、紙面にプリントされ、不揮発性の電子又は磁気記憶媒体に記憶され、ローカルエリアネットワーク又はインターネット上で伝送され、あるいはその他の方法で処理される。好ましくは、放射線医又はその他の操作者が入力手段60を介してコンピュータ断層撮影スキャナ12を制御し、それにより、撮影セッションを設定し、撮影セッションを変更し、撮影セッションを実行し、撮影セッションを監視し、あるいはその他の方法でスキャナ12を制御するように、スキャン制御器62をプログラムする。
図1に加えて図2をも参照するに、時間的重み付けプロファイルプロセッサ又はアルゴリズム64は、ボクセルに依存せず且つ投影に依存する時間的重み付けプロファイル66を計算する。この時間的重み付けプロファイル66は位相点ウィンドウの複製から成る。すなわち、これは各心サイクル内で予め特定された位相点34に中心を有する関数:
(c)=Λ(c/c’)
から成る。ただし、Λは1をピーク値とする三角形の形状をした関数であり、cは測定点と該測定点に最も近い位相点との間の時間的なずれであり、c’は再構成に寄与する位相点ウィンドウの半値幅である。
最高の分解能は、完全なデータセットを実現するc’の最小幅を求めることによって得られる。位相点ウィンドウは全ボクセルの再構成に十分なデータを提供する幅を用いて計算される。
図1に加えて図3をも参照するに、Z領域分割(若しくはZ領域重み付け関数決定)手段、プロセッサ又はアルゴリズム68は、Zプロファイル70を軸方向Zに、交互になった第1領域及び第2の領域、すなわち固定加重領域及び可変加重領域72、74に分割する。各々の2つの隣接する第1及び第2の領域72、74を合わせた軸方向の長さは:
=ν・RR
と定義される。ただし、νは寝台の速度、RRは患者のECGから抽出された心周期である。
第2の領域74の各々は対応する位相点34を囲むように定められる。一実施形態において、縦線で図示される各位相点34は関連する第2の領域74の実質的に中心に位置している。第2の領域74の軸方向の長さLは:
=γ・ν・RR
と定義される。ただし、νは寝台の速度、RRは患者のECGから抽出された心周期、γは心周期内の第2の領域の長さを規定するパラメータである。
上述のように選択されるパラメータαはパラメータγ:
γ=α+mod((WWmin−RT/2)/RRmax−α,1)
の最小値を設定する。ただし、WWminは再構成される全ボクセルを照射するウィンドウの組からの最小の照射ウィンドウ幅、RTはガントリーの回転時間、RRmaxはスキャン中の最大心周期である。
図1及び3に加えて再び図2を参照するに、時間ウィンドウ決定アルゴリズム又はプロセッサ88は、対応する再構成されるボクセルA、B、C、・・・の各々に関する時間ウィンドウ90を決定する。より具体的には、時間ウィンドウ90は、Nを整数として、所与のボクセルの照射ウィンドウ94内に完全に含まれるN個の連続する位相点34の位相点ウィンドウ92から収集された投影のみを再構成において含むように選定される。所与のボクセルの照射ウィンドウ94に完全には含まれない位相点ウィンドウは、直線96にて図示されるように切り捨てられる。各位相点34の周りの狭い位相点ウィンドウは一層良好な分解能をもたらすが、少ない数のデータ点は再構成画質を低下させることになる。
連続する位相点34の数Nは:
N=floor((WWmin−RT/2)/RRmax−α)
によって決定される。ただし、WWminは再構成される全ボクセルを照射するウィンドウの組からの最小の照射ウィンドウ幅、RTはガントリーの回転時間、RRmaxはスキャン中の最大心周期、αは前もって選定され、0.3−0.5の典型値を想定するパラメータである。
図3の例においては、連続する位相点の数Nは2であり、例えば、照射ウィンドウ94は2つの位相点ウィンドウへと短く切り詰められている。斯くして、第1の領域72内の各々の関連するボクセルが放射線を受ける動きの範囲を表す照射ウィンドウ94の各々は、正確な数の心サイクルに対応する一定の時間ウィンドウ90に置き換えられる。時間照射ウィンドウ90の何れかの側への時間的重み付けプロファイルは切り捨てられる。
重み付けプロセッサ100は、Z座標に関して全ボクセルに、滑らかに変化する加重を適用する。より具体的には、固定重み付け(又は時間ウィンドウ)プロセッサ又はアルゴリズム102が、例えばボクセルA、C等の第1の領域72内に位置するボクセルのために、正規化された逆投影加重を計算する。固定重み付けプロセッサ102は、時間的重み付けプロファイル66のうちのボクセルに依存する時間ウィンドウ90に入る区画を選定する。選定された区画は、同一のθ∈[0,π)に折り畳まれる全投影に与えられる加重の合計が1に等しくなるように正規化される。各測定値に与えられる正規化された加重の総和は:
W=W(c)
である。
可変加重プロセッサ110は、例えばボクセルB等の第2の領域74内に位置する全ボクセルに、滑らかに変化する重み付けを適用する。第2の領域74内のボクセルBに適用される正規化された逆投影加重が、隣接する第1の領域の正規化された重み付けプロファイルの間への線形補間によって計算される。この補間の加重は、第2の領域74とそれに隣接する先導及び後続の第1の領域116、118との第1及び第2の境界線112、114と、ボクセルBのZ座標との間の距離d1、d2に従って計算される。
逆投影器120(又は逆投影プロセッサ又はアルゴリズム)は正規化された投影を画像メモリに逆投影する。
斯くして、ボクセルに寄与する位相点の数を、第1の領域内のボクセルの位相点の正確な一定数に丸めること、及び第2の領域内のボクセルの位置に関する連続した滑らかに変化する正規化された逆投影重み付けを提供することにより、Z方向での或る1つのボクセルから別の1つのボクセルへと急激にCT数が変化する現象が平滑化され、あるいは実質的に排除され、それにより、アーチファクトが抑制される。
好適な実施形態を参照しながら本発明を説明してきた。以上の詳細な説明を読み、理解した者によって改良及び改変が想到されることは明らかである。本発明は、添付の請求項又はその均等物の範囲内に入る限りにおいて、そのような全ての改良及び改変を含むとして解釈されるものである。
コンピュータ断層撮影システムを示す図である。 ボクセルに依存しない重み付けプロファイルを示す図である。 交互にされた重み付けプロファイル領域に分割されるZ区間の一部を示す図である。

Claims (20)

  1. 被検体が検査領域内で動くときに1つ又は複数の連続する心サイクル内に現れる予め選択された位相点で撮影を行う診断用撮影システムであって:
    前記予め選択された位相点を取り囲む位相点ウィンドウ内で前記被検体の再構成されるボクセルの各々を照射する照射ウィンドウ;
    Nを整数として、前記位相点ウィンドウが前記照射ウィンドウ内に完全に含まれるN個の連続する心サイクルを含む時間ウィンドウを決定し、且つ決定された時間ウィンドウに従って前記照射ウィンドウを短く切り詰める、時間ウィンドウ決定プロセッサ又はアルゴリズム;
    軸方向に沿って互いに交互にされた第1及び第2の領域を決定するZ領域分割プロセッサ又はアルゴリズム;
    前記時間ウィンドウの前記第1及び第2の領域にある前記再構成されるボクセルの測定値に、正規化された加重を適用する重み付けプロセッサ;及び
    正規化された加重を適用された測定値を立体的な画像表現へと3次元的に逆投影する逆投影プロセッサ又はアルゴリズム;
    を有するシステム。
  2. 投影に依存する時間的重み付けプロファイルを決定する時間的重み付けプロファイルプロセッサ;
    を更に有し、且つ
    前記重み付けプロセッサは:
    前記時間的重み付けプロファイルのうちの前記時間ウィンドウに対応する区画を選定し、且つ関連する第1の領域内にある各ボクセルに、選定された区画の正規化された一定の加重を適用する、時間ウィンドウプロセッサ
    を含む、
    請求項1に記載のシステム。
  3. 前記重み付けプロセッサは更に:
    関連する第2の領域内にある各ボクセルの測定値に、連続的且つ滑らかな重み付けを施す可変加重プロセッサ
    を含む、請求項2に記載のシステム。
  4. 前記関連する第2の領域内にある前記ボクセルに与えられる加重は、該第2の領域に隣接する先導及び後続の第1の領域内にあるボクセルに与えられる前記一定の加重の間で補間される、請求項3に記載のシステム。
  5. 回転式ガントリー;
    前記回転式ガントリーに配置された、前記検査領域を横切る円錐ビーム放射線の放射線源;及び
    前記放射線源に対向するように前記回転式ガントリーに配置され、前記放射線が前記検査領域を通過した後で該放射線を検出し、且つ検出された放射線を電子的な投影データ様式に変換する、放射線検出器;
    を更に有する請求項1に記載のシステム。
  6. 前記時間ウィンドウの大きさを決定する連続する心サイクルの数Nは:
    WWminは前記再構成される全ボクセルを照射する照射ウィンドウの組からの最小の照射ウィンドウ幅、RTはガントリーの回転時間、RRmaxはスキャン中の最大心周期、αは前もって選定されたパラメータであるとして、
    N=floor((WWmin−RT/2)/RRmax−α)
    である、請求項5に記載のシステム。
  7. 各々の隣接し合う前記第1及び第2の領域を合わせた長さLは:
    νは前記軸方向での前記被検体の速度、RRは前記心サイクルの周期であるとして、
    =ν・RR
    である、請求項1に記載のシステム。
  8. 前記第2の領域の長さLは:
    νは前記軸方向での前記被検体の速度、RRは前記心サイクルの周期、γは前記心サイクル内の該第2の領域の長さLを規定するパラメータであるとして、
    =γ・ν・RR
    である、請求項1に記載のシステム。
  9. 少なくとも複数の照射ウィンドウ内で投影データを収集するCTスキャナであり:
    回転式ガントリー;
    前記検査領域を横切る円錐ビーム放射線の放射線源;
    前記放射線が前記検査領域を通過した後で該放射線を検出し、それを投影データ様式に変換する放射線検出器であり、前記放射線源及び該放射線検出器は、前記検査領域の周りを連続的に回転するように前記回転式ガントリー上に対向し合うように配置されている、放射線検出器;及び
    前記立体的な画像表現を表示するディスプレー;
    を有するCTスキャナ;
    を更に有する請求項1に記載のシステム。
  10. 予め選択された位相点を取り囲む位相点ウィンドウ内で被検体の再構成されるボクセルの各々を照射する照射段階であり、該位相点の各々は前記被検体が検査領域内で動くときに1つ又は複数の連続する心サイクル内に現れる、照射段階;
    Nを整数として、前記位相点ウィンドウが照射ウィンドウ内に完全に含まれるN個の連続する心サイクルを含む時間ウィンドウを決定し、且つ決定された時間ウィンドウに従って前記照射ウィンドウを短く切り詰める段階;
    軸方向に沿って前記時間ウィンドウ内で第1及び第2の領域を決定する決定段階であり、該第1及び第2の領域は互いに交互にされる、決定段階;
    前記時間ウィンドウの前記第1及び第2の領域にある前記再構成されるボクセルの測定値に、正規化された加重を適用する段階;及び
    正規化された加重を適用された測定値を立体的な画像表現へと3次元的に逆投影する段階;
    を有する撮影方法。
  11. 投影に依存する時間的重み付けプロファイルを決定する段階;
    前記時間的重み付けプロファイルのうちの前記時間ウィンドウに対応する区画を選定する段階;及び
    関連する第1の領域内にある各ボクセルに、選定された区画の正規化された共通の加重を適用する段階;
    を更に有する請求項10に記載の方法。
  12. 前記加重を適用する段階は更に:
    関連する第2の領域内にある各ボクセルの測定値に、連続的且つ滑らかな加重を適用する段階、
    を含む、請求項11に記載の方法。
  13. 前記関連する第2の領域内にある前記ボクセルに与えられる加重は、該第2の領域に隣接する第1の領域群にあるボクセルに与えられる加重の間で補間される、請求項12に記載の方法。
  14. 前記検査領域の周りで円錐ビーム放射線を回転させる段階;及び
    前記被検体を横切った放射線を検出する段階;
    を更に有する請求項10に記載の方法。
  15. 前記時間ウィンドウに対応する連続する心サイクルの数Nは:
    WWminは前記再構成される全ボクセルを照射する照射ウィンドウの組からの最小の照射ウィンドウ幅、RTはガントリーの回転時間、RRmaxはスキャン中の最大心周期、αは前もって選定されたパラメータであるとして、
    N=floor((WWmin−RT/2)/RRmax−α)
    である、請求項14に記載の方法。
  16. 隣接し合う前記第1及び第2の領域を合わせた長さLは:
    νは前記被検体の速度、RRは前記心サイクルの周期であるとして、
    =ν・RR
    である、請求項10に記載の方法。
  17. 前記第2の領域の長さLは:
    νは前記被検体の速度、RRは前記心サイクルの周期、γは前記心サイクル内の該第2の領域の長さLを規定するパラメータであるとして、
    =γ・ν・RR
    である、請求項10に記載の方法。
  18. 請求項10に記載の段階群を実行するCTスキャナ。
  19. 検査領域を横切る円錐ビーム放射線の放射線源;
    前記放射線が前記検査領域を通過した後で該放射線を検出し、それを投影データ様式に変換する放射線検出器;
    円錐ビーム投影データを3次元再構成画像に再構成する画像プロセッサであり:
    前記投影データを、選択された心サイクルの各々の間に収集されたデータセットに仕分けする段階;
    収集されたデータを平行な光線の様式にリビニングする段階;
    前記平行な光線の様式のデータをフィルタリングする段階;
    フィルタリングされたデータを畳み込む段階;
    連続した心サイクルを整数個含む時間ウィンドウを決定する段階;
    投影に依存する時間的重み付けプロファイルを決定する段階;
    軸方向に沿って前記時間ウィンドウ内で第1及び第2の領域を決定する段階であり、該第1及び第2の領域は交互にされる、決定する段階;
    前記時間的重み付けプロファイルのうちの前記時間ウィンドウに対応する区画を選定する段階;
    関連する第1の領域内にある各ボクセルの測定値に、前記時間的重み付けプロファイルのうちの選定された区画の正規化された共通の加重を適用する段階;及び
    関連する第2の領域内にある各ボクセルの測定値に、連続的且つ滑らかな重み付けを施す段階;
    を実行するようにプログラムされた画像プロセッサ;並びに
    畳み込まれ、正規化され、重み付けられた逆投影データを人が視認可能な画像様式で表示するディスプレー;
    を有する診断用撮影システム。
  20. 円錐ビーム放射線の放射線源を用いて検査領域を横切らせる段階;
    前記放射線が前記検査領域を通過した後で前記放射線を検出し、該放射線を電子的な投影データ様式に変換する段階;
    前記円錐ビームの投影データを3次元再構成画像へと再構成する段階であり:
    前記投影データを、選択された複数の心サイクルの各々の間に収集されたデータセットに仕分けし;
    前記収集されたデータを平行な光線の様式にリビニングし;
    前記平行な光線の様式のデータをフィルタリングし;
    フィルタリングされたデータを畳み込み;
    整数個の連続した心サイクルを含む時間ウィンドウを決定し;且つ
    投影に依存する時間的重み付けプロファイルを決定する;
    ことを含む再構成する段階;
    軸方向に沿って前記時間ウィンドウ内で第1及び第2の領域を決定する段階であり、該第1及び第2の領域は互いに交互にされる、決定する段階;
    前記時間的重み付けプロファイルのうちの前記時間ウィンドウに対応する区画を選定する段階;
    関連する第1の領域内にある各ボクセルの測定値に、前記時間的重み付けプロファイルのうちの選定された区画の正規化された共通の加重を適用する段階;
    関連する第2の領域内にある各ボクセルの測定値に連続的な重み付けを施す段階;及び
    畳み込まれ、正規化され、加み付けられた逆投影データを、人が視認可能な画像様式で表示する段階;
    を有する心臓撮影方法。
JP2008517643A 2005-06-22 2006-06-13 心臓ct撮影のバンドアーチファクトの抑制 Expired - Fee Related JP5199081B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US69275305P 2005-06-22 2005-06-22
US60/692,753 2005-06-22
PCT/IB2006/051889 WO2006136973A2 (en) 2005-06-22 2006-06-13 Bands artifact reduction for cardiac ct imaging

Publications (2)

Publication Number Publication Date
JP2008546464A JP2008546464A (ja) 2008-12-25
JP5199081B2 true JP5199081B2 (ja) 2013-05-15

Family

ID=37570815

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008517643A Expired - Fee Related JP5199081B2 (ja) 2005-06-22 2006-06-13 心臓ct撮影のバンドアーチファクトの抑制

Country Status (6)

Country Link
US (1) US8121374B2 (ja)
EP (1) EP1895906B1 (ja)
JP (1) JP5199081B2 (ja)
CN (1) CN101203181B (ja)
AT (1) ATE542480T1 (ja)
WO (1) WO2006136973A2 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE602007011309D1 (de) * 2006-10-18 2011-01-27 Philips Intellectual Property Bildgebungssystem, bildgebungsverfahren und computerprogramm zur abbildung eines interessensbereiches
JP5887057B2 (ja) * 2008-02-07 2016-03-16 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 患者の呼吸安定を測定及び予測する装置
CN102306239B (zh) * 2011-07-22 2014-11-05 上海赛提菲克医疗器械有限公司 基于锥形束ct图像ct值校正技术的放疗剂量评估和优化方法
DE102013201822B4 (de) * 2013-02-05 2022-12-29 Siemens Healthcare Gmbh Verfahren zur Erzeugung eines PET- oder SPECT-Bilddatensatzes und Hybridbildgebungsmodalität hierfür
CN104331914B (zh) * 2014-09-18 2017-06-30 东莞南方医大松山湖科技园有限公司 重建心脏相位图像的方法和装置
US10245000B2 (en) 2014-12-12 2019-04-02 General Electric Company Method and system for defining a volume of interest in a physiological image
US11308662B2 (en) 2016-08-02 2022-04-19 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
CN106296765B (zh) * 2016-08-02 2020-06-02 上海联影医疗科技有限公司 图像重建方法及系统
US10347014B2 (en) 2016-08-02 2019-07-09 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
CN106548473B (zh) * 2016-11-07 2019-03-08 赛诺联合医疗科技(北京)有限公司 一种构建相位图像的方法及装置
CN107157505B (zh) * 2017-06-09 2019-10-18 北京理工大学 一种提高锥形束spect成像分辨率的方法
US20190180481A1 (en) * 2017-12-13 2019-06-13 General Electric Company Tomographic reconstruction with weights
KR102137765B1 (ko) * 2018-08-20 2020-07-24 연세대학교 산학협력단 사전 영상을 이용하는 고품질의 4차원 콘빔 전산화 단층 촬영 시스템
EP3756458A1 (en) * 2019-06-26 2020-12-30 Viking Genetics FmbA Weight determination of an animal based on 3d imaging
US11160523B2 (en) * 2019-09-30 2021-11-02 GE Precision Healthcare LLC Systems and methods for cardiac imaging

Family Cites Families (13)

* 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
US6597803B1 (en) * 1999-10-29 2003-07-22 Ge Medical Systems Global Technology Company, Llc Hybrid reconstruction for high pitch multi-slice helical cardiac imaging
US6522712B1 (en) * 1999-11-19 2003-02-18 General Electric Company Reconstruction of computed tomographic images using interpolation between projection views
DE10129631A1 (de) * 2001-06-20 2003-01-02 Philips Corp Intellectual Pty Verfahren zur Rekonstruktion eines hoch aufgelösten 3D-Bildes
DE10162768A1 (de) * 2001-12-20 2003-07-03 Philips Intellectual Property Computertomograph
DE10207623B4 (de) * 2002-02-22 2004-05-06 Siemens Ag Verfahren für die Computertomographie sowie Computertomographie (CT)-Gerät
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
DE10244180B4 (de) * 2002-09-23 2009-08-27 Siemens Ag Verfahren zur Bilderstellung in der Computertomographie eines periodisch bewegten Untersuchungsobjektes und CT-Gerät zur Durchführung des Verfahrens
US7289841B2 (en) * 2002-10-25 2007-10-30 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
DE10251448A1 (de) * 2002-11-05 2004-05-19 Siemens Ag Verfahren für die Computertomographie eines periodisch sich bewegenden Untersuchungsobjektes, sowie ein CT-Gerät zur Durchführung dieses Verfahrens
EP1615561B1 (en) 2003-04-10 2009-06-24 Philips Intellectual Property & Standards GmbH Computer tomography method for a periodically moving object
JP2004337515A (ja) * 2003-05-19 2004-12-02 Hitachi Medical Corp X線ct装置
CA2555201A1 (en) * 2004-02-09 2005-08-18 Institut De Cardiologie De Montreal Computation of a geometric parameter of a cardiac chamber from a cardiac tomography data set

Also Published As

Publication number Publication date
ATE542480T1 (de) 2012-02-15
CN101203181A (zh) 2008-06-18
JP2008546464A (ja) 2008-12-25
WO2006136973A3 (en) 2007-07-26
WO2006136973A2 (en) 2006-12-28
US8121374B2 (en) 2012-02-21
EP1895906B1 (en) 2012-01-25
EP1895906A2 (en) 2008-03-12
CN101203181B (zh) 2010-06-16
US20080193003A1 (en) 2008-08-14

Similar Documents

Publication Publication Date Title
JP5199081B2 (ja) 心臓ct撮影のバンドアーチファクトの抑制
US7558362B2 (en) Streak artifact reduction in cardiac cone beam CT reconstruction
US7340027B2 (en) Metal artifact correction in computed tomography
JP5221394B2 (ja) ラドンデータから画像関数を再構成する方法
JP5248648B2 (ja) コンピュータ断層撮影システムおよび方法
EP1800264B1 (en) Image reconstruction with voxel dependent interpolation
JP4606414B2 (ja) 円錐形状光線束を用いるコンピュータ断層撮影方法
US20080219527A1 (en) Cardiac Region Detection From Motion Analysis of Small Scale Reconstruction
JP2008519636A (ja) 周期運動物体の検査のためのコンピュータ断層撮影方法
US20120294417A1 (en) Image domain based noise reduction for low dose computed tomography fluoroscopy
JP2002034970A (ja) マルチ・スライスct走査の螺旋再構成の方法及び装置
Hsieh Tomographic reconstruction for tilted helical multislice CT
US8315351B2 (en) System and method for tomographic reconstruction utilizing circular trajectory and scanogram to reduce artifacts
US6999550B2 (en) Method and apparatus for obtaining data for reconstructing images of an object
Zbijewski Model-based image reconstruction in X-ray Computed Tomography
Fränkel et al. Total variation regularization in digital breast tomosynthesis
Buzug Three-dimensional fourier-based reconstruction methods
Parallel et al. Principles of Computerized Tomographic Imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090611

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110830

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20111125

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20111202

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120224

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130207

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

Free format text: PAYMENT UNTIL: 20160215

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees