JP2000083948A - 断層写真画像を作る装置と方法 - Google Patents
断層写真画像を作る装置と方法Info
- Publication number
- JP2000083948A JP2000083948A JP11237663A JP23766399A JP2000083948A JP 2000083948 A JP2000083948 A JP 2000083948A JP 11237663 A JP11237663 A JP 11237663A JP 23766399 A JP23766399 A JP 23766399A JP 2000083948 A JP2000083948 A JP 2000083948A
- Authority
- JP
- Japan
- Prior art keywords
- data
- image
- projection data
- data array
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims description 8
- 238000003491 array Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims 1
- 230000002776 aggregation Effects 0.000 abstract 1
- 238000004220 aggregation Methods 0.000 abstract 1
- 238000002591 computed tomography Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013519 translation Methods 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- 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 Processing (AREA)
Abstract
(57)【要約】
【課題】人為効果を低レベルにして、品質の高い画像が
得られるような形で、４ビーム螺旋走査で得られたデー
タから、計算機式断層写真（ＣＴ）画像を作る装置と方
法を提供する。 【解決手段】 この発明の１形式では、４ビーム螺旋走
査によって得られるデータを使って、画像の再構成を実
施する装置を提供する。画像を再構成するとき、投影デ
ータ配列を発生する。次にこの投影データに加重係数を
加重して、加重済み投影データ配列を発生する。加重済
み投影データ配列をフィルタ作用及び逆投影にかけて、
画像データ配列を発生する。次に各ビームに対する画像
データ配列を加算して、スライス画像データ配列を発生
する。
得られるような形で、４ビーム螺旋走査で得られたデー
タから、計算機式断層写真（ＣＴ）画像を作る装置と方
法を提供する。 【解決手段】 この発明の１形式では、４ビーム螺旋走
査によって得られるデータを使って、画像の再構成を実
施する装置を提供する。画像を再構成するとき、投影デ
ータ配列を発生する。次にこの投影データに加重係数を
加重して、加重済み投影データ配列を発生する。加重済
み投影データ配列をフィルタ作用及び逆投影にかけて、
画像データ配列を発生する。次に各ビームに対する画像
データ配列を加算して、スライス画像データ配列を発生
する。
Description
【０００１】
【発明の背景】この発明は全体的に医療用作像、更に具
体的に言えば、多重スライス・スキャナによる螺旋走査
から収集された投影データから、画像を再構成すること
に関する。典型的には計算機式断層写真（ＣＴ）装置と
呼ばれる少なくとも１つの公知の医療用作像装置では、
ｘ線源が扇形ビームを投射し、このビームは、直角座標
系のＸ−Ｙ平面内にあるようにコリメートされており、
この平面が一般的に「作像平面」と呼ばれる。ｘ線ビー
ムが患者のような作像される物体を通過する。ビーム
は、物体による減衰を受けた後、放射検出器の配列に入
射する。検出器配列が受取った減衰したビーム放射の強
度は、物体によるｘ線ビームの減衰に関係する。配列の
各々の検出素子が、検出器がある場所でのビームの減衰
の測定値である別々の電気信号を発生する。すべての検
出器からの減衰測定値を別々に収集して、透過分布を作
る。
体的に言えば、多重スライス・スキャナによる螺旋走査
から収集された投影データから、画像を再構成すること
に関する。典型的には計算機式断層写真（ＣＴ）装置と
呼ばれる少なくとも１つの公知の医療用作像装置では、
ｘ線源が扇形ビームを投射し、このビームは、直角座標
系のＸ−Ｙ平面内にあるようにコリメートされており、
この平面が一般的に「作像平面」と呼ばれる。ｘ線ビー
ムが患者のような作像される物体を通過する。ビーム
は、物体による減衰を受けた後、放射検出器の配列に入
射する。検出器配列が受取った減衰したビーム放射の強
度は、物体によるｘ線ビームの減衰に関係する。配列の
各々の検出素子が、検出器がある場所でのビームの減衰
の測定値である別々の電気信号を発生する。すべての検
出器からの減衰測定値を別々に収集して、透過分布を作
る。
【０００２】公知の第３世代ＣＴ装置では、ｘ線源及び
検出器配列が、作像平面内で、作像しようとする物体の
周りをガントリと共に回転し、この為、ｘ線ビームが物
体と交わる角度が絶えず変化する。１つのガントリ角度
で検出器配列から得られるｘ線減衰測定値の群、即ち、
投影データが「図」と呼ばれる。物体の「走査」は、ｘ
線源及び検出器が１回転する間の相異なるガントリ角度
で求められた１組の図で構成される。軸方向走査では、
投影データが、物体を通る２次元スライスに対応する画
像を構成するように処理される。１組の投影データから
画像を再構成する１つの方法が、業界では、フィルタ式
逆投影方式と呼ばれている。この過程は走査で得られた
減衰測定値を、「ＣＴ数」又は「ハウンズフィールド単
位」と呼ばれる整数に変換し、こういう整数を使って、
陰極線管表示装置の対応する画素の輝度を制御する。
検出器配列が、作像平面内で、作像しようとする物体の
周りをガントリと共に回転し、この為、ｘ線ビームが物
体と交わる角度が絶えず変化する。１つのガントリ角度
で検出器配列から得られるｘ線減衰測定値の群、即ち、
投影データが「図」と呼ばれる。物体の「走査」は、ｘ
線源及び検出器が１回転する間の相異なるガントリ角度
で求められた１組の図で構成される。軸方向走査では、
投影データが、物体を通る２次元スライスに対応する画
像を構成するように処理される。１組の投影データから
画像を再構成する１つの方法が、業界では、フィルタ式
逆投影方式と呼ばれている。この過程は走査で得られた
減衰測定値を、「ＣＴ数」又は「ハウンズフィールド単
位」と呼ばれる整数に変換し、こういう整数を使って、
陰極線管表示装置の対応する画素の輝度を制御する。
【０００３】多数のスライスに必要な合計の走査時間を
短縮する為、「螺旋」走査を実施することが出来る。
「螺旋」走査を実施するには、所定数のスライスに対す
るデータを収集する間、患者を移動させる。螺旋走査で
得られたデータから、画像を再構成する為に利用するこ
とが出来る画像再構成アルゴリズムが、１９９５年５月
９日に出願され、この出願の被譲渡人に譲渡された継続
中の米国特許出願通し番号０８／４３６，１７６号に記
載されている。
短縮する為、「螺旋」走査を実施することが出来る。
「螺旋」走査を実施するには、所定数のスライスに対す
るデータを収集する間、患者を移動させる。螺旋走査で
得られたデータから、画像を再構成する為に利用するこ
とが出来る画像再構成アルゴリズムが、１９９５年５月
９日に出願され、この出願の被譲渡人に譲渡された継続
中の米国特許出願通し番号０８／４３６，１７６号に記
載されている。
【０００４】扇形ビーム螺旋走査で収集される投影デー
タは、Ｐ（θ，γ，ｚ）で表すことが出来る。ここでθ
はある基準（例えばｙ軸）に対する扇形ビームの中心線
の角度、γは扇形ビーム内にある特定のｘ線の中心のｘ
線に対する角度、ｚは測定が行われる時点でのガントリ
の軸方向位置である。実際の投影データが得られない各
々の場所ｚ0では、普通に使われている公知の螺旋再構
成アルゴリズムが、ｚ方向の１次補間を使って、位置ｚ
0に於けるスライスに対する生データを発生する。具体
的に言うと、Ｐ（θ，γ，ｚ0）を発生するには、同じ
θ及びγで、ｚ0に出来るだけ近いが、その両側にある
ｚに於ける投影データを用いる。例えば、ｚ1及びｚ2が
Ｐ（θ，γ，ｚ）を利用出来るｚの値であって、ｚ1≦
ｚ0≦ｚ2である場合、下記の式を使って、１次補間によ
り、Ｐ（θ，γ，ｚ1）及びＰ（θ，γ，ｚ2）からＰ
（θ，γ，ｚ0）を推定することが出来る。
タは、Ｐ（θ，γ，ｚ）で表すことが出来る。ここでθ
はある基準（例えばｙ軸）に対する扇形ビームの中心線
の角度、γは扇形ビーム内にある特定のｘ線の中心のｘ
線に対する角度、ｚは測定が行われる時点でのガントリ
の軸方向位置である。実際の投影データが得られない各
々の場所ｚ0では、普通に使われている公知の螺旋再構
成アルゴリズムが、ｚ方向の１次補間を使って、位置ｚ
0に於けるスライスに対する生データを発生する。具体
的に言うと、Ｐ（θ，γ，ｚ0）を発生するには、同じ
θ及びγで、ｚ0に出来るだけ近いが、その両側にある
ｚに於ける投影データを用いる。例えば、ｚ1及びｚ2が
Ｐ（θ，γ，ｚ）を利用出来るｚの値であって、ｚ1≦
ｚ0≦ｚ2である場合、下記の式を使って、１次補間によ
り、Ｐ（θ，γ，ｚ1）及びＰ（θ，γ，ｚ2）からＰ
（θ，γ，ｚ0）を推定することが出来る。
【０００５】
【数１１】
【０００６】螺旋走査では、各々の３６０゜の回転で同
じｘ線が２回測定される、即ち、Ｐ（θ，γ，ｚ）＝Ｐ
（θ＋２γ＋１８０゜，−γ，ｚ）であるから、ｚの標
本化が実効的に２倍になる。このように標本化が増える
ことにより、合計走査時間を短縮することが出来る。勿
論、人為効果を低レベルにして、又はその数を少なくし
て、品質の高い画像が得られるような形で、４ビーム螺
旋走査で得られたデータから、画像を再構成することが
望ましい。更に、この様な画像を再構成する為に必要な
合計時間を短縮することが望ましい。更に、悉くの軸方
向の位置に対してデータが必ずしも得られないことがあ
るから、高い品質の画像を発生することができる様な形
で、この様な投影データを推定するアルゴリズムを提供
することも望ましい。
じｘ線が２回測定される、即ち、Ｐ（θ，γ，ｚ）＝Ｐ
（θ＋２γ＋１８０゜，−γ，ｚ）であるから、ｚの標
本化が実効的に２倍になる。このように標本化が増える
ことにより、合計走査時間を短縮することが出来る。勿
論、人為効果を低レベルにして、又はその数を少なくし
て、品質の高い画像が得られるような形で、４ビーム螺
旋走査で得られたデータから、画像を再構成することが
望ましい。更に、この様な画像を再構成する為に必要な
合計時間を短縮することが望ましい。更に、悉くの軸方
向の位置に対してデータが必ずしも得られないことがあ
るから、高い品質の画像を発生することができる様な形
で、この様な投影データを推定するアルゴリズムを提供
することも望ましい。
【０００７】
【発明の要約】上記並びにその他の目的を、１実施例で
は４扇形ビーム螺旋走査で各々の扇形ビームによって収
集された投影データから、投影空間データ配列を発生す
る装置によって達成することが出来る。この時、各々の
配列内データに装置によって加重して、患者の並進運動
に対する補正をすると共に、データの冗長効果を打消
す。その後、加重データを使って、画像を再構成する。
は４扇形ビーム螺旋走査で各々の扇形ビームによって収
集された投影データから、投影空間データ配列を発生す
る装置によって達成することが出来る。この時、各々の
配列内データに装置によって加重して、患者の並進運動
に対する補正をすると共に、データの冗長効果を打消
す。その後、加重データを使って、画像を再構成する。
【０００８】更に具体的に言うと、画像を再構成すると
き、装置が作像しようとするスライスに関連するデータ
平面に対応する投影データ配列を発生する。次に、装置
によってこのデータ配列に加重係数が適用され、各々
の、特定のデータ要素に対して重みを割当てる。その
後、加重済み投影データ配列をフィルタ作用及び逆投影
にかけて、画像データ配列を発生する。次に画像データ
配列を加算してスライス画像データ配列を発生する。
き、装置が作像しようとするスライスに関連するデータ
平面に対応する投影データ配列を発生する。次に、装置
によってこのデータ配列に加重係数が適用され、各々
の、特定のデータ要素に対して重みを割当てる。その
後、加重済み投影データ配列をフィルタ作用及び逆投影
にかけて、画像データ配列を発生する。次に画像データ
配列を加算してスライス画像データ配列を発生する。
【０００９】実際には投影データをそこで測定しなかっ
た特定の場所ｚ0に於ける特定のスライスに対して、ス
ライスを再構成する場合、１実施例では、特定のｚ0に
最も近いが、その両側にあるｚ位置に対する投影データ
を確認する。その後、このスライスに対する投影データ
は、こうして確認されたｚ位置で収集された投影データ
を使って推定する。次にスライス画像をこの推定された
投影データを用いて再構成することが出来る。
た特定の場所ｚ0に於ける特定のスライスに対して、ス
ライスを再構成する場合、１実施例では、特定のｚ0に
最も近いが、その両側にあるｚ位置に対する投影データ
を確認する。その後、このスライスに対する投影データ
は、こうして確認されたｚ位置で収集された投影データ
を使って推定する。次にスライス画像をこの推定された
投影データを用いて再構成することが出来る。
【００１０】患者の４ビーム螺旋走査は、合計の患者走
査時間が短縮するという利点がある。更に、前に述べた
画像再構成アルゴリズムは、患者テーブルの並進速度を
高くしても、高い品質の画像スライスを発生することが
出来るという利点をもたらす。
査時間が短縮するという利点がある。更に、前に述べた
画像再構成アルゴリズムは、患者テーブルの並進速度を
高くしても、高い品質の画像スライスを発生することが
出来るという利点をもたらす。
【００１１】
【図面の詳しい説明】次に一例の多重スライスＣＴ装置
を説明する。１実施例の装置をこれから詳しく説明する
が、この発明はこの代りの多くの多重スライス装置に関
連して利用することが出来ることを承知されたい。例え
ば、特定の１つの検出器を説明するが、この発明は他の
検出器に関連して用いることが出来るし、この発明を実
施するには何ら特定の１形式の多重スライス検出器に制
限されない。具体的に言うと、これから説明する検出器
は複数個のモジュールを持ち、各々のモジュールが複数
個の検出セルを含む。これから説明する特定の検出器の
代りに、ｘ軸並びに／又はｚ軸に沿った多重素子を何れ
かの方向に一緒に結合した多数のモジュールを持ってい
て、多重スライス走査データを同時に収集する検出器を
利用することが出来る。全般的に言うと、この装置は多
重スライスモードで、１つ又は更に多くのデータ・スラ
イスを収集する様に作用し得る。この装置を用いて軸方
向走査及び螺旋走査を実施することができ、走査された
物体の断面画像を処理し、再構成し、表示し並びに／又
は記録することが出来る。
を説明する。１実施例の装置をこれから詳しく説明する
が、この発明はこの代りの多くの多重スライス装置に関
連して利用することが出来ることを承知されたい。例え
ば、特定の１つの検出器を説明するが、この発明は他の
検出器に関連して用いることが出来るし、この発明を実
施するには何ら特定の１形式の多重スライス検出器に制
限されない。具体的に言うと、これから説明する検出器
は複数個のモジュールを持ち、各々のモジュールが複数
個の検出セルを含む。これから説明する特定の検出器の
代りに、ｘ軸並びに／又はｚ軸に沿った多重素子を何れ
かの方向に一緒に結合した多数のモジュールを持ってい
て、多重スライス走査データを同時に収集する検出器を
利用することが出来る。全般的に言うと、この装置は多
重スライスモードで、１つ又は更に多くのデータ・スラ
イスを収集する様に作用し得る。この装置を用いて軸方
向走査及び螺旋走査を実施することができ、走査された
物体の断面画像を処理し、再構成し、表示し並びに／又
は記録することが出来る。
【００１２】図１及び２について説明すると、計算機式
断層写真（ＣＴ）作像装置１０が、「第３世代」ＣＴス
キャナを表すガントリ１２を含むことが示されている。
ガントリ１２がｚ線源１４を持ち、このｘ線源が、ｘ線
ビーム１６を、ガントリ１２の反対側にある検出器配列
１８に向かって投射する。検出器配列１８は検出素子２
０によって形成され、これらの検出素子が一緒になっ
て、医療の患者２２を通過した、投射されたｘ線を感知
する。各々の検出素子２０が、入射するｘ線ビームの強
度、従って患者２２を通過したときのビームの減衰を表
す電気信号を発生する。ｘ線投影データを収集する為の
走査の間、ガントリ１２及びその上に装着された部品が
回転中心２４の周りを回転する。
断層写真（ＣＴ）作像装置１０が、「第３世代」ＣＴス
キャナを表すガントリ１２を含むことが示されている。
ガントリ１２がｚ線源１４を持ち、このｘ線源が、ｘ線
ビーム１６を、ガントリ１２の反対側にある検出器配列
１８に向かって投射する。検出器配列１８は検出素子２
０によって形成され、これらの検出素子が一緒になっ
て、医療の患者２２を通過した、投射されたｘ線を感知
する。各々の検出素子２０が、入射するｘ線ビームの強
度、従って患者２２を通過したときのビームの減衰を表
す電気信号を発生する。ｘ線投影データを収集する為の
走査の間、ガントリ１２及びその上に装着された部品が
回転中心２４の周りを回転する。
【００１３】ガントリ１２の回転並びにｘ線源１４の動
作がＣＴ装置１０の制御機構２６によって制御される。
制御機構２６が、ｘ線源１４に対する電力及びタイミン
グ信号を供給するｘ線制御器２８と、ガントリ１２の回
転速度及び位置を制御するガントリ・モータ制御器３０
を含む。制御機構２６にあるデータ収集装置（ＤＡＳ）
３２が、検出素子２０からのアナログ・データを標本化
し、そのデータをこの後の処理の為にディジタル信号に
変換する。画像再構成装置３４は、ＤＡＳ ３２からの
標本化されてディジタル化されたｘ線データを受取り、
高速画像再構成を行う。再構成された画像がコンピュー
タ３６に対する入力として印加され、このコンピュータ
が大量記憶装置３８に画像を記憶する。
作がＣＴ装置１０の制御機構２６によって制御される。
制御機構２６が、ｘ線源１４に対する電力及びタイミン
グ信号を供給するｘ線制御器２８と、ガントリ１２の回
転速度及び位置を制御するガントリ・モータ制御器３０
を含む。制御機構２６にあるデータ収集装置（ＤＡＳ）
３２が、検出素子２０からのアナログ・データを標本化
し、そのデータをこの後の処理の為にディジタル信号に
変換する。画像再構成装置３４は、ＤＡＳ ３２からの
標本化されてディジタル化されたｘ線データを受取り、
高速画像再構成を行う。再構成された画像がコンピュー
タ３６に対する入力として印加され、このコンピュータ
が大量記憶装置３８に画像を記憶する。
【００１４】コンピュータ３６は、ユーザ・インターフ
ェイス又はグラフ式ユーザ・インターフェイス（ＧＵ
Ｉ）を通じて信号を受取ると共に供給する。具体的に言
うと、コンピュータが、キーボード及びマウス（図に示
していない）を持つコンソール４０を介して、オペレー
タからの指令及び走査パラメータを受取る。関連した陰
極線管表示装置４２が、オペレータが、コンピュータ３
６からの再構成された画像及びその他のデータを観察す
ることが出来る様にする。オペレータから供給された指
令及びパラメータが、ｘ線制御器２８、ガントリ・モー
タ制御器３０、ＤＡＳ ３２及びテーブル・モータ制御
器４４に対する制御信号及び情報を供給する為に、コン
ピュータ３６によって使われる。
ェイス又はグラフ式ユーザ・インターフェイス（ＧＵ
Ｉ）を通じて信号を受取ると共に供給する。具体的に言
うと、コンピュータが、キーボード及びマウス（図に示
していない）を持つコンソール４０を介して、オペレー
タからの指令及び走査パラメータを受取る。関連した陰
極線管表示装置４２が、オペレータが、コンピュータ３
６からの再構成された画像及びその他のデータを観察す
ることが出来る様にする。オペレータから供給された指
令及びパラメータが、ｘ線制御器２８、ガントリ・モー
タ制御器３０、ＤＡＳ ３２及びテーブル・モータ制御
器４４に対する制御信号及び情報を供給する為に、コン
ピュータ３６によって使われる。
【００１５】図３に示す様に、４扇形ビーム装置では、
４列の検出器が構成されている。ｘ線扇形ビームが、事
実上、回転ｚ軸に沿って変位した４つの扇形ビームに分
割される。隣り合ったビームの中心の間の距離は、ガン
トリの回転軸線のところで測定したとき、Ｄである。画
像再構成装置３４が図４に更に詳しく示されている。具
体的に言うと、各々の扇形ビームに対して、ＤＡＳ ３
２から得られるデータの各図を夫々のプリプロセッサ５
２Ａ−Ｄで受取り、そこで、夫々のビームは、ビーム硬
化、オフセット及び検出器及びチャンネル利得の変動の
様な種々の周知の誤差を補正する為に前処理される。更
に、負の対数を求めて、投影データを作り、それが投影
データ配列５４Ａ−Ｄに記憶される。
４列の検出器が構成されている。ｘ線扇形ビームが、事
実上、回転ｚ軸に沿って変位した４つの扇形ビームに分
割される。隣り合ったビームの中心の間の距離は、ガン
トリの回転軸線のところで測定したとき、Ｄである。画
像再構成装置３４が図４に更に詳しく示されている。具
体的に言うと、各々の扇形ビームに対して、ＤＡＳ ３
２から得られるデータの各図を夫々のプリプロセッサ５
２Ａ−Ｄで受取り、そこで、夫々のビームは、ビーム硬
化、オフセット及び検出器及びチャンネル利得の変動の
様な種々の周知の誤差を補正する為に前処理される。更
に、負の対数を求めて、投影データを作り、それが投影
データ配列５４Ａ−Ｄに記憶される。
【００１６】各々のビームに対する投影データ配列５４
Ａ−Ｄを読出して、掛け算器５８Ａ−Ｄで対応する加重
関数５６Ａ−Ｄを適用する。加重済み投影データが加重
済み投影データ配列６０Ａ−Ｄの対応する場所に書き込
まれ、この加重済み投影データを６２Ａ−Ｄのところで
フィルタ作用にかけると共に逆投影して、ビーム画像デ
ータ配列６４Ａ−Ｄを作る。
Ａ−Ｄを読出して、掛け算器５８Ａ−Ｄで対応する加重
関数５６Ａ−Ｄを適用する。加重済み投影データが加重
済み投影データ配列６０Ａ−Ｄの対応する場所に書き込
まれ、この加重済み投影データを６２Ａ−Ｄのところで
フィルタ作用にかけると共に逆投影して、ビーム画像デ
ータ配列６４Ａ−Ｄを作る。
【００１７】その後、６６のところで画像データ配列６
４Ａ−Ｄを加算してスライス画像データ配列６８を作
る。具体的に言うと、ビーム１配列にある各々の画素の
大きさが、ビーム２、ビーム３及びビーム４配列にある
対応する画素の大きさと加算される。この結果得られた
スライス画像配列６８は、後で使う為に記憶してもよい
し、或いはオペレータに表示することも出来る。画像デ
ータ配列６４Ａ−Ｄを発生した後に、データを加算する
代りに、同じガントリ（又は図の）角度であるが、異な
る検出器列からの投影を、データのフィルタ作用及び逆
投影の前に組合わせることが出来る。この様な組合せ
は、処理の負担を軽減することが出来る。
４Ａ−Ｄを加算してスライス画像データ配列６８を作
る。具体的に言うと、ビーム１配列にある各々の画素の
大きさが、ビーム２、ビーム３及びビーム４配列にある
対応する画素の大きさと加算される。この結果得られた
スライス画像配列６８は、後で使う為に記憶してもよい
し、或いはオペレータに表示することも出来る。画像デ
ータ配列６４Ａ−Ｄを発生した後に、データを加算する
代りに、同じガントリ（又は図の）角度であるが、異な
る検出器列からの投影を、データのフィルタ作用及び逆
投影の前に組合わせることが出来る。この様な組合せ
は、処理の負担を軽減することが出来る。
【００１８】この発明の１実施例は、特に、ある予定の
条件の下で４ビーム走査が実施されたとき、加重済み投
影データ配列６０Ａ−Ｄを作ることに関係する。以下の
説明で、ｄはガントリの回転軸線のところで測った検出
器列の（ｚ）間隔、ｓはテーブル送り速度（回転毎）、
ｐはｄ及びｓの比である。即ち、 ｐ＝ｄ／ｓ （２） 図５に示す様に、データ平面Ｐ1、Ｐ2、Ｐ3及びＰ4が、
線Ｌ1Ｍ1、Ｌ2Ｍ2、Ｌ 3Ｍ3及びＬ4Ｍ4で、再構成しよう
とするスライスＰと交わる。これらの線関数は次の様に
表すことが出来る。
条件の下で４ビーム走査が実施されたとき、加重済み投
影データ配列６０Ａ−Ｄを作ることに関係する。以下の
説明で、ｄはガントリの回転軸線のところで測った検出
器列の（ｚ）間隔、ｓはテーブル送り速度（回転毎）、
ｐはｄ及びｓの比である。即ち、 ｐ＝ｄ／ｓ （２） 図５に示す様に、データ平面Ｐ1、Ｐ2、Ｐ3及びＰ4が、
線Ｌ1Ｍ1、Ｌ2Ｍ2、Ｌ 3Ｍ3及びＬ4Ｍ4で、再構成しよう
とするスライスＰと交わる。これらの線関数は次の様に
表すことが出来る。
【００１９】 Ｚ[1-4]＝［−１．５，−０．５，０．５，１．５］ｄ βn＝２Ｚnｐπ 即ち、Ｌ1Ｍ1：β1＝−３ｐπ Ｌ2Ｍ2：β2＝−３ｐπ （３） Ｌ3Ｍ3：β3＝ｐπ Ｌ4Ｍ4：β4＝３ｐπ ここでβはガントリ角度に等しい。線Ｌ1Ｍ1、Ｌ2Ｍ2、
Ｌ3Ｍ3及びＬ4Ｍ4は、下記の＋及び−の集合で表す「鏡
像」線を持っている。
Ｌ3Ｍ3及びＬ4Ｍ4は、下記の＋及び−の集合で表す「鏡
像」線を持っている。
【００２０】
【数１２】
【００２１】テーブル送り速度ｓ及び検出器ｚの間隔ｄ
が、γmを扇形の角度の半分と定義して、（２π／（π
−２γm））ｄ＜ｓ＜（４π／（π＋２γm））ｄという
関係を満たすとき、Ｗ1（β，γ）、Ｗ2（β，γ）、Ｗ
3（β，γ）及びＷ4（β，γ）と記す、各々のデータ集
合に適用すべき螺旋加重係数は、次の様になる。
が、γmを扇形の角度の半分と定義して、（２π／（π
−２γm））ｄ＜ｓ＜（４π／（π＋２γm））ｄという
関係を満たすとき、Ｗ1（β，γ）、Ｗ2（β，γ）、Ｗ
3（β，γ）及びＷ4（β，γ）と記す、各々のデータ集
合に適用すべき螺旋加重係数は、次の様になる。
【００２２】
【数１３】
【００２３】ここで
【００２４】
【数１４】
【００２５】であり、
【００２６】
【数１５】
【００２７】ここで
【００２８】
【数１６】
【００２９】であり、
【００３０】
【数１７】
【００３１】式（５）−（８）で、βM及びβmは、夫々
β1+及びβ4-の内の大きい方並びに小さい方を表す。式
（５）−（８）に示す加重関数は連続している。しか
し、この各々の式の１次微分は境界線で不連続である。
必要であれば、境界に互る数チャンネル（〜１０チャン
ネル）のフェザリングによって、この不連続性をなくす
ことが出来る。
β1+及びβ4-の内の大きい方並びに小さい方を表す。式
（５）−（８）に示す加重関数は連続している。しか
し、この各々の式の１次微分は境界線で不連続である。
必要であれば、境界に互る数チャンネル（〜１０チャン
ネル）のフェザリングによって、この不連続性をなくす
ことが出来る。
【００３２】一旦加重済み投影データ配列６０Ａ−Ｄを
作る様に加重したら、処理時間を短縮する為、フィルタ
作用及び逆投影にかける前に、同じガントリ（又は図
の）の角度であるが、異なる検出器列からの投影を組合
わせることが出来る。データ列１にある或る投影図は、
データ列４にある対応する投影図から３６０゜離れてい
る。これらの図の対をフィルタ作用にかける前に更に組
合わせて、処理の負担の不必要な増加をなくすことが出
来る。
作る様に加重したら、処理時間を短縮する為、フィルタ
作用及び逆投影にかける前に、同じガントリ（又は図
の）の角度であるが、異なる検出器列からの投影を組合
わせることが出来る。データ列１にある或る投影図は、
データ列４にある対応する投影図から３６０゜離れてい
る。これらの図の対をフィルタ作用にかける前に更に組
合わせて、処理の負担の不必要な増加をなくすことが出
来る。
【００３３】隣接するスライスを再構成する為に必要な
投影データは、図の角度βの原点を垂直方向にシフトさ
せて、再構成しようとする新しいスライスと整合させる
ことにより、確認することが出来る。大抵の場合、各々
のデータ集合にある或るスライスと隣接するスライスに
対するデータの間にはかなりの重なりがある。加重の
前、前処理はスライスの位置に関係しない。この為（螺
旋加重をしない）前処理の結果を、将来再び使う為にバ
ッファに記憶することが出来る。これは多くの場合、特
に遡及形再構成の場合、前処理の負担を大幅に減らすこ
とがある。
投影データは、図の角度βの原点を垂直方向にシフトさ
せて、再構成しようとする新しいスライスと整合させる
ことにより、確認することが出来る。大抵の場合、各々
のデータ集合にある或るスライスと隣接するスライスに
対するデータの間にはかなりの重なりがある。加重の
前、前処理はスライスの位置に関係しない。この為（螺
旋加重をしない）前処理の結果を、将来再び使う為にバ
ッファに記憶することが出来る。これは多くの場合、特
に遡及形再構成の場合、前処理の負担を大幅に減らすこ
とがある。
【００３４】所望のスライスの輪郭が、データ並びに前
述の再構成アルゴリズムによって支援される輪郭より厚
手である場合、所望のスライスの輪郭内にある幾つかの
薄いスライスを加算することにより、一層厚手のスライ
スを導き出すことが出来る。幾つかの薄いスライス自体
に関心がなければ、投影領域の早期に対応する加算を実
施することにより、幾つかの薄いスライスを再構成する
中間工程を回避することが出来る。こうすると、計算の
負担並びに画像を記憶する負担が減少する。その結果得
られる加重関数は、データ平面の対応するシフトした形
を加算することによって導き出すことが出来る。
述の再構成アルゴリズムによって支援される輪郭より厚
手である場合、所望のスライスの輪郭内にある幾つかの
薄いスライスを加算することにより、一層厚手のスライ
スを導き出すことが出来る。幾つかの薄いスライス自体
に関心がなければ、投影領域の早期に対応する加算を実
施することにより、幾つかの薄いスライスを再構成する
中間工程を回避することが出来る。こうすると、計算の
負担並びに画像を記憶する負担が減少する。その結果得
られる加重関数は、データ平面の対応するシフトした形
を加算することによって導き出すことが出来る。
【００３５】この発明の種々の実施例についてこれまで
説明したところから、この発明の目的が達成されたこと
は明らかである。この発明を詳しく説明し、図面に示し
たが、これは例示の為であって例に過ぎず、この発明を
制約するものと解してはならないことを承知されたい。
例えば、ここで説明したＣＴ装置は、ｘ線源及び検出器
の両方がガントリと共に回転する「第３世代」装置であ
る。検出器がリング全体に互る不動の検出器であって、
ｘ線源だけがガントリと共に回転する「第４世代」装置
を含む他の多くのＣＴ装置も使うことが出来る。更にこ
の発明は、４つの検出器列を持つ装置だけに制限されな
い。従って、この発明の範囲は、特許請求の範囲の記載
のみによって限定されることを承知されたい。
説明したところから、この発明の目的が達成されたこと
は明らかである。この発明を詳しく説明し、図面に示し
たが、これは例示の為であって例に過ぎず、この発明を
制約するものと解してはならないことを承知されたい。
例えば、ここで説明したＣＴ装置は、ｘ線源及び検出器
の両方がガントリと共に回転する「第３世代」装置であ
る。検出器がリング全体に互る不動の検出器であって、
ｘ線源だけがガントリと共に回転する「第４世代」装置
を含む他の多くのＣＴ装置も使うことが出来る。更にこ
の発明は、４つの検出器列を持つ装置だけに制限されな
い。従って、この発明の範囲は、特許請求の範囲の記載
のみによって限定されることを承知されたい。
【図１】ＣＴ作像装置の見取図。
【図２】図１に示した装置の簡略ブロック図。
【図３】ｚ軸に沿った断面で示す４ビームｘ線の略図。
【図４】図２のＣＴ作像装置の一部分を構成する画像再
構成装置のブロック図。
構成装置のブロック図。
１２ ガントリ １８ 検出器配列 ３４ 画像再構成装置
Claims (12)
- 【請求項１】 ４列の検出器配列を持つと共に画像再構
成装置を持ち、螺旋走査で収集された投影データから物
体の断層写真画像を作る装置に於て、前記画像再構成装
置が、 （イ）ｘ線扇形ビームの各々から得られたデータに対応
する投影データ配列を作り、 （ロ）工程（イ）で発生された各々の投影データ配列に
加重関数を適用して、夫々の加重データ配列を発生し、
Ｗ1（β，γ）、Ｗ2（β，γ）、Ｗ3（β，γ）及びＷ4
（β，γ）と記す各々のデータ集合に適用すべき加重関
数が下記の式によって与えられ、 【数１】 ここで 【数２】 であり、 【数３】 ここで 【数４】 であり、 【数５】 （ハ）工程（ロ）で発生された投影データ配列から、ス
ライス画像を構成する為に使う画像データ配列を発生す
るように構成されている装置。 - 【請求項２】 画像データ配列を発生することが、各々
の加重投影データ配列にフィルタ作用及び逆投影を加え
る工程で構成されている請求項１記載の装置。 - 【請求項３】 前記フィルタ作用及び逆投影を行う前
に、同じガントリ角度であるが、異なる検出器列からの
データ配列が組合わされる請求項２記載の装置。 - 【請求項４】 第１のデータ列にある投影図が第４のデ
ータ列にある投影図から３６０゜である場合、フィルタ
作用及び逆投影にかける前に、これらの図を組合わせる
請求項３記載の装置。 - 【請求項５】 各々の投影データ配列に加重関数を適用
する前に、相次ぐスライスを再構成する為にデータが装
置のメモリに記憶される請求項１記載の装置。 - 【請求項６】 前記画像再構成装置が、更に、所望のス
ライスの輪郭がデータ配列によって支援される最も薄い
輪郭より厚手である場合、所望のスライス輪郭内にある
多数の薄いスライスを加算するように構成されている請
求項１記載の装置。 - 【請求項７】 螺旋走査で４列の検出器配列から収集さ
れた投影データから物体の断層写真画像を作る方法に於
て、 （イ）各々のｘ線扇形ビームから得られたデータに対応
する投影データ配列を作り、 （ロ）工程（イ）で発生された各々の投影データ配列に
加重関数を適用して、夫々の加重投影データ配列を発生
し、Ｗ1（β，γ）、Ｗ2（β，γ）、Ｗ3（β，γ）及
びＷ4（β，γ）と記す各々のデータ集合に適用される
加重係数が下記の式で与えられ、 【数６】 ここで 【数７】 であり、 【数８】 ここで 【数９】 であり、 【数１０】 （ハ）工程（ロ）で発生された投影データ配列から、ス
ライス画像を構成する為に使われる画像データ配列を発
生する工程を含む方法。 - 【請求項８】 画像データ配列を発生することが、各々
の加重投影データ配列にフィルタ作用及び逆投影を行う
工程で構成される請求項７記載の方法。 - 【請求項９】 フィルタ作用及び逆投影を行う前に、同
じガントリ角度であるが異なる検出器列からのデータ配
列が組合わされる請求項８記載の方法。 - 【請求項１０】 第１のデータ列の投影図が第４のデー
タ列の投影図から３６０゜であれば、フィルタ作用及び
逆投影を行う前にこれらの図を組合わせる請求項９記載
の方法。 - 【請求項１１】 各々の投影データ配列に加重関数を適
用する前に、相次ぐスライスを再構成する為に装置のメ
モリにデータを記憶する請求項８記載の方法。 - 【請求項１２】 所望のスライス輪郭がデータ配列によ
って支援される最も薄い輪郭よりも厚手である場合、所
望のスライス輪郭内にある多数の薄いスライスが加算さ
れる請求項８記載の方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/140106 | 1998-08-25 | ||
| US09/140,106 US6055291A (en) | 1998-08-25 | 1998-08-25 | Systems, methods and apparatus for reconstructing images |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2000083948A true JP2000083948A (ja) | 2000-03-28 |
Family
ID=22489785
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11237663A Withdrawn JP2000083948A (ja) | 1998-08-25 | 1999-08-25 | 断層写真画像を作る装置と方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (2) | US6055291A (ja) |
| EP (1) | EP0982680A3 (ja) |
| JP (1) | JP2000083948A (ja) |
| IL (1) | IL131377A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002153453A (ja) * | 2000-08-15 | 2002-05-28 | Ge Medical Systems Global Technology Co Llc | マルチスライス型ｃｔスキャンのハイピッチ再構成 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| US6055291A (en) * | 1998-08-25 | 2000-04-25 | General Electric Company | Systems, methods and apparatus for reconstructing images |
| US6373919B1 (en) * | 1999-08-10 | 2002-04-16 | Ge Yokogawa Medical Systems, Limited | Projection data correction method and apparatus, and radiation tomographic imaging method and apparatus |
| US6418184B1 (en) * | 2000-12-29 | 2002-07-09 | Ge Medical Systems Global Technology Company, Llc | Helical rowwise view weighting of computed tomographic images |
| JP4363834B2 (ja) * | 2002-02-19 | 2009-11-11 | 株式会社東芝 | Ｘ線コンピュータ断層撮影装置 |
| DE102012217163B4 (de) * | 2012-09-24 | 2022-06-02 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von CT-Bilddaten mit gewichteter Rückprojektion, einschließlich Recheneinheit und CT-System für dieses Verfahren |
| US9508164B2 (en) * | 2013-01-23 | 2016-11-29 | Czestochowa University Of Technology | Fast iterative image reconstruction method for 3D computed tomography |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5541970A (en) * | 1995-05-09 | 1996-07-30 | General Electric Company | Image reconstruction for a CT system implementing a four fan beam helical scan |
| US5559847A (en) * | 1995-12-06 | 1996-09-24 | General Electric Company | Systems, methods and apparatus for reconstructing images in a CT system implementing a helical scan |
| US5974110A (en) * | 1997-11-26 | 1999-10-26 | General Electric Company | Helical reconstruction algorithm |
| US6055291A (en) * | 1998-08-25 | 2000-04-25 | General Electric Company | Systems, methods and apparatus for reconstructing images |
-
1998
- 1998-08-25 US US09/140,106 patent/US6055291A/en not_active Expired - Lifetime
-
1999
- 1999-08-12 IL IL13137799A patent/IL131377A/en not_active IP Right Cessation
- 1999-08-20 EP EP99306608A patent/EP0982680A3/en not_active Withdrawn
- 1999-08-25 JP JP11237663A patent/JP2000083948A/ja not_active Withdrawn
-
2000
- 2000-03-28 US US09/537,018 patent/US6278762B1/en not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002153453A (ja) * | 2000-08-15 | 2002-05-28 | Ge Medical Systems Global Technology Co Llc | マルチスライス型ｃｔスキャンのハイピッチ再構成 |
Also Published As
| Publication number | Publication date |
|---|---|
| US6278762B1 (en) | 2001-08-21 |
| IL131377A0 (en) | 2001-01-28 |
| IL131377A (en) | 2004-07-25 |
| EP0982680A3 (en) | 2003-05-14 |
| US6055291A (en) | 2000-04-25 |
| EP0982680A2 (en) | 2000-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6141398A (en) | Protocol driven image reconstruction, display, and processing in a multislice imaging system | |
| US6421412B1 (en) | Dual cardiac CT scanner | |
| US5663995A (en) | Systems and methods for reconstructing an image in a CT system performing a cone beam helical scan | |
| US5400255A (en) | Reconstruction of images from cone beam data | |
| US6888914B2 (en) | Methods and apparatus for computing volumetric perfusion | |
| US6937690B2 (en) | Method for computed tomography of a periodically moving object to be examined, and a CT unit for carrying out this method | |
| EP1324699B1 (en) | Cardiac helical half scan reconstructions for multiple detector row ct | |
| US5559847A (en) | Systems, methods and apparatus for reconstructing images in a CT system implementing a helical scan | |
| US6493416B1 (en) | Method and apparatus for noise reduction in computed tomographic systems | |
| US6490333B1 (en) | Methods and apparatus for cone-tilted parallel sampling and reconstruction | |
| US6324247B1 (en) | Partial scan weighting for multislice CT imaging with arbitrary pitch | |
| US5606585A (en) | Methods and apparatus for multislice helical image reconstruction in a computer tomography system | |
| JP4175791B2 (ja) | 画像生成方法およびｘ線ｃｔ装置 | |
| US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
| JPH08263638A (ja) | 物体の断層写真像を作成するシステム及び方法 | |
| EP1372115B1 (en) | Methods and apparatus for reconstructing an image of an object | |
| JP2004160218A (ja) | Ｘ線コンピュータ断層装置、ｘ線コンピュータ断層装置制御方法、及びｘ線コンピュータ断層撮影プログラム | |
| US5974110A (en) | Helical reconstruction algorithm | |
| JP3917684B2 (ja) | 物体の断層写真像を作成する方法及び装置 | |
| JPH07250832A (ja) | 断層写真像を作成する方法 | |
| EP0989521B1 (en) | Fluoroscopy image reconstruction | |
| JP2000083948A (ja) | 断層写真画像を作る装置と方法 | |
| US6490334B1 (en) | Methods and apparatus for high pitch helical computed tomography image reconstruction | |
| US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
| US5764720A (en) | Methods and apparatus for simplified pre-processing of data in a computed tomography system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060815 |
|
| A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20080730 |