JP3917684B2 - 物体の断層写真像を作成する方法及び装置 - Google Patents

物体の断層写真像を作成する方法及び装置 Download PDF

Info

Publication number
JP3917684B2
JP3917684B2 JP11404396A JP11404396A JP3917684B2 JP 3917684 B2 JP3917684 B2 JP 3917684B2 JP 11404396 A JP11404396 A JP 11404396A JP 11404396 A JP11404396 A JP 11404396A JP 3917684 B2 JP3917684 B2 JP 3917684B2
Authority
JP
Japan
Prior art keywords
projection data
data array
data
image
generating
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
JP11404396A
Other languages
English (en)
Other versions
JPH0998968A (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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of JPH0998968A publication Critical patent/JPH0998968A/ja
Application granted granted Critical
Publication of JP3917684B2 publication Critical patent/JP3917684B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • 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/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
    • 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

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

Description

【０００１】
【産業上の利用分野】
本発明は、一般的には計算機断層写真（ＣＴ）作像に関し、更に具体的に言えば、４つの扇形ビームを用いた螺旋走査から収集されたデータからの像の再生に関する。
【０００２】
【従来の技術】
ＣＴシステムでは、Ｘ線源が扇形ビームを投射し、この扇形ビームは、「作像平面」と呼ばれる直角座標系のＸ−Ｙ平面内に来るようにコリメートされる。Ｘ線ビームは患者のような作像されるべき物体を通過し、放射線検出器の線形配列に入射する。透過した放射線の強度は、物体によるＸ線ビームの減衰に関係する。線形配列の各々の検出器は、ビームの減衰量の測定値である別々の電気信号を発生する。すべての検出器からの減衰測定値を別々に収集して、透過輪郭を形成する。
【０００３】
ＣＴシステムのＸ線源及び線形検出器配列は、ガントリと共に作像平面内で物体の周りを回転し、Ｘ線ビームが物体と交わる角度が絶えず変化するようにする。１つのガントリ角度における検出器配列からの一群のＸ線減衰測定値が、「ビュー（図）」と呼ばれる。物体の「走査（スキャン）」は、Ｘ線源及び検出器が１回転する間の相異なるガントリ角度で求められた１組のビューで構成される。軸方向走査では、物体を通るように切られた２次元スライスに対応する像を構成するようにデータが処理される。１組のデータから像を再生する１つの方法が、この分野ではフィルタ付き逆投影（バック・プロジェクション）方式と呼ばれている。この方法は、走査からの減衰測定値を「ＣＴ数」又は「ハウンズフィールド単位」と呼ばれる整数に変換し、このような整数を用いて、陰極線管表示装置における対応する画素の輝度を制御する。
【０００４】
多数のスライスに必要な合計走査時間を短縮するために、「螺旋（ヘリカル）」走査を実施することができる。「螺旋」走査を実施するためには、所定数のスライスに対するデータを収集する間に、患者を移動させる。４ビーム螺旋走査を用いて、合計走査時間を更に短縮することができる。このような走査は、４列の検出器が相異なる軸方向の位置における投影測定値を同時に収集するＣＴシステムを利用して行われる。４つの検出器列は４つの扇形ビームを画定する。このようなシステムは、普通の扇形ビーム１つの螺旋走査からの単一の螺旋とは対照的に、互いに織り込まれた四重の螺旋を発生する。４つの扇形ビームによって写像（マップ）された互いに織り込まれた螺旋が、投影データを発生し、この投影データから各々の所定のスライス内の像を再生することができる。
【０００５】
勿論、４ビーム螺旋走査で得られたデータから、アーティファクトを最低レベルにして品質の高い像が得られるような形で像を再生することが望ましい。このような像を再生するのに要する合計時間を短縮することも望ましい。
【０００６】
【発明の要約】
本発明の一形式では、４ビーム螺旋走査によって得られたデータを用いて、像の再生を実行する装置を提供する。更に具体的に言うと、本発明の一形式によれば、各々の扇形ビームによって収集された投影データから、投影空間データ配列が選択される。次に、各々の配列内のデータに加重して、患者の並進運動を補正し、データ冗長効果を打ち消す。その後、加重データを用いて像を再生する。
【０００７】
更に具体的に言うと、像を再生するときに、投影データ配列を発生する。投影データ配列は、作像されるべきスライスに関連したデータ平面に対応する。次に、データ配列に加重係数を適用して、各々の特定のデータ要素に重みを割り当てる。この後、加重投影データ配列をフィルタ作用にかけて逆投影し、像データ配列を発生する。その後、像データ配列を加算して、スライス像データ配列を発生する。
【０００８】
患者の４ビーム螺旋走査を用いると、合計の患者走査時間が短くなるという利点がある。更に、上に述べた像再生アルゴリズムにより、患者テーブル並進速度が高くなっても、品質の高い像スライスを発生することができるという利点が得られる。
【０００９】
【実施例】
図１及び図２について説明すると、計算機断層写真（ＣＴ）作像システム１０が、「第三世代」ＣＴスキャナを代表するガントリ１２を含んでいる。ガントリ１２はＸ線源１４を有しており、Ｘ線源１４は、ガントリ１２の反対側にある検出器配列１８に向かってＸ線ビーム１６を投射する。検出器配列１８は検出器素子２０によって形成されており、これらの検出器素子は、医療の患者２２を通過する投射されたＸ線を一緒に感知する。検出器素子２０の各々は、入射するＸ線ビームの強度、従って患者２２を通過するときのビームの減衰を表す電気信号を発生する。Ｘ線投影データを収集する走査の間に、ガントリ１２及びその上に取り付けられた部品が回転中心２４の周りを回転する。
【００１０】
ガントリ１２の回転及びＸ線源１４の動作は、ＣＴシステム１０の制御機構２６によって制御される。制御機構２６は、Ｘ線源１４に対する電力及びタイミング信号を供給するＸ線制御装置２８と、ガントリ１２の回転速度及び位置を制御するガントリ・モータ制御装置３０とを含んでいる。制御機構２６にあるデータ収集システム（ＤＡＳ）３２が、検出器素子２０からのアナログ・データをサンプリングし、この後の処理のために、このデータをディジタル信号に変換する。像再生装置３４が、ＤＡＳ３２からサンプリングされてディジタル化されたＸ線データを受け取り、高速の像再生を行う。再生された像はコンピュータ３６に対する入力として印加され、このコンピュータは、大容量記憶装置３８に像を記憶する。
【００１１】
コンピュータ３６は又、キーボートを有しているコンソール４０を介して、オペレータからの指令及び走査パラメータを受け取る。関連した陰極線管表示装置４２が、オペレータが再生された像及びコンピュータ３６からのその他のデータを観察することができるようにする。オペレータから供給された指令及びパラメータはコンピュータ３６で用いられて、ＤＡＳ３２、Ｘ線制御装置２８及びガントリ・モータ制御装置３０に対する制御信号及び情報を供給する。更に、コンピュータ３６は、ガントリ１２内で患者２２を位置決めするためにモータ式テーブル４６を制御するテーブル・モータ制御装置４４を動作させる。特に、テーブル４６は患者２２の一部をガントリ開口４８内に通す。
【００１２】
図３に示すように、４扇形ビーム・システムでは４列の検出器が画定されている。Ｘ線扇形ビームは、事実上、回転のＺ軸に沿って変位した４つの扇形ビームに分割される。隣接したビームの中心の間の距離は、ガントリの回転軸線の所で測定したときに、Ｄである。
像再生装置３４は図４に更に詳しく示されている。具体的に言うと、各々の扇形ビームに対するＤＡＳ３２からのデータの各々のビューをそれぞれの前処理装置５２Ａ〜５２Ｄで受け取り、そこで、それぞれのビームは、ビーム硬化、検出器及びチャンネル利得のオフセット及び変動のような周知の種々の誤差を補正するために前処理される。更に、負の対数を求めて投影データとし、この投影データは投影データ配列５４Ａ〜５４Ｄに記憶される。
【００１３】
各々のビーム配列５４Ａ〜５４Ｄの投影データを読み出し、対応する加重関数５６Ａ〜５６Ｄが掛算器５８Ａ〜５８Ｄで適用される。加重投影データは、加重投影データ配列６０Ａ〜６０Ｄ内の対応する位置に書き込まれ、この加重投影データをフィルタ作用及び逆投影にかけて（６２Ａ〜６２Ｄ）、ビーム像データ配列６４Ａ〜６４Ｄを発生する。
【００１４】
この後、像データ配列６４Ａ〜６４Ｄを参照番号６６の所で加算して、スライス像データ配列６８を発生する。具体的に言うと、ビーム１の配列にある各々の画素の大きさを、ビーム２の配列、ビーム３の配列及びビーム４の配列にある対応する画素の大きさと加算する。その結果得られたスライス像配列６８は、後で用いるために記憶され、又はオペレータに対して表示され得る。像データ配列６４Ａ〜６４Ｄを発生した後にデータを加算する代わりに、同じガントリ（ビュー）角度であるが、異なる検出器列からの投影を、データのフィルタ作用及び逆投影の前に組み合わせることができる。このような組み合わせは、処理の負担を軽減させる。
【００１５】
本発明は一形式では、ある所定の条件の下で４ビーム走査が実行されたときに、加重投影データ配列６０Ａ〜６０Ｄを発生することに特に関連する。以下の説明では、ｄはガントリの回転軸線の所で測定された検出器列（Ｚ方向）の間隔を表し、ｓは（１回転当たりの）テーブル送り速度を表し、ｐはｄとｓとの比を表す。即ち、
ｐ＝ｄ／ｓ （１）
図５に示すように、データ平面Ｐ1 、Ｐ2 、Ｐ3 及びＰ4 が、再生しようとするスライスＰと線Ｌ1 Ｍ1 、Ｌ2 Ｍ2 、Ｌ3 Ｍ3 及びＬ4 Ｍ4 で交わっている。これらの線の関数は、次のように表すことができる。
【００１６】
Ｌ1 Ｍ1 ： β１＝−３ｐπ
Ｌ2 Ｍ2 ： β２＝−ｐπ
Ｌ3 Ｍ3 ： β３＝ｐπ
Ｌ4 Ｍ4 ： β４＝３ｐπ （２）
ここで、βはガントリ角度に等しい。線Ｌ1 Ｍ1 、Ｌ2 Ｍ2 、Ｌ3 Ｍ3 及びＬ4 Ｍ4 は、下に示すように＋及び−の組で表す「鏡像（ミラー）」線を有している。
【００１７】
β１± ＝−３ｐπ±π−２γ
β２± ＝−ｐπ±π−２γ
β３± ＝ｐπ±π−２γ
β４± ＝３ｐπ±π−２γ （３）
ここで、γは検出器角度に等しい。
【００１８】
テーブル送り速度ｓ及び検出器のＺ軸間隔ｄが、式（４π／（π−２γm ））ｄ＜ｓ＜（６π／（π−２γm ））ｄの関係を満たすときに（ここで、γm は扇形角度の半分と定義する。）、Ｗ１（β，γ）、Ｗ２（β，γ）、Ｗ３（β，γ）及びＷ４（β，γ）で表す各々のデータの組に適用されるべき螺旋加重係数は、次のようになる。
【００１９】
β≦β３-のとき、 Ｗ１（β，γ）＝０
β３-＜β≦β１のとき、Ｗ１（β，γ）＝（１／２）・（β−β３-）／（β１−β３-）
β１＜β＜β２のとき、Ｗ１（β，γ）＝（１／２）・（β−β２）／（β１−β２）
β≧β２のとき、 Ｗ１（β，γ）＝０ （４）
β≦β３-のとき、 Ｗ２（β，γ）＝０
β３-＜β≦β２のとき、Ｗ２（β，γ）＝ｗ２1（β，γ）＋ｗ２2（β，γ）
β２＜β＜β３のとき、 Ｗ２（β，γ）＝（β−β３）／（β２−β３）
β≧β３のとき、 Ｗ２（β，γ）＝０ （５）
ここで、
β≦β１のとき、 ｗ２1（β，γ）＝０
β１＜β≦β２のとき、 ｗ２1（β，γ）＝（１／２）・（β−β１）／（β２−β１）
β≧β２のとき、 ｗ２1（β，γ）＝０
β≦β４-のとき、 ｗ２2（β，γ）＝０
β４-＜β≦β２のとき、ｗ２2（β，γ）＝（１／２）・（β−β４-）／（β２−β４-）
β＞β２のとき、 ｗ２2（β，γ）＝０
β≦β２のとき、 Ｗ３（β，γ）＝０
β２＜β＜β３のとき、 Ｗ３（β，γ）＝（β−β３）／（β２−β３）
β３≦β＜β２+のとき、Ｗ３（β，γ）＝ｗ３1（β，γ）＋ｗ３2（β，γ）
β≧β２+のとき、 Ｗ３（β，γ）＝０ （６）
ここで、
β＜β３のとき、 ｗ３1（β，γ）＝０
β３≦β＜β１+のとき、ｗ３1（β，γ）＝（１／２）・（β−β１+）／（β３−β１+）
β≧β１+ のとき、 ｗ３1（β，γ）＝０
β＜β３のとき、 ｗ３2（β，γ）＝０
β３≦β＜β４のとき、 ｗ３2（β，γ）＝（１／２）・（β−β４）／（β３−β４）
β≧β４のとき、 ｗ３2（β，γ）＝０
β≦β３のとき、 Ｗ４（β，γ）＝０
β３＜β≦β４のとき、 Ｗ４（β，γ）＝（１／２）・（β−β３）／（β４−β３）
β４＜β＜β２+のとき、Ｗ４（β，γ）＝（１／２）・（β−β２+）／（β４−β２+）
β≧β２+のとき、 Ｗ４（β，γ）＝０ （７）
加重関数によって定められたデータ平面内の領域が図６及び図７に示されている。図７の陰影線を施した区域は、像スライスを再生するために用いられる各々のデータ平面からのデータを表す。扇形角度２γm ＝π／４のときに、この像再生アルゴリズムの作用し得る範囲は、１６／３ｄ＜ｓ＜８ｄである。
【００２０】
式（４）から式（７）に示した加重関数は連続している。しかしながら、この各々の式の１次微分は境界線で不連続である。必要であれば、境界にわたって数チャンネル（〜１０チャンネル）にフェザリングを施すことによって、この不連続性を除去することができる。
一旦加重して加重投影データ配列６０Ａ〜６０Ｄを作成したら、処理時間を短縮するために、フィルタ作用及び逆投影の前に、同じガントリ（ビュー）角度であるが、異なる検出器列からの投影を組み合わせることができる。データ列１にある投影ビューは、データ列４にある対応する投影ビューから３６０°離れている。これらのビューの対をフィルタ作用の前に更に組み合わせて、処理負担の不必要な増加を除くことができる。
【００２１】
前に述べたアルゴリズムについて述べると、概略的に、１つのスライスを再生するために３．６７π相当の投影データを前処理する必要がある。従来の走査（軸方向及び螺旋の両方）で像を再生するのに必要なデータ量に比べると、４ビーム螺旋走査では、前処理の負担が１．８倍に増加する。
隣接するスライスを再生するのに必要な投影データは、ビューの角度βの原点を垂直にずらして、再生しようとする新しいスライスと整合させることにより確認することができる。大抵の場合、各々のデータの組にある１つのスライスに対するデータと、隣接するスライスに対するデータとの間にはかなりの重なりがある。加重する前には、前処理はスライスの位置には関係しない。このため、（螺旋加重をしない）前処理の結果を、将来再び用いるためにバッファに記憶することができる。これにより、大抵の場合には、前処理の負担が大幅に軽減される。遡及形の再生では特に、前処理の負担が大幅に軽減される。
【００２２】
所望のスライス輪郭が、上に述べたデータ及び再生アルゴリズムによって形成される輪郭よりも一層厚手である場合には、所望のスライス輪郭内にある多数の薄いスライスを加算することにより、一層厚手のスライスを導き出すことができる。多数の薄いスライス自体に関心がない場合には、多数の薄いスライスを再生する中間工程を側路して、投影領域の早い段階で、対応する加算を実行することができる。こうすると、計算の負担及び像を記憶する負担が減少する。このときの加重関数は、データ平面を対応してずらしたものを加算することにより導き出すことができる。
【００２３】
本発明の種々の実施例についてこれまで説明したところから、本発明の目的が達成されたことは明らかである。本発明を詳細に説明すると共に図面に例示したが、これは例示に過ぎず、本発明を制約するものと解してはならないことは言うまでもない。例えば、ここで説明したＣＴシステムは、Ｘ線源及び検出器の両方がガントリと共に回転する「第三世代」システムである。検出器が完全リングの不動の検出器であって、Ｘ線源のみがガントリと共に回転する「第四世代」システムを含めて、他の多くのＣＴシステムを用いることができる。従って、本発明の要旨は、特許請求の範囲のみによって限定されることを承知されたい。
【図面の簡単な説明】
【図１】ＣＴ作像システムの見取図である。
【図２】図１に示すシステムの概略ブロック図である。
【図３】Ｚ軸に沿った４ビームＸ線の概略断面図である。
【図４】図２のＣＴ作像システムの一部を形成している像再生装置のブロック図である。
【図５】像スライスと４つのデータ平面との交差を示す図である。
【図６】データ平面及び種々のデータ領域を示す図である。
【図７】像スライスを構成するために利用される選択されたデータ領域及び４つのデータ平面を示す図である。
【符号の説明】
１０ 計算機断層写真（ＣＴ）作像システム
１２ ガントリ
１４ Ｘ線源
１６ Ｘ線ビーム
１８ 検出器配列
２０ 検出器素子
２４ 回転中心
２６ 制御機構
２８ Ｘ線制御装置
３０ ガントリ・モータ制御装置３０
３２ データ収集システム（ＤＡＳ）
３４ 像再生装置
３６ コンピュータ
３８ 大容量記憶装置
４０ コンソール
４２ 陰極線管表示装置
４６ モータ式テーブル
４８ ガントリ開口
５４Ａ〜５４Ｄ 投影データ配列
５６Ａ〜５６Ｄ 加重関数
５８Ａ〜５８Ｄ 掛算器
６０Ａ〜６０Ｄ 加重投影データ
６２Ａ〜６２Ｄ フィルタ作用及び逆投影
６４Ａ〜６４Ｄ 像データ配列
６８ スライス像データ配列

Claims (9)

  1. ４つのＸ線扇形ビームを用いた螺旋走査で収集された投影データから物体の断層写真像を作成する方法であって、
    （ａ） 各々のＸ線扇形ビームから得られたデータに対応する投影データ配列を発生する工程と、
    （ｂ） 工程（ａ）で発生された前記投影データ配列の各々に加重関数を適用して、それぞれの加重投影データ配列を発生する工程であって、各々のデータの組に適応される加重係数は、Ｗ１（β，γ）、Ｗ２（β，γ）、Ｗ３（β，γ）及びＷ４（β，γ）として表すと、
    テーブル送り速度ｓ及び検出器のＺ軸間隔ｄが、式（４π／（π−２γ m ））ｄ＜ｓ＜（６π／（π−２γ m ））ｄの関係を満たすとき（γ m は扇形角度の半分）であって、
    β≦β３-のとき、 Ｗ１（β，γ）＝０
    β３-＜β≦β１のとき、Ｗ１（β，γ）＝（１／２）・（β−β３-）／（β１−β３-）
    β１＜β＜β２のとき、Ｗ１（β，γ）＝（１／２）・（β−β２）／（β１−β２）
    β≧β２のとき、 Ｗ１（β，γ）＝０
    β≦β３-のとき、 Ｗ２（β，γ）＝０
    β３-＜β≦β２のとき、Ｗ２（β，γ）＝ｗ２1（β，γ）＋ｗ２2（β，γ）
    β２＜β＜β３のとき、 Ｗ２（β，γ）＝（β−β３）／（β２−β３）
    β≧β３のとき、 Ｗ２（β，γ）＝０
    ここで、
    β≦β１のとき、 ｗ２1（β，γ）＝０
    β１＜β≦β２のとき、 ｗ２1（β，γ）＝（１／２）・（β−β１）／（β２−β１）
    β≧β２のとき、 ｗ２1（β，γ）＝０
    β≦β４-のとき、 ｗ２2（β，γ）＝０
    β４-＜β≦β２のとき、ｗ２2（β，γ）＝（１／２）・（β−β４-）／（β２−β４-）
    β＞β２のとき、 ｗ２2（β，γ）＝０
    β≦β２のとき、 Ｗ３（β，γ）＝０
    β２＜β＜β３のとき、 Ｗ３（β，γ）＝（β−β３）／（β２−β３）
    β３≦β＜β２+のとき、Ｗ３（β，γ）＝ｗ３1（β，γ）＋ｗ３2（β，γ）
    β≧β２+のとき、 Ｗ３（β，γ）＝０
    ここで、
    β＜β３のとき、 ｗ３1（β，γ）＝０
    β３≦β＜β１+のとき、ｗ３1（β，γ）＝（１／２）・（β−β１+）／（β３−β１+）
    β≧β１+ のとき、 ｗ３1（β，γ）＝０
    β＜β３のとき、 ｗ３2（β，γ）＝０
    β３≦β＜β４のとき、 ｗ３2（β，γ）＝（１／２）・（β−β４）／（β３−β４）
    β≧β４のとき、 ｗ３2（β，γ）＝０
    β≦β３のとき、 Ｗ４（β，γ）＝０
    β３＜β≦β４のとき、 Ｗ４（β，γ）＝（１／２）・（β−β３）／（β４−β３）
    β４＜β＜β２+のとき、Ｗ４（β，γ）＝（１／２）・（β−β２+）／（β４−β２+）
    β≧β２+のとき、 Ｗ４（β，γ）＝０
    ここで、
    β１± ＝−３ｐπ±π−２γ
    β２± ＝−ｐπ±π−２γ
    β３± ＝ｐπ±π−２γ
    β４± ＝３ｐπ±π−２γ であり、
    βはガントリ角度であり、γは検出器角度である、
    である、加重投影データ配列を発生する工程と、
    （ｃ） 工程（ｂ）で発生された投影データ配列から、スライス像を構成するために用いられる像データ配列を発生する工程とを備えた物体の断層写真像を作成する方法。
  2. 前記像データ配列を発生する工程は、各々の加重投影データ配列にフィルタ作用及び逆投影を行う工程を含んでいる請求項１に記載の方法。
  3. フィルタ作用及び逆投影を行う前に、同じガントリ角度であるが、相異なる検出器列からのデータ配列が組み合わされる請求項２に記載の方法。
  4. 第１のデータ列にある投影ビューが第４のデータ列にある投影ビューから３６０°である場合、フィルタ作用及び逆投影を行う前に、これらのビューを組み合わせる請求項３に記載の方法。
  5. 前記投影データ配列の各々に加重関数を適用する前に、該データを相次ぐスライスを再生するためのシステム・メモリに記憶する請求項２に記載の方法。
  6. 所望のスライス輪郭がデータ配列によって形成される輪郭よりも一層厚手である場合に、所望のスライス輪郭内にある多数の薄いスライスが加算される請求項２に記載の方法。
  7. 螺旋走査で収集された投影データから物体の断層写真像を作成する装置であって、各々のＸ線扇形ビームから得られたデータに対応する投影データ配列を発生する手段と、該投影データ配列を発生する手段により発生された前記投影データ配列の各々に加重関数を適用して、それぞれの加重投影データ配列を発生する手段であって、各々のデータの組に適用される加重関数は、Ｗ１（β，γ）、Ｗ２（β，γ）、Ｗ３（β，γ）及びＷ４（β，γ）として表すと、
    テーブル送り速度ｓ及び検出器のＺ軸間隔ｄが、式（４π／（π−２γ m ））ｄ＜ｓ＜（６π／（π−２γ m ））ｄの関係を満たすとき（γ m は扇形角度の半分）であって、
    β≦β３-のとき、 Ｗ１（β，γ）＝０
    β３-＜β≦β１のとき、Ｗ１（β，γ）＝（１／２）・（β−β３-）／（β１−β３-）
    β１＜β＜β２のとき、Ｗ１（β，γ）＝（１／２）・（β−β２）／（β１−β２）
    β≧β２のとき、 Ｗ１（β，γ）＝０
    β≦β３-のとき、 Ｗ２（β，γ）＝０
    β３-＜β≦β２のとき、Ｗ２（β，γ）＝ｗ２1（β，γ）＋ｗ２2（β，γ）
    β２＜β＜β３のとき、 Ｗ２（β，γ）＝（β−β３）／（β２−β３）
    β≧β３のとき、 Ｗ２（β，γ）＝０
    ここで、
    β≦β１のとき、 ｗ２1（β，γ）＝０
    β１＜β≦β２のとき、 ｗ２1（β，γ）＝（１／２）・（β−β１）／（β２−β１）
    β≧β２のとき、 ｗ２1（β，γ）＝０
    β≦β４-のとき、 ｗ２2（β，γ）＝０
    β４-＜β≦β２のとき、ｗ２2（β，γ）＝（１／２）・（β−β４-）／（β２−β４-）
    β＞β２のとき、 ｗ２2（β，γ）＝０
    β≦β２のとき、 Ｗ３（β，γ）＝０
    β２＜β＜β３のとき、 Ｗ３（β，γ）＝（β−β３）／（β２−β３）
    β３≦β＜β２+のとき、Ｗ３（β，γ）＝ｗ３1（β，γ）＋ｗ３2（β，γ）
    β≧β２+のとき、 Ｗ３（β，γ）＝０
    ここで、
    β＜β３のとき、 ｗ３1（β，γ）＝０
    β３≦β＜β１+のとき、ｗ３1（β，γ）＝（１／２）・（β−β１+）／（β３−β１+）
    β≧β１+ のとき、 ｗ３1（β，γ）＝０
    β＜β３のとき、 ｗ３2（β，γ）＝０
    β３≦β＜β４のとき、 ｗ３2（β，γ）＝（１／２）・（β−β４）／（β３−β４）
    β≧β４のとき、 ｗ３2（β，γ）＝０
    β≦β３のとき、 Ｗ４（β，γ）＝０
    β３＜β≦β４のとき、 Ｗ４（β，γ）＝（１／２）・（β−β３）／（β４−β３）
    β４＜β＜β２+のとき、Ｗ４（β，γ）＝（１／２）・（β−β２+）／（β４−β２+）
    β≧β２+のとき、 Ｗ４（β，γ）＝０
    ここで、
    β１± ＝−３ｐπ±π−２γ
    β２± ＝−ｐπ±π−２γ
    β３± ＝ｐπ±π−２γ
    β４± ＝３ｐπ±π−２γ であり、
    βはガントリ角度であり、γは検出器角度である、
    加重投影データ配列を発生する手段と、該投影データ配列からスライス像を構成する手段とを備えた物体の断層写真像を作成する装置。
  8. 前記像データ配列を発生する手段は、各々の加重投影データ配列にフィルタ作用及び逆投影を行う手段を含んでいる請求項７に記載の装置。
  9. 所望のスライス輪郭内にある多数の薄いスライスを加算する手段を更に含んでいる請求項７に記載の装置。
JP11404396A 1995-05-09 1996-05-09 物体の断層写真像を作成する方法及び装置 Expired - Fee Related JP3917684B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/436,176 US5541970A (en) 1995-05-09 1995-05-09 Image reconstruction for a CT system implementing a four fan beam helical scan
US08/436176 1995-05-09

Publications (2)

Publication Number Publication Date
JPH0998968A JPH0998968A (ja) 1997-04-15
JP3917684B2 true JP3917684B2 (ja) 2007-05-23

Family

ID=23731422

Family Applications (1)

Application Number Title Priority Date Filing Date
JP11404396A Expired - Fee Related JP3917684B2 (ja) 1995-05-09 1996-05-09 物体の断層写真像を作成する方法及び装置

Country Status (4)

Country Link
US (1) US5541970A (ja)
JP (1) JP3917684B2 (ja)
DE (1) DE19618547A1 (ja)
IL (1) IL118131A0 (ja)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5974110A (en) * 1997-11-26 1999-10-26 General Electric Company Helical reconstruction algorithm
US6275562B1 (en) 1998-04-28 2001-08-14 General Electric Company Apparatus and methods for performing scalable multislice computed tomography scan
US6055291A (en) * 1998-08-25 2000-04-25 General Electric Company Systems, methods and apparatus for reconstructing images
US6198791B1 (en) 1998-08-25 2001-03-06 General Electric Company Scalable multislice imaging system
US6141398A (en) 1998-08-25 2000-10-31 General Electric Company Protocol driven image reconstruction, display, and processing in a multislice imaging system
EP1059879B1 (en) 1998-12-30 2007-08-08 General Electric Company Image thickness selection for multislice imaging system
US6185271B1 (en) 1999-02-16 2001-02-06 Richard Estyn Kinsinger Helical computed tomography with feedback scan control
US6466645B1 (en) * 2000-11-22 2002-10-15 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for tube-spit correction
US6522714B1 (en) * 2001-11-09 2003-02-18 Ge Medical Systems Global Technologies Company Llc Row-wise full helical view weighting method and apparatus for CT scanners
US7013034B2 (en) * 2002-10-10 2006-03-14 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for reconstructing an image of an object
WO2008133765A2 (en) * 2007-02-13 2008-11-06 Sentinel Scanning Corporation Ct scanning and contraband detection
US8340245B2 (en) * 2009-06-05 2012-12-25 Sentinel Scanning Corporation Transportation container inspection system and method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
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
US4821210A (en) * 1987-04-02 1989-04-11 General Electric Co. Fast display of three-dimensional images
US5047931A (en) * 1989-05-30 1991-09-10 General Electric Company Parallel computation of fan-beam back-projection reconstruction algorithm in computed tomography
US5270923A (en) * 1989-11-02 1993-12-14 General Electric Company Computed tomographic image reconstruction method for helical scanning using interpolation of partial scans for image construction
US5233518A (en) * 1989-11-13 1993-08-03 General Electric Company Extrapolative reconstruction method for helical scanning
US5253171A (en) * 1990-09-21 1993-10-12 General Electric Company Parallel processing method and apparatus based on the algebra reconstruction technique for reconstructing a three-dimensional computerized tomography (CT) image from cone beam projection data
US5265142A (en) * 1992-05-08 1993-11-23 General Electric Company Image reconstruction technique for a computer tomography system

Also Published As

Publication number Publication date
IL118131A0 (en) 1996-09-12
JPH0998968A (ja) 1997-04-15
DE19618547A1 (de) 1996-11-14
US5541970A (en) 1996-07-30

Similar Documents

Publication Publication Date Title
US7903779B2 (en) Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system
US6937690B2 (en) Method for computed tomography of a periodically moving object to be examined, and a CT unit for carrying out this method
US5559847A (en) Systems, methods and apparatus for reconstructing images in a CT system implementing a helical scan
US5663995A (en) Systems and methods for reconstructing an image in a CT system performing a cone beam helical scan
US6256366B1 (en) Apparatus and method for reconstruction of volumetric images in a computed tomography system using sementation of slices
US6256365B1 (en) Apparatus and method for reconstruction of images in a computed tomography system using oblique slices
US6072851A (en) Helical computerized tomography system and method providing quality images with faster scans
US6490333B1 (en) Methods and apparatus for cone-tilted parallel sampling and reconstruction
US5513236A (en) Image reconstruction for a CT system implementing a dual fan beam helical scan
JPH08509408A (ja) 円すい状ビームデータからの画像の再構成
EP1061476B1 (en) Volumetric image reconstruction
EP1077430A2 (en) Apparatus and method for reconstruction of volumetric images in a computed tomography system
US6904117B2 (en) Tilted gantry helical cone-beam Feldkamp reconstruction for multislice CT
US7529335B2 (en) Voxel-driven spiral reconstruction for cone-beam computer tomography
JPH05192322A (ja) 平面と交差する輻射線を用いた断層像再生方式
JPH08308824A (ja) 物体の断層写真像を作成するシステム
JPH07194590A (ja) 計算機式断層写真装置及び撮像されるべき物体の像を発生する方法
JP3917684B2 (ja) 物体の断層写真像を作成する方法及び装置
EP1372115A2 (en) Methods and apparatus for reconstructing an image of an object
EP0991022A2 (en) Methods and apparatus for image reconstruction
JPH07250832A (ja) 断層写真像を作成する方法
US6061423A (en) Fluoroscopy image reconstruction
US5974110A (en) Helical reconstruction algorithm
US5546439A (en) Systems, methods and apparatus for incrementally reconstructing overlapped images in a CT system implementing a helical scan
US6055291A (en) Systems, methods and apparatus for reconstructing images

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060418

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20060712

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20060718

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060831

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061003

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061215

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070209

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