JPH01141646A - フーリエ変換形ｘ線ｃｔ - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明は、フーリエ変換形ＸＩＩＣＴの改良に関し、更
に詳しくは、空間分解能が高くアーティファクトが少な
い再構成画像が得られるフーリエ変換形Ｘ線ＣＴに関す
る。

（従来の技術） フーリエ変換形Ｘ線ＣＴは、各方向毎に平行ビームＸ線
に騒づく被検体の断面のプロジエクションデータを求め
、各プロジェクションデータ又は値Ｏのデータを該デー
タの両側又は片側に付加したデータを１次元フーリエ変
換し、フーリエ変換された各プロジェクションデータを
極座ＩＩＡ／直交座標変換を施した後、これに２次元逆
フーリエ変換を施して、被検体の断面の像を再構成する
ように構成されている。このようなフーリエ変換形Ｘ線
ＣＴは、フーリエ変換法の特質に由来して高速な画像再
構成が行える利点があり、又、ｉｉ像再構成の手段とし
て汎用のフーリエ変換装置を用いることができ、他の形
式のＸ線ＣＴ１例えばフィルタート・バックプロジェク
ション法によるＸａＣ丁のように専用のバックプロジェ
クタ等を持つ必要がないという利点がある。

ところで、フィルタート・バックプロジェクション法の
Ｘ線ＣＴ等においては、再構成画像の空間分解能を高め
、且つアーティファクトを低減するために、はとんどの
場合、所■クォーターオフセット法によるサンプルデー
タ収集を行う。このサンプルデータ収集法は、例えば第
３世代のＸ線ＣＴでは被検体を間に置いて互いに対向す
るＸ線源と多チャネルの検出器を被検体の周りで回転さ
せて多方向のサンプルデータを収集する場合、×６源か
ら回転の中心を通って多チヤネル検出器の中央チャネル
に照射されるＸ線が、チャネルの中心からチャネル間隔
即ちサンプル間隔の１／４だけずれた点に入射するよう
に検出器を位置決めしてデータを収集するものである。

このようなデータ収集法によって得られた多方向のプロ
ジェクションデータにおいては、プロジェクションの方
向が反対なもの同士は、Ｘ線ビームの経路が検出器のチ
ャネル間隔即ちサンプル間隔の１／２だけずれたものと
なるので、そのような関係にあるデータを組合わせると
、検出器のチャネル間隔即ちサンプル間隔が１／２に細
かくなったのと等価なプロジェクションデータが得られ
、そのようなデータに基づいて画像再構成することによ
り、空間分解能が高くアーティファクトが少ない画像を
得ることができる。

（発明が解決しようとする問題点） このクォーターオフセット法によるサンプルデータ収集
を単純にフーリエ変換形Ｘ線ＣＴに適用しても、検出器
のオフセット配置によるデータのサンプリング点のずれ
を活用した解像力が高くアーティファクトが少い良好な
画像は得られない。

本発明は上記問題点に鑑みてさなれたもので、その目的
は、クォーターオフセット法等のオフセット検出測定法
でサンプルデータを収集し、且つこのクォーターオフセ
ット法等のオフセット検出測定法の効果を生かし、空間
分解能が高くアーティファクトが少ない再構成画像を得
られるフーリエ変換形Ｘ線ＣＴを提供することにある。

（問題点を解決するための手段） 上記問題点を解決するための本発明は、放射線発生源よ
り放射され被検体を透過した放が紗亡検出器で検出し、
各方向毎にサンプルデータ群として収集し、該サンプル
データ群に前処理等の処理を施して平行ビームの各１ビ
ュー分のプロジェクションデータを求め、該データ又は
値Ｏを該データの両側又は片側に付加したデータを１次
元フーリエ変換し、周波数空間で極座標→直交座標交換
を施した後、これに２次元フーリエ逆変換を施して被検
体の断面の再構成画像を求めるフーリエ変換形ＸＩＩＣ
Ｔにおいて、前記サンプルデータの収集をクォーターオ
フセット法等のオフセット検出測定法に基づいて行い、
且つフーリエ変換後の周波数成分にオフセット検出処理
を施して得た周波数成分を新たな周波数成分として用い
ることを特徴とするものである。

（作用） 本発明に係るフーリエ変換形Ｘｌ１ｌＣＴでは、クォー
ターオフセット法等のオフセット検出測定法を用いてサ
ンプルデータの収集を行い、該サンプルデータに基づき
得られたプロジェクションデータに対しては、フーリエ
変換後の周波数空間でオフセラトル検出処理を行う。

（実施例） 以下、図面を用いて本発明の詳細な説明する。

第１図は本発明の一実施例を示すブロック図である。図
において、Ｇはガントリ、ＴＡは被検体Ｂを載せるテー
ブルで、これらはテーブル・ガン１−り制御装置！ＴＧ
Ｃによって位置制御が行われる。

ＸＧＣはＸ線発生ｉＩ＋ＩＪｌｌＩ］装置で、高圧・Ｘ
線管Ｉ！Ｉ１１罪部ＸＲを介してＸＳｔ管〈Ｘ線源）Ｘ
に接続され、Ｘ線発生の制御を行う。Ｓは断層像の再構
成領域ＰＡ（被検体Ｂ）を透過したＸ線を受ける多チャ
ネルの検出器である。本発叫ではクォーターオフセット
法に基づきサンプルデータの収集を行うので、Ｘ線管Ｘ
から再構成領域ＰＡの中心Ｏを通って多チヤネル検出器
Ｓの中央チャネルに照射されるＸＩＱが、チャネルの中
心からチャネル間隔即ちサンプル間隔の１／４だけずれ
た点に入射するように、検出器Ｓの位置決めがなされて
いる。ＤＡＳは検出器Ｓの出力を収集するデータ収集装
置で、検出器Ｓの出力を増幅、積分し、Ａ／Ｄ変換した
後、データ記憶装［ＤＳｌに格納する。

ＰＰＲは前処］！ｌ！ｖｉ＊−ｃ’、データ収集装置１
ＤＡｓで収集された被検体ＢのＸ線透過データに対し、
データ収集装ｆｉ！ＦＤＡＳのオフセット補正、対数変
換、Ｘ線ビーム硬化補正等々の各種補正演算を施し、Ｘ
ＩＪのプロジェクションデータを求める部分である。こ
の結果はデータ記憶装置ＤＳ２に格納される。

ＦＦＴはフーリエ変換装置で、平行ビームの１ビュー分
のプロジェクションデータ毎に離散フーリエ変換（ＤＦ
Ｔ）を施し、結果の周波数成分をデータ記憶１ＨｆｆＤ
ｓ２に格納する。ここでの演算式は、Ｎ個の実数データ
ａＯ＊　ａｌ　＊　Ｃ２＊　”’　。

ａＮ−ｔ　　（Ｎ−２’　　、γは整数）に対し、次の
ものを用いる（ｊ−Ｅ］）。

Ａ　　（ｎ　　）−Ａｒ　　（ｎ　　＞＋ｊ　　−Ａｉ
　　（ｎ　　＞・・・　（２） 但し　ｎ−０，１，２，・・・　ｉ Ｒ（４−）　−Ｒ（０）　。

１　　（−”−）　−Ｉ　　（０）　。

Ｙ　　（８）　　−Ｙ　　（０） △　（ｎ＋−） ＝Ａｒ　　（−Ｔ−ｎ　）　−ｊ　　−Ａｉ　　（２ｎ
　）・・・　（２）′ ９３．　　旦 ここで、ｎ−１，２，、−１ 式（１）の計算は所謂高速フーリエ変換（Ｆ　ａｓｔＦ
Ｏｕｒｉ６ｒＴｒａｎｓｒｏｒｍ）により高速に演算す
る。

ａＯ＊　８１　＊　”’＋　ａＮ−１のフーリエ変換の
結果Ａｏ。

八！、・・・、ＡＮ−１は式〈２）と（２）′に示す通
りとなる。Ａｒ（ｎ）とＡｔ（ｎ＞はＡ（ｎ）の各々実
数部と虚数部に対応する。式（１）は所謂圧縮フーリエ
変換の式で、Ｎ個の実データをＮ／２組の複素データと
して高速演算を計り、Ｒ柊的にアンバック処理（２−１
）、（２−２）により、求めるフーリエ変換結果Ａ　＜
ｎ　）　　（ｎ　−０，１゜・・＋　、　Ｊｉ、　）を
得る。

０ＦＰはオフセット検出処理装置で、Ｘ線源（Ｘ線発生
焦点）ｘ、検出器群Ｓの中心（中央に位置する検出器の
中心）Ｃ２再構成酒域ＰＡの中心Ｏに対し、直線ＸＣを
直線Ｘｏから故意にずらしてクォーターオフセット法に
より収集し且つ前述の処理が施され、データ記憶装置１
ｆｆＤＳ２に格納されたデータに対し、対応する処理を
施す装置である。ここでの処理の演算式は次のものであ
る。

Ｆ　（ｎ　）−Ｆｒ　　（ｎ　）＋ｊ　−Ｆｉ　　（ｎ
　）−ノγ −Ａ　（ｎ　）・ｅ　　　　　　　・・・（３）Ｆｒ　
　（ｎ　）　−Ａｒ　　（ｎ　）　−ｃｏｓ　７＋Ａｉ
　　（ｎ　）　−５ｉｎ　７　−（３−１）Ｆｉ（ｎ　
）　＝Ａｉ　　（ｎ　）−ｃｏｓγ−Ａｒ　　（ｎ　）
　−５ｉｎ　７−（３−２）但し　ｎ　−０，１，２，
・・・、Ｊｉ−Ｆ（ｎ十工） −Ｆｒ　（８−ｎ　）　−ｊ　−Ｆｉ　（子−ｎ　＞・
・・（３）′ Ｆｒ　　（ｎ十亙）−Ｆｒ（−−ｎ） ・・・（３−１）’ ＦＨ（ｎ　＋ＪＬ）　−−ＦＨ（＋−ｎ）・・・　（３
−２＞’ 但し　ｎ−１，２，・・・、ニー１ ４πｎ＆ γ−−１．−　　　　　　　　　　　　・・・（４）Ｎ
は入力１２素周波数成分の対の数（ＦＦＴの出力周波数
成分の数）であり、ｄはＰＡの中心Ｏでの平行ビームの
サンプル間隔（■）、δはオフセット検出でのオフセッ
ト１１（線）で、基準座標に対して右方向（正方向）の
オフセットを十とする。

クォーターオフセット法では、δ−−ｄ／４又はｄ／４
である。Ａ　（ｎ　）は入力周波数成分で、Ａｒ（ｎ）
は実数部、Ａｉ（ｎ）は虚数部である。

Ｆ　（ｎ　）は出力の周波数成分で、Ｆｒ（ｎ）はその
実数部、Ｆｉ（ｎ）は虚数部である。フーリエ変換され
オフセット処理を施された周波数成分は、周波数空間（
ξ−η平面）上で第２図（ｂ）のように直線ω上に並ん
でいる。これと１８００ずれた周波数成分は丁度光の成
分と逆の向きに並んでいる。元の周波数成分Ｐ（ω０．
θＩ）。

１８０°ずれた方向の周波数成分Ｑ（ωｎ、θｌ＋π）
に対し、合成された周波数成分Ｆ（ωｎ。

θ−）を次のように演算により求める。

Ｆｒ（ωｎ、θｍ　）　−＋（Ｐｒ　（ωｎ　、６ｍ）
＋Ｑｒ　　（ωｎ、θ−＋π）） ・・−（５−１） Ｆｉ（ωｎ、θｌ）−子（ｐｉ（ωｎ、θ層）−Ｑｉ　
　（ωｎ、θ霞＋π）） ・・・（５−２） １：　　（ωｎ　　、　　０膳　）　　−１Ｆｒ　　　
（ωｎ、　　６ｍ）＋ｊ−Ｆｉ（ωｎ、θ閣　） ・・・（５） ここで、ｎ−０，１，２，・・・、−！−Ｆｒ　（ω、
＋ｇ、θ偏）−Ｆｒ（ωト７．θ履）・・・＜５−１）
’ Ｆｉ　（ω、ｌ＋トθ−）−−Ｆｉ（ωトｎ、θ糟）・
・・（５−２＞’ Ｆ（ω７＋＋五、θ鵬　）−Ｆｒ（ω、÷凡、θｍ　）
＋ｊ−Ｆｉ（ω７＋号、θ−） ・・・（５）′ ここで、ｌ’ｌ　−０，１、２，−、ＪＬ−’１式（３
）、（３）’　、（５）、（５）’は、高解像性、低偽
像性、高画質性（低ＳＤ、低ノイズ等）等々から重要で
ある。ＯＦＰの出力はデータ記憶＠［ＤＳ３に格納され
る。ＯＦＦは本発明の最も特徴的な部分である。

ＰＲＣは極座標／直交座標変換装置である。実空間での
プロジェクションデータと再構成領域の関係を第２図（
ａ）に示す。Ｘ線源はＸ側にあり、ＯＸはＸ軸止方向と
θ−十丁の角をなす。ｘｙ座標系とθ腫をなす直角座標
系Ｕｖに対し、平行ビームのプロジェクションデータｂ
Ｏ＊　ｂ　！　＋・・・。

ｂＭ−１はＵ軸上に図のように配置される。実際のプロ
ジェクションデータｂｂ　　は、ｘｍの方向が直線に１
”Ｋｚに向い、直線ＫＩＫＺ上での放射線吸収係数の線
積分であるが、これをＵ軸上の点Ｕｋにあるとみなす。

サンプルデータｂＯ＊ｂｌ＊・・・。

ｂｏの両端に適当な数の値Ｏをつけ且つ配列変更したデ
ータａＯ＋ａｌｅ・・・−’Ｎ−１に対して離散フーリ
エ変換を施し、オフセット検出に対する処理を施した周
波数成分データＦｋ　　（ｋ−０，１，２゜・・・、Ｎ
−１）は、周波数空間上（ξ−η平面上）で第２図（ｂ
）のようにξ軸とθ−の角をなすω軸上の離散点に対応
する。極座標／直交座標変換は、この極座標ω−θの各
離散データより近似計算により直交座標ξ−ηの各離散
点のデータを求める変換処理である。（ｂ）図のξ−η
座標上の点Ｐの近傍でω−θ座標上の４点はＡ、Ｂ、Ａ
’　。

Ｂ′となる。点Ｐ、Ａの座標を各々Ｐ（ξに、ηｎ）ｓ
Ｐ（ｋ、Ａ’）、Ａ（ωｎ、　　６ｍ　＞　！！Ａ　（
ｎ。

■）とし、点Ｐ、Ａの周波数成分を各々Ｇ（ξｋ。

７７７　）ｍＧ　（ｋ　、　／）、　Ｆ　（ωｎ　、θ
ｌ）！Ｆ（０，■）とすると、本変換の演算式として次
のものをあげることができる。

Ｇｒ　　（ξｋ　、　７７４　）　ａ　Ｆ　ｒ’　（（
１）　ｎ　＊　θ霞）十トエ（（ω−ωｎ　）　・立＋
（θ−θ−）　・立）２ｐ、、　Ｐ／　　　　　　　　
　　　　　　　　　　６θＯψ Ｇｉ　　（ξに＋　　”ｆ）ｘ　　）＝Ｆｉ　　（ωｎ
、　　θ請　）十ト和（ω−ωｎ）・斎＋（θ−θｍ）
・☆）２Ｇ（ξに、　　７７Ｊ　　）＝Ｇｒ　　（ξに
＋’７ｎ）＋ｊ　　−Ｇｉ　　（ξｋ　、η１　）・・
・　（６） Ｇｒ　　（ξに、７７＊）、Ｇｉ（ξｋ　、　”ｆ）ｓ
　）は各々Ｇ（ξに、η１）の実数部、虚数部（点Ｐに
おけるＧ（ξ、η）の実数部、虚数部）であり、Ｆｒ（
ωｎ、θ−）、Ｆｉ（ωｎ、θ醜）は各々Ｆ（ωｎ、θ
ｍ）の実数部、虚数部（点ＡにおけるＦ（ω、θ）の実
数部、虚数部）である。ｐ、しは整数で、Ｌは近似の正
確さにより例えばＬ−３意味で（６−２＞式についても
同様である。

回部微分し、更に、θでｑ回偏微分した関数のω−ωｎ
、θ−θ−における値である）、本近似式は、前述のフ
ーリエ変換前の０ｆｉｌｌの数と共に、近似の精度に重
要な影響を有するので、式の次数（Ｌ）をどう選ぶか、
高速で高精度の近似式として何を選ぶか等慎重に検討さ
れる。この変換結果はデータ記憶装置ＤＳ４に格納され
る。

ＩＦＦＴは２次元フーリエ逆変換装置で、第２図（ｂ）
のξ−η座標（直交座標）の各点（各格子点）に対応す
る周波数成分Ｇ（ξ、η）から、離散２次元フーリエ逆
変換により実空間イメージデータを求める部分である。

演算式は以下の通りＵ　　（Ｘ　ｎ　＊　Ｖｍ　　） 一μ（Ｘｎ、Ｖｍ）＋ｊ　・ν（Ｘ、ｎ、Ｖｒｎ）・・
・（７−１） ｎ、Ｉｌ−一工　−−Ｌ−＋　１、−、０　、　１　、
　・。

２　′　　　　２ μ（Ｘｎ、％））が求めるものである。即ち、２次元フ
ーリエ逆変換（７）の実数部が求めるものである。画像
データとしては更に次式のＣ（ＸＴ＋。

ん　）を使用する（ａ　、　ｂは定数）。

Ｃ（Ｘｎ、　Ｖｍ　　）＝ａ　・μ（Ｘｔ＋、Ｖｍ）＋
ｂ・・・（８） （ア）式は離散高速フーリエ変換として高速に演算する
。

ＩＭは映像データ記憶装置で、２次元フーリエ逆変換装
置ＩＦＦＴにより得られた映像データを格納する。この
映鍮データは画像表示装置ＧＯＣにて表示され、像写真
Ｗｉ像装置ＭＦＣにて写真撮影されるようになっている
。尚、前述の各装置は庵影制ｔＩＩ装置ＳＣＣにて、制
御されるように構成されている。

次に上記実施例の動作を第３図を参照しながら説明する
。この例では、各方向の１ビューデータは、収集方向の
位置の制御がなされ、複数方向の平行ビームの成分とな
るデータの集合であり、見かけ上フアン状に発散するデ
ータの集合である（通常のファンビームデータと異なる
）。

被検体Ｂの所要の断面は再構成領１ｉ！ＰＡとしてテー
ブルＴＡとガン１−リＧの傾斜の適宜な位Ｉｔ選択によ
り得られる。この断面ＰＡに対して各方向からＸＩＤの
発生と検出器Ｓによる被検体Ｂの透過Ｘｔａの収集をデ
ータ収集装［ＤＡＳにより行う。

データ収集装置！ｔＤＡｓでは入力の積分、増幅、Ａ／
Ｄ変換等がなされ、ファン状の１ビュー分のデータとし
て１ビュー毎にデータ記憶装＠ＤＳ１に転送格納する。

）７ンビームに対して前処理を施し、結果をデータ記憶
装置ｆＤＳ１に格納する（ステップ１）。

ファンビームの並べ換えにより平行ビームを得るが、１
ビュー分の平行ビームが揃うと（ステップ２）、前記の
前処理を施しくステップ３）、結果をデータ記憶装置Ｄ
Ｓ２に格納し、全ビュー・１／２＋１ビューが完了する
と（ステップ４）、次のステップに移る。

１ビュー分の平行ビームデータは第２図（ａ　）のよう
な位置でサンプルされたＭ個（ＭＬｔ奇数）のデータｂ
Ｏ＋　ｂ　！　＋　ｂ２　＊　””　、Ｉ））４−１で
ある。この両端に°０値データをつけ、 Ｎ／２個 Ｎ／２個 夕を基にして、これを配列変更したデータ列（ａ　ｉ　
　；ｉ　−Ｑ、　１．−、　Ｎ−１）ａ　Ｏ−ｂＬ−ｉ
、　ａ　１−ｂ旦址、−、ａｋｌ−ｂ、−１、０゜０、
−−１ＨＨ＋ＨＨ，ａＮ−Ｑ−ｂ　ｏ　、−、ａＮ−２
−ｂ早ｔａＮ、　−ｂＬ１ に対し、前記の離散フーリエ変換ＤＦＴ、即ち式％式％
（２） ′を施しくステップ５）、結果をデータ記憶！！装置Ｄ
Ｓ２へ格納する。

ステップ５の演算は原理的に次式に相当する。

ｂ（ｒ　）は実測値（離散値）ｂｉと次の関係をもつ。

以下余白 ・・・　（１１−１） この離散フーリエ変換の結果に対して、オフセット検出
処理、即ち（３）、（３〜１）、（３−２）、　＜３）
’　、　（３−１）’　、　＜３−２）’式を施した後
（ステップ６）、対向ビューデータとの間で周波数成分
の合成、即ち（５）、（５−１）、（５−２）、（５）
’　、（５−１＞’　、（５−２）′式を施すくステッ
プ７）。

全ビューの平行ビームに対しステップ１〜ステツプ７の
処理が終了したら（ステップ８）、極座標／直交座標変
換を式（６）、（６−１）、（６−２）に従って行う（
ステップ９）。

次に得られたＬ　（Ｃｏｍｐｉｅｘ）　Ｘ　Ｌの直交座
標データより１式（７）、＜７−１）による２次元離散
フーリエ逆変換を行う（ステップ１０）。そして、この
ステップ１０で得られた結果の実数部を使用し、イメー
ジ処理（式（８）等各種）を施すくステップ１１）。ス
キャンと再構成処理の統轄制御は撮影制御装置１ｓｅｃ
によりなされる。各装置は、バイブライン構造等により
並列動作が可能で並列処理を行う。

ステップ１０の演算は原理的に次式に相当する。

Ｇ（ξ、η）は離散データＧ（ξｐ、ηｑ）と次の間係
にある。

Ｇ（ξ、η） 尚、本発明は上記の実施例に限るものではなく、例えば
、以下に示すような種々の変形が可能である。

■オフセット検出処理装置０ＦＰの後（極座標／直交座
標変換装置ＰＲＣの前）にフィルタ装置ＦＩＬＴを設け
たもの。

フィルタ処理の例としては、例えばハニングフィルタ（
）（ａｎｎｉｎｇ　　Ｆｉｌｔｅｒ）等がある。

■２次元フーリエ逆変換＠ＩＩＦＦＴの後に侵処理装置
ＰＯ８Ｔを設けたもの。

後処Ｉ！！！装けＰＯ８Ｔでは、式（８）を含む各種の
処理を行う。この場合、２次元フーリエ逆変換装置！Ｉ
ＦＦＴは式（７）、（７−１）の純２次元フーリエ逆変
換を行う。

■オフセット検出処理装置０ＦＰの後に周波数成分生成
装ｆｆｆＧＦＲＱを設けたもの。

周波数成分生成装置！ＧＦＲＱでは式（５）。

（５−１）、（５−２）、（５）’　。

（５−１）’　、　　（５−２＞’　（７）ｖ４＄１［
処理を行う。

この場合、オフセット検出処理装置ＯＦＦは式％式％） ＜３−１）’　、（３−２）’の演算処理を行う。

■装置の分刈 （ｉ）前処理装置ｆＰＰＰを２つ以上の前処理装置に分
離して並列動作を可能としたもの （ｉｉ）撮影線御装置１ｆＳＣＣをマイクロプロセッサ
。

マイクロプログラムメモリ、デコーダ等で構成したもの
。

■装置の統合 （ｉ　）フーリエ変換装置ＦＦＴと２次元フーリエ逆変
換装置［ＦＦＴを同一の装置とし、時分割使用するもの
。

（ｉｉ）複数の装置を統合して構成したもの。

■記憶装置の分割又は統合 （ｉ）１つのデータ記憶装置（例えばＤＳ２）を複数に
分割したもの。

（ｉｉ）　２つ以上のデータ記憶装置を統合したもの。

■バス、制御ライン、信号ライン等の分離又は統合。

■外部装置の付加 ・磁気テープ装ＷＭＴ、フロッピーディスク装置ＦＤＤ
、専用オペコン等の付加。

・データ処理用に汎用の情報処理装置ＣＰＵや専用のア
レイプロセッサＡＰ等を設けたもの。

■同一の演算処理（例えば式（３）、（３−１）。

（３−２＞、（３）’　、（３−１＞’　、（３−２）
’　、（５）、（５−１＞、（５−２）。

（５）’　、（５−１＞’　、（５−２＞’等）を別の
構成（例えば情報処理装置ｃｐｕ、アレイプロセッサＡ
Ｐ等の装置とソフトウェア）で実現したもの。

（発明の効果） 以上説明したように、本発明では、クォーターオフセッ
ト法等のオフセット検出測定法によりデータを収集する
と共に、１次元フーリエ変換後の周波数成分にオフセッ
ト検出処理を行っているため、クォーターオフセット法
等のオフセット検出測定法の効果をフーリエ変換形ＸＩ
ＣＴにも生かすことができる。即ち、高速フーリエ変換
による高速な処理性を保持しつつ、空間分解能が高く、
アーティアクトが少ないフーリエ変換形ＸＩＩＣＴを実
現できる。又、本発明によれば、逆投影装置のような専
用装置は不要で、汎用のフーリエ演算装置を主体とした
シンプルな構成をもつ、経済性の高いＸＩ！１ｌｃＴを
実現可能である。更に、Ｆ（ωｎ、θ信）が原点に関し
て複素共役対称である性質によりＦ〈ω、θ）のための
データメモリの容量（サイズ）を通常の場合の１／２に
成し得る等多くの長所をもつ。

【図面の簡単な説明】

第１図は本発明の一実施例の構成を示すブロック図、第
２図は実空間での射影や周波数空間での周波数成分に関
する説明図、第３図は第１図の実施例の動作説明図であ
る。 Ｇ・・・ガントリ　　　　ＴＡ・・・テーブルＢ・・・
被検体　　　　　ＰＡ・・・再構成領域Ｏ・・・再構成
領域の中心 ＴＧＣ・・・テーブル・ガントリ制御′ＶＡａＸＧＣ・
・・Ｘ？！Ａ発生ＩＩ線ｔｌ装置Ｘ・・・Ｘ１ｓ管（Ｘ
線発生焦点） ＸＲ・・・高圧、Ｘ線管υＵ部 Ｓ・・・検出器群　　　　ＤＡＳ・・・データ収集装置
ＰＰＰ・・・前処理装置 ＦＦＴ・・・フーリエ変換ｖ装置 ＯＦＦ・・・オフセット検出処理装置 ＰＲＣ・・・極座標／直交座標変換Ｖ装置ＩＦＦＴ・・
・２次元フーリエ逆変換装置Ｏ８Ｉ〜ＤＳ４・・・デー
タ記憶装置 ＩＭ・・・映像データ記憶装置 ＧＤＣ・・・画像表示装置 ＭＦＣ・・・像写真搬影装置 ＳＣＣ・・・搬影制御装置

Claims (2)

【特許請求の範囲】
1. （１）放射線発生源より放射され被検体を透過した放射
   線を検出器で検出し、各方向毎にサンプルデータ群とし
   て収集し、該サンプルデータ群に前処理等の処理を施し
   て平行ビームの各１ビュー分のプロジェクションデータ
   を求め、該データ又は値０を該データの両側又は片側に
   付加したデータを１次元フーリエ変換し、周波数空間で
   極座標→直交座標交換を施した後、これに２次元フーリ
   エ逆変換を施して被検体の断面の再構成画像を求めるフ
   ーリエ変換形Ｘ線ＣＴにおいて、前記サンプルデータの
   収集をクォーターオフセット法等のオフセット検出測定
   法に基づいて行い、且つフーリエ変換後の周波数成分に
   オフセット検出処理を施して得た周波数成分を新たな周
   波数成分として用いることを特徴とするフーリエ変換形
   Ｘ線ＣＴ。
2. （２）各方向のサンプルデータ群である各方向の１ビュ
   ーのデータ群が、収集方向の位置制御がなされ１又は複
   数の方向の平行ビームの成分の集合であり、並べ換えに
   より各方向の平行ビームになし得るデータ群であるか、
   若しくは収集方向の位置制御がなされ純ファンビームと
   して放射線発生源から発散する等の形態をもち補間演算
   等の計算により平行ビームになし得るデータ群である特
   許請求の範囲第１項記載のフーリエ変換形Ｘ線ＣＴ。