JP2008505676A - 心臓の円錐ビームct再生における筋状アーチファクト低減 - Google Patents
心臓の円錐ビームct再生における筋状アーチファクト低減 Download PDFInfo
- Publication number
- JP2008505676A JP2008505676A JP2007519927A JP2007519927A JP2008505676A JP 2008505676 A JP2008505676 A JP 2008505676A JP 2007519927 A JP2007519927 A JP 2007519927A JP 2007519927 A JP2007519927 A JP 2007519927A JP 2008505676 A JP2008505676 A JP 2008505676A
- Authority
- JP
- Japan
- Prior art keywords
- data
- noise reduction
- filter
- reduction factor
- radiation
- 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.)
- Pending
Links
- 230000009467 reduction Effects 0.000 title claims abstract description 47
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 16
- 210000003205 muscle Anatomy 0.000 title 1
- 230000005855 radiation Effects 0.000 claims abstract description 42
- 230000003044 adaptive effect Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 28
- 238000001914 filtration Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 6
- 238000013170 computed tomography imaging Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 210000003484 anatomy Anatomy 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008929 regeneration Effects 0.000 description 4
- 238000011069 regeneration method Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007408 cone-beam computed tomography Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 239000013585 weight reducing agent Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/503—Apparatus 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (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)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- Cardiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Pulmonology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
診断用画像形成システム(10)では2次元の投射データがデータメモリ(30)に収集される。このデータは、選択された心臓のフェーズで収集されたデータセットにソートされる。再構成プロセッサ(38)は、投射データを並列放射線フォーマットに再構成する。適応フィルタ(70)は、各読取値の計算された光子ノイズに基づき任意のノイズ低減係数が各フィルタチャネルで異なることを前提に、異なる複数フィルタチャネルの各々により並列放射線フォーマットデータをフィルタ処理する。畳込器(78)は、各フィルタチャネルによりフィルタ処理後のデータを畳み込む。ノイズ低減係数プロセッサ(84)は、角θ∈[0,π]で各ボクセルvに寄与する読取値の加重平均のために生じる実際のノイズ低減係数を判定する。重み付けプロセッサ(90)は、当該チャネルのノイズ低減係数及び当該実際のノイズ低減係数に基づき各フィルタチャネルにつき畳込データを重み付ける。背面投射器(98)は、人が読取可能なフォーマットで表示可能な呈示画像を形成するため投射情報の合算された加重和を背面投射する。
Description
本願は、診断用の画像形成技術に関する。また、本発明は、対象物の心臓コンピュータ断層撮影画像形成における特定の用途を見出すものであり、これに特に関連して後述するものである。但し、本発明は、他のタイプのコンピュータ断層撮影画像形成、単光子放射型コンピュータ断層撮影法(SPECT)、陽電子放射型断層撮影法(PET)、3次元x線画像形成その他における用途を見出しうるものでもある。
一般に、コンピュータ断層撮影システムは、x線源と検査すべきオブジェクトの周りを回転するx線検出器とを有する。複数の方向から、当該オブジェクトはそのx線源からのx線ビームが照射される。これらそれぞれの方向において、x線検出器は、当該オブジェクトを通過したx線放射を受け、問題となっている方向の減衰特性を形成する。この減衰特性は、x線の拡散及び当該問題となっている方向における当該オブジェクトを通じるx線の経路に沿うx線の吸収による当該オブジェクトにおける入射x線の減衰を示す。
このCT画像の品質は、しばしば、過剰なx線量子ノイズに起因する筋状のアーチファクトにより悪化する。フィルタ処理の施された背面投射再生の処理において、全てのボクセルは、当該照射窓において利用可能な光線により再生される。概して円錐形ビームによる再生アルゴリズムの構成においては、投射空間における各サンプルが画像領域において直線にマップ化される。近接したラインの中の正及び負の寄与の度合いは組み合わせられ、最終的画像に直線は現われない。この寄与の値がノイズを含むとき、当該再生処理はその正及び負の寄与の値を正しく組み合わせることができず、線又は線分例えば筋のようなものとなる。この筋により大幅に画質が低下し、当該領域の構造を消してしまう。
心臓の円錐ビーム画像は、フェーズ選択型アルゴリズムを用いて再生される。普通は、当該心臓の特定のフェーズが心臓画像発生のために選ばれる。当該選択されたフェーズに時間的に近いところで、すなわち同じ心臓フェーズに対応するものの心臓の拍動サイクルが異なる時間的ポイントで得られたデータだけが、再生処理に用いられる。スキャンパラメータや患者の心拍数及びその変化、心臓の通門の窓の幅や位置に基づいて、ボクセルの各々の再生のためにサイクルの変数が用いられる。当該心臓の通門の窓及び各ボクセルの照射窓の範囲内にある時間期間により、再生ボリュームにおける各ボクセルにつき利用可能な冗長データの量が決まる。これら時間期間において、得られるデータは、不均一な角度分布を呈する可能性がある。例えば0から180度まで、唯1つのデータポイントがボクセルの再生に寄与する角度がある可能性がある一方で、他の角度については、3つのデータポイントによる寄与があると目される。このような不均等な分布は、不均等なノイズをもたらし、パターンノイズ筋をもたらす。
過剰な光子ノイズに起因する筋状アーチファクトを抑制するための幾つかの対策が提案されている。或る対策は、測定された減衰特性における幾つかのポイントに適応フィルタを適用するものである。このようなフィルタは、専ら、患者の解剖学的構造及び照射線量の変動によるノイズ寄与の差に適応可能なものである。このフィルタは、高い吸収角度において比較的高いx線光子ノイズについて補償するように変形される。フェーズ選択型アルゴリズムにおいて用いられるデータの均質でない分布についても、また再生の処理についても、考慮に入れられていない。
読取値(測定値)のノイズや再生アルゴリズムのフェーズ選択型の性質を考慮した心臓の円錐ビーム画像形成における筋状アーチファクトを抑制する技術の必要性がある。本発明は、前述した制約その他を克服する方法及び装置を狙いとするものである。
本発明の一態様によれば、診断用画像形成システムが開示される。或る手段は、投射データを並列放射線フォーマットに構成する。或る手段は、この並列放射線フォーマットのデータを、異なる複数のフィルタチャネルの各々でフィルタ処理し、各チャネルが異なるノイズ削減係数を持つものとしている。或る手段は、当該フィルタチャネルの各々でフィルタ処理されたデータを畳み込む。或る手段は、角度θ∈[0,π]において各ボクセルに寄与する読取値の加重平均を行う。或る手段は、当該加重平均に基づいて実際のノイズ低減係数を判定する。或る手段は、当該チャネルのノイズ低減係数及び当該判定された実際のノイズ低減係数に基づいてフィルタチャネルの各々について当該畳み込まれたデータの重み付けをする。
本発明の他の態様によれば、円錐ビームCT画像形成システムにおける筋状アーチファクトを補正する方法が開示される。投射データは、並列放射線フォーマットに構成される。この並列放射線フォーマットのデータは、各々が異なるノイズ低減係数を有する異なる複数のフィルタチャネルによりフィルタ処理される。このフィルタチャネルの各々によりフィルタ処理されたデータは、畳み込まれる。角度θ∈[0,π]において各ボクセルに寄与する読取値は、加重平均化される。実際のノイズ低減係数は、この加重平均に基づいて判定される。当該フィルタチャネルの各々についての当該畳み込まれたデータは、チャネルのノイズ低減係数及び当該判定された実際のノイズ低減係数に基づいて重み付けされる。
本願の1つの効果は、筋状のアーチファクトが低減することにある。
もう1つの効果は、ボクセルごとのノイズ補正を回避する技術の簡単化にある。
他の効果は、心臓その他のゲート型画像形成技術による画像の改善にある。
数多くの付加的な効果及び利点は、以下の好適実施例の詳細な説明を読むことにより当業者に明らかとなる。
本発明は、様々な構成要素及び構成要素の様々な構成に具体化することが可能であり、また、様々な処理動作及び処理動作の様々な構成に具体化することが可能である。図面は、好適実施例を例示することを目的としているに過ぎず、本発明を限定するものと解釈してはならない。
図1を参照すると、画像形成システム10は、検査領域16に向けられる放射線ビーム、好ましくは円錐形又は楔形ビームを生成する放射線源14を有するコンピュータ断層撮影スキャナ12を含む。この放射線ビームは、相互干渉し、検査領域16に配される画像形成対象物の関心領域を横断するので部分的に吸収され、当該検査領域を通過することにより当該放射線の空間的に変化する吸収が生じることとなる。放射線検出器18、好ましくは2次元ディテクタは、放射線が検査領域16を通過した後に吸収ないしは減衰された放射線を検出する。線源14と検出器18の放射線検出素子の各々との間の経路は、放射線として示される。
好ましくは、放射線源14はx線の円錐ビームを生成する。放射線源14及び検出器18は、回転台20上に対向する形で実装され、検出器18が放射線源14からのx線を連続的に受けるようにしている。線源14及び検出器18は、回転台20において検査領域16に対し連続的に回転するので、複数の回転状態において視界(像)が得られる。各視界又は2次元アレイのデータは、検出器18の検出素子の同時サンプリングにより収集される線源14に頂点を有する円錐の放射線を表す。代表的なものとして、対象支持部26は、モータ手段28により軸方向すなわちz方向において線形に移動可能なものとされる。
オプションとして、円錐ビームコンピュータ断層撮影投射データは、(i)対象支持部26が各軸方向スキャンにおいて静止し軸方向スキャン状態の一方と他方との間では線形に進むようにして、或いは(ii)当該対象支持部が螺旋状の軌道を規定するように連続的に移動するようにして、複数の回転状態において得られる。放射線検出器18の検出素子の出力は、データメモリ30に記憶される電気的に得られる積分された減衰投射値μd0に変換される。各投射データμd0は、放射線源14から検出器18の検出素子の対応の1つまでのラインに沿った減衰の線積分に対応する。
代表的な円錐ビームの形状寸法について、線積分インデックスは、当該読取値を測定するために用いられる検出素子に対応するのが普通である。但し、線積分インデックスは検出素子番号との直接の対応関係を欠く可能性もあることを意図している。かかる直接の対応関係の欠如は、例えば、構成される投射の一方と他方との補間に起因する可能性がある。
心臓モニタ32は、患者の心臓周期を監視する。選別手段34は、選択された心臓フェーズの各々において収集されたデータセットに当該減衰データをソート(選別又は並べ替え)する。これらデータセットは、フェーズメモリ36に記憶される。好ましくは、再構成プロセッサ(re-binning processor)38は、円錐からのデータを平行ビームの形状寸法に構成し直して1組の平行な画像にする。各画像は、等距離のπラインを含む。このπラインは、軸平面、すなわちスキャンFOVに交差する回転軸に垂直な面に含まれ、基準座標θπ(伝搬角度∈[0,π])と1(アイソセンタからの距離)とにより特徴づけられる線積分として規定される。特に短い時間的窓により規定される心臓フェーズに対しては、1つの心臓フェーズのデータは、複数の回転状態の各々における短い円弧部分において収集されたデータに対応する。ここでの円弧部分のデータは、個々に、小さくフルのデータセットとすべきものである。複数の回転にわたる各選択された心臓フェーズのデータが連結される場合、各πラインに対応する放射線の数は変わる。或る幾つかのπラインは、3つ又は4つの冗長性減衰データ値を有することができ、かかる値は、加重平均後に冗長性ノイズ低減のために低いノイズレベルを呈しうるものとなる。他の幾つかのπラインは、唯1つのデータ値を有することができ、かかる値は、当該加重平均のために冗長性ノイズ低減が生じうるのでその元のノイズレベルを保つものと想定される。このノイズ低減は、一方のボクセルから他方のボクセルにおいて変化するので、ボクセルに応じたフィルタ処理が必要とされ、計算上複雑となる。以下により詳しく説明するように、このような複雑さは、ボクセルに依存しないNチャネルフィルタ処理を適用し、ボクセルに応じたノイズ低減の影響を、全てのチャネルの寄与程度がボクセルを単位として重み付けされる背面投射ステップに後回しすることによって、回避される。
画像プロセッサ50は、各選択された心臓フェーズのデータからの画像データを処理して画像メモリ52に記憶される対応の3次元画像とする。好適な再生に関する詳細な数学的解説は、Tuyの発行された米国特許第6,104,775号に提示されている。減衰データは、冗長性の差及び他のノイズの差のために調整するよう再生の間に調整される。ビデオプロセッサ54は、画像メモリ52の内容の幾つか又は全てを処理し、3次元レンダリング、選択された画像スライス、最大強度投影、CINEアニメーション又はこれらに同類のものの如き人が視認可能な画像再生情報を生成する。当該視認可能な画像再生情報は、ユーザインターフェース58のディスプレイ56に表示される。かかるインターフェースとしては、パーソナルコンピュータ、ワークステーション、ラップトップコンピュータなどが好ましい。オプションとして、画像メモリ52の選択された内容は、紙に印刷され、不揮発性の電子的又は磁気的記憶媒体に記憶され、ローカルエリアネットワーク又はインターネットにおいて伝送され、又は処理される。好ましくは、放射線専門医又は他の操作者は、入力手段60を介してコンピュータ断層撮影画像形成スキャナ12を制御して、画像化セッションをセットアップし、画像化セッションを変更し、画像化セッションを実行し、画像化セッションを監視し、又はスキャナ12を操作するようにスキャンコントローラ62をプログラムする。
引き続き図1を参照しさらに図2を参照すると、適応フィルタ70は、適応フィルタ補正のタイプ及び程度により互いに異なるN個の独立フィルタチャネルにおいて各選択された心臓フェーズのデータ値をフィルタ処理する。さらに、各チャネルは、データノイズが、例えば冗長な読取値に対する加重平均のために異なる程度で低減されたものとみなすことにより、同様に当該データノイズの関数としての当該データをフィルタ処理する。より詳しくは、光ノイズ算出ルーチン又は手段72は、当該放射線に沿った光子ノイズを計算する。この光ノイズは、パス(すなわち患者の解剖学的構造)に沿った減衰線積分、管電流その他同類のもののようないくつかの要因により影響を受けるのが普通である。オプションとして、患者の解剖学的構造に関するアプリオリがフィルタ処理関数を調整することを補助するために組み込まれる。これは、或る経路に対してはx線光子が通過して当該検出器に到達することが比較的に難しいことが知られていることによる。
特定の角度θ∈[0,π]において、当該角θに寄与する種々の投射からの可変な数のデータポイントが得られる可能性がある。ボクセルvのハイブリッド読取値は、ボクセルvに交差し角θに収められた種々の投射の読取値の加重和である。簡単化された例として、角θにおけるボクセルvに対しては、ボクセルvに寄与する単一のポイントを各々が有する2つの投射がある可能性がある。当該計算をより簡単にするため、冗長データポイントの線積分が概ね同じであると仮定する。さらに、検出器に到達する冗長データポイントに対応する光子の数が概ね同じであると仮定する。例えば、ノイズ概算は、当該冗長データポイントが同じ量の固有ノイズを有することを前提にすることにより行われる。フィルタチャネル選択手段76は、角θにおけるボクセルvの冗長性に対応するノイズ低減係数NRn∈[1,…,N]∈[0,1]を有する適応フィルタチャネルn∈[1,2,…,N]を選択する。上述した例においては、2つのポイントが同じ心臓加重をもって角θにおけるボクセルvに対して寄与するものであるが、ノイズ低減は、√2の割合で行われる。例えば、フィルタチャネル選択手段76は、1/√2のノイズ低減係数NRn∈[1,…,N]を有する適応フィルタチャネルnを選択する。ノイズ低減係数NRn∈[1,…,N]は、各適応フィルタチャネルに特有のものであり、予め各チャネルに割り当てられる。適応フィルタ70は、既知の適応フィルタの1つであり、かかるフィルタは、各チャネルに対する異なるノイズ低減係数の前提を除き、当該フィルタチャネルの各々に対して同様に用いられるものである。この適応フィルタ処理技術は、患者の解剖学的構造及び照射線量により生じるノイズの他、ハイブリッド読取値の概算されたノイズを考慮する。ノイズの抑制は、好ましくはフローティングの閾値を適用することにより過剰ノイズを持ったデータポイントに対して行われる。
データがフィルタ処理された後、畳込器78は、チャネル1ないしNの各々においてランプカーネル(ramp kernel)により1次元の畳み込みを行う。この畳み込みは、検出器の行部分の各々の読取値の並列セットに沿って行われる。データは、2次元データセットを完成させるために角度ごとに畳み込まれ、θ∈[0,π]の角度範囲をカバーする。
図1及び図2を引き続き参照すると、背面投射手段80は、当該畳み込まれたデータの重み付けされた背面投射を行い3次元画像の呈示とする。より詳しくは、加重平均手段82は、各チャネルn∈[1,2,…,N]のハイブリッド読取値に寄与する複数の読取値、例えば角θ∈[0,π]における各ボクセルに寄与する読取値を平均化する。ノイズ低減係数判定手段84は、その加重平均に因るノイズ低減係数を判定する。この加重平均は、普通の加重平均とは異なるものとすることができる。何故なら、同じハイブリッド読取値に寄与する異なる読取値に対して与えられる異なる心臓の重みがあるからである。各読取値に与えられる心臓の重みは、既知の技術の1つによって計算され、普通は、患者の心電図の各インターバル内のゲーティング窓内の位置及び対応する窓幅に応じたものとなる。各読取値につき判定されたノイズ低減係数は、低減係数メモリ86に記憶される。
加重手段90は、判定されたノイズ低減係数に応じて各畳み込まれたデータチャネルからの寄与程度に重み付けする。検索手段92は、角θにおけるボクセルvに寄与する種々の投射の読取値に関連した冗長性の程度すなわちノイズ低減係数を検索する。かかる冗長性の程度に基づいて、加重選択手段94は、特定のボクセル又はボクセルのサブブロックに対する各チャネルに与えるべき重みWn v,pを判定する。すなわち、
Wn v,p=WNAF v,p・Wn(NRv,θ)
である。
Wn v,p=WNAF v,p・Wn(NRv,θ)
である。
ここで、Wn v,pはチャネルnによりボクセルvへ背面投射された投射pの読取値に与えられる重みである。
WNAF v,pは投射pの読取値に適用され当該分野において既知の方法の1つにより計算されたボクセル位置に依存して合成された照射心臓重みである。
Wn(NRv,θ)は、ボクセルに依存した重み関数であり、ここでのNRv,θは、予め計算されたノイズ低減係数である。
合算手段96は、各チャネルからの重み付けされた寄与量の和をとる。背面投射器98は、この加重和を画像メモリ52に背面投射する。各放射線につき、加重選択手段94は、チャネルのうち最も近くに割り当てられたノイズ低減を呈する単一のチャネルを選択することができ、或いは最も重く重み付けされた最も近いノイズ低減と、複数又は全てのチャネルからの値の一部とを組み合わせることができる。
以上、本発明を好適実施例に基づいて説明したが、これまでの詳細な説明を読み理解することにより種々の変形や変更が可能であることは明らかである。本発明は、添付の請求項の範囲又はその等価な範囲内にある限り、このような変形及び変更の全てを含むものと解釈すべきことを意図するものである。
Claims (22)
- 診断用画像形成システムであって、
投射データを並列放射線フォーマットに構成する手段と、
異なる複数のフィルタチャネルの各々により前記並列放射線フォーマットのデータをフィルタ処理する手段と、
前記フィルタチャネルの各々によりフィルタ処理されたデータを畳み込む手段と、
角θ∈[0,π]における各ボクセルに寄与する読取値を加重平均化する手段と、
前記加重平均により生成されるべき実際のノイズ低減係数を判定する手段と、
当該チャネルのノイズ低減係数及び当該判定された実際のノイズ低減係数に基づいて前記フィルタチャネルの各々について当該畳み込まれたデータを重み付ける手段と、
を有するシステム。 - 請求項1に記載のシステムであって、前記フィルタ処理手段は、独立した複数のフィルタチャネルの各々において異なる初期ノイズ低減係数を適用し、当該初期ノイズ低減係数は、予め定められ各フィルタチャネルに割り当てられる、システム。
- 請求項1に記載のシステムであって、前記フィルタチャネルに対して適応フィルタ関数を選択するために各データ値の光子ノイズを計算する手段をさらに含むシステム。
- 請求項1に記載のシステムであって、当該判定された実際のノイズ低減係数に基づいて当該別々にフィルタ処理され畳み込まれたデータ値に与えるべき重みを判定する手段をさらに含むシステム。
- 請求項4に記載のシステムであって、角θにおけるボクセルvに寄与する種々の投射の読取値の当該判定された実際のノイズ低減係数を検索する手段をさらに含むシステム。
- 請求項5に記載のシステムであって、背面投射の前に当該別々にフィルタ処理され畳み込まれたデータ値を合算する手段をさらに含むシステム。
- 請求項1に記載のシステムであって、2次元投射データを収集する手段をさらに含み、当該収集されたデータの各データ値は、各放射線に沿ってなされる再生対象オブジェクトの線積分に関係している、システム。
- 請求項1に記載のシステムであって、当該重み付けされ畳み込まれたデータを3次元的に背面投射する手段をさらに含むシステム。
- 請求項1に記載のシステムであって、選択された心臓フェーズの各々において収集されたデータセットに当該照射データをソートする手段をさらに含むシステム。
- 請求項1に記載のシステムであって、少なくとも選択された複数の時間窓の内にある投射データを得るCTスキャナをさらに含み、このスキャナは、
回転台と、
検査領域を横断する円錐ビーム放射源と、
前記検査領域を通過した後の放射線を検出し、これを前記投射データフォーマットに変換するものであって、前記検査領域に関しての連続的な回転のために前記回転台において前記放射源と相対して位置づけられた放射線検出器と、
当該2次元放射線検出器により受け取られたx線を体積的な呈示画像に再生する背面投射器と、
前記体積的な画像を表示するディスプレイと、
を含む、
システム。 - 請求項1に記載のシステムであって、当該フィルタ処理は、独立した並列処理チャネルにおいて行われる、システム。
- 円錐ビームCT画像形成システムにおける筋状アーチファクトを補正する方法であって、
投射データを平行放射線フォーマットに構成するステップと、
異なる複数のフィルタチャネルの各々により前記平行放射線フォーマットのデータをフィルタ処理するステップと、
前記フィルタチャネルの各々によりフィルタ処理されたデータを畳み込むステップと、
角θ∈[0,π]における各ボクセルに寄与する読取値を加重平均するステップと、
加重平均により生成されるべき実際のノイズ低減係数を判定するステップと、
当該チャネルのノイズ低減係数及び当該判定された実際のノイズ低減係数に基づいて前記フィルタチャネルの各々について当該畳み込まれたデータを重み付けるステップと、
を有する方法。 - 請求項12に記載の方法であって、独立した各フィルタチャネルにおいて異なる初期ノイズ低減補正係数を適用し、当該初期ノイズ低減補正係数が、予め定められ各フィルタチャネルに割り当てられるようにするステップをさらに含む方法。
- 請求項12に記載の方法であって、前記フィルタチャネルに対して適応フィルタ関数を選択するために当該データポイントの光子ノイズを計算するステップをさらに含む方法。
- 請求項12に記載の方法であって、各独立したフィルタチャネルの当該畳み込まれたデータに付与すべき重みを判定するステップをさらに含む方法。
- 請求項15に記載の方法であって、当該重みは、0から1の数である、方法。
- 請求項15に記載の方法であって、
当該判定された実際のノイズ低減係数を取得するステップと、
前記実際のノイズ低減係数に基づいて角θにおけるボクセルvについて当該畳み込まれたデータの重みを判定するステップと、
をさらに含む方法。 - 請求項15に記載の方法であって、最終的画像を再生するために用いられるフィルタチャネルに係る投射を選択するステップをさらに含む方法。
- 請求項12に記載の方法であって、2次元データを収集するステップをさらに含み、当該収集されたデータの各データ値は、各放射線に沿ってなされる再生対象オブジェクトの線積分に関係している、方法。
- 請求項12に記載の方法であって、当該重み付けされた畳み込みデータを3次元的に背面投射するステップをさらに含む方法。
- 請求項12に記載のステップを行うようにプログラムされた画像処理器。
- 診断用画像形成システムであって、
検査領域を横断する円錐ビーム放射線源と、
前記検査領域を通過した後の放射線を検出しこれを電子的投射データフォーマットに変換する放射線検出器と、
当該円錐ビーム投射データを3次元再生画像に再生する画像処理器と、
を含み、
前記画像処理器は、
当該収集されたデータを並列放射線フォーマットに構成するステップと、
異なる複数のフィルタチャネルの各々により前記並列放射線フォーマットのデータをフィルタ処理するステップと、
前記フィルタチャネルの各々によりフィルタ処理されたデータを畳み込むステップと、
角θ∈[0,π)における各ボクセルvに寄与する読取値を加重平均するステップと、
前記加重平均により生成される各寄与読取値のノイズ低減係数を判定するステップと、
各寄与読取値に与えられる心臓の重み及び前記ノイズ低減係数に基づいて前記フィルタチャネルの各々について当該畳み込まれたデータを重み付けるステップと、
を行うようにプログラムされ、
当該システムは、人間が視認可能な画像フォーマットにて当該背面投射されるデータを表示するディスプレイを有する、
システム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US58596804P | 2004-07-07 | 2004-07-07 | |
PCT/IB2005/052079 WO2006006090A1 (en) | 2004-07-07 | 2005-06-23 | Streak artifact reduction in cardiac cone beam ct reconstruction |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008505676A true JP2008505676A (ja) | 2008-02-28 |
Family
ID=34977027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007519927A Pending JP2008505676A (ja) | 2004-07-07 | 2005-06-23 | 心臓の円錐ビームct再生における筋状アーチファクト低減 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7558362B2 (ja) |
EP (1) | EP1765174A1 (ja) |
JP (1) | JP2008505676A (ja) |
CN (1) | CN100563570C (ja) |
WO (1) | WO2006006090A1 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012250043A (ja) * | 2011-06-06 | 2012-12-20 | Toshiba Corp | ノイズ低減方法 |
JP2014023936A (ja) * | 2012-07-30 | 2014-02-06 | Toshiba Corp | X線コンピュータ断層撮影装置及び画像再構成方法 |
US9872661B2 (en) | 2014-11-19 | 2018-01-23 | Toshiba Medical Systems Corporation | X-ray CT apparatus, and image processing apparatus |
JP2021502842A (ja) * | 2017-11-14 | 2021-02-04 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 複数の異なるタイプの投影データのボクセル毎に1回のジオメトリ演算を行う単一のct逆投影器 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5123192B2 (ja) * | 2005-09-22 | 2013-01-16 | ウイスコンシン アラムナイ リサーチ ファウンデーシヨン | 機能的磁気共鳴イメージング用の、画像の取得及び再構成の方法 |
JP4509971B2 (ja) * | 2006-06-09 | 2010-07-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
JP2008006032A (ja) * | 2006-06-29 | 2008-01-17 | Ge Medical Systems Global Technology Co Llc | X線ct装置およびx線ct撮影方法 |
US8378310B2 (en) * | 2009-02-11 | 2013-02-19 | Prismatic Sensors Ab | Image quality in photon counting-mode detector systems |
DE102010013360B4 (de) * | 2010-03-30 | 2017-01-19 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von Bilddaten eines zyklisch sich bewegenden Untersuchungsobjektes |
EP2583241A1 (en) * | 2010-06-21 | 2013-04-24 | Koninklijke Philips Electronics N.V. | Method and system for noise reduction in low dose computed tomography |
US8270561B2 (en) * | 2010-10-13 | 2012-09-18 | Kabushiki Kaisha Toshiba | Motion weighting in computed tomography (CT) with cone angle |
RU2586968C2 (ru) * | 2010-10-27 | 2016-06-10 | Конинклейке Филипс Электроникс Н.В. | Уменьшение уровня шума в низкодозной компьютерной томографии |
CN103299345B (zh) * | 2011-01-10 | 2016-10-19 | 皇家飞利浦电子股份有限公司 | 双能量断层摄影成像系统 |
JP6122269B2 (ja) * | 2011-12-16 | 2017-04-26 | キヤノン株式会社 | 画像処理装置、画像処理方法、及びプログラム |
US10388000B2 (en) * | 2014-09-15 | 2019-08-20 | Analogic Corporation | Noise reduction in radiation image |
US10013780B2 (en) * | 2016-02-29 | 2018-07-03 | General Electric Company | Systems and methods for artifact removal for computed tomography imaging |
CN109543763B (zh) * | 2018-11-28 | 2022-10-21 | 重庆大学 | 一种基于卷积神经网络的拉曼光谱分析方法 |
CN109920020B (zh) * | 2019-02-27 | 2022-10-18 | 西北工业大学 | 一种锥束ct病态投影重建伪影抑制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1033521A (ja) * | 1996-07-23 | 1998-02-10 | Hitachi Medical Corp | X線ct装置 |
JPH11502757A (ja) * | 1996-01-17 | 1999-03-09 | アナロジック コーポレーション | コンピュータ断層撮影用条痕抑制フィルタ |
JP2000135214A (ja) * | 1998-10-29 | 2000-05-16 | Picker Internatl Inc | イメ―ジ再構成 |
JP2003180675A (ja) * | 2001-11-21 | 2003-07-02 | Ge Medical Systems Global Technology Co Llc | 計算機式断層写真法システムの雑音低減の方法及び装置 |
WO2004010383A2 (en) * | 2002-07-23 | 2004-01-29 | General Electric Company | Method and apparatus for deriving motion information from projection data |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4580219A (en) | 1983-05-02 | 1986-04-01 | General Electric Company | Method for reducing image artifacts due to projection measurement inconsistencies |
US5416815A (en) * | 1993-07-02 | 1995-05-16 | General Electric Company | Adaptive filter for reducing streaking artifacts in x-ray tomographic images |
US5625660A (en) * | 1995-06-30 | 1997-04-29 | Picker International, Inc. | Image reconstruction from helical partial cone-beam data |
US5818896A (en) * | 1996-11-18 | 1998-10-06 | General Electric Company | Methods and apparatus for three-dimensional and maximum intensity projection image reconstruction in a computed tomography system |
FR2810141B1 (fr) * | 2000-06-07 | 2002-08-23 | Commissariat Energie Atomique | Procede de reconstruction accelere d'une image tridimentionnelle |
JP4330881B2 (ja) * | 2001-03-29 | 2009-09-16 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 対象画像におけるアーチファクトを減少する方法 |
WO2002093490A1 (en) * | 2001-05-16 | 2002-11-21 | Koninklijke Philips Electronics N.V. | Computed-tomography system with filtering |
US7058440B2 (en) | 2001-06-28 | 2006-06-06 | Koninklijke Philips Electronics N.V. | Dynamic computed tomography imaging using positional state modeling |
DE10143484A1 (de) * | 2001-09-05 | 2003-04-03 | Siemens Ag | Adaptives Filter |
US7187794B2 (en) | 2001-10-18 | 2007-03-06 | Research Foundation Of State University Of New York | Noise treatment of low-dose computed tomography projections and images |
US6775346B2 (en) | 2002-10-21 | 2004-08-10 | Koninklijke Philips Electronics N.V. | Conebeam computed tomography imaging |
-
2005
- 2005-06-23 CN CN200580022988.0A patent/CN100563570C/zh active Active
- 2005-06-23 WO PCT/IB2005/052079 patent/WO2006006090A1/en not_active Application Discontinuation
- 2005-06-23 EP EP05758566A patent/EP1765174A1/en not_active Ceased
- 2005-06-23 JP JP2007519927A patent/JP2008505676A/ja active Pending
- 2005-06-23 US US11/571,660 patent/US7558362B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11502757A (ja) * | 1996-01-17 | 1999-03-09 | アナロジック コーポレーション | コンピュータ断層撮影用条痕抑制フィルタ |
JPH1033521A (ja) * | 1996-07-23 | 1998-02-10 | Hitachi Medical Corp | X線ct装置 |
JP2000135214A (ja) * | 1998-10-29 | 2000-05-16 | Picker Internatl Inc | イメ―ジ再構成 |
JP2003180675A (ja) * | 2001-11-21 | 2003-07-02 | Ge Medical Systems Global Technology Co Llc | 計算機式断層写真法システムの雑音低減の方法及び装置 |
WO2004010383A2 (en) * | 2002-07-23 | 2004-01-29 | General Electric Company | Method and apparatus for deriving motion information from projection data |
JP2006513734A (ja) * | 2002-07-23 | 2006-04-27 | ゼネラル・エレクトリック・カンパニイ | 投影データから運動情報を導く方法及び装置 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012250043A (ja) * | 2011-06-06 | 2012-12-20 | Toshiba Corp | ノイズ低減方法 |
JP2014023936A (ja) * | 2012-07-30 | 2014-02-06 | Toshiba Corp | X線コンピュータ断層撮影装置及び画像再構成方法 |
US9872661B2 (en) | 2014-11-19 | 2018-01-23 | Toshiba Medical Systems Corporation | X-ray CT apparatus, and image processing apparatus |
JP2021502842A (ja) * | 2017-11-14 | 2021-02-04 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 複数の異なるタイプの投影データのボクセル毎に1回のジオメトリ演算を行う単一のct逆投影器 |
Also Published As
Publication number | Publication date |
---|---|
WO2006006090A1 (en) | 2006-01-19 |
US20070195925A1 (en) | 2007-08-23 |
CN1980604A (zh) | 2007-06-13 |
US7558362B2 (en) | 2009-07-07 |
CN100563570C (zh) | 2009-12-02 |
EP1765174A1 (en) | 2007-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008505676A (ja) | 心臓の円錐ビームct再生における筋状アーチファクト低減 | |
JP5198443B2 (ja) | 画像の分解能を高めるシステム及び方法 | |
RU2565507C2 (ru) | Система и способ для улучшения качества изображения | |
JP5199081B2 (ja) | 心臓ct撮影のバンドアーチファクトの抑制 | |
US8666137B2 (en) | Apparatus and method for processing projection data | |
Taguchi | Temporal resolution and the evaluation of candidate algorithms for four‐dimensional CT | |
RU2557466C2 (ru) | Устройство компьютерной томографии | |
JP5133690B2 (ja) | ボクセルに依存する補間を用いる画像再構成 | |
EP2512345B1 (en) | Computed tomography apparatus | |
JP4360817B2 (ja) | 放射線断層撮影装置 | |
US20070248253A1 (en) | Artifact Reduction | |
US20110150305A1 (en) | Method and system for correcting artifacts in image reconstruction | |
US7928727B2 (en) | Adapting acquisition time in nuclear imaging | |
CN103717137B (zh) | 图像生成装置 | |
JP2016198504A (ja) | 画像生成装置、x線コンピュータ断層撮影装置及び画像生成方法 | |
US20090213985A1 (en) | System and method for tomographic reconstruction utilizing circular trajectory and scanogram to reduce artifacts | |
Pien et al. | Conventional and Newer Reconstruction Techniques in CT | |
Manzke et al. | Image reconstruction with multi-slice CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080620 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20101215 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110120 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110630 |