JP4409671B2 - 断層写真画像を作る装置と方法 - Google Patents
断層写真画像を作る装置と方法 Download PDFInfo
- Publication number
- JP4409671B2 JP4409671B2 JP23487499A JP23487499A JP4409671B2 JP 4409671 B2 JP4409671 B2 JP 4409671B2 JP 23487499 A JP23487499 A JP 23487499A JP 23487499 A JP23487499 A JP 23487499A JP 4409671 B2 JP4409671 B2 JP 4409671B2
- Authority
- JP
- Japan
- Prior art keywords
- correction
- factor
- coefficient
- detector
- projection data
- 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 - Lifetime
Links
- 238000004519 manufacturing process Methods 0.000 title 1
- 238000012937 correction Methods 0.000 claims description 178
- 238000001514 detection method Methods 0.000 claims description 65
- 239000011159 matrix material Substances 0.000 claims description 49
- 238000000034 method Methods 0.000 claims description 36
- 239000013598 vector Substances 0.000 claims description 25
- 238000002591 computed tomography Methods 0.000 claims description 23
- 238000005259 measurement Methods 0.000 description 20
- 230000008569 process Effects 0.000 description 15
- 238000004364 calculation method Methods 0.000 description 11
- 230000004907 flux Effects 0.000 description 11
- 238000012545 processing Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 7
- 210000000988 bone and bone Anatomy 0.000 description 6
- 210000003484 anatomy Anatomy 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- OVSKIKFHRZPJSS-UHFFFAOYSA-N 2,4-D Chemical compound OC(=O)COC1=CC=C(Cl)C=C1Cl OVSKIKFHRZPJSS-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004164 analytical calibration Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000004175 ponceau 4R Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 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
-
- Y—GENERAL 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
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Description
【発明の分野】
この発明は全般的に計算機式断層写真(CT)作像、更に具体的に言えば、異なる個別の利得を持つx線検出セルの出力信号を組合わせることによってそのデータに導入される誤差に対して画像データを補正することに関する。
【0002】
【背景】
少なくとも1つの公知のCTの装置の形式では、x線源が扇形ビームを投射し、このビームは、直角座標系のX−Y平面内にある様にコリメートされ、この平面が一般的に「作像平面」と呼ばれる。x線ビームが患者の様な作像される物体を通過する。ビームは、物体による減衰を受けた後、放射検出器の配列に入射する。検出器配列が受け取る減衰したビーム放射の強度は、物体によるx線ビームの減衰に関係する。配列の各々の検出素子又はセルが、この検出器の場所でのビームの減衰の測定値である別々の電気信号を発生する。すべての検出セルからの減衰測定値を別々に収集して、透過分布を作る。
【0003】
公知のCT装置では、x線源及び検出器配列が作像平面の中で、作像しようとする物体の周りをガントリと一緒に回転し、この為、x線ビームが物体と交わるときの回転角度が絶えず変化する。ガントリの1つの回転角度に於ける検出器配列からの一群のx線検出減衰測定値、即ち、投影データが「図」と呼ばれる。物体の「走査」は、x線源及び検出器が1回転する間のガントリの相異なる回転速度又は図の角度で求められた1組の図で構成される。軸方向走査では、投影データが、物体を通り抜ける2次元スライスに対応する画像を構成するように処理される。1組の投影データから画像を再構成する1つの方法は、業界ではフィルタ式逆投影方式と呼ばれている。この方法は、走査からの減衰測定値を「CT数」又は「ハウンズフィールド単位」と呼ばれる整数に変換され、こういう整数を使って、陰極線管表示装置の対応する画素の輝度を制御し、こうして画像を発生する。
【0004】
前に述べた様に、画像は、物体の2次元スライス又は断面に対応する。理論的には、画像が2次元であるから、厚さを持たず、従ってz方向の測定線束及び検出器利得は存在せず、その為、一様である。しかし、現実には、z方向に有限の厚さを持つ検出器から供給される減衰測定値から構成される。検出器の厚さが有限であることにより、一様でない測定線束並びに一様でない検出器利得の惧れを生ずる。検出器の固有の厚さが、正確な2次元画像を構成しようとするときに問題を生ずる。これは、走査される物体の組成(従って測定線束)並びに検出器の利得が、検出器の厚さに互って変化し得るからである。言い換えると、ある場所に於ける減衰測定値が、一部分が骨、そして一部分が空気で構成される場合、実際の減衰測定値に基づいて構成される画像は、骨と空気の組合せを表す。同様に、骨が検出器の厚さ全体に互っているが、検出器利得がこの厚さに互って変化する場合、構成された画像は、骨の組成を正確に反映しない。この為、検出器からの減衰測定値は、画像の厚さに互る検出器利得及び測定線束の変動を考慮して補正しなければならない。
【0005】
一般的には2−D検出器の名前で知られている検出器を使う第3世代CT装置では、こういう問題が更に著しくなる。2−D検出器は複数個の列及び行の検出セルで構成され、異なるx位置で整合した検出セルが行を形成し、異なるz位置で整合した検出セルが列を形成する。この様な2−D検出器を持つCT装置は、多重スライス装置と呼ばれる場合もあるが、こういう装置では、多数の行並びに/又は列の検出セルの測定値を組合わせることにより、画像が形成されることがある。1列のセルの加算により「マクロ・セル」ができ、同じ平面内のマクロ・セルが「マクロ行」を作る。検出セルの測定値がデータ収集装置に対する入力として供給される。組合わされる検出セルの測定値が、相異なる個別の利得を持つ検出セル並びに異なる線束から得られたものである場合、組合わされた測定値は、個別の検出セル測定値の加重和を表し、異なる検出セルの利得及び線束によって、加重が異なる。検出セルの利得及び線束の違いによって導入される合成の誤差は、物体に依存していて、標準的な較正によって取除くことが出来ない。
【0006】
この様なデータから画像を更に正確に作る為に、異なる個別の利得を持つx線検出セルからのデータを組合わせることにより誤差を推定することが知られている。この誤差を推定する装置が、この出願の被譲渡人に譲渡された米国特許第5,734,691号、発明の名称「計算機式断層写真装置に於ける検出器Z軸利得非一様性の補正」及び1995年1月23日に出願され、この出願の被譲渡人に譲渡された継続中の米国特許通し番号08/376,813号、発明の名称「CT装置に対する検出器Z軸利得の補正」に記載されている。しかし、こういう補正装置は多数の点で改良することが出来る。こういう装置は、z方向の検出器利得及び測定線束の変動を補正することが出来るが、複雑度並びに要求される計算資源が補正の次数、即ち、補正しなければならないzに於ける変動の複雑度と共に増加する。この関係により、画像の品質を最適にすることと、再構成時間及びコストを最適にすることとの間に競合が起こることがある。
【0007】
更に、こういう装置は、多数のz軸補正係数を持っていて、各々の補正係数が、補正の次数、並びに補正を適用すべき検出器セグメントの群に関係する。具体的に言うと、こういう装置は、検出器をx方向に重なり合う多数のセグメントに分割して、検出器の各セグメントに対して異なる補正係数を定める。その後、各セグメントに互って、x方向にz軸補正を正規化し、各セグメントに対する補正が、そのセグメント内にある悉くの検出セルに適用される。この様にして、重なり合うセグメントは2回処理され、それを平滑して、セグメントからセグメントへの移り変わりが一様になる様にする。しかし、検出器の大きなセグメントの一様な補正は、セグメント内にある個別の検出セル又はセルの列の各々に対して適用される補正にスキューを生ずる傾向がある。更に、複数個の重なり合う列を含めることは、計算時間及び労力の重複投資を必要とする。
【0008】
従って、必要とされるz軸補正の次数に関係なく、データから画像を更に正確に且つ効率よく作ることが好ましい。又、装置のコストをあまり増加せずに、この様な作像を行うことが望ましい。
【0009】
【発明の要約】
この発明は、一様でない患者の解剖学的な部分が、検出セルの利得の差の変動と相互作用することによって生ずる投影強度のz軸勾配から起こる画像の人為効果(アーティファクト)の補正の計算効率並びに画像の品質を改善する装置と方法を提供する。この発明では、x線を放出するx線源と、複数個の検出セルを持ち、このx線を受取ってx線の減衰を表す投影データを出力として発生する検出器と、投影データの間の非一様性から生ずる画像の人為効果を補正する為に、投影データから抽出された誤差係数に補正係数を適用する画像プロセッサとを有する、投影データから物体の断層写真画像を作る装置を提供する。補正係数は、一様でない投影データの複雑度に関係なく、略同じままである計算の複雑度を持つことが有利である。投影データがべクトル式であって、補正係数が、一様でない投影データの複雑度に関係する可変係数を持つ行列乗数であることが好ましい。更に検出器は複数個のチャンネルを持っていてよく、1列の検出セルがこの複数個のチャンネルの各々を形成する。補正係数が複数個のチャンネル補正係数で構成され、複数個のチャンネルの各々が、各々のチャンネル補正係数を夫々決定する為に使われる対応するセグメントの中心チャンネルである。補正係数は画像を構成する際に適用されることが好ましい。
【0010】
この発明の別の一面では、複数個のx線減衰測定値から物体の断層写真画像を作る装置が、x線を放出するx線源と、複数個の検出セルを持ち、複数個の検出セルの各々がx線を受取ってx線の減衰を表す出力信号を発生する検出器と、出力信号の間のzー勾配の変動を考慮する最終補正を、初期誤差係数及び最終補正係数に基づいて計算する画像プロセッサとを有する。前に述べた装置と同様に、最終補正係数の計算の複雑度は実質的に次数に無関係である。画像プロセッサが、計算が複雑で次数に関係する初期補正係数を計算し、この初期及び最終補正係数は装置の較正に基づく。最終補正は、出力信号及び最終補正係数の組合せに基づいて決定される。最終補正が、最終補正係数と初期誤差係数を組合わせる行列操作として計算される。
【0011】
更に検出器が複数個のチャンネルを持ち、1列の検出セルが、複数個のチャンネルの各々を形成する。最終補正係数は複数個のチャンネル補正係数で構成され、複数個のチャンネルの各々が、各々のチャンネル補正係数を夫々決定する為に使われる対応するセグメントの中心チャンネルである。最終補正係数VFは次の式によって定義することが出来る。
【0012】
【数5】
【0013】
ここで
VFi,j=V(i)N'+1,j; j=1,…,N
N=2N’+1
であって、V(i)NxNは次の式で表される補正行列である。
【0014】
【数6】
【0015】
VFiがこのとき、行列VFから行iを抽出することによって得られる1×Nベクトルである。
この発明の更に別の一面では、計算機式断層写真装置で検出セルの利得誤差に対して投影データを補正する方法を開示する。装置は、複数個の検出セルを持つ検出器に向かってx線ビームを放出するx線源を持ち、検出セルの出力信号から投影データが発生される。この方法は、一様でない投影データから生ずる画像の人為効果を補正する為に、投影データから抽出された初期誤差係数に補正係数を適用することを含み、補正係数の複雑度は、一様でない投影データの複雑度に関係なく、略同じままである。投影データがベクトル式であり、補正係数が、一様でない投影データの複雑度に依存する可変係数を持つ行列乗数であることが望ましい。更にこの方法は、補正の対象とする複数個のチャンネルを決定することを含んでいて良く、1列の検出セルが複数個チャンネルの各々を形成する。更に補正係数は複数個のチャンネル補正係数で構成すことができ、複数個のチャンネルの各々が、各々のチャンネル補正係数を夫々決定する為に使われる対応するセグメントに対する中心チャンネルである。補正係数が、画像を構成する際に適用されることが好ましい。
【0016】
この発明の別の一面では、計算機式断写真装置で検出セルの利得誤差に対して投影データを補正する方法を開示する。装置が、複数個の検出セルを持つ検出器に向かってx線ビームを投射するx線源を持ち、検出セルの出力信号から投影データが発生される。この方法は、実質的に次数に無関係な計算の複雑度を持つ最終補正係数に基づいて、出力信号中のzー勾配の変動を考慮に入れる最終誤差係数を決定することを含む。最終誤差係数が、装置の較正の間に決定された初期誤差係数と最終補正係数の組合せに基づいて決定される。
【0017】
更にこの方法は、補正の対象にする複数個のチャンネルを決定することを含んでいて良く、1列の検出セルが複数個のチャンネルの各々を形成する。更に、最終補正係数が補正の対象とされる複数個のチャンネルの各々に対して適用され、最終補正係数は複数個のチャンネル補正係数で構成される。最終補正係数内の複数個のチャンネルの各々は、各々のチャンネル補正係数を夫々決定する為に使われる対応するセグメントの中心チャンネルである。最終補正係数eiは次の関数によって定義することが出来る。
【0018】
【数7】
【0019】
ここで
VFi,j=V(i)N'+1,j; j=1,…,N
N=2N’+1
であって、V(i)NxNは次の式で表される補正行列である。
【0020】
【数8】
【0021】
【発明の詳しい説明】
図1及び2について説明すると、計算機式断層写真(CT)作像装置10が、「第3世代」CTスキャナを表すガントリ12を含むことが示されている。ガントリ12がx線源14を持ち、これがガントリ12の反対側にある検出器配列18に向かってx線ビーム16を投射する。検出器配列18は、検出素子20によって形成され、これらが一緒になって医療の患者22を通過する投射されたx線を感知する。各々の検出素子20が、入射したx線ビームが患者22を通過したときのその強度、従ってビームの減衰を表す電気信号を発生する。x線投影データを収集する走査の間、ガントリ12及びその上に装着された部品が回転中心24の周りを回転する。
【0022】
ガントリ12の回転及びx線源14の動作がCT装置10の制御機構26によって制御される。制御機構26が、x線源14に対する電力及びタイミング信号を供給するx線制御器28と、ガントリ12の回転速度及び位置を制御するガントリ・モータ制御器30を含む。制御機構26にあるデータ収集装置(DAS)32が、検出素子20からのアナログ・データを標本化し、そのデータをこの後の処理の為にディジタル信号に変換する。画像再構成装置34が、DAS 32からの標本化されてディジタル化されたx線投影データを受取り、高速画像再構成を実施する。再構成された画像がコンピュータ36に対する入力として印加され、このコンピュータが大量記憶装置38に画像を記憶する。
【0023】
コンピュータ36は、キーボードを持つコンソール40を介して、オペレータからの指令及び走査パラメータをも受取る。関連する陰極線管表示装置42が、コンピュータ36からの再構成された画像及びその他のデータをオペレータが観察することが出来る様にする。オペレータから供給された指令及びパラメータが、DAS 32、x線制御器28及びガントリ・モータ制御器30に対する制御信号及び情報を供給する為に、コンピュータ36によって使われる。更に、コンピュータ36が、患者22をガントリ12内で位置ぎめする為にモータ式テーブル46を制御するテーブル・モータ制御器44を作動する。
【0024】
図3について説明すると、x線ビーム16が源14(図3には示していない)の焦点スポット50から出てくる。x線ビーム16が患者の手前にあるコリメータ52によってコリメートされ、検出器配列18の検出セル20に向かって投射される。一般的に「扇形ビーム平面」と呼ばれる平面54が、焦点スポット50の中心線及びビーム16の中心線を含んでいる。
【0025】
図4は、スイッチ(例えば電界効果トランジスタ(FET))102に結合された1列の検出セル100を示している。検出器の列100は、1列に配置された複数個の検出セルで構成されている。図面に示していないが、完全な検出器は、z軸に沿って検出セルの行を形成する複数個の検出器列で構成されている。前に説明した様に、各々の検出セルが、入射するx線ビームが患者を通過するときのその強度、従ってビームの減衰を表す電気信号を発生する。各セルの出力がFET 102を介してプリアンプ104に供給され、このプリアンプが増幅された信号を、DAS 32(図2)の一部分であるアナログ・ディジタル変換器106に供給する。ディジタル化された信号が、更に処理して画像を再構成する為に、画像再構成装置34に供給される。
【0026】
動作について説明すると、FET 102が各々の検出セルの列からプリアンプ104への出力信号の供給を制御する。例えば、FET 102がスイッチ制御集成体(図面に示していない)の制御の下に「開き」並びに「閉じる」。特定のFETが閉じると、対応する検出セルからの出力信号がプリアンプ104に供給される。FETが開くとき、このセルからプリアンプ104には信号が供給されない。
【0027】
FET 102は、特定のサンプル時間の間、1つ又は1つより多くのセルを付能することが出来る。例えば、1列内の1つの検出セルが、各々のサンプル時間の間に付能されることがある。各々のサンプル時間の間、2つのセルが付能されてもよい。プリアンプ104は、こういう信号の増幅された出力をA/D変換器106に供給する。
【0028】
各々のサンプル時間の間に各チャンネルで作動されるセルの数が、再構成しようとする画像のスライス寸法によって決定される。例えば、図5に示す様に、1列の中に16個の検出セルが配置される。図5では水平に示してあるが、図5のセルは図4に示した列に対応することを承知されたい。更に、1列の中に16個より多くの又は少ないセルを配置することが出来ることを承知されたい。1番上の列110は、寸法が4×1.25 mmである画像スライスに対するセル出力に対応する。1番下の列116は、寸法が4×5.00 mm である画像スライスに対するセルの組合せに対応する。
【0029】
薄いセル(例えば4×1.25 mmのスライス)では、検出セルの加算は実施しない。言い換えると、図5で、4つの付能された検出セルの各々の出力が、DAS 32に送られる別々の信号を構成する。一層厚手のスライス(例えば4×2.50 mmのスライス)では、検出セルの加算が実施される。具体的に言うと、4×2.50 mmのスライスでは、2つの隣接セル(陰影線で示す列2(112))の加算出力を表す4つの別々の信号が、DAS 32に送られる。4×3.75 mmのスライス(列3(116))では、3つのセルが加算され、4×5.00 mmのスライス(列4(118))では、4つのセルが加算される。こういう加算は一層厚手のスライスのデータを収集するときに利用される。これは、一層厚手のスライスでは、前に述べた様に、検出セルの出力、従って投影データを加算することにより、適切なカバーが得られ、処理時間を短縮することが出来るからである。
【0030】
しかし、検出セルの出力又は投影データを加算するとき、各々の検出セルが異なる利得を持つことにより、加算された信号に誤差が導入される。検出セルの利得が違うことを考慮せずに、検出セルの出力を加算すると、その結果得られる信号(例えばDAS 32のA/D変換器106から出力されるディジタル化信号)には、利得の違いによる誤差が含まれる。補正をしないと、この誤差により、再構成された画像に人為効果が生ずることがある。更に、走査される物体が変化する解剖学的な部分を持つとき、検出セルの出力、即ち、測定されるx線束に変動が起こる。こういう線束の変動は、それらが走査に関係していて、変化し得るセルの利得と組合わさったとき、利得の較正によって補正することが出来ない誤差を生ずるので、特にz方向では問題になる。従って、検出セルの利得及び線束の変動が、組合わさったとき、再構成される画像に人為効果を招く原因になり得る2つの重要な因子である。
【0031】
信号又は投影データを組合わせることから生ずる誤差に対して投影データを補正する装置が、この出願の被譲渡人に譲渡された米国特許第5,734,691号、発明の名称「計算機式断層写真装置に於ける検出器z軸利得の非一様性の補正」並びに1995年1月23日に出願され、同じくこの出願の被譲渡人に譲渡された係属中の米国特許出願通し番号08/376,813号、発明の名称「CT装置に対する検出器z軸利得の補正」に記載されている。こういう装置を使うと、z方向の個別のセル利得の変動を考慮に入れる初期補正係数を決定する較正が実施される。初期補正係数は希望する補正の次数に関係する。その為、夫々異なる補正の次数に対して1つずつの多数の初期補正係数が出てくることがある。更に、検出セルの利得の変動の組合せ、従って初期補正係数と、そのスライス内での線束とに基づいて、各々のスライスに対して最終補正係数が決定される。初期補正係数と同様に、最終補正係数が次数に関係し、その為、希望する補正の次数毎に1つずつ、多数の最終補正係数が出てくることがある。その後、最終補正係数が画像を再構成する際に、投影データに適用され、画像の人為効果を減らす。更に、こういう方法では、補正の計算の複雑度が、補正の次数の増加と共に増加する。
【0032】
更に、初期補正係数及び最終補正係数の両方が、複数個のセグメント補正係数で構成される。セグメント補正係数は、検出器を多数の重なり合うセグメントに分割することにより、x方向で、初期補正係数及び最終補正係数を平滑する。セグメント補正係数は、セグメント内の各チャンネルに対して決定され、重なり合う検出チャンネルは、あるセグメントと次のセグメントとのフェザリングによるセグメント補正係数を受取る。こうして、最終補正係数は、検出器に互って滑らかな補正をしようとする。
【0033】
こういう装置では、補正の次数を変えることが出来ることが重要である。これは、補正の次数を高くすると、更に複雑な誤差を補正することが出来るが、計算の複雑度と、画像を再構成する為に必要な時間と費用が増加するという欠点がある。この為、画像の品質と計算の時間及び費用との間の折り合いを必要とする画像再構成の条件が付けられることがある。
【0034】
他方、この発明は、画像の品質又は計算の時間と費用を犠牲にせずに、次数に無関係な最終補正係数を作る装置を有利に提供する。最終補正係数の複雑度が、従来の補正技術に比べて著しく低下し、画像再構成の時点で完全に透明になる。この発明の装置は、最終補正係数に次数に関係する計算の複雑度を必要とする従来技術と対照的に、最終補正係数を流線形にする。この発明の装置は、従来の様に再構成過程に複雑度を入れるのとは対照的に、最終補正係数を決定するのに利用される較正過程に補正の次数の複雑度を有利に移し替える。典型的には比較的長い期間、例えば6ヶ月に互って1回だけ較正が実施されるから、較正過程及び最終補正係数に関係する高くなった複雑度が、夫々の画像の最終補正及び再構成の複雑度に悪影響を与えない。この為、この発明の装置は、最終補正及び画像再構成の効率を最大にする。
【0035】
この発明の装置は、フェザリングの為に必要とされる計算の様な、最終補正のすべての重複計算を有利に取除き、それらの計算を較正に入れる。言い換えると、初期補正係数を決定する較正過程が、検出器利得及び検出器の線束の変動を補正するすべての計算の複雑度を包括する。これによって、各々の個別の走査の投影データ、従ってその線束の変動を、最終補正係数に詰め込んで、投影データを補正する為に使われる最終補正係数を決定することが出来る。この発明は、較正過程で初期補正係数を決定するときに、最終補正係数の補正の次数を考慮に入れることが出来ることを認識する。この為、最終補正係数に基づいて補正を決定する計算は、補正の次数に関係なく、同じであり、計算のある定数だけが次数によって変化する。従って、投影データの最終補正の計算の速度及び効率が、画像の品質を保ちながら高くなる。
【0036】
更に、この発明は、最終補正係数を定めるのに移動するセグメントを利用することにより、改良された画像の品質を提供する。移動するセグメントは、そのチャンネルに対する補正を決定する為に、最終補正係数に含まれるチャンネル毎に、それを中心とする。こうすることにより、この発明はチャンネル毎の補正を同一の形で取扱い、これは画像の品質を改善することにつながる。その補正を決定する為に各々のチャンネルの中心合せをすることにより、この発明は従来の一定の重なり合うセグメントに関連するチャンネル毎に取扱いを変えることによる不連続の問題を避ける。この為、最終補正係数が、移動するセグメントを利用することによって、各々のチャンネルに対する補正を最適にする。
【0037】
更に具体的に言うと、この発明の装置は、較正工程を変更して、前に述べた初期及び最終補正係数に対応する行列の畳込み又は行列乗数の形で、z−勾配の形を書き直すプロセッサ・アーキテクチュアを活用する。行列の寸法は、最小自乗計算でずっと多数のデータ点から少数のパラメータの推定値を計算するときの数値安定性の理由で、zの補正の次数の影響を受けるが、行列の寸法は画像の品質の検討から前以って決定することが出来る。一旦寸法が選ばれれば、行列の係数だけがzの次数に関係する。
【0038】
更に、この発明は、固定の重なり合うセグメントとは対照的に、滑りセグメントを利用し、セグメントの中心チャンネルだけが最終補正行列に利用される。所定のチャンネルに対する行列補正を定める為に関心のあるセグメントを滑らせることにより、セグメントの重なりをなくす様な、1つの最終的な補正行列が定まり、最終補正の計算の回数が約20%減少し、その構成を著しく流線形にする。
【0039】
滑りセグメント方式に基づいて1つの最終補正行列を構成すると共に、次数に関係する複雑度を較正過程に一体化したことの組合せの結果として、zー勾配補正のすべての複雑度が較正過程に移し替えられる。その結果、走査時間に於ける補正処理は非常に流線形になって効率がよくなる。
この発明の装置は、当業者であれば容易に判る様に、色々な実施例を持つことが出来る。特定の装置の一例を次に述べるが、この例はこの発明を制約するものと解してはならない。4つの検出器がz方向に組合わされて、5 mmのスライスを構成すると仮定すると、アルゴリズムは次に述べる様に作用する。組合わされる4つの検出器が、個別の利得gk(k=1,2,3及び4)を持つ。各々の個別の検出器から見たx線強度はIk(ここでk=1,2,3,4)である。Yと記した測定データは次の様なモデルにすることが出来る。
【0040】
【数9】
【0041】
検出器の利得を正規化すると、利得が正規化されたデータImは次の式で与えられる。
【0042】
【数10】
【0043】
ここでGは対象とするマクロ行の平均利得である。即ち、
【0044】
【数11】
【0045】
求めようとする測定値をIと記すと、これは次の様になる。
【0046】
【数12】
【0047】
各々の個別の検出器の利得は次の様に表される。
gk=G+δgk (5)
ここでδgkがgkの残りの部分である。δgkの物理的な意味は、検出器の利得変動である。式(2)と共に式(5)を用い、Gを因数として使い、式(4)を利用すると、次の様になる。
【0048】
【数13】
【0049】
式(6)は、真の信号I、測定されたデータから導き出された信号Im及び検出器のz軸利得変動による誤差を関係付けている。
理想的な投影Pi及び測定投影Pを次の式
【0050】
【数14】
【0051】
但し
【0052】
【数15】
【0053】
この式が、1つの形式としても、別の形式でも、zー勾配補正アルゴリズムの実質的な出発点である。信号Ikに対する次数を制限した展開を特定することにより、直接的に構成することが出来る補正式が得られる。z方向の1次では、
【0054】
【数16】
【0055】
他方
【0056】
【数17】
【0057】
従って1次では、次の様になる。
【0058】
【数18】
【0059】
この式は、1次zー勾配補正を構成する根拠になる。
更に、Nチャンネルで構成される検出器バンクのセグメントを考え、[mi,j]pxqを素子mijを持つp行q列の行列を表すとする。装置較正過程が、正規化公称利得の行列[GNchannel x 16]を定める。ここでNchannelは、zー勾配補正の対象となる検出器に互るチャンネルの総数、16は検出器にある検出セルの行の数(zに於ける行の数)である。こういう公称利得は、前に示した式(10)のδg/Gと密接な関係にあるが、それとは幾分異なる。その違いは、正規化過程が、各セルの利得をx方向で計算したセル利得の平均値によって除算することによるものである。次に、所定のセグメントに関連するGの部分行列並びに一緒に「加算」されるK個の所定の行(16の内の)を考えると、勾配利得係数(SGF)は、Nチャンネルに対し、次の式によって定義される。
【0060】
【数19】
【0061】
ここで[Z]Kx1は、マクロ行を定義する為に一緒に平均を取るK個のセル行に関係するzー基本関数である。SGF項は、式(13)の下記の項
【0062】
【数20】
【0063】
に対応する。例えば、4行の平均をすると仮定すると
【0064】
【数21】
【0065】
であり、ここでは距離の単位をzのセル長にとる。
「真」の線積分サンプル、検出器及びその後の処理によって得られる線積分サンプルの推定値、線積分のz(実際には円錐角)に対する微分、及び線積分サンプルの補正済み推定値を夫々次のベクトル
【0066】
【数22】
【0067】
で表すと、検出の理論の解析から、1次では、次の様になることが判る。
【0068】
【数23】
【0069】
これは、式(10)及び(13)をベクトル形式で書き直したものである。式(18)で、星印”*”は、下記の点毎のベクトル乗数を表す。
【0070】
【数24】
【0071】
指数i(チャンネル方向)に作用し、”x”に沿って適用される広域フィルタをHPで表し、真の線積分データがxに沿って急速に変化しないと仮定すると、
【0072】
【数25】
【0073】
であり、従って下記の「誤差ベクトル」Eが得られる。
【0074】
【数26】
【0075】
適当な境界処理をすると、
【0076】
【数27】
【0077】
上に示した式(21)では、線積分のzー微分がゆっくりと変化し、広域オペレータから引き出すことが出来ると仮定している。データから推定すべき独立パラメータの数を減らすことが望ましい。従って、zー勾配ベクトル[l’]を(p+1)個の基本ベクトル{[B0],…,[Bp]}で”x”に沿って展開すると、
【0078】
【数28】
【0079】
推定しなければならないパラメータの数が、N個のlから、pをこの方法のパラメータとして、(p+1)個の勾配展開係数{Sk},k=0,…pに減少し、結果は、p=1(単純な相関)から9までのpまでの最小自乗を用いて得られる。基本ベクトルの2つの例を下に示す。
1)単項多項式展開
【0080】
【数29】
【0081】
2)チェビシェフ多項式展開
【0082】
【数30】
【0083】
上に記した式で、ベクトルに作用する関数は、ベクトルの各点にこの関数が適用されたときのベクトルを表す短針記法である。式(21)に代入すると、誤差ベクトルに対する下記の表式が得られる。
【0084】
【数31】
【0085】
N行、(p+1)列の行列AAが次の式によって定義される。
【0086】
【数32】
【0087】
従って、その要素は次の式で表される。
【0088】
【数33】
【0089】
上に概略的に説明した方法は、すべての計算された誤差ベクトル成分が、zー勾配の人為効果によるものであることを仮定している(式(21)及び式(23)参照)ことに注意することが重要である。当然ながら、広域フィルタHPの結果は、変化する解剖学的な部分(肋骨の縁のような)による投影の変化の速い成分による寄与をも含むので、近似にすぎない。2つの工程によって、zー勾配の誤差を本来の解剖学的な変化から分離しようとする。第1に、誤差ベクトルをクリップして、zー勾配誤差項目に対する寄与が制限される様にする。稠密な骨の縁に予想される様な一層大きな変化は、zー勾配が原因で起こったと誤って解釈してはならない。第2に、誤差ベクトルを多数の図に互って平均する。実際、勾配の人為効果が(主にリング、円弧、帯並びに種々の「汚れ」の形で)画像の中に表れる為には、3つの条件が必要である。(1)z及びxの利得が一様でないこと、(2)zに解剖学的な部分の勾配があること(又は更に一般的に言えばzに信号の勾配があること)、(3)この勾配が多数の図に互ってゆっくりと変化することである。典型的には、利得の変動は、検出器によって定められ、図から図へと変化しないから、誤差ベクトルを幾つかの図に互って平均すると、(骨の縁等による)過渡的な寄与を平均して除くことが出来る。
【0090】
このとき、勾配のみ(1次)の補正をする装置が、次に述べる工程で得られる。
1)誤差推定値を計算する。
【0091】
【数34】
【0092】
2)行列AAの逆行列を介して勾配展開係数を計算する。
【0093】
【数35】
【0094】
次に勾配及び誤差推定値を計算する。
【0095】
【数36】
【0096】
次にもとの線積分推定値から減算して下記を求める。
【0097】
【数37】
【0098】
ここで波形の記号は、線積分サンプルの補正済み推定値を表す。逆行列AA-1は、特異値分解(SVD)方法を使って(較正で)計算される。SVDは、未知数より多い式を用いて、線形システム問題の最小自乗解を与える。この過程は、セグメント毎に行われ、その後、半走査加重から最初に導き出されたフェザリング関数を使って、セグメントの重なりに互ってこの補正が混合される。
【0099】
上に述べた様な補正をする装置は、実際の「zー勾配データの挙動」に関して多数の仮定に頼っており、まだ相当な数のパラメータを含んでいる。更に、3つの基本関数でzー信号を展開することにより、改善された結果が得られる。これは、3次で次の様に書くことによって行われる。
【0100】
【数38】
【0101】
ここでl”及びl"'は、夫々zに対する信号の2次及び3次微分を夫々表す。CGF及びTFGは夫々曲率利得係数及び3次利得係数であり、(1つのセグメントに対し)下記の式によって定義される。
【0102】
【数39】
【0103】
l’に対してした様に、l”及びl"'をxで展開すると、
【0104】
【数40】
【0105】
となり、ここでckはzー曲率”x”展開係数であり、
【0106】
【数41】
【0107】
となり、tkはzー3次”x”展開係数である。従って、行列AAは(1つのセグメントに対し)次の様に定義し直すことが出来る。
【0108】
【数42】
【0109】
このとき、勾配、曲率及び3次展開係数は、走査データ補正工程のときに、次の式によって与えられる。
【0110】
【数43】
【0111】
上の式で、”x”展開の根拠となる次数pは、勾配、曲率及び3次項に対して同じであると仮定している。しかし、そうである必要はない。式(40)は逆行列AA-1を介して、初期誤差ベクトルEを勾配、曲率及び3次展開係数(si)i=0,...,p,(ci)i=0,...,p,(ti)i=0,...,p,に関係付けている。このとき、最終誤差ベクトル推定値は次の式で表される。
【0112】
【数44】
【0113】
ここでBは基本関数の行列、SGF、CGF及びTGFは、夫々勾配利得係数、曲率利得係数及び3次利得係数である。
この発明は、式(14)及び式(35)−(38)を式(41)に代入すると共に、式(40)を使ってEを因数として取り出すことにより、e及びEを有利に関係付ける。その結果は、行列形式にすると、次の様になる。
【0114】
【数45】
【0115】
ここでVNxNは次の式で表される。
【0116】
【数46】
【0117】
記号は次に定義する。Diag()は対角線行列、α、β及びγは、行列AA-1から、夫々最初の(p+1)行、2番目の(p+1)行及び最後の(p+1)行を抽出することによって得られた行列、fiは、セグメント混合の重み又はフェザリングの重みである。勿論、Vはzー次数が必要とする数の行列の和として表すことが出来る。しかし、補正の時点では、1つの行列だけに初期誤差ベクトルを乗算して、最終誤差推定値を取出すことに注意されたい。
【0118】
しかし、上に述べた様にして得られた行列Vを用いて処理しても、依然として検出器チャンネルの(”x”に於ける)セグメント分割が必要であり、長さNの各々のセグメントがその隣りと(例えばN/3チャンネルだけ)重なり合う。この重なりがフェザリングの重みfによって考慮される。重なり領域にある各チャンネルは、2倍の量の処理を必要とすることが明らかである。更に、データの流れは継ぎ目がないわけではなく、メモリに頻繁にロードすることを必要とするが、これは処理を遅くする。
【0119】
図6について説明すると、この発明は、寸法(Nチャンネル×N)の1つの最終補正行列VF 120を有利に定義することにより、こういう問題を解決する。ここでNチャンネルは、zー勾配誤差に対して補正すべきチャンネルの総数であり、Nは逆行列AA-1の計算に必要なセグメント長である。これは、チャンネルI 126を中心として、奇数の寸法N=2N’+1を持つ初期誤差ベクトルE 124の所定のセグメントS(i) 122を考慮することによって達成される。セグメントS(i) 122が奇数の寸法であることにより、関心が持たれているチャンネルI 126が、セグメントの中心にあることが保証される。計算が上に述べた様に進められて、V(i)N×N 128を定める。V(i) 128の中から、1つの行130、即ち、行N’+1だけを取って、最終行列VF 120の1つの行132を定める。ここでVFNチャンネル×Nは次の式
【0120】
【数47】
【0121】
によって定義される。
次に、セグメントS(i) 122を1チャンネルだけシフトして、S(i+1) 134及び対応するV(i+1) 136を定める。行i+1である、VFの1つの行138が、V(i+1) 136から、行N’+1である1つの行140を抽出することによって定められ、以下同様になる。
【0122】
この過程は、各チャンネルに対し、Nチャンネルを考慮して逆行列を計算するので、従来の較正方法よりもずっと多くの計算を必要とする。しかし、この処理は、zー勾配の較正に必要となる数ヶ月毎に1回だけ行われ、速くする必要がない。これは、最適化Cとは対照的に、PV−WAVEの様な第4世代言語を使うことによって行うことも出来る。しかし、次に述べる利点が得られる。補正処理の重なりがなくなり(フェザリングの重みは1.0に設定する)、こうして補正処理時間のかなりの部分(重なりがN/3であった上に述べた例では1/3)だけ節約され、重要なことは、装置が著しく流線形になることである。補正行列は、1度に1つのチャンネルを更新することが出来る。補正すべきチャンネルの範囲全体を単一コード・ループ内で、割込み又はブランチ飛越しを用いないで処理することが出来る。
【0123】
このとき最終補正の式は次の様になる。
【0124】
【数48】
【0125】
上に述べた装置は大規模に試験を行ったが、従来と比べて、補正効率の低下はないことが判った。2次補正を用いたこの発明の画像の品質は、従来の2次の補正の場合と等しいか又はそれよりもよく、従来の1次の場合よりも更に速く進められる。
画像の品質の観点からすると、この発明の滑りセグメントは、セグメントを選ぶのに最大限の融通性が得られる。逆行列の計算の安定性という様なIQの観点によって決定される所定のセグメントの寸法に対し、各チャンネルから「見た」補正係数(行列VFの1行)は、そのチャンネルがセグメントの中心にあったかの様に計算される。更に、ターゲット・アーキテクチュアに最終ソフトウエア構成を使って、計算の利点を定量化した。再構成ハードウエアでは、この発明の2次補正は、従来の2次補正より6.1倍速く進められ、従来の1次補正よりも3.3倍速く進められる。
【0126】
この発明の1実施例では、典型的には、コンピュータ36(図2)に供給されるデータは、最初に、ビーム硬化の様な周知の種々の誤差を補正する為に、画像再構成装置34によって予備処理される。この発明の補正アルゴリズムは、ビーム硬化補正より後であるが、PCAL補正より前のこの様な予備処理の一部分となる様に構成することが出来る。この代りに、この発明は、画像再構成の前の任意の時点で実施することが出来る。
【0127】
要約すると、この発明は、一様でない患者の解剖学的な構造が、検出セル毎の利得の差の変動と相互作用することによって引き起こされる投影強度のzー軸勾配から生ずる画像の人為効果を補正して、画像の品質並びに計算効率を改善する。前に述べた様に、zー軸勾配利得較正を実施して、個別の検出セルの利得を決定し、希望する補正の次数に基づいて、初期zー軸勾配補正係数を発生する。最終補正のすべての複雑度が、利得の較正に移し替えられ、最終補正の速度及び効率が改善される。更に、この発明の装置は、最終補正行列VFを定める。各々のチャンネルiに対し、初期誤差推定値Eから抽出されたセグメントS(i)に対応する逆行列V(i)を定める。V(i)の内の中心の行だけを残してVFの行iを定める。セグメントS(i)を1チャンネルだけシフトして、S(i+1)を定め、この過程を繰り返す。この為、最終補正行列VFの係数だけが、所望の補正の次数に応じて変化する。従って、この発明は、補正の次数に関係なく、最終補正係数を適用する計算の複雑度を低下させ、最終補正から生ずる画像の品質を更に改善する。
【0128】
上に述べたこの発明は、z方向の複数個の検出セルに用いることが出来る。例えば、上に述べた装置は、双スライスCT装置、即ち、z方向に2列の検出セルを持つ装置を使って構成することが出来る。同様に、この発明は単一スライスCT装置を用いて構成することが出来る。この代りに、この発明は、z方向に2列より多くの列の検出セルを持つ多重スライスCT装置で実施することが出来る。
【0129】
この発明を詳しく説明して図面に示したが、これは例としてであって、例に過ぎず、この発明を制約するものと解してはならないことを明瞭に承知されたい。例えば、ここで説明したCT装置は、x線源及び検出器がガントリと共に回転する「第3世代」装置である。しかし、この発明は、検出器がリング全体の不動の検出器であって、x線源だけがガントリと共に回転する「第4世代」装置を含むこの他の多くのCT装置に用いることが出来る。この発明は、螺旋走査形CT装置と同じく、ステップ・アンド・シュート方式に関連しても利用することが出来る。更に、ここで説明した補正アルゴリズムが、単一スライス及び多重スライス形の両方のCT装置に関連して使うことが出来る。更に、この発明の1形式をビーム硬化の補正後のデータに対して実施されるものとして説明したが、この発明は、データ補正/処理の中の種々の点で実施することが出来る。従って、この発明の変形並びに変更は当業者に明らかであろうし、特許請求の範囲は、これらのすべての変更及び均等物をカバーするものであることを承知されたい。
【図面の簡単な説明】
【図1】CT作像装置の見取図。
【図2】図1に示した装置の簡略ブロック図。
【図3】ビームの焦点スポット、コリメータ及び検出器の略図。
【図4】検出器及び関連した制御装置の1列の検出セルのブロック図形式の見取図。
【図5】種々の厚さの画像スライスに対する検出セル・データの組合せを例示する図。
【図6】この発明の最終補正行列の定義を示す略図。
【符号の説明】
14 x線源
20 検出セル
34 画像再構成装置
Claims (10)
- 投影データから物体の断層写真画像を作る装置に於て、
x線を放出するx線源と、
複数個の検出セルを持ち、前記x線を受け取って、該x線の減衰を表す投影データを出力として発生する検出器と、
前記投影データの間の非一様性から生ずる画像アーティファクトを補正する為に、異なる検出セルの利得によって導入される誤差を補正するためのベクトル式の初期誤差係数であって、前記投影データから抽出した前記初期誤差係数に補正係数を適用する画像プロセッサとを有し、
該補正係数が、基本関数行列B、勾配利得係数SGF、曲率利得係数CGF及び3次利得係数TGFの関数であり、
1つの行列に初期誤差係数を乗算することにより前記補正係数が導出される、装置。 - 前記投影データがベクトル式であって、前記補正係数が、可変係数を持つ行列乗数である請求項1記載の装置。
- 前記検出器が複数個のチャンネルを持ち、1列の検出セルが、前記複数個のチャンネルの各々を形成しており、前記補正係数が複数個のチャンネル補正係数で構成され、前記複数個のチャンネルの各々が、各々の前記チャンネル補正係数を決定する為に夫々用いられる対応するセグメントの中心チャンネルである請求項1記載の装置。
- 前記補正係数が画像を構成する間に適用される請求項1記載の装置。
- x線ビームを複数個の検出セルを持つ検出器に向かって放出するx線源を持ち、該検出セルからの出力信号から投影データが発生される様な計算機式断層写真装置で異なる検出セルの利得によって導入される誤差に対して投影データを補正する方法に於て、非一様な投影データから生ずる画像アーティファクトを補正する為に、投影データから抽出されたベクトル式の初期誤差係数に補正係数を適用する工程を含み、前記補正係数が、基本関数行列B、勾配利得係数SGF、曲率利得係数CGF及び3次利得係数TGFの関数であり、
1つの行列に初期誤差係数を乗算することにより前記補正係数が導出される、方法。 - 前記投影データがベクトル式であり、前記補正係数が可変係数を持つ行列乗数を持つ請求項6記載の方法。
- 更に、補正の対象とする複数個のチャンネルを決定し、1列の検出セルが該複数個のチャンネルの各々を形成しており、前記補正係数は複数個のチャンネル補正係数で構成されていて、前記複数個のチャンネルの各々が、各々のチャンネル補正係数を夫々決定する為に使われる対応するセグメントの中心チャンネルである請求項6記載の方法。
- 前記補正係数が画像を構成する間に適用される請求項6記載の方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/139,992 US6173029B1 (en) | 1998-08-25 | 1998-08-25 | Higher order detector z-slope correction for a multislice computed tomography system |
US09/139992 | 1998-08-25 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000083945A JP2000083945A (ja) | 2000-03-28 |
JP4409671B2 true JP4409671B2 (ja) | 2010-02-03 |
Family
ID=22489233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP23487499A Expired - Lifetime JP4409671B2 (ja) | 1998-08-25 | 1999-08-23 | 断層写真画像を作る装置と方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6173029B1 (ja) |
EP (1) | EP0982682A3 (ja) |
JP (1) | JP4409671B2 (ja) |
IL (1) | IL131375A (ja) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6418185B1 (en) | 1999-08-18 | 2002-07-09 | General Electric Company | Methods and apparatus for time-multiplexing data acquisition |
US7031993B1 (en) * | 2000-02-18 | 2006-04-18 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for fast natural log(X) calculation |
CA2305249A1 (en) * | 2000-04-14 | 2001-10-14 | Branko Sarcanin | Virtual safe |
JP3852746B2 (ja) * | 2001-03-08 | 2006-12-06 | インターナショナル・ビジネス・マシーンズ・コーポレーション | データ補正装置及びデータ補正方法 |
EP1244058A1 (de) * | 2001-03-22 | 2002-09-25 | VAMP Verfahren und Apparate der Medizinischen Physik GmbH | Computertomograph mit merkmalsgewichteter Bild-zu-Volumen Interpolation in Echtzeit |
US6956977B2 (en) * | 2001-08-08 | 2005-10-18 | General Electric Company | Methods for improving contrast based dynamic range management |
US6885764B2 (en) * | 2001-11-21 | 2005-04-26 | Ge Medical Systems Global Technology Company, Llc | High Speed Z-smoothing method and apparatus for CT imaging system |
US7379524B2 (en) * | 2002-09-13 | 2008-05-27 | Xoran Technologies, Inc. | Computer tomography scanner |
JP4062232B2 (ja) * | 2003-10-20 | 2008-03-19 | 株式会社日立製作所 | X線ct装置及びx線ct装置による撮像方法 |
US7235790B2 (en) * | 2004-02-17 | 2007-06-26 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for radiation detection |
JP2006034785A (ja) * | 2004-07-29 | 2006-02-09 | Ge Medical Systems Global Technology Co Llc | X線ct画像処理方法およびx線ct装置 |
US7734079B2 (en) * | 2004-09-28 | 2010-06-08 | General Electric Company | Methods and apparatus for image reconstruction |
US7379527B2 (en) * | 2005-12-22 | 2008-05-27 | General Electric Company | Methods and apparatus for CT calibration |
FR3015035B1 (fr) * | 2013-12-13 | 2016-01-22 | Safran | Mesure non intrusive de la densite volumique d'une phase dans une piece |
US10373092B2 (en) * | 2015-09-08 | 2019-08-06 | Oracle International Corporation | System and method for providing a tradeoff between time and quality for a task |
CN105678711B (zh) * | 2016-01-29 | 2018-08-21 | 中国科学院高能物理研究所 | 一种基于图像分割的衰减校正方法 |
US10741247B1 (en) * | 2019-06-21 | 2020-08-11 | Macronix International Co., Ltd. | 3D memory array device and method for multiply-accumulate |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4284950A (en) * | 1978-08-05 | 1981-08-18 | E M I Limited | Imaging systems |
US5301108A (en) | 1993-02-01 | 1994-04-05 | General Electric Company | Computed tomography system with z-axis correction |
US5579359A (en) | 1995-12-21 | 1996-11-26 | General Electric Company | Methods and apparatus for calibrating detector cell output signals |
US5610963A (en) * | 1996-02-06 | 1997-03-11 | General Electric Company | Methods and systems for determining the z-axis profile of a detector in a CT system |
-
1998
- 1998-08-25 US US09/139,992 patent/US6173029B1/en not_active Expired - Lifetime
-
1999
- 1999-08-12 IL IL13137599A patent/IL131375A/en not_active IP Right Cessation
- 1999-08-20 EP EP99306632A patent/EP0982682A3/en not_active Withdrawn
- 1999-08-23 JP JP23487499A patent/JP4409671B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
EP0982682A3 (en) | 2000-05-17 |
IL131375A (en) | 2002-11-10 |
IL131375A0 (en) | 2001-01-28 |
US6173029B1 (en) | 2001-01-09 |
JP2000083945A (ja) | 2000-03-28 |
EP0982682A2 (en) | 2000-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4409671B2 (ja) | 断層写真画像を作る装置と方法 | |
US9245320B2 (en) | Method and system for correcting artifacts in image reconstruction | |
JP3708434B2 (ja) | 反復円錐ビームct画像再構成 | |
US7391844B2 (en) | Method and apparatus for correcting for beam hardening in CT images | |
US7379575B2 (en) | Method for post- reconstructive correction of images of a computer tomograph | |
US5165100A (en) | Over-range image artifact reduction in tomographic imaging | |
US7747057B2 (en) | Methods and apparatus for BIS correction | |
US20070217568A1 (en) | Method and apparatus for correcting motion in image reconstruction | |
JP2002320613A (ja) | 計算機式断層写真法イメージング・システムのための逆投影方法及び装置 | |
JP2002345808A (ja) | Ctスカウト画像処理のための方法及び装置 | |
US6944258B2 (en) | Beam hardening post-processing method and X-ray CT apparatus | |
JP2004188163A (ja) | 断層撮影装置 | |
US20060193429A1 (en) | Calculation of additional projection data from projection data acquired with a divergent beam | |
JPH0838467A (ja) | 検出器チャンネル利得較正係数を求める方法 | |
JP2013085955A (ja) | 連続マルチスケール再構成において詳細画像を補うx線コンピュータ断層撮像装置(x線ct装置)、医用画像処理装置及び医用画像処理方法 | |
JP6446361B2 (ja) | X線ct装置および補正処理装置 | |
CN103054599B (zh) | X射线ct装置及x射线ct装置的动作方法 | |
JP2002541896A (ja) | 半撮影域のみをカバーする縮小サイズ検出器を利用するコンピュータ断層撮影システムに用いる装置及び方法 | |
JPH08252249A (ja) | 像データから検出器の劣化によるアーティファクトを除去する方法及び物体の断層写真像を作成するシステム | |
JPH10295685A (ja) | 投影データを補正する方法及び物体の断層画像を形成するシステム | |
US7379525B2 (en) | Method and system for efficient helical cone-beam reconstruction | |
JP4582997B2 (ja) | 高速コンピュータ断層撮影方法 | |
US9508164B2 (en) | Fast iterative image reconstruction method for 3D computed tomography | |
JP2000083948A (ja) | 断層写真画像を作る装置と方法 | |
JP4387758B2 (ja) | Spect装置及びspect画像再構成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060818 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090331 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090603 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20090603 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20090603 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090623 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090629 |
|
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: 20091020 |
|
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: 20091112 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4409671 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121120 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131120 Year of fee payment: 4 |
|
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 |
|
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 |
|
EXPY | Cancellation because of completion of term |