JPH0669452B2 - 入れなおしを用いた扇状ビームらせん走査法 - Google Patents
入れなおしを用いた扇状ビームらせん走査法Info
- Publication number
- JPH0669452B2 JPH0669452B2 JP2315692A JP31569290A JPH0669452B2 JP H0669452 B2 JPH0669452 B2 JP H0669452B2 JP 2315692 A JP2315692 A JP 2315692A JP 31569290 A JP31569290 A JP 31569290A JP H0669452 B2 JPH0669452 B2 JP H0669452B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- projection set
- scan
- parallel beam
- beam projection
- 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
- 238000000034 method Methods 0.000 title claims description 35
- 238000003780 insertion Methods 0.000 title 1
- 238000003384 imaging method Methods 0.000 claims description 27
- 238000013213 extrapolation Methods 0.000 claims description 7
- 230000008707 rearrangement Effects 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000002591 computed tomography Methods 0.000 description 5
- 238000003325 tomography Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000013170 computed tomography imaging Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000013519 translation Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009429 distress Effects 0.000 description 1
- 238000002546 full scan Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
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/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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
に関するものである。更に詳しく述べると、本発明はら
せん走査で断層撮影投影データを取得することによって
生じる像アーチファクトを少なくするための像再構成方
法に関するものである。
コリメーションされて、規定された扇状ビーム角で扇状
ビームが形成される。扇状ビームは「イメージング平
面」と呼ばれるデカルト座標系のx−y平面内にあるよ
うに、またイメージング対象を透過してイメージング平
面内に配向されたx線検出器列に達するように配向され
る。検出器列は多数の検出素子で構成される。各検出素
子はx線源からその特定の検出素子に投射される射線に
沿って透過した放射線の強度を測定する。これらの検出
素子はそれぞれ扇状ビームの異なる射線に沿ったx線源
からのx線を遮えぎるように円弧状に配列することがで
きる。透過する放射線の強度はイメージング対象による
射線に沿ったx線ビームの減衰によってきまる。
イメージング平面内でガントリ上で回転させることがで
きる。これにより扇状ビームは異なる角度でイメージン
グ対象を横切る。各角度で、各検出素子からの強度信号
で構成される投影が取得される。次にガントリを新しい
角度まで回転して、上述の過程を反復することにより様
々の角度での多数の投射を収集して、1つの断層撮影投
影組を形成する。
る。これをコンピュータで処理して、当業者には既知の
再構成アルゴリズムに従ってスライス像を「再構成」す
ることができる。再構成されたスライス像は従来の陰極
線管にディスプレーしてもよいし、コンピュータ制御の
カメラによってフイルム記録に変換してもよい。
象の一連のスライスのイメージングが行なわれ、この一
連のスライスはx軸およびy軸に垂直なz軸に沿って増
分的に位置がずれている。これにより第3空間次元の情
報が得られる。放射線医はz軸に沿った位置の順にスラ
イス像を見ることによってこの第3次元を重い浮かべる
ことができる。あるいは再構成されたスライスの組を構
成する数値データをコンピュータ・プログラムで編集し
て、イメージング対象の三次元の陰影付き斜視図を作成
することもできる。
z次元で付加的なスライスが必要となる。断層撮影検査
の時間および費用は必要なスライス数がふえるにつれて
増大する。また、走査時間が長くなると、断層撮影像再
構成の忠実度を維持するためにほぼ不動でなければなら
ない患者の苦痛が増大する。したがって、一連のスライ
スを得るために必要な時間を減らすことにかなり関心が
集まっている。
時間は部分的に次の4つの構成要素によってきまる。す
なわちa)ガントリを走査速度まで加速するために必要
な時間、b)完全な1つの断層撮影投影組を得るために
必要な時間、c)ガントリを減速するために必要な時
間、およびd)次のスライスのためにz軸方向に患者を
再位置ぎめするために必要な時間によってきまる。全ス
ライス列を得るために必要な時間の短縮はこの4つのス
テップのいずれかを完了するために必要な時間を短縮す
ることによって行なうことができる。
信するケーブルではなくてスリップリングを使用する断
層撮影システムでは避けることができる。スリップリン
グによって、ガントリを連続的に回転することができ
る。以下に説明するCTシステムではスリップリングまた
は同等のものをそなえることにより360゜を超えて連続
的に回転することができるものとする。
ることが難しい。現在のCTスキャナでは1つのスライス
に対する投影組を取得するのに約1秒乃至2秒必要であ
る。この走査時間はガントリをより早い速度で回転させ
ることによって短縮することができる。一般に、ガント
リ速度が早くなると、取得したデータの信号対雑音比は
回転速度上昇率の平方根だけ小さくなる。これは透過形
断層撮影装置ではx線管の放射線出力を大きくすること
によりある程度は克服することができるが、このような
装置ではパワーに限界がある。
てz軸方向に患者を並進させることによって達成するこ
とができる。ガントリの回転中にz軸に沿って患者を一
定速度で並進させながら投影データを取得する方式は
「らせん走査」と呼ばれ、イメージング対象の物体上の
基準点に対するガントリ上の一点の見掛けの径路を表わ
している。ここで使用されているように、「らせん走
査」は一般に断層撮影イメージング・データの取得中に
患者またはイメージング対象の連続的な並進を使用する
ことを指す。また「一定z軸走査」は取得期間中に患者
またはイメージング対象を並進させることなく断層撮影
データ組を取得することを指す。
査相互の合間で患者を再位置ぎめするために通常必要と
される時間長がなくなり、与えられた数のスライスの取
得に必要とされる総走査時間が短縮される。しかし、ら
せん走査では取得された断層撮影投影組のデータについ
てあるエラーが生じる。断層撮影再構成の数学では一定
z軸スライス平面に沿って断層撮影投影組が取得される
と仮定している。らせん走査径路は明らかにこの条件か
らずれており、このずれの結果、z軸方向に対象に著し
い変化がある場合には再構成されたスライス像に像アー
チファクトが生じる。像アーチファクトのひどさは一般
に、走査データのテーブル位置と所望のスライス平面の
z軸値との差として測定された投影データの「らせんオ
フセット」によってきまる。らせん走査によって生じる
誤差はまとめて「スキュー」エラーと呼ばれる。
法が使用されてきた。1989年6月26日出願の米国特許出
願第371,332号(特願平2−163057号)「らせん投影走
査でスキュー像アーチファクトを減らすための方法」に
開示された第1の手法では、非一様なテーブルの動きを
使用することにより患者に加わる加速力を制限すると共
にらせん状に取得される投影をスライス平面の近くに集
中させている。
平2−295600号)「らせん走査のためのコンピュータ断
層撮影像再構成法」では180゜と扇状ビーム角度との和
の角度にわたるガントリ回転のみをそれぞれ必要とする
2つの半走査データの間で補間することによってスキュ
ーアーチファクトを減らしている。半走査に必要なガン
トリ回転が少なくなるので、テーブルの動きが少なくな
り、これにより投影データの全体のらせんオフセットが
少なくなる。
平2−304189号)「らせん走査用の補外式再構成方法」
に述べられている第3の手法では、180゜のみのガント
リ回転の2つの部分投影組の間で補間および補外を行な
うことによってスキューアーチファクトを少なくする。
2つの部分投影組は上記の半走査手法よりも更に少ない
ガントリ回転しか必要としないので、投影データの全体
のらせんオフセットは更に少なくなる。
転で取得される投影データから断層撮影像を作成するこ
とができる。一般に、この結果は180゜離れたガントリ
角度で取得された投影ではある射線の減衰が等しいこと
によって生じる。断層撮影像を再構成するこの方法は
「半走査」再構成と呼ばれる。半走査データ組からの像
の重み付けと再構成についてはメデイカル・フィジック
ス誌、9(2)、1982年3/4月号所載のデニス・エル
・パーカによる論文「扇状ビームに対する最適短走査コ
ンボリューション再構成」に述べられている。
ーム半走査からのらせんオフセットを小さくした投影組
の補間および補外を行なうことによりらせん状に取得し
たデータのスキューアーチファクトを小さくする。半走
査は接合(splicing)手順により2πのガントリ回転の
みで取得された扇状ビーム投影から作成される。
リ回転中に取得され、対応する平行ビーム投影組に入れ
なおされる。2つの半走査は入れなおされた平行ビーム
投影組から分割される。これらの半走査からのデータは
接合されて全2πの平行ビーム投影が作成される。半走
査は重み付けされてスライス平面に対する補間および補
外を行なえるようにした像、再構成されて像を形成す
る。
ス像に対する投影データを取得できるようにすることで
ある。接合過程により、平行ビーム半走査を360゜で取
得することができる。与えられた走査ピッチに対して72
0゜で取得された2つの全走査ではなくて360゜で取得さ
れた2つの平行ビーム半走査を使うことにより、らせん
走査で必要とされるz軸行程が短かくなる。これによ
り、取得される投影がスライス平面により近い点に集中
するので、補間と補外の正確さが向上し、部分的な容積
アーチファクトが少なくなる。
像の投影データを取得できるようにすることである。断
層撮影投影組の投影データの取得中の患者の動きによっ
て像アーチファクトが生じ得る。与えられたガントリ速
度に対して、360゜のみガントリ回転で取得される平行
半走査を使用することによって、モーション・アーチフ
ァクトすなわち動きによるアーチファクトの起りにくい
像の再構成が可能にになる。
ことである。360゜のガントリ回転で取得されたデータ
を接合して2つの平行ビーム半走査を形成することによ
り、患者のx先総被曝量を減らすことができる。
ら明らかとなる。以下の説明で参照する付図は本発明の
一部を形成するものであり、本発明の実施例を図示して
いる。しかし、このような実施例はかならずしも本発明
の全範囲を表わすものではないので、発明の範囲の解釈
にあたっては請求の範囲を参照しなければならない。
すCTガントリ16は、イメージング対象12を通して検出器
列18に扇状x線ビームの24を投影するように配向された
x線源10を含んでいる。扇状ビーム24はデカルト座標系
のx−y平面すなわち「イメージング平面」に沿った方
向を向いており、イメージング平面に沿って測った「扇
状角度」を形成する。検出器列18は多数の検出素子26で
構成される。多数の検出素子26はx線がイメージング対
象12を透過することによって生じる投影像を受けて、そ
の大きさに比例した値を検出する。
ントリに付設された制御モジュール48に結合されている
ので、360゜より大きい角度にわたって自由に連続的に
回転して投影データを取得することができる。
メージング過程に対する妨害を最小限にするようにテー
ブル22は放射線に対して半透明になっている。
切って、イメージング対象12に対して規定されたスライ
ス平面14を動かすことによって、x−yイメージング平
面に垂直なz軸に沿ってテーブル22の上表面が並進する
ようにテーブル22を制御することができる。簡単のため
に以後、テーブル22は一定速度で動き、したがってテー
ブル22のz軸位置はガントリ16の角度位置θに比例する
ものと仮定する。したがって、取得される断層撮影投影
はzまたはθによって規定することができる。
およびイメージング対象に対するイメージング平面のz
軸位置はそれぞれ一定z軸走査およびらせん走査に対し
て投影矢印20で示される。第2a図に示される一定z軸走
査では、各断層撮影投影組は一定z軸位置で取得され、
イメージング対象はこのような取得の合間にz軸に沿っ
て次のスライス平面へと動かされる。
場合には、イメージング平面に対するイメージング対象
のz軸位置は各断層撮影投影組の取得中に絶えず変化す
る。したがって矢印20はz軸に沿ってイメージング対象
の中でらせんを描く。らせんのピッチを走査ピッチと呼
ぶ。
ジング装置の制御システムはガントリに結合された制御
モジュール48をそなえている。制御モジュール48には電
力信号およびタイミング信号をx線源10に供給するx線
制御器54、ガントリ16の回転速度および位置を制御して
情報をコンピュータ60に供給するガントリ電動機制御器
56、ガントリ位置に関するデータ取得システム62、およ
びデータ取得システム62を介して検出器列18からサンプ
ルおよびディジタル化された信号を受けて当業者には知
られている方法に従って高速像再構成を行なう像再構成
器68が含まれている。上記の各々はスリップリング50を
介してガントリ16上のそれに対応する素子に接続するこ
とができ、コンピュータ60の種々のガントリ機能に対す
るインタフェースの役目を果す。
動機制御器52を介してコンピュータ60に伝えられ、コン
ピュータ60によって制御される。コンピュータ60は操作
卓64を介して指令および走査パラメータを受ける。操作
卓は一般にCRTディスプレーおよびキーボードであり、
これにより操作者は走査用のパラメータを入力したり、
コンピュータ60からの再構成された像等の情報を表示す
ることができる。大容量記憶装置66はCTイメージング装
置のためのオペレーティング・プログラムおよび操作者
が将来参照するための像デーアを記憶する手段を提供す
る。
データ要素は角度θおよびφによって表わすことができ
る。第4a図に示すように、角度φは第1図に示す扇状ビ
ーム24の真ん中の射線20から測ったものであり、扇状ビ
ーム24の中の射線21とそれに対応する検出器26を表わ
す。φは扇状ビーム角度と呼ばれる。θは(第1図に示
される)ガントリ16の角度位置であり、扇状ビームの真
ん中の射線20が垂直で下向きになっているとき任意に基
準値0とされる。線源10とガントリ16の回転の中心との
間の距離はDと表わされ、後で参照する。
影データが取得され、スライス像に再構成される。第5a
図に示されるように、投影組70に対するデータは垂直軸
測定引数θおよび水平軸測定引数φを持つデカルト「扇
状ビーム」空間の中の長方形の領域を充たす。θが一定
の水平線はガントリ位置θで取得された単一の投影を表
わし、−φmax<φ<+φmaxと表わされる角度からの検
出器信号を含んでいる。最も低い投影のガントリ角度θ
が任意に0に割り当てられ、投影組70の第1の投影であ
る。上記のらせん走査技術により、テーブル22がz軸に
沿って進んでいる間にガントリ角度θがθ=2πラジア
ンまで増大していく状態で相次ぐ投影が取得される。
0からπまで進められて、第1の部分的な扇状ビーム投
影組72が取得される。この取得の終りに、(第1図に示
される)イメージング対象12のスライス平面14がイメー
ジング平面とそろっている。次に第2の部分的な扇状ビ
ーム投影組74が開始され、ガントリ角度θ=πで始まっ
てガントリ角度θ=2πまで続く。当業者には既知の扇
状ビーム再構成技術によって2πラジアンの完全な扇状
ビーム投影組70を像に変換することができる。
投影組72および74は「平行ビーム」投影に入れなおすこ
とができる。このような入れなおし(rebinning)は米
国特許第4,852,132号「x線断層撮影のためのデータ収
集方法」に述べられている。名前が示すように、平行ビ
ーム投影は各投影が平行射線のみを持っているような投
影である。
び検出器26′によって得られる。平行投影組のデータ要
素は変数βおよびtによって表わすことができる。平行
ビーム24′の真ん中の射線20′から測った距離tは平行
ビーム24′の各射線21′およびそれに対応する検出器2
6′を表わし、平行ビーム・オフセットと呼ばれる。β
はガントリ16(図示されていない)の角度位置であり、
各射線21′の角度を規定する。扇状ビーム方式のθと同
様、平行ビームの真ん中の射線20′が垂直で下向きのと
き任意に基準値0とされる。
際、第5a図の扇状ビーム投影の各射線21が分離され、新
しい平行投影に分類される。分類は扇状ビーム装置によ
って取得された投影データと平行ビーム方式によって取
得された投影データとの間の以下の関係によって支配さ
れる。扇状ビーム投影組と平行ビーム投影組のそれぞれ
の任意の2つのデータ要素P1およびP2に対して P1(θ,φ)=P2(β,t) (1) 但し、 β=θ+φ (2) t=Dsin(φ) (3) とする。
記のようにβを測る垂直軸およびγを測る水平軸を持
つ、第5b図のデカルト「平行ビーム」空間に示されるよ
うな平行ビーム投影組76が得られる。
フセットを小さくした新しい投影組を補間および補外す
ることにより像再構成の前に第5a図の入れなおされた平
行ビーム投影に対してスキューアーチファクトの補正を
行なうことができる。
の走査データである必要はない。180゜の投影データし
かない平行ビーム投影組を像に再構成した後、補外およ
び補間の過程のために使用することができる。このよう
に削減された投影組は「半走査」と呼ばれる。
360゜の全平行ビーム投影組の中のデータの冗長性によ
って生じる。この冗長性の起源は第4b図を点検すること
によって明らかとなる。非らせん走査の場合、すなわち
走査中にイメージング対象12が動かない場合には、ガン
トリ角度βの任意の投影の中の射線24はガントリ角度β
+πラジアンで取得される投影の射線21に対して丁度逆
の180゜となる。イメージング対象12による射線21の減
衰は射線21の方向に無関係であるので、2つの同時に生
じるが逆向きの射線21に対して得られるデータ要素は同
じとなる。そしてこの2つのガントリ角度に対する投影
データは同じになる。但し、データの順序は逆になる。
更に詳しく述べると、任意の2つのデータ要素Paおよび
Pbについて次式が成り立つ。
ージング対照12はガントリ16の回転につれて動くので、
対向する角度の2つの射線21に対して得られる投影デー
タは異なる。それにも拘わらず、上記の式(8)は他の
データ要素対に比べて非常に相関度が高くなると予想さ
れ得る投影相互の間のデータ要素対を記述している。ら
せん走査で得られるデータに対する式(8)の関係は
「冗長性」と呼ばれる。
必要とし、2πの平行ビーム投影で2つの半走査を得る
ことができる。
状ビーム投影組70を入れなおしたデータから2πの完全
な平行ビーム投影データを得ることはできないというこ
とがわかる。詳しく述べると、2π+γ<β<2πであ
る領域82および0<β<γである領域84はある角度βに
対する投影の欠落部分を表わす。
データを得るために、平行ビーム空間の別のところから
これらの領域82および84にデータを接合しなければなら
ない。このようなデータは上記の冗長性の式(8)によ
って領域82および84の欠落データに関連付けることが好
ましい。信号対雑音比に配慮すると、再構成過程で処理
されない領域からのデータを接合することが好ましい。
2π<β<2π+γである領域86およびγ<β<0であ
る領域88はこれらの必要条件を満たす。したがって、式
(8)の関係に従って領域88のデータが領域82に接合さ
れ、領域86のデータが領域84に接合される。
る。第1の半走査は0<β<πの範囲のデータから構成
され、第2の半走査はπ<β<2πの範囲のデータから
構成される。次にこの2つの半走査は適宜重み付けして
加算することによりスライス平面に補間および補外する
ことができる。あるいは代案として、そしてより効率的
には、像再構成過程の陰の加算によって重み付けと加算
を行なってもよいことは当業者には理解されよう。
な補間および補外の重みは、上記の式(8)による対応
する冗長データ要素のスライス平面からの距離に対す
る、データ要素のスライス平面からの距離によってきま
る。重み付けは冗長データ要素の値にそれぞれの重みを
乗算することによって行なわれる。
けるP2(β,γ)の任意の2つの冗長データ要素に対し
て、ZSPスライス平面に対する線形の補間または補外の
ための点P1に対する重みW1は次式で表わされる。
データ要素を決定する必要がある。第6a図に示すよう
に、扇状ビーム・データ70から入れなおされた元の平行
ビーム・データ76は式(8)によりβ>π+γの組90お
よびβ<π+γの組92に分割される。領域90および92は
冗長データの組を表わすので、上記の式(9)および
(10)による別々の重み関数を必要とする。更に接合動
作により組92および90のデータ要素のいくつかが転置さ
れ、転位を補償する付加的な独自の重み関数を必要とす
るような転位されたデータの付加的な領域が生じる。
4つの領域が作られ、各領域は異なる重みを必要とす
る。
わされている組80および78の一部分としてのそれらの起
源を反映している。
れらの領域のデータ要素のzの値を決定しなければなら
ない。各データ要素のzの値は扇状ビーム投影組の対応
するデータ要素に対するθの値に比例する。したがっ
て、 z(β,γ)=k(θ) (11) =k(β−γ)[式(7)および(8)による](12) スライス平面のzの値は前に定めたようにk(π)であ
る。
に決定することができる。
同じであるが、接合手順の結果として2πだけ移され
る。したがって 領域2′に対する重み係数W2′は次式で表わされる。
るため、領域1,1′,2および2′の間の境界は不連続に
なる。これらの不連続により最終像に縞状の像アーチフ
ァクトが生じ得る。それらの領域の界面の近くでW1,
W1′,W2およびW2′のフェザリング(feathering)を行
なうことによって不連続をなくすことができる。高さω
の領域相互の間の領域でフェザリングが行なわれる。10
個の検出素子26が張る角度に等しいωの値が充分である
ことがわかる。
フェザリング関数f1(β,γ),f1′(β,γ),f
2(β,γ),f2′(β,γ)を乗算し、積を投影組全体
のデータに印加する。ここで 但し、 但し、 但し、 但し、 本発明の趣旨と範囲に入る実施例の多数の変形および変
更は当業者には明らかであろう。たとえば、線形補間以
外の補間方法を使うことができる。第1および第2の半
走査の前後の付加的な半走査からのデータを使う補間方
法、より高次の補間方法等である。更に、線形関係のフ
ェザリング関数のような他もフェザリング関数を使うこ
ともできる。
を含むCT装置の簡略斜視図であり、相対角度とそれと結
合された軸を示す。第2a図および第2b図は、第1図のイ
メージング対象を簡略に表わした斜視図であり、それぞ
れ一定z軸走査とらせん走査の場合のイメージング対象
に対するガントリとイメージング平面の相対配向を示
し、らせん走査のピッチはわかりやすくするため誇張し
て示してある。第3図は、第1図のCT装置に使用し得る
CT制御システムのブロック図である。第4a図は、x線扇
状ビームCT装置の幾何学的配置、ならびに扇状ビーム投
影組の各データ要素を規定する引数θおよびφの関係を
示す簡略平面図である。第4b図は、x線平行ビームCT装
置の幾何学的配置、ならびに平行ビーム投影組の各デー
タ要素を規定する変数βおよびtの関係を示す簡略平面
図である。第5a図は、第1図のCT装置に於いてらせん走
査で取得された扇状ビーム投影組の投影データに対応す
る引数θおよびφを表わしたグラフである。第5b図は、
第5a図の扇状ビーム投影組を入れなおすことによって作
成される平行ビーム投影組の投影データに対応する引数
βおよびγを表わしたグラフである。第6a図は、入れな
おした投影組の中の冗長データを表わした第5b図と同様
なグラフである。第6b図は、2つの平行ビーム半走査を
作るために接合データとともに入れなおされた平行ビー
ム投影組を表わした第5b図および第6a図と同様なグラフ
である。 [主な符号の説明] 10……x線源、 12……イメージング対象、 14……スライス平面、 70……扇状ビーム投影組、 76……平行ビーム投影組、 90,92……平行ビーム・データを分割した組。
Claims (5)
- 【請求項1】らせん走査で取得されたデータからイメー
ジング対象の断層撮影像を作成する方法であって、該デ
ータがz軸を中心とした複数のガントリ角度θに於ける
イメージ平面内の一連の扇状ビーム投影として取得さ
れ、該投影が扇状ビーム角度φに於ける複数のデータを
含んでいる断層撮影像作成方法に於て a)イメージ平面と平行な、イメージング対象に対する
スライス平面ZSPを定めるステップ、 b)2πの線源の回転にわたって扇状ビーム投影組のデ
ータを取得するステップ、 c)上記の線源の回転を行いながらz軸に沿ってイメー
ジング対象を動かすことにより、扇状ビーム投影組の取
得中にイメージング平面がスライス平面を横切るように
するステップ、 d)扇状ビーム投影組を、平行ビーム・ガントリ角度β
および扇状ビーム・オフセットを有する平行ビーム投影
組に入れなおすステップであって、この入れなおしによ
る該平行ビーム投影組には2πの平行ビーム・ガントリ
角度にわたる完全な平行ビーム投影組と比べて冗長なデ
ータおよび欠落データがあるステップ、 e)上記の入れなおしにより得られた平行ビーム投影組
を2つの半走査に分割するステップ、 f)少なくとも冗長なデータについての平行ビーム・ガ
ントリ角度を変えて、一方の半走査の冗長なデータを他
方の半走査の欠落データとして用いることによって、2
つの半走査相互の間でデータを接合することにより、2
πの平行ビーム・ガントリ角度にわたる完全な平行ビー
ム投影組を作成するステップ、 g)半走査のデータに対して補外および補間を行なって
スライス平面の平行ビーム投影組を得るステップ、およ
び h)スライス平面の平行ビーム投影組を再構成してスラ
イス像を得るステップ、 を含むことを特徴とする断層撮影像作成方法。 - 【請求項2】半走査のデータに対して補外および補間を
行なうため、半走査に重み関数を印加し、ステップ
(h)で該重み付けされた半走査を再構成してスライス
像を得る請求項1記載の断層撮影像作成方法。 - 【請求項3】各半走査の中の冗長データ対に対する重み
関数は加算すると定数となり、任意の冗長データに対す
る重みがθの関数である請求項2記載の断層撮影像作成
方法。 - 【請求項4】半走査に対してフェザリング重みを印加す
るステップを含む請求項1記載の断層撮影像作成方法。 - 【請求項5】扇状ビーム投影組の取得の中間にスライス
平面がイメージング平面を横切る請求項1記載の断層撮
影像作成方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US440,530 | 1989-11-22 | ||
US07/440,530 US5216601A (en) | 1989-11-22 | 1989-11-22 | Method for fan beam helical scanning using rebinning |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH03186250A JPH03186250A (ja) | 1991-08-14 |
JPH0669452B2 true JPH0669452B2 (ja) | 1994-09-07 |
Family
ID=23749125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2315692A Expired - Lifetime JPH0669452B2 (ja) | 1989-11-22 | 1990-11-22 | 入れなおしを用いた扇状ビームらせん走査法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US5216601A (ja) |
EP (1) | EP0430549B1 (ja) |
JP (1) | JPH0669452B2 (ja) |
CA (1) | CA2021623A1 (ja) |
DE (1) | DE69033923T2 (ja) |
IL (1) | IL96320A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003530928A (ja) * | 2000-04-14 | 2003-10-21 | ゼネラル・エレクトリック・カンパニイ | 投影ビューの間の補間を用いた計算機式断層画像の再構成 |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5454019A (en) * | 1991-03-15 | 1995-09-26 | Hitachi, Ltd. | Computed tomography system |
US5315665A (en) * | 1991-03-20 | 1994-05-24 | Kabushiki Kaisha Toshiba | X-ray computerized tomography apparatus for obtaining consecutive tomographic images in either a forward or backward direction without performing interpolation for views over an entire 360 |
US5396528A (en) * | 1991-06-28 | 1995-03-07 | General Electric Company | Tomographic image reconstruction using cross-plane rays |
US5469486A (en) * | 1992-08-07 | 1995-11-21 | General Electric Company | Projection domain reconstruction method for helical scanning computed tomography apparatus with multi-column detector array employing overlapping beams |
US5524130A (en) * | 1992-08-28 | 1996-06-04 | Kabushiki Kaisha Toshiba | Computed tomography apparatus |
US5406479A (en) * | 1993-12-20 | 1995-04-11 | Imatron, Inc. | Method for rebinning and for correcting cone beam error in a fan beam computed tomographic scanner system |
US5469487A (en) * | 1993-12-30 | 1995-11-21 | General Electric Company | CT system with twin fan beam helical scan |
US5515409A (en) * | 1994-12-22 | 1996-05-07 | General Electric Company | Helical interpolative algorithm for image reconstruction in a CT system |
US5491735A (en) * | 1994-12-22 | 1996-02-13 | General Electric Company | Image reconstruction apparatus and method for helical scanning |
US6091840A (en) * | 1995-12-21 | 2000-07-18 | General Electric Company | Methods and apparatus for single slice helical image reconstruction in a computed tomography system |
US6408088B1 (en) * | 1995-12-21 | 2002-06-18 | General Electric Company | Methods and apparatus for single slice helical image reconstruction in a computed tomography system |
DE19711963C2 (de) * | 1996-05-02 | 2000-04-13 | Siemens Ag | Verfahren zur Bildrekonstruktion bei einem Computertomographen |
CN1107290C (zh) * | 1996-05-02 | 2003-04-30 | 西门子公司 | 计算机层析x射线摄影机的图像再现方法 |
US6061422A (en) * | 1997-04-30 | 2000-05-09 | Hitachi Medical Corporation | X-ray ct apparatus and an image reconstructing method for obtaining a normal reconstructed image by adding weights to a plurality of partial reconstructed images |
US6324242B1 (en) * | 1999-04-09 | 2001-11-27 | Arch Development Corporation | Fast reconstruction with uniform noise properties in half-scan tomography |
US6639965B1 (en) * | 1999-09-30 | 2003-10-28 | General Electric Company | Methods and apparatus for cardiac imaging with conventional computed tomography |
US6411670B1 (en) | 1999-11-17 | 2002-06-25 | General Electric Company | Data rebinning to increase resolution in CT image reconstruction |
US6522712B1 (en) | 1999-11-19 | 2003-02-18 | General Electric Company | Reconstruction of computed tomographic images using interpolation between projection views |
US6452996B1 (en) | 2001-03-16 | 2002-09-17 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus utilizing generalized helical interpolation algorithm |
US6999550B2 (en) * | 2004-02-09 | 2006-02-14 | Ge Medical Systems Global Technology | Method and apparatus for obtaining data for reconstructing images of an object |
JP4498023B2 (ja) | 2004-06-15 | 2010-07-07 | キヤノン株式会社 | X線ct装置 |
US20060020200A1 (en) * | 2004-07-08 | 2006-01-26 | Medow Joshua E | Artifact-free CT angiogram |
US7583777B2 (en) | 2004-07-21 | 2009-09-01 | General Electric Company | Method and apparatus for 3D reconstruction of images |
US7570733B2 (en) * | 2005-06-10 | 2009-08-04 | General Electric Company | Step-and-shoot cardiac CT imaging |
US8218720B2 (en) * | 2007-03-12 | 2012-07-10 | Varian Medical Systems, Inc. | Method and apparatus to facilitate reconstructing an image using fan-beam data |
JP6125148B2 (ja) | 2012-03-14 | 2017-05-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成方法、画像生成装置および放射線断層撮影装置並びにプログラム |
JP5696963B2 (ja) | 2012-06-21 | 2015-04-08 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成方法、画像生成装置および放射線断層撮影装置並びにプログラム |
CN103961122B (zh) | 2013-01-31 | 2018-07-31 | 通用电气公司 | 用于非等γ角CT系统中数据转换的方法和装置 |
CN110533738B (zh) * | 2019-09-02 | 2021-06-18 | 上海联影医疗科技股份有限公司 | 重建数据处理方法、装置、医学成像系统及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4284896A (en) * | 1979-08-24 | 1981-08-18 | General Electric Company | Computerized tomographic reconstruction method and system utilizing reflection |
JPS59111738A (ja) * | 1982-12-16 | 1984-06-28 | 株式会社東芝 | X線断層撮影装置 |
US4707822A (en) * | 1985-05-09 | 1987-11-17 | Kabushiki Kaisha Toshiba | Tomographic apparatus |
US4789929A (en) * | 1987-05-14 | 1988-12-06 | Hitachi Medical Corporation | CT system for spirally scanning subject on a movable bed synchronized to X-ray tube revolution |
NL8800321A (nl) * | 1988-02-10 | 1989-09-01 | Philips Nv | Computertomografie-inrichting voor spiraalsgewijze aftasting. |
JPH0728862B2 (ja) * | 1989-02-13 | 1995-04-05 | 株式会社東芝 | Ct装置 |
-
1989
- 1989-11-22 US US07/440,530 patent/US5216601A/en not_active Expired - Lifetime
-
1990
- 1990-07-19 CA CA002021623A patent/CA2021623A1/en not_active Abandoned
- 1990-11-13 IL IL9632090A patent/IL96320A/en not_active IP Right Cessation
- 1990-11-21 EP EP90312638A patent/EP0430549B1/en not_active Expired - Lifetime
- 1990-11-21 DE DE69033923T patent/DE69033923T2/de not_active Expired - Fee Related
- 1990-11-22 JP JP2315692A patent/JPH0669452B2/ja not_active Expired - Lifetime
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003530928A (ja) * | 2000-04-14 | 2003-10-21 | ゼネラル・エレクトリック・カンパニイ | 投影ビューの間の補間を用いた計算機式断層画像の再構成 |
JP4773667B2 (ja) * | 2000-04-14 | 2011-09-14 | ゼネラル・エレクトリック・カンパニイ | 投影ビューの間の補間を用いた計算機式断層画像の再構成 |
Also Published As
Publication number | Publication date |
---|---|
DE69033923D1 (de) | 2002-04-11 |
IL96320A0 (en) | 1991-08-16 |
EP0430549A3 (en) | 1992-05-20 |
US5216601A (en) | 1993-06-01 |
JPH03186250A (ja) | 1991-08-14 |
EP0430549B1 (en) | 2002-03-06 |
EP0430549A2 (en) | 1991-06-05 |
DE69033923T2 (de) | 2002-10-24 |
CA2021623A1 (en) | 1991-05-23 |
IL96320A (en) | 1995-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0669452B2 (ja) | 入れなおしを用いた扇状ビームらせん走査法 | |
US5270923A (en) | Computed tomographic image reconstruction method for helical scanning using interpolation of partial scans for image construction | |
US5233518A (en) | Extrapolative reconstruction method for helical scanning | |
US5046003A (en) | Method for reducing skew image artifacts in helical projection imaging | |
JP4360817B2 (ja) | 放射線断層撮影装置 | |
US6256365B1 (en) | Apparatus and method for reconstruction of images in a computed tomography system using oblique slices | |
EP0520778B1 (en) | Tomographic image reconstruction using cross-plane rays | |
JPH08509408A (ja) | 円すい状ビームデータからの画像の再構成 | |
JP4612294B2 (ja) | X線コンピュータ断層装置、x線コンピュータ断層装置制御方法、及びx線コンピュータ断層撮影プログラム | |
JP2001087260A (ja) | ヘリカル・スキャン円錐形ビームによるコンピュータ断層撮影システムにおける三次元画像を再構成するための装置および方法 | |
US5208746A (en) | Method for helical scanning with a stationary detector using rebinning and splicing to create detector vertex projection sets | |
JPH07194590A (ja) | 計算機式断層写真装置及び撮像されるべき物体の像を発生する方法 | |
JPH10146331A (ja) | 対象物の関心ある領域の3次元コンピュータトモグラフ撮影のための走査及びデータ取得方法および装置 | |
US5170346A (en) | Method for reducing patient translation artifacts in tomographic imaging | |
US7215734B2 (en) | Method and system for three-dimensional reconstruction of images | |
JPH0998968A (ja) | 物体の断層写真像を作成する方法及び装置 | |
JP2006527618A (ja) | 余剰な測定値を使用するコンピュータ断層撮影法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080907 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090907 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100907 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110907 Year of fee payment: 17 |
|
EXPY | Cancellation because of completion of term | ||
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110907 Year of fee payment: 17 |