(57)【要約】本発明は、検査ゾーンが円錐形の放射ビームによってらせん状の軌道に沿って照射されるコンピュータ断層撮影方法に関する。検査ゾーンによって透過される放射は、検出器ユニットによって測定され、それから、検査ゾーンにおける吸収分布が、近似なしで再構成される。再構成は、異なる投影の平行な放射線の測定値の微分、Ｋラインに沿ったこれらの値の積分、これらの値の重みづけ及び逆投影を含む。
2006505337.tif 000002
【特許請求の範囲】
【請求項１】
放射線源を使用して、検査ゾーン又は前記検査ゾーンに存在するオブジェクトを通過する円錐形の放射ビームを生成するステップと、
前記放射線源と、前記検査ゾーン又は前記オブジェクトとの間に相対運動を生じさせるステップであって、該相対運動が、回転軸を中心とする回転と、前記回転軸に平行な移動とを含み、らせんとして形作られる、ステップと、
前記相対運動の間、検出器ユニットによって、前記検査ゾーンの向こう側で、前記放射ビームの強度に依存する測定値を取得するステップと、
前記測定値から前記検査ゾーンのＣＴ画像を再構成するステップと、
を含み、前記再構成において、
前記放射線源の角度位置に関して、異なる放射線源位置による平行な放射線の測定値の偏導関数を決定するステップと、
Ｋラインに沿って前記微分された測定値の重み付き積分を行うステップと、
すべての測定値に、当該測定値に関する放射線の円錐角のコサインに対応する重み係数を乗じるステップと、
すべての測定値に、当該測定値に関するビームの扇角のコサインの逆数の値に対応する重み係数を乗じるステップと、
前記測定値の逆投影によって、各オブジェクトポイントの吸収を再構成するステップと、
を含む正確な３次元逆投影が実行される、コンピュータ断層撮影方法。
【請求項２】
前記再構成するステップにおいて、前記逆投影の前に、前記測定値のリビニングが、各グループが複数の平面を含む複数のグループを形成するように実施され、前記複数の平面は、互いに及び前記回転軸に平行に延在し、前記平面内には、個々の扇形ビームが位置する、請求項１に記載のコンピュータ断層撮影方法。
【請求項３】
Ｋラインに沿った前記測定値の前記積分は、
それぞれの前記放射線源位置及び前記検査ゾーン内の再構成されるべきそれぞれの位置について、Ｋ平面を決定するステップと、
前記Ｋ平面と、前記検出器ユニットの検出器表面との間の交線であるＫラインを決定するステップと、
各Ｋライン上の測定値に、角度Ｋのサインの逆数の値に対応する重み係数を乗じるステップと、
前記Ｋラインに沿って前記測定値を積分するステップと、
を含む、請求項１に記載のコンピュータ断層撮影方法。
【請求項４】
請求項１に記載の方法を実行するためのコンピュータ断層撮影装置であって、
放射線源と、検査ゾーン又は前記検査ゾーンに存在するオブジェクトを通過する放射ビームを生成するために前記検査ゾーンと前記放射線源との間に位置する絞り構造と、
前記放射線源に結合される検出器ユニットと、
前記検査ゾーンに存在するオブジェクト及び前記放射線源を、回転軸を中心として及び／又は前記回転軸に平行に、相対的に移動させる役割を果たす駆動装置と、
前記検出器ユニットによって取得される測定値から、前記検査ゾーン内の吸収の空間分布を再構成する再構成ユニットと、
請求項１に開示されるステップに従って、前記放射線源、前記検出器ユニット、前記駆動装置及び前記再構成ユニットを制御する制御ユニットと、
を有するコンピュータ断層撮影装置。
【請求項５】
請求項１に開示されるステップを実行するように、コンピュータ断層撮影装置の放射線源、絞り構造、検出器ユニット、駆動装置及び再構成ユニットを制御する制御ユニット用のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、検査ゾーンが円錐形のＸ線ビームによってらせん状の軌道に沿って照射される、コンピュータ断層撮影方法に関する。更に、本発明は、コンピュータ断層撮影装置及びコンピュータ断層撮影装置を制御するためのコンピュータプログラムに関する。
【背景技術】
【０００２】
上述した種類の知られている方法では、近似を利用して、検査ゾーンにおける放射の吸収又は減衰の空間内の変化が、検出器ユニットによって取得される測定値から再構成されることができる。しかしながら、このような近似は、再構成される画像にアーチファクトを生じさせ、放射ビームのアパーチャの角度が、回転軸の方向において大きくなるにつれて、前記アーチファクトは、より顕著になる（Artifact Analysis of Approximative Cone-Beam CT Algorithms, Medical Physics, Vol. 29, pp. 51-64, 2002）。
【０００３】
知られている正確な方法は、通常、逆ラドン変換（Radon inversion）に基づく。これらの方法は、大量の計算作業を必要とし、再構成される画像に離散化誤差を生じさせる。
【０００４】
更に、フィルタ補正逆投影法（filtered back projection）を利用する正確な方法が、「Analysis of an Exact Inversion Algorithm for Spiral Cone-Beam CT」（Physics Medicine and Biology, Vol. 47, pp. 2583-2597）(以下、文献Ｅ１と呼ぶ)から知られている。この方法は、同様に、大量の計算作業を必要とし、従って長い再構成時間を生じさせる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
従って、本発明の目的は、検査ゾーンにおける吸収分布の、より高速であり、正確な再構成を実現する方法を提供することである。
【課題を解決するための手段】
【０００６】
本発明の方法に関して、この目的は、
−放射線源を使用して、検査ゾーン又は検査ゾーンに存在するオブジェクトを通過する（traverse、横断する）円錐形の放射ビームを生成するステップと、
−放射線源と、検査ゾーン又はオブジェクトとの間に相対運動を生じさせるステップであって、該相対運動が、回転軸を中心とする回転及び回転軸に平行な移動を含み、らせんとして形作られる、ステップと、
−相対運動の間、検出器ユニットによって、検査ゾーンの向こう側で、放射ビームの強度に依存する測定値を取得するステップと、
−測定値から検査ゾーンのＣＴ画像を再構成するステップと、
を含み、前記再構成において、
−放射線源の角度位置に関して、異なる放射線源位置による平行な放射線の測定値の偏導関数を決定するステップと、
−Ｋラインに沿って微分された測定値の重み付き積分を行うステップと、
−すべての測定値に、当該測定値に関するビームの円錐角のコサインに対応する重み係数を乗じるステップと、
−すべての測定値に、当該測定値に関するビームの扇角（fan angle）のコサインの逆数値に対応する重み係数を乗じるステップと、
−測定値の逆投影によって、それぞれのオブジェクトポイントの吸収を再構成するステップと、
を含む正確な３次元逆投影が実施される、コンピュータ断層撮影方法によって達成される。
【０００７】
文献Ｅ１から知られている方法によれば、逆投影の前に、測定値は、検査ゾーンにおける再構成されるべきオブジェクトポイントの位置に依存する重み係数を乗じられなければならない。これとは対照的に、本発明によれば、測定値は、逆投影の前に、検出器ユニット上の測定値の位置に依存する重み係数を乗じられる。再構成されるべきオブジェクトポイントの数は、概して、検出器素子の数よりずっと少ないので、正確な再構成が可能であるとともに、かなり少ない量の計算作業のみが必要とされる。更に、文献Ｅ１から知られている方法と異なり、本発明による逆投影の積分間隔は、オブジェクトポイントに依存せず、従って、逆投影の間、それぞれのオブジェクトポイントについて積分間隔を決定する必要がなく、要求される計算作業の量の更なる低減をもたらす。
【０００８】
請求項２は、他の方法と比較してより少ない計算作業量を含むとともに、良好な画像品質をもたらす好適な再構成方法を記述する。
【０００９】
請求項３は、Ｋラインに沿った積分の好適なバージョンを開示し、特に高い画像品質をもたらす重み係数による乗算を開示する。
【００１０】
本発明の方法を実施するためのコンピュータ断層撮影装置が、請求項４に開示されている。請求項５は、請求項４に記載のコンピュータ断層撮影装置を制御するためのコンピュータプログラムを規定する。
【発明を実施するための最良の形態】
【００１１】
以下、本発明について図面を参照して詳しく説明する。
【００１２】
図１に示されるコンピュータ断層撮影装置は、図１に示される座標系のｚ方向に平行に延びる回転軸１４を中心に回転可能なガントリ１を有する。このために、ガントリ１は、好適には一定であるが調整可能な角速度で、モータ２によって駆動される。例えばＸ線源のような放射線源Ｓが、ガントリ１に取り付けられる。放射線源は、コリメータ構造３を備える。コリメータ構造３は、放射線源Ｓによって生成される放射から、円錐形の放射ビーム４、すなわち、ｚ方向及びｚ方向に垂直な方向（すなわち回転軸に対して垂直な平面）においてゼロ以外の有限なディメンジョンを有する放射ビームを形成する。
【００１３】
放射ビーム４は、例えば患者テーブル上の患者（両方とも不図示）又は技術的なオブジェクトのようなオブジェクトが位置しうる、柱状の検査ゾーン１３を通過する。検査ゾーン１３を通過したのち、放射ビーム４は、ガントリ１に取り付けられた検出器ユニット１６に入射する。検出器ユニット１６は、本実施例において、マトリックスの形の行及び列に配置された複数の検出器素子を含む検出器表面を有する。検出器の列は、回転軸１４に平行に延在する。検出器の行は、回転軸に対して垂直に延在する平面に位置し、好適には放射線源Ｓを中心とする円弧上に位置する。しかしながら、検出器の行はまた、別の形状を有することもでき、例えば回転軸１４を中心とする円弧の形でもよく、又は直線状であってもよい。一般に、検出器の行は、検出器の列（例えば１６）より多くの検出器素子（例えば１０００）を含む。放射ビーム４が当たるそれぞれの検出器素子は、放射線源のそれぞれの位置において、放射ビーム４の放射線の測定値を生成する。
【００１４】
参照符号αｍａｘによって示される放射ビーム４のアパーチャの角度は、検査されるべきオブジェクトが測定値の取得中に位置するオブジェクトシリンダの直径を決定する。アパーチャの角度は、放射線源Ｓ及び回転軸１４によって規定される平面に対し、回転軸に垂直な平面において、放射ビーム４の端に位置する放射線によって囲まれる角度として規定される。検査ゾーン１３又はオブジェクト若しくは患者テーブルは、モータ５によって、回転軸１４又はｚ軸に平行に移動されることができる。しかしながら、等価なものとして、ガントリが、この方向において移動されることもできる。
【００１５】
患者の代わりに、技術的なオブジェクトが関係する場合、検査中、オブジェクトが回転されることができ、放射線源Ｓ及び検出器ユニット１６は静止したままである。
【００１６】
モータ２及び５が同時に動作するとき、放射線源Ｓ及び検出器ユニット１６は、検査ゾーン１３に対してらせん状の軌道を示す。しかしながら、回転軸１４の方向への移動用のモータ５が静止したままであり、モータ２がガントリを回転させる場合、放射線源及び検出器ユニット１６について、円形の軌道が、検査ゾーン１３に対して得られる。以下、らせん状の軌道のみが考慮される。
【００１７】
検出器ユニット１６によって取得される測定値は、画像処理コンピュータ１０に与えられる。画像処理コンピュータ１０は、例えばワイヤレスデータ伝送（図示せず）を介して、検出器ユニット１６に接続される。画像処理コンピュータ１０は、検査ゾーン１３における吸収分布を再構成し、それを、例えばモニタ１１上に再生する。２つのモータ２及び５、画像処理コンピュータ１０並びに放射線源Ｓ及び検出器ユニット１６から画像処理コンピュータ１０への測定値の伝送は、制御ユニット７によって制御される。
【００１８】
他の実施例において、取得された測定値は、まず、再構成のために、１つ又は複数の再構成コンピュータに与えられることができ、それらのコンピュータは、例えば光ファイバケーブルを介して、再構成されたデータを画像処理コンピュータに与える。
【００１９】
図２は、図１に示されるコンピュータ断層撮影装置によって実行されることができる測定及び再構成方法の１つのバージョンを説明するフローチャートを示している。
【００２０】
ステップ１０１における初期化ののち、ガントリは、本実施例において一定である角速度で回転する。しかしながら、角速度は、例えば時間又は放射線源位置に依存して、変化してもよい。ステップ１０３において、検査ゾーン又はオブジェクト若しくは患者テーブルが、回転軸に平行に移動され、放射線源Ｓの放射のスイッチが入れられ、それによって、検出器ユニット１６は、複数の角度位置からの放射を検出することができる。
【００２１】
次のステップを理解するために、「Analysis of an Exact Inversion Algorithm for Spiral Cone-Beam CT（Physics Medicine and Biology, Vol. 47, pp. 2583-2597）から次の方程式を参照する。
【数１】
2006505337.tif 000003
【００２２】
この方程式は、測定値の逆投影による吸収の正確な再構成を記述する。上式において、ｆ（ｘ）は、検査ゾーン内の位置ｘにおける空間的な吸収分布を示し、ＩＰＩ（ｘ）は、ＰＩライン３１によって囲まれるらせんの一部を表す。
【００２３】
検査ゾーン内の位置ｘにおけるオブジェクトポイント３５のＰＩライン３１が、図３及び図４に示されており、以下に詳しく記述される。放射線源は、らせん状のパス１７をたどって、オブジェクトポイント３５の周囲を、検査ゾーンと相対的に移動する。ＰＩライン３１は、２つの位置でらせんと交差するとともに、オブジェクトポイント３５と交差するラインであり、このラインによって囲まれるらせん状のセグメントＩＰＩ（ｘ）は、２πより小さい角度をカバーする。
【００２４】
更に、方程式（１）において、ｓは、任意であるが固定の基準角度位置に関連するらせん上の放射線源Ｓの角度位置であり、ｙ（ｓ）は、３次元空間における放射線源の位置である。
【００２５】
測定値Ｄｆ（ｙ，Θ）は、以下の線積分（２）によって記述されることができる。
【数２】
2006505337.tif 000004
【００２６】
上式の単位ファクタ（unity factor）Θは、測定値に関連する放射線の方向を示す。
【００２７】
ステップ１０５において、方程式（１）に従って、測定値は、ｑに関して、すなわち位置ｑ＝ｓにおける放射線源の角度位置に関して、偏微分される。この点で、ｙだけがΘではなくｑに依存し、従って、平行な放射線の測定値が、微分のために考慮されなければならないことに留意すべきである。平行な放射線は、同じ円錐角を有し、放射線の円錐角は、回転軸を通るとともに回転軸に対して垂直に延びる放射線に対し、図１に示される座標系のｘｚ平面への放射線の投影によって囲まれる角度である。図５に示すように、同じ円錐角を有する放射線は、フォーカス中心の検出器の場合には同じ検出器の行に入り、従って、偏微分のために、同じ行にあるが異なる投影からの測定値が考慮される。微分は、例えば有限差分法によって行われることができる。
【００２８】
単位ファクタΘは、いわゆるＫ平面５１によって表されることができるＫ角度γに依存する。以下、Ｋ平面５１について詳しく説明する。
【００２９】
Ｋ平面５１を決定するために、関数、
【数３】
2006505337.tif 000005
が取り入れられ、この関数は、非負の整数値ｎ及びｍに依存する。ここで、ｎ＞ｍである。この実施例において、ｎ＝２及びｍ＝１である。しかしながら、ｎ、ｍの他の値が選択されることもできる。それにもかかわらず、方程式（１）は正確なままであり、Ｋ平面５１の位置のみが変化する。更に、ベクトル関数、
【数４】
2006505337.tif 000006
及び単位ベクトル、
【数５】
2006505337.tif 000007
が規定される。ベクトルβは、放射線源位置ｙ（ｓ）から位置ｘまでを指す。Ｋ平面を決定するために、値ｓ２εＩＰＩ（ｘ）は、ｙ（ｓ）、ｙ（ｓ１（ｓ、ｓ２））ｙ（ｓ２）及びｘが１つの平面内に位置するように選択される。この平面は、Ｋ平面５１と呼ばれ、Ｋ平面５１と検出器表面との間の交線は、Ｋライン５３と呼ばれる。図６は、Ｋ平面の扇様の一部を示す。扇形の辺は、放射線源の位置で交わる。Ｋ平面５１のこの定義は、ｓ２に関する方程式、
【数６】
2006505337.tif 000008
の解に等しい。こうして、ｕは、Ｋ平面５１の法線ベクトルである。ベクトル関数Θ（ｓ，ｘ，ｙ）を決定するために、ベクトル、
【数７】
2006505337.tif 000009
が規定される。β及びｅに関する定義を使用して、ベクトル関数Θ（ｓ，ｘ，ｙ）は、次のように表現されることができる。
【数８】
2006505337.tif 000010
ベクトルβ及びｅは両方ともｕに垂直な方向であるので、Ｋ角度γは、ベクトルΘの方向を示し、それゆえＫ平面内の放射線の方向を示す。
【００３０】
Ｋ平面及びＫラインは、本願明細書において明確に参照される文献Ｅ１に詳しく記述されている。
【００３１】
ステップ１０７において、Ｋライン沿いに得られた測定値は、Ｋ角度γのサインの逆数に対応する重み係数を乗じられ、方程式（１）に従って積分される。このために、検査ゾーンの各位置ｘ及び各投影角について、Ｋラインが決定される。上述したように、値ｓ２εＩＰＩ（ｘ）は、ｙ（ｓ）、ｙ（ｓ１（ｓ，ｓ２））、ｙ（ｓ２）及びｘが１つの平面、すなわちＫ平面に位置するように、選択される。Ｋラインは、Ｋ平面と検出器表面との間の交線として決定される。重み係数による乗算及び積分は、例えばフーリエフィルタリングによって実施されることができる。
【００３２】
微分され、積分された測定値は、次の方程式によって表されることができる。
【数９】
2006505337.tif 000011
上式で、ｐ（ｙ（ｓ）、Φ（ｓ，ｘ））は、微分され、積分された測定値を示し、Φ（ｓ，ｘ）は、放射線源位置ｙ（ｓ）から、検査ゾーンの位置ｘの方向を指す単位ファクタである。
【００３３】
方程式（１）のまだ欠けている積分ステップ又は測定値の逆投影は、次の方程式によって記述されることができる。
【数１０】
2006505337.tif 000012
【００３４】
この方程式に従って、それぞれの測定値は、検査ゾーンにおける空間的な吸収分布を再構成するために、係数１／｜ｘ−ｙ（ｓ）｜を乗じられなければならない。この係数は、位置ｘに依存し、従って、放射線源位置ｙ（ｓ）及び位置ｘのそれぞれの組み合わせについて改めて計算されなければならない。更に、ｓに関する積分は、らせんのセグメントＩＰＩ（ｘ）に沿って行われる。このように、積分間隔は、吸収が決定されるべき位置ｘに依存し、従って、積分間隔は、それぞれの位置ｘについて決定されなければならない。方程式（１０）に従う積分は、これらの理由により、大量の計算作業を必要とするので、積分変数ｓは、以後、投影角φと置き換えられる。投影角φは、回転軸に対して垂直な平面（以下、ｘｙ平面と称される）に投影されるオブジェクトポイントｘと、放射線源から発せられ、位置ｘを通る放射線のｘｙ平面への投影とのＰＩラインによって囲まれる角度である。
【００３５】
方程式（１０）における積分は、らせんセグメントＩＰＩ（ｘ）に沿って実施される。このセグメントは、ＰＩラインによって囲まれており、従って、積分は、積分変数の置換ののち、それぞれの測定値について０からπまで実施されなければならない。こうして、積分間隔は、検査ゾーンのそれぞれの位置ｘについて同じである。
【００３６】
積分変数ｄｓ及びｄφ間の関係は、図７から導き出されることができる。この図は、らせん１７、位置ｘにおけるオブジェクトポイント３５、及び放射線源位置ｙ（ｓ）及びｙ（ｓ＋ｄｓ）のｘｙ平面への投影を示す。次の方程式は、図７から与えられる。
【数１１】
2006505337.tif 000013
上式で、Ｐｘｙは、ベクトルのｘｙ平面への投影のための射影演算子を示し、Ｒは、らせん１７の半径を示す。扇角εは、放射線源位置から回転軸までの垂線と、放射線源位置から発せられ、位置ｘを通る放射線のｘｙ平面への投影と、によって囲まれる角度である。添え字ｘ及びｙは、ベクトルのｘ及びｙ成分を示す。成分は、図１に示されるデカルト座標系に関する。
【００３７】
こうして、次の方程式が、円錐角λについて得られる。
【数１２】
2006505337.tif 000014
【００３８】
方程式（１０）、（１１）及び（１２）は、次式を与える。
【数１３】
2006505337.tif 000015
【００３９】
この方程式に従って、ステップ１０９において、測定値は、円錐角λのコサインに対応する第１の重み係数と、扇角εのコサインの逆数の値に対応する第２の重み係数とを乗じられる。更に、測定値は、半径Ｒの逆数を乗じられることができる。半径は、取得中、一定であるので、後者の乗算は、逆投影ののちに実施されることもできる。
【００４０】
小さい角度λ及びεの場合、これらの角度のコサインはほぼ１であるので、重み係数ｃｏｓ（λ）及び１／ｃｏｓ（ε）による乗算は、無視されることができる。
【００４１】
ステップ１０９における重み係数は、円錐角λ及び扇角εに依存する。こうして、重み係数は、同じ放射線によって通過される検査ゾーン内のすべての位置ｘについて同じであり、このことは、これらの位置について、重み係数が一度だけ計算されればよいことを意味する。方程式（１）の重み係数１／｜ｘ−ｙ（ｓ）｜による知られている重みづけと比較して、前述のものは、必要な計算作業の量を大幅に低減させる。
【００４２】
逆投影の前に、測定値のリビニング（rebinning）が、ステップ１１１において実施されることができる。リビニング処理の結果として、測定値は、あたかもそれらが別の放射線源（らせんの一部に配置され、互いに平行な扇形ビームをその都度放出することができる細長い放射線源）及び別の検出器（回転軸１４を含む平坦な矩形の「仮想」の検出器）によって測定されたかのように、ソートしなおされ、補間しなおされる。
【００４３】
これについて、図８を参照して詳細に説明する。図８において、参照数字１７は、らせん状の軌道を示し、放射線源は、そこから検査ゾーンを照射する。参照数字４３は、放射線源位置Ｓ０から発する扇形の放射ビームを示し、その放射線は、回転軸１４を含む平面において広がる。位置Ｓ０における放射線源によって放出される円錐形の放射ビームは、回転軸１４と平行な平面に位置し、放射線源位置Ｓ０において交わる複数の平坦な扇形ビームからなるとみなされることができる。図８は、これらの扇形ビームの１つ、すなわち扇形ビーム４３だけを示している。
【００４４】
更に、図８は、扇形ビーム４３と平行に延在し、互いに及び回転軸１４に平行である平面に位置する他の扇形ビーム４１、４２及び４４、４５を示している。関連する放射線源位置Ｓ−２、Ｓ−１及びＳ１、Ｓ２は、それぞれ、放射線源位置Ｓ０に達する前後に放射線源Ｓによって占められる。扇形ビーム４１乃至４５のすべての放射線は、同じ投影角を有する。
【００４５】
扇形ビーム４１乃至４５は、テント様の形状をもつ放射ビーム７０を規定する。扇形ビームのそれぞれのグループについて、矩形の仮想検出器１６０が規定され、この仮想検出器１６０は、回転軸１４を含む平面に位置し、あるグループの平行な扇形ビームに対して垂直である。仮想検出器１６０のコーナーポイントは、外側の放射線源位置から、この平面を通って、反対に位置するらせんのセグメントに入射する放射線のパンクチャ（puncture、突き刺し）ポイントを構成する。図８の扇形ビーム７０の場合、これらは、扇形ビーム４１及び４５と、らせんとの交点である。矩形の検出器１６０上には、デカルト形態で、すなわち行及び列に並べられる構成要素が規定され、これに関して、測定値が、補間しなおされる。
【００４６】
リビニング後に得られる測定値は、そののち、本実施例において方程式（１３）に従う逆投影によって、検査ゾーンにおける吸収分布を再構成するために使用される。
【００４７】
ステップ１１３において、ボクセルＶ（ｘ，ｙ，ｚ）が、選択可能な（ｘ，ｙ，ｚ）ゾーン（視野又はＦＯＶ）内で決定される。そののち、ステップ１１５において、投影角が、レンジ［φ０，φ０＋π］内で選択される。ここで、φ０は、ボクセルＶ（ｘ，ｙ，ｚ）が放射ビームに入る角度である。ステップ１１７において、投影の放射線が、ボクセルＶ（ｘ，ｙ，ｚ）の中心を通るかどうか調べられる。投影のどの放射線もボクセルの中心を通らない場合、関連する値が、隣り合う放射線の測定値の補間によって決定されなければならない。ボクセルを通る放射線に関連しうる測定値又は補間によって得られる測定値は、ステップ１１９において、ボクセルＶ（ｘ，ｙ，ｚ）に関して累積される。ステップ１２１において、投影角φ０乃至φ０＋πのすべての投影が考慮されたかどうか調べられる。まだすべての投影が調べられていない場合、フローチャートは、ステップ１１５に分岐する。他の場合、ステップ１２３において、ＦＯＶ内のすべてのボクセルＶ（ｘ，ｙ，ｚ）が処理されたかどうか調べられる。すべてのボクセルＶ（ｘ，ｙ，ｚ）が処理されていない場合、プロシージャはステップ１１３から続く。しかしながら、ＦＯＶ内のすべてのボクセルＶ（ｘ，ｙ，ｚ）が処理されている場合、吸収が、ＦＯＶ全体で決定され、再構成方法が終了される（ステップ１２５）。
【図面の簡単な説明】
【００４８】
【図１】本発明による方法を実施するのに適したコンピュータ断層撮影装置を示す図。
【図２】本発明による方法のフローチャート。
【図３】検査ゾーンのあるポイントに関する、ＰＩライン及び積分経路を示す図。
【図４】回転軸に対して垂直な平面に投影される検査ゾーンのポイントに関する、ＰＩライン及び積分経路を示す図。
【図５】それぞれ異なるビーム位置をもつ平行ビームを示す図。
【図６】Ｋ平面及びＫラインを示す図。
【図７】回転軸に対して垂直な平面に投影される、放射線源の２つの異なる放射線源位置、検査ゾーン内のポイント及び軌道の概略図。
【図８】平行な面にリビニングすることによって形成される扇形ビームを示す図。
【符号の説明】
【００４９】
１  ガントリ
２、５  モータ
３  コリメータ構造
４  円錐形の放射ビーム
７  制御ユニット
１０  画像処理コンピュータ
１１  モニタ
１３  検査ゾーン
１４  回転軸
１６  検出器ユニット
１７  らせん
３１  ＰＩライン
３５  オブジェクトポイント
４１．．．４５  扇形の放射ビーム
５１  Ｋ平面
５３  Ｋライン
１６０  仮想検出器
Ｓ  放射線源
Ｓ−２．．．Ｓ２  放射線源位置
【図１】
2006505337.tif 000016
【図２】
2006505337.tif 000017
【図３】
2006505337.tif 000018
【図４】
2006505337.tif 000019
【図５】
2006505337.tif 000020
【図６】
2006505337.tif 000021
【図７】
2006505337.tif 000022
【図８】
2006505337.tif 000023
【国際調査報告】
2006505337.tif 000024
2006505337.tif 000025