(57)【要約】      （修正有）【課題】２次元Ｘ線検出器による回転撮影データ、特にＣ型アームによる回転撮影データから、鮮明な３次元的Ｘ線ＣＴ像を再構成するための、幾何学パラメータの決定方法を提供する。【解決手段】Ｘ線源と、Ｘ線画像データを出力する２次元Ｘ線検出器と、前記Ｘ線源及び前記２次元Ｘ線検出器を、回転移動させるための回転手段と、前記Ｘ線画像データに基づいて画像再構成演算を行なう画像処理装置と、を備えたＸ線ＣＴ装置において、前記画像処理装置は、基準回転中心軸の位置と実際の回転移動を示す回転中心軸の位置との相対変位量を、前記基準回転中心軸及び前記２次元Ｘ線検出器の回転軌道面の近傍に位置したファントムの投影像の特徴量に基づいて算出する回転中心軸投影位置計算手段を備え、前記逆投影手段は、前記相対変位量に基づいて前記Ｘ線画像データの修正を行い、画像を再構成することを特徴とするＸ線ＣＴ装置である。【選択図】      図１
2006130060.tif 000002
【特許請求の範囲】
【請求項１】
被検体にＸ線を照射するＸ線源と、
前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、
前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、
前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、
を備えたＸ線ＣＴ装置において、
前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動するように設計され、
前記画像処理装置は、
前記基準回転中心軸の位置と実際の回転移動を示す回転中心軸の位置との相対変位量を、前記基準回転中心軸及び前記２次元Ｘ線検出器の回転軌道面の近傍に位置したファントムの投影像の特徴量に基づいて算出する回転中心軸投影位置計算手段を備え、
前記逆投影手段は、前記相対変位量に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する、
ことを特徴とするＸ線ＣＴ装置。
【請求項２】
被検体にＸ線を照射するＸ線源と、
前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、
前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、
前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、
を備えたＸ線ＣＴ装置において、
前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動し、かつ、前記２次元Ｘ線検出器は、前記基準回転中心軸に対して所定の基準角度をなして取り付けられると設計され、
前記画像処理装置は、
前記基準回転中心軸に対する前記２次元Ｘ線検出器の実際の取り付け角度と前記基準角度との角度差を修正するために必要な補正角度を、前記基準回転中心軸近傍かつ前記２次元Ｘ線検出器の回転軌道面から所定距離離して位置したファントムの投影像の特徴量に基づいて算出する検出器取り付け補正角計算手段を備え、
前記逆投影手段は、前記補正角度に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する、
ことを特徴とするＸ線ＣＴ装置。
【請求項３】
被検体にＸ線を照射するＸ線源と、
前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、
前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、
前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、
を備えたＸ線ＣＴ装置において、
前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動し、かつ、前記２次元Ｘ線検出器は、前記基準回転中心軸に対して所定の基準角度をなして取り付けられると設計され、
前記画像処理装置は、
前記基準回転中心軸及び前記２次元Ｘ線検出器の回転軌道面の近傍に位置したファントムの投影像の特徴量と、前記基準回転中心軸近傍かつ前記２次元Ｘ線検出器の回転軌道面から所定距離離して位置したファントムの投影像の特徴量と、前記所定距離と、に基づいて、前記基準回転中心軸の位置と実際の回転移動を示す回転中心軸の位置との相対変位量と、前記基準回転中心軸に対する前記２次元Ｘ線検出器の実際の取り付け角度と前記基準角度との角度差を修正するために必要な補正角度と、を算出する幾何学パラメータ計算手段を備え、
前記逆投影手段は、前記相対変位量及び前記補正角度に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する、
ことを特徴とするＸ線ＣＴ装置。
【請求項４】
前記被検体を載置した寝台を、更に備え、
前記ファントムは、前記寝台に載置される前記被検体の体軸方向に沿って配置された円柱状又は中空のパイプ状のファントムとして構成される、
ことを特徴とする請求項１乃至３の一つに記載のＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ装置に係り、特に、本発明は、２次元Ｘ線検出器による回転撮影データから、３次元的Ｘ線ＣＴ像を再構成するのに有効な技術に関する。
【背景技術】
【０００２】
従来のコーンビームＸ線ＣＴ装置は、Ｘ線源とこのＸ線源に対向して配置される２次元Ｘ線検出器（Ｘ線イメージインテンシファイアとテレビカメラとを組み合わせたものを含む。）とを同一回転中心の円軌道面上で回転移動させながら、回転中心軸上に位置する被検体のＸ線透視像を２次元Ｘ線検出器で撮影し、更に前記Ｘ線透視像を基に画像再構成演算を行うことによって、３次元的Ｘ線ＣＴ像を得るものである。
【０００３】
このような２次元Ｘ線検出器を用いたコーンビームＸ線ＣＴ装置の画像再構成演算アルゴリズムとして、非特許文献１に記載のフェルドカンプの方法が代表的である。
【０００４】
また、コーンビームＸ線ＣＴ装置では機械的製作誤差に伴うアーチファクトのない３次元的Ｘ線ＣＴ像を得るために、回転中心軸の設計位置とのずれや、Ｘ線源の焦点の軌道を含む平面である回転軌道面等をパラメータ（以下、幾何学パラメータと称す）として求め、画像再構成演算において補正しながら演算することがある。幾何学パラメータを求める方法として特許文献１及び特許文献２、またガントリタイプのコーンビームＸ線ＣＴ装置の回転中心軸を精度良く求める方法として特許文献３に記載の従来技術がある。
【０００５】
特許文献１は、Ｘ線吸収係数の大きな小金属球を撮影し、その投影像のサイノグラムからのずれを用い、回転中心の幾何学補正を行う方法を開示する。
【０００６】
特許文献２は、棒状の支持体にＸ線吸収係数の大きい小金属球を複数個入れたファントムを回転中心軸に平行に配置して、３６０度回転撮影し、その金属球が２次元Ｘ線透過像上に描く楕円軌道を利用して、回転軌道面、回転中心軸、検出器取付け補正角などの幾何学パラメータを求める方法を開示する。
【０００７】
特許文献３は、Ｘ線吸収係数の大きなワイヤー線を撮影し、評価関数によりその３次元的Ｘ線ＣＴ像の鮮明度を評価しながら再構成演算で使用する幾何学パラメータ、回転中心軸投影を精度良く求める方法を開示する。
【非特許文献１】「Practical Cone-Beam Algorithm; L.A.Feldkamp, et al.; J.Optical Society of America, A/Vol. 1(6), (1984), pp.612-619」
【特許文献１】特開平７-１６２２０号公報
【特許文献２】特開平９-１７３３３０号公報
【特許文献３】特開２０００-２０１９１８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
本発明者は、前記従来技術による撮影系幾何学パラメータの決定方法を検討した結果、以下の問題点を見いだした。
【０００９】
特許文献１に記載の幾何学パラメータ補正方法は、投影像上での補間演算を用いた幾何学パラメータ調整方法である。そのためＣＴ再構成に必要な精度が出ない場合があるという問題があった。
【００１０】
特許文献２に記載の幾何学パラメータ補正方法も、投影像上での補間演算を用いた幾何学パラメータ調整方法であるため、ＣＴ再構成に必要な精度が出ないという問題がある。
【００１１】
また特許文献３に記載の回転中心軸算出方法は、ワイヤー断面の再構成像をそのアーチファクトを含め、評価関数でその鮮明度を評価しながら幾何学パラメータを求めていく手法であるため、ガントリタイプのように回転軌道の真円性の高いコーンビームＸ線ＣＴ装置では精度良く回転中心軸を算出することができるが、Ｃ型アームのようにたわみやぶれがある回転軌道の場合には正しい幾何学パラメータが得られないという問題がある。
【００１２】
本発明はこのような事情に鑑みてなされたもので、２次元Ｘ線検出器による回転撮影データ、特にＣ型アームによる１８０°＋検出器開き角の回転撮影データから、鮮明な３次元的Ｘ線ＣＴ像を再構成するための、幾何学パラメータの決定方法を提供することを目的とする。
【課題を解決するための手段】
【００１３】
前記目的を達成するために、本発明に係るＸ線ＣＴ装置は、被検体にＸ線を照射するＸ線源と、前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、を備えたＸ線ＣＴ装置において、前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動するように設計され、前記画像処理装置は、前記基準回転中心軸の位置と実際の回転移動を示す回転中心軸の位置との相対変位量を、前記基準回転中心軸及び前記２次元Ｘ線検出器の回転軌道面の近傍に位置したファントムの投影像の特徴量に基づいて算出する回転中心軸投影位置計算手段を備え、前記逆投影手段は、前記相対変位量に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する。
【００１４】
また、本発明に係るＸ線ＣＴ装置は、被検体にＸ線を照射するＸ線源と、前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、を備えたＸ線ＣＴ装置において、前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動し、かつ、前記２次元Ｘ線検出器は、前記基準回転中心軸に対して所定の基準角度をなして取り付けられると設計され、前記画像処理装置は、前記基準回転中心軸に対する前記２次元Ｘ線検出器の実際の取り付け角度と前記基準角度との角度差を修正するために必要な補正角度を、前記基準回転中心軸近傍かつ前記２次元Ｘ線検出器の回転軌道面から所定距離離して位置したファントムの投影像の特徴量に基づいて算出する検出器取り付け補正角計算手段を備え、前記逆投影手段は、前記補正角度に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する。
【００１５】
また、本発明に係るＸ線ＣＴ装置は、被検体にＸ線を照射するＸ線源と、前記Ｘ線源に対向して配置され、前記被検体を透過した前記Ｘ線を検出して前記被検体のＸ線画像データを出力する２次元Ｘ線検出器と、前記Ｘ線源及び前記２次元Ｘ線検出器を、同一回転中心軸を中心に回転移動させるための回転手段と、前記Ｘ線画像データに基づいて画像再構成演算を行なう逆投影手段を備えた画像処理装置と、を備えたＸ線ＣＴ装置において、前記Ｘ線源及び前記２次元Ｘ線検出器は、基準回転中心軸を中心に回転移動し、かつ、前記２次元Ｘ線検出器は、前記基準回転中心軸に対して所定の基準角度をなして取り付けられると設計され、前記画像処理装置は、前記基準回転中心軸及び前記２次元Ｘ線検出器の回転軌道面の近傍に位置したファントムの投影像の特徴量と、前記基準回転中心軸近傍かつ前記２次元Ｘ線検出器の回転軌道面から所定距離離して位置したファントムの投影像の特徴量と、前記所定距離と、に基づいて、前記基準回転中心軸の位置と実際の回転移動を示す回転中心軸の位置との相対変位量と、前記基準回転中心軸に対する前記２次元Ｘ線検出器の実際の取り付け角度と前記基準角度との角度差を修正するために必要な補正角度と、を算出する幾何学パラメータ計算手段を備え、前記逆投影手段は、前記相対変位量及び前記補正角度に基づいて前記Ｘ線画像データの修正を行い、画像を再構成する。
【００１６】
また、本発明に係るＸ線ＣＴ装置は、前記被検体を載置した寝台を、更に備え、前記ファントムは、前記寝台に載置される前記被検体の体軸方向に沿って配置された円柱状又は中空のパイプ状のファントムとして構成される。
【発明の効果】
【００１７】
本発明によれば、２次元Ｘ線検出器による回転撮影データ、特にＣ型アームによる１８０°＋検出器開き角の回転撮影データから、鮮明な３次元的Ｘ線ＣＴ像を再構成するための、幾何学パラメータの決定方法を提供することができる。
【発明を実施するための最良の形態】
【００１８】
以下、添付図面に従って本発明に係るＸ線ＣＴ装置の好ましい実施の形態について詳説する。
【００１９】
図１は、本発明を適用したＣ型アーム方式のコーンビームＸ線ＣＴ装置の概略構成を示すブロック図である。図２は、画像歪み補正テーブルの生成に使用されるホールチャート１８を示す概念図である。図３は、前記ホールチャート１８を２次元Ｘ線検出器のＸ線入射面に取り付ける際に生じる検出器取付け角を説明するための図である。
【００２０】
図１のコーンビームＸ線ＣＴ装置１は、被検体４０に対してＸ線を照射し、被検体４０のＸ線透過画像を撮影してＸ線画像データを得る撮影部１０と、撮影部１０の各構成要素を制御したり、Ｘ線画像データに基づいて被検体４０の３次元的Ｘ線ＣＴ像を再構成したりする制御演算部２０とを備える。
【００２１】
（撮影部１０）
撮影部１０は、被検体４０を載せる寝台１７と、寝台１７に載せられた被検体４０にＸ線を照射するＸ線源１１、Ｘ線源１１に対向する位置に設置され、被検体４０を透過したＸ線を検出することによりＸ線画像データを出力する２次元Ｘ線検出器１２と、Ｘ線源１１及び２次元Ｘ線検出器１２を機械的に接続するＣ型アーム１３とを備える。また撮影部１０は、Ｃ型アーム１３を保持するＣ型アーム保持体１４と、Ｃ型アーム保持体１４を天井に取り付ける天井支持体１５と、天井支持体１５を図１の状態で、前後左右の２次元方向に移動可能に支持する天井レール１６とを備える。
【００２２】
Ｘ線源１１は、Ｘ線を発生するＸ線管１１ｔと、Ｘ線管１１ｔからのＸ線照射の方向を円錐または四角錐状に制御するコリメータ１１ｃとを備える。
【００２３】
２次元Ｘ線検出器１２は、Ｘ線透過像を可視光像に変換するＸ線イメージインテンシファイア１２ｉと、Ｘ線イメージインテンシファイア１２ｉによる可視光像を撮影するテレビカメラ１２ｃとを備える。なお本実施の形態では、２次元Ｘ線検出器１２は、Ｘ線イメージインテンシファイア１２ｉ及びテレビカメラ１２ｃを備えるが、２次元Ｘ線検出器１２は、フラットパネルディテクター（ＦＰＤ）であってもよい。またその形状は円形、方形等のいかなる形状であってもよい。そして、２次元Ｘ線検出器１２はその検出器素子列が回転中心軸３１に平行（０°）あるいは９０°の角度をなして設置される。例えば、２次元Ｘ線検出器１２として長方形のフラットパネルディテクター（ＦＰＤ）を用いる場合、その長辺を回転中心軸と９０°の角度をなして設置すると胸部、腹部等、大視野の断面像の撮影に適合するし、長辺を回転中心軸と平行（０°）の角度をなして設置すると、頭頸部、四肢等の撮影に有用である。２次元Ｘ線検出器１２は、その検出器素子列が、回転中心軸３１に対し所定の基準角度だけ手動もしくは電動で回転できるようになっていてもよい。
【００２４】
上記Ｃ型アーム１３は、被検体４０の撮影に際して、所定の投影角度毎に回転中心軸３１を中心として回転移動する。これにより、上記Ｘ線源１１及び２次元Ｘ線検出器１２は、ほぼ同一の円軌道上で回転移動しながら、Ｘ線撮影を行う。この回転移動については、画像再構成演算に使用される幾何学パラメータが存在する。すなわち、Ｃ型アーム１３が回転移動することにより、Ｘ線源１１と２次元Ｘ線検出器１２とが描く円軌道を含む面である回転軌道面（ミッドプレーン）３０と、回転中心軸３１、及び検出器取付け角の基準角度（０°あるいは９０°）からのずれである。
【００２５】
これらの幾何学パラメータは、原理的には、撮影部１０の設計データから定まるものであるが、実際的には撮影部１０の製作誤差や部材の変形に起因して、個々のコーンビームＸ線ＣＴ装置１に固有の値となる。
【００２６】
（制御演算部２０）
制御演算部２０は、撮影部１０を制御する撮影部制御手段１００と、撮影部１０が出力したＸ線画像データを収集して格納する画像収集手段１１０と、収集されたＸ線画像データに基づいて３次元的Ｘ線ＣＴ像を再構成する再構成手段２００と、撮影部１０の機械的製作上の誤差を数値的に表わし、再構成手段２００における３次元再構成の際に補正データとして用いる幾何学パラメータを求めるための幾何学パラメータ計算手段３００とを備える。更に、再構成手段２００で生成した３次元的Ｘ線ＣＴ像を表示する画像表示手段２１０を備える。
【００２７】
撮影部制御手段１００は、Ｃ型アーム１３が、回転中心軸３１の回りを回転する（以下、「プロペラ回転」という。）回転移動を制御する撮影系回転制御手段１０１と、天井支持体１５の天井レール１６上での位置を制御してＣ型アーム１３の被検体４０に対する位置を２次元的に制御する撮影系位置制御手段１０２とを備える。更に撮影部制御手段１００は、Ｘ線管１１ｔに流す管電流のＯＮ、ＯＦＦなどを制御するＸ線照射制御手段１０３と、寝台１７の位置を制御して被検体４０の位置を調整するための寝台制御手段１０４と、２次元Ｘ線検出器１２によるＸ線透過像の撮影を制御する検出系制御手段１０５とを備える。
【００２８】
（再構成手段２００）
再構成手段２００は、前処理手段２０１と、画像歪み補正手段２０２と、フィルタリング手段２０３と、逆投影手段２０４とを備える。
【００２９】
前処理手段２０１は、画像収集手段１１０が収集したＸ線画像データをＸ線吸収係数の分布像に変換するための手段である。本実施の形態では先ず、被検体４０及び寝台１７を撮影視野内に配置しない状態で予め撮影された空気のＸ線透過像の各画素データに対して自然対数変換演算を施す。次に被検体４０を寝台１７に載せた状態で撮影したＸ線透過画像の各画素データに対して自然対数変換演算を施す。上記２つのＸ線透過画像の差分を取ることにより、被検体４０及び寝台１７のＸ線吸収係数の分布像を得る。
【００３０】
画像歪み補正手段２０２は、前処理手段２０１が生成したＸ線吸収係数の分布像の画像歪みを補正する。この画像歪みは、Ｘ線イメージインテンシファイア１２ｉによってＸ線透過像を可視光像に変換する際に生ずる画像歪みであり、後述する画像歪み補正テーブル格納手段３３０により格納されている画像歪み補正テーブルを用いて、前処理手段２０１で得られたＸ線吸収係数の分布像の画像歪みを補正する。
【００３１】
フィルタリング手段２０３は、Ｘ線ＣＴ画像再構成におけるフィルタリング処理を行う。
【００３２】
逆投影手段２０４は、フィルタリング処理後のＸ線画像データに基づいて逆投影演算を行う手段であり、３次元的Ｘ線ＣＴ像を生成する。
【００３３】
（幾何学パラメータ計算手段３００）
幾何学パラメータ計算手段３００は、再構成手段２００で画像再構成を行う際に必要となる幾何学パラメータを算出する手段である。幾何学パラメータ計算手段３００が扱う幾何学パラメータは、回転軌道面３０と、回転中心軸３１と、Ｘ線イメージインテンシファイア１２ｉによる画像歪み、及び２次元Ｘ線検出器１２の取付け角である。図１における回転軌道３０と回転中心軸３１、Ｘ線イメージインテンシファイア１２ｉによる画像の歪み、それと後述する検出器取付け補正角である。これらの幾何学パラメータを算出するために、幾何学パラメータ計算手段３００は、画像歪み補正テーブル生成手段３２０と、画像歪み補正テーブル格納手段３３０と、回転軌道面算出手段３５０と、回転中心軸投影位置計算手段３７０と、検出器取付け補正角計算手段３８０とを備える。上述した幾何学パラメータ計算手段３００の各構成要素については、後述する。
【００３４】
画像歪み補正テーブル生成手段３２０は、上述の画像歪み補正手段２０２で使用する画像歪み補正テーブルを生成する。画像歪み補正テーブル生成手段３２０は、図２に示すホールチャート１８を用いて画像歪み補正テーブルを生成する。図２のホールチャート１８は、Ｘ線に対する吸収の大きな材料で形成した板材に、格子状配列で多数の小さなホール１８ｈを穿設して形成されている。画像歪み補正テーブルを生成するには、まず、Ｘ線イメージインテンシファイア１２ｉの前面にホールチャート１８を固定し、このホールチャート１８のＸ線投影像（以下「ホールチャート歪み投影像」という。）を撮影する。次に、前処理手段２０１が、ホールチャート歪み投影像に前処理を行う。そして、ホールチャート歪み投影像における各ホール１８ｈの投影位置を検出する。次に、ホールチャート１８が歪みなくＸ線入射面に投影された場合を仮想した各ホール１８ｈの仮想投影位置を演算する。そして、検出した各ホール１８ｈの投影位置が、仮想投影位置に一致するように変換する計算をホール１８ｈ毎に行い、画像歪み補正テーブルを生成する。
【００３５】
画像歪み補正テーブル格納手段３３０は、画像歪み補正テーブル生成手段３２０で生成した画像歪み補正テーブルを磁気ディスク等に格納する。そして、画像歪み補正手段２０２がＸ線透過像の歪みを補正する際に、格納した画像歪み補正テーブルを読み出す。
【００３６】
回転軌道面算出手段３５０は、Ｘ線源の焦点の軌道を含む平面である回転軌道面３０の２次元Ｘ線透過像上の座標を算出する手段である。これは例えば、棒状の支持体にＸ線吸収係数の大きい小金属球を複数個入れたファントムを回転中心軸に平行に配置して、回転撮影を行い、その金属球が２次元Ｘ線透過像上に描く楕円軌道を利用して求めることができる。
【００３７】
コーンビームＸ線ＣＴ装置１は、図４及び図５に示すように、円柱形もしくは中空のパイプ状ファントム５０のＸ線透過画像（以下「ファントム投影像」という。）を撮影し、ファントムＸ線画像データを出力する。
【００３８】
回転中心軸投影位置計算手段３７０は、ファントム投影像の再構成演算から、回転中心軸の２次元Ｘ線透過像上への投影座標（以下「回転中心軸投影位置」という。）を算出する手段である。
【００３９】
検出器取付け補正角計算手段３８０は、２次元Ｘ線検出器１２の画像上での座標軸と、実際の回転中心軸との角度の差（検出器取付け補正角）を算出する手段である。
【００４０】
図３を用い検出器取付け角について説明する。３２は回転中心軸３１の投影、３０は回転軌道面である。２２はホールチャート１８に固定された縦軸、２１はホールチャート１８の横軸であり、各々、回転中心軸３１の投影３２、回転軌道面３０から角度ｂ（２４）だけ傾いて取り付いている。２４を検出器取付け角とよび、角度ｂ（２４）がゼロでないと、画像歪み補正手段２０２で画像歪み補正後の画像は角度ｂだけ回転する。この回転角ｂを検出器取付け補正角計算手段３８０で求め、Ｘ線画像を検出器取付け角ｂとは逆方向に回転してこれを補正する。ただし、本実施の形態ではＸ線撮像系としてＸ線イメージインテンシファイア１２ｉ及びテレビカメラ１２ｃを用いたが、ＴＦＴ素子等を用いた２次元Ｘ線検出器に置き換える場合には、Ｘ線検出器の取付け角度のずれが角度ｂになる。
【００４１】
上記のＣ型アーム方式コーンビームＸ線ＣＴ装置１の仕様例は次のとおりである。Ｘ線管１１ｔと回転中心軸３１との距離は８００ｍｍ、回転中心軸３１と２次元Ｘ線検出器１２のＸ線入力面、すなわちＸ線イメージインテンシファイア１２ｉのＸ線入射面との距離は４００ｍｍ、Ｘ線イメージインテンシファイア１２ｉのＸ線入射面の大きさは４００ｍｍ、画像サイズは１０２４×１０２４（走査線数）である。そして２次元Ｘ線検出器１２のピッチは０．４ｍｍである。
【００４２】
撮影系回転制御手段１０１は、被検体４０の頭部をＣ型アーム１３側に向けて寝台１７上に載置した状態で、２次元Ｘ線検出器１２を、被検体４０の左手の方向（−１００°）から天井方向（０°）を通過し、被検体４０の右手方向（＋１００°）まで移動させる（この時、Ｘ線源１１は被検体４０の右手方向（＋１００°）から寝台下方向を通過し、被検体４０の左手方向（−１００°）まで移動する）ことにより、２００度の投影角度にわたって被検体４０の２次元Ｘ線透過画像が撮影される。Ｃ型アーム１３の回転速度の代表例は１秒当たり４０度で、スキャン時間は５秒である。なお、被検体４０の足部をＣ型アーム１３側に向けて寝台１７上に載置した状態で、２次元Ｘ線検出器１２を、被検体４０の右手の方向（−１００°）から天井方向（０°）を通過し、被検体４０の左手方向（＋１００°）まで移動させる場合もある。
【００４３】
次に、本実施形態のＣ型アーム方式コーンビームＸ線ＣＴ装置１による撮影における動作の概要について説明する。
【００４４】
先ず撮影系回転制御手段１０１はＣ型アーム１３のプロペラ回転を開始する。回転加速期間を経たのち、Ｘ線照射制御手段１０３は、Ｘ線管１１ｔにＸ線を照射させる。検出系制御手段１０５はテレビカメラ１２ｃによる可視光像の撮像を開始する。Ｘ線管１１ｔから照射されたＸ線は、被検体４０を透過した後、Ｘ線イメージインテンシファイア１２ｉへ入射する。Ｘ線透過像はＸ線イメージインテンシファイア１２ｉで可視光像に変換され、テレビカメラ１２ｃに取り込まれる。テレビカメラ１２ｃは可視光像をビデオ信号に変換し、そのビデオ信号がＡ／Ｄ変換を経た後、デジタル信号からなる２次元のＸ線画像データとして画像収集手段１１０に記録される。テレビカメラ１２ｃの撮影における標準走査モードは毎秒３０フレーム、走査線数１０２４本である。回転角度ピッチは１．３３度で、５秒間に１５０枚のＸ線透過画像を取得する。２００度の回転撮影が完了すると、Ｘ線照射制御手段１０３はＸ線管１１ｔのＸ線照射を終了し、撮影系回転制御手段１０１は回転を停止する。
【００４５】
再構成手段２００は、以上のような撮影に並行し、あるいは撮影終了後に画像収集手段１１０から２次元のＸ線画像データを読み出し、このＸ線画像データに基づいて画像再構成演算を行い、被検体４０の３次元的Ｘ線ＣＴ像の再構成演算を行う。画像表示手段２１０は、３次元的Ｘ線ＣＴ像を、ＣＲＴ装置や液晶ディスプレイ装置等からなる表示装置８０に表示する。なお画像表示手段２１０は、画像収集手段１１０に記録されたＸ線画像データに基づき、２次元画像を表示してもよい。
【００４６】
次に、図４及至６に基づいて、幾何学パラメータ計算手段３００における処理の内容を説明する。図４は、円柱ファントム５０が回転撮影される状態を説明するための概念図である。図５は、図４の状態を、寝台１７からＣ型アーム１３に向かって見た状態を示す概念図である。図６は、幾何学パラメータ計算手段３００が、幾何学パラメータ計算処理を行う手順を示すフローチャートである。以下、図６のステップ順に説明をする。
【００４７】
（ステップＳ３００）
幾何学パラメータ計算処理を開始する（Ｓ３００）。
【００４８】
（ステップＳ３１０）
ホールチャート１８をＸ線イメージインテンシファイア１２ｉのＸ線入射面に固定して回転撮影を行い、画像収集手段１１０により各投影角度におけるホールチャート歪み投影像の収集を行う。収集するホールチャート歪み投影像は、本実施の形態の場合であれば上述のように１５０枚になる。
【００４９】
（ステップＳ３２０）
画像歪み補正テーブル生成手段３２０は、ステップＳ３１０で得られた１５０枚のホールチャート歪み投影像に基づいて、画像歪み補正テーブルを生成する。
【００５０】
（ステップＳ３３０）
画像歪み補正テーブル格納手段３３０は、ステップＳ３２０で生成計算された画像歪み補正テーブルを格納する。
【００５１】
（ステップＳ３４０）
小金属球が複数個入った棒状ファントムの回転撮影を行い、画像収集手段１１０により金属球ファントムのＸ線透過画像を収集する。
【００５２】
（ステップＳ３５０）
ステップＳ３４０で収集された金属球ファントムのＸ線透過画像から、Ｘ線透過画像内の小金属球が描く楕円軌道を求め、回転軌道面の２次元Ｘ線透過像上の座標を決定する。
【００５３】
（ステップＳ３６０）
円柱ファントム５０の回転撮影を行う。図４に示すように、ファントム保持体５０ａを寝台１７上に載置する。そして、ファントム保持体５０ａは、図４及び図５に示すように、円柱ファントム５０を回転中心軸３１にできるだけ近い位置に、寝台１７から突出させて載置させる。そして、Ｘ線源１１及び２次元Ｘ線検出器１２は、円柱ファントム５０の回転撮影を行い、円柱ファントム５０のＸ線透過像を撮影し、ファントムＸ線画像データを出力する。この円柱ファントム５０のＸ線透過像には、円柱ファントム５０の透過像であるファントム投影像が含まれる。画像収集手段１１０は、ファントムＸ線画像データを収集する。円柱ファントム５０の撮影が終了すると、円柱形ファントム５０及び円柱ファントム保持体５０ａは、寝台１７上から取り除かれる。
【００５４】
（ステップＳ３７０）
回転中心軸投影位置計算手段３７０は、ステップＳ３６０で収集されたファントムＸ線画像データに含まれるファントム投影像を用い、３次元的Ｘ線ＣＴ像を結像するために基準となるような回転中心軸投影パラメータを決定する。計算の具体的な処理の詳細は、図７乃至図１２に基づいて後述する。なお、ステップＳ３７０では、回転軌道面（ミッドプレーン）３０に限定したファントムＸ線画像データを用いた処理を行えば十分であり、こうすることで演算量を少なくすることができる。
【００５５】
（ステップＳ３８０）
ステップＳ３６０で収集された円柱ファントム５０のＸ線透過像から、検出器取付け補正角計算手段３８０により、検出器取付け補正角を算出する。なお、ステップＳ３８０はステップＳ３７０と異なり、回転軌道面から離れた面で円柱ファントム５０の再構成像を生成する。
【００５６】
ステップＳ３５０、Ｓ３７０、及びＳ３８０において算出された幾何学パラメータは、逆投影手段２０４において３次元的Ｘ線ＣＴ像を結像させるのに用いられる。
【００５７】
次に、図７乃至図１２に基づいて、回転中心軸投影位置計算手段３７０と、検出器取付け補正角計算手段３８０の詳細について説明する。図７は、回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の前半の投影データから再構成される円柱ファントムの再構成断面について示した図である。図８は、回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の後半の投影データから再構成される円柱ファントムの再構成断面について示した図である。図９は、回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面について示した図である。図１０は、回転中心軸の投影パラメータが−ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面について示した図である。図１１は、本発明の実施の形態における回転中心軸投影位置計算手段３７０が、回転中心軸投影位置計算処理を行う手順を示すフローチャートである。図１２は、本発明の実施の形態における検出器取付け補正角計算手段３８０が、検出器取付け補正角計算処理を行う手順を示すフローチャートである。
【００５８】
まず、図７乃至図１０に基づいて、回転中心軸３１の初期位置が異なっている場合に得られるファントム再構成像５３について説明する。
【００５９】
分かりやすくするため、図７に回転の前半の投影データから再構成像が生成される様子、図８に回転の後半の投影データから再構成像が生成される様子を示す。図７及び図８は、回転中心軸５１の投影位置を幾何学的に正しい値５２よりもｕ軸方向に大きい値で再構成を行った場合を示す。また図７及び図８は、被検体４０の足方向（すなわち寝台１７からＣ型アーム１３に向かって）からみたもので、簡単のためＸ線イメージインテンシファイア１２ｉ及びテレビカメラ１２ｃのみを示し、他は省略した。Ｘ線イメージインテンシファイア１２ｉ及びテレビカメラ１２ｃは被検体４０の左手の方向から、天井の方向を通過し、被検体４０の右手の方向まで移動する。
【００６０】
図７に示すように、回転の前半の投影データからの再構成断面像５４は、正しい大きさにならず、再構成断面像５４が示すファントム半径は小さくなる。そして縮小した再構成ファントムの半径は、幾何学的に正しい回転中心軸で計測されるファントム端の投影位置５１と、現在使用している回転中心軸で計測されるファントム端の投影位置５２との差に等しくなる。
【００６１】
一方、図８に示す回転の後半の投影データからの再構成断面像５５のファントム半径は大きくなる。この場合、増大した再構成ファントムの半径は、幾何学的に正しい回転中心軸で計測されるファントム端の投影位置５１と、現在使用している回転中心軸で計測されるファントム端の投影位置５２との差に等しくなる。
【００６２】
図９に、回転中心軸の投影位置を幾何学的に正しい値よりもｕ軸方向に大きい値で再構成を行った場合に、実際に得られる再構成像５６を、図１０に、回転中心軸の投影位置を幾何学的に正しい値よりもｕ軸方向に小さい値で再構成を行った場合に、得られる再構成像５７を示す。ファントム再構成断面は円形とならず、左右にアンバランスなものとなる。なお図９及び図１０は分かりやすくするために模式的にしたものであり、実際は左右データの境界が突き出ることはなく、滑らかにつながった再構成像が得られる。
【００６３】
次に、図１１に基づいて、回転中心軸投影位置計算手段３７０が、回転中心軸投影位置を決定する手順を詳細に説明する。なお上述したように、以下の処理では、回転中心軸投影位置計算手段３７０では、回転軌道面（ミッドプレーン）３０に限定した投影データを用い、演算処理を行えば十分である。
【００６４】
（ステップＳ３７０）
回転中心軸投影位置計算処理を開始する（Ｓ３７０）。
【００６５】
（ステップＳ３７１）
回転中心軸投影位置（以下、centerと称す）を初期位置（設計位置）に設定する。
【００６６】
（ステップＳ３７２）
前処理手段２０１は、ステップＳ３６０で収集されたファントムＸ線画像データの前処理を行う。
【００６７】
（ステップＳ３７３）
画像歪み補正手段２０２は、ステップＳ３７２で前処理を行ったファントムＸ線画像データに対して画像歪み補正処理を行う。
【００６８】
（ステップＳ３７４）
フィルタリング手段２０３は、ステップＳ３７３で画像歪み補正処理を行ったファントムＸ線画像データについて、フィルタリング処理を行う。
【００６９】
（ステップＳ３７５）
逆投影手段２０４は、ステップＳ３７４でフィルタリング処理を行ったファントムＸ線画像データについて、回転中心軸投影位置（center）の値で、回転軌道面３０上の逆投影演算を行い、ファントム再構成像５３を生成する。
【００７０】
（ステップＳ３７６）
ステップＳ３７５で生成したファントム再構成像５３の、左半分のファントム再構成半径（rA）、右半分のファントム再構成半径（rB）を算出する。ここで、再構成半径の算出には、例えば円曲線の方程式へのパラメータフィッテングの方法などを使用することができる。
【００７１】
（ステップＳ３７７）
ステップＳ３７６で算出した左右のファントム再構成半径（rA、rB）が等しいと見なせるかを判定する。再構成半径Raと再構成半径rBとの差がゼロと見なせれば現在のcenterの値を出力して終了する（Ｓ３７９）。ゼロと見なせなければステップＳ３７８へ移行する。
【００７２】
（ステップＳ３７８）
ステップＳ３７７で弧の半径がゼロと見なせなかった場合、図９の場合か図１０の場合であるかを判断し、それに応じてcenterの値が正しくなるように加減算して補正する。そして、ステップＳ３７５からの処理を再度行う。
【００７３】
次に図１２を用いて、検出器取付け補正角計算手段３８０が、検出器取付け補正角を決定する手順を詳細に説明する。
【００７４】
（ステップＳ３８０）
検出器取付け補正角計算処理を開始する（Ｓ３８０）。
【００７５】
（ステップＳ３８１）
検出器取付け補正角（以下、twistと称す）を初期位置（設計位置、例えば０°）に設定する。
【００７６】
（ステップＳ３８２）
前処理手段２０１は、ステップＳ３７２と同様、ステップＳ３６０で収集されたファントムＸ線画像データの前処理を行う。
【００７７】
（ステップＳ３８３）
画像歪み補正手段２０２は、ステップＳ３７３と同様、ステップＳ３８２で前処理を行ったファントムＸ線画像データに対して画像歪み補正処理を行う。
【００７８】
（ステップＳ３８４）
フィルタリング手段２０３は、ステップＳ３７４と同様、ステップＳ３８３で画像歪み補正処理を行ったファントムＸ線画像データについて、フィルタリング処理を行う。
【００７９】
（ステップＳ３８５）
逆投影手段２０４は、ステップＳ３８４でフィルタリング処理を行ったファントムＸ線画像データについて、検出器取付け補正角（twist）、及びステップＳ３８０で決定した回転中心軸投影位置（center）の値で、回転軌道面３０から離れた１つまたは複数の面での逆投影演算を行い、ファントム再構成像５３を生成する。
【００８０】
（ステップＳ３８６）
ステップＳ３８５で生成したファントム再構成像５３の、左半分のファントム再構成半径（rA）、右半分のファントム再構成半径（rB）を算出する。
【００８１】
（ステップＳ３８７）
ステップＳ３８７で算出した左右のファントム再構成半径（rA、rB）が等しいと見なせるかを判定する。ファントム再構成半径rAと、ファントム再構成半径rBとの差がゼロと見なせれば現在のtwistの値を出力して終了する（Ｓ３８９）。ゼロと見なせなければステップＳ３８８へ移行する。
【００８２】
（ステップＳ３８８）
ステップＳ３８７で弧の半径がゼロと見なせなかった場合、図９の場合か図１０の場合であるかを判断し、それに応じてcenterの値が正しくなるように加減算して補正する。そして、ステップＳ３８５からの処理を再度行う。
【００８３】
次に、本発明の第二の実施形態について図１３乃至図１５を用いて説明する。上記、第一の実施形態では、回転中心軸投影位置計算手段３７０が回転中心軸の投影位置を算出した後に、検出器取付け補正角計算手段３８０による補正角度の算出を実施するが、第二の実施形態は、回転中心軸投影位置計算手段３７０及び検出器取付け補正角計算手段３８０と同等の処理を一度に施すものである。図１３は、第二の実施の形態におけるＣ型アーム方式のコーンビームＸ線ＣＴ装置１の概略構成を示すブロック図である。図１４は、第二の実施の形態を説明するためのフローチャート、図１５は、図１４の処理から回転中心軸投影位置と検出器取付け補正角が求まる様子を示す図である。図１５左側の楕円図は、回転中心軸投影位置と検出器取付け角、両方の補正が正しく行われないで円柱ファントム５０を再構成したときに得られる再構成断面像、及び検出器取付け角ｂがゼロでないことにより再構成断面がスライス方向に変化していく様子を模式的に示す斜図である。この図では、３スライスの再構成断面像の斜図を示すが、回転中心軸投影位置と検出器取付け角、両方の補正が正しく行われない場合、左右アンバランスがスライス方向に連続的に変化していく円柱体として再構成される様子を示している。図１５の特別な場合として、検出器取付け角の補正のみが正しく行われている場合には、左右アンバランスがスライス位置によらず一定な再構成断面像が得られる。
【００８４】
図１３に示すように、第二の実施の形態に係るＸ線ＣＴ装置１の幾何学パラメータ計算手段３００は、第一の実施の形態の幾何学パラメータ３００における回転中心軸投影位置計算手段３７０、及び検出器取付け補正角計算手段３８０に代えて、回転中心軸投影位置及び検出器取付け補正角計算手段４７０を備える。その他の構成は、第一の実施の形態に係るＸ線ＣＴ装置と同様である。
【００８５】
図１４に基づいて、第二の実施形態における、回転中心軸投影位置及び検出器取付け補正角計算処理手段４７０について説明する。
【００８６】
（ステップＳ４７０）
回転中心軸投影位置及び検出器取付け補正角計算処理を開始する（Ｓ４７０）。
【００８７】
（ステップＳ４７１）
ステップＳ３７１と同様、回転軌道面での回転中心軸投影位置ＣＭを初期位置（設計位置）に設定する。
【００８８】
（ステップＳ４７２）
前処理手段２０１は、ステップＳ３７２と同様、ステップＳ３６０で収集されたファントムＸ線画像データの前処理を行う。
【００８９】
（ステップＳ４７３）
画像歪み補正手段２０２は、ステップＳ３７３と同様、ステップＳ４７２で前処理を行ったファントムＸ線画像データに対して画像歪み補正処理を行う。
【００９０】
（ステップＳ４７４）
フィルタリング手段２０３は、ステップＳ３７４と同様、ステップＳ４７３で画像歪み補正処理を行ったファントムＸ線画像データについて、フィルタリング処理を行う。
【００９１】
（ステップＳ４７５）
逆投影手段２０４は、ステップＳ４７４でフィルタリング処理を行ったファントムＸ線画像データについて、回転軌道面に平行な複数個の断面（回転軌道面からの距離をＤとする）について逆投影演算を行い、ファントム再構成像５３を生成する。
【００９２】
（ステップＳ４７６）
ステップＳ４７５で生成した複数個のファントム再構成像５３の、左半分のファントム再構成半径（rA）、右半分のファントム再構成半径（rB）を算出する。再構成半径の算出には、第一の実施形態と同様、円曲線の方程式へのパラメータフィッテングの方法などを使用することができる。
【００９３】
（ステップＳ４７７）
ステップＳ４７６で算出した左右のファントム再構成半径（rA、rB）から、各再構成断面について回転中心軸投影位置の補正量Ｃ（Ｄ）を算出する。
【００９４】
（ステップＳ４７８）
ステップＳ４７７で算出した回転中心投影位置の補正量Ｃ（Ｄ）を用い、図１５に示すように、回転中心軸投影位置ＣＭと検出器取付け角ｂを関係式：Ｃ（Ｄ）＝ＣＭ＋Ｄ・tanｂへのフィッティングにより決定する。決定した回転中心軸投影位置ＣＭと検出器取付け角ｂを出力して終了する（Ｓ４７９）。
【００９５】
以上説明した幾何学パラメータ処理計算手段３００により、３次元的Ｘ線ＣＴ像を生成し表示するコーンビームＸ線ＣＴ装置１において、鮮明な３次元的Ｘ線ＣＴ像を生成するための幾何学パラメータ、回転中心軸３１の投影位置、検出器取付け角２４を求めることができる。
【００９６】
以上説明したように、本実施の形態の幾何学パラメータ計算手段を搭載する事により、Ｃ型アーム等の回転撮影から収集される２次元Ｘ線像から３次元的Ｘ線ＣＴ像を生成するコーンビームＸ線ＣＴ装置において、回転中心軸及び検出器取付け角が正しい設計値から変位した場合にも、鮮明な３次元的Ｘ線ＣＴ像を生成し表示することができる。これにより、頭部、腹部等の造影撮影、並びに歯顎、腰椎、四肢等の整形分野の診断性能を向上させることができる。
【００９７】
なお本発明は上記実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々に変形することができる。例えば、上記実施の形態では、２次元Ｘ線検出器としてＸ線イメージインテンシファイアとテレビカメラの組み合わせで形成するようにしていたが、これに代えて、例えばフラットパネルディテクター（ＦＰＤ）を用いて２次元Ｘ線検出器を構成するようにしてもよい。また本発明はＣ型アーム方式以外のコーンビームＸ線ＣＴ装置にも適用可能であるし、回転撮影の全角度が２００度以外の場合にも勿論、適用可能である。
【図面の簡単な説明】
【００９８】
【図１】本発明を適用したＣ型アーム方式のコーンビームＸ線ＣＴ装置１の概略構成を示すブロック図である。
【図２】画像歪み補正テーブルの生成に使用されるホールチャート１８を示す概念図である。
【図３】ホールチャート１８を２次元Ｘ線検出器のＸ線入射面に取り付ける際に生じる検出器取付け角を説明するための図である。
【図４】円柱ファントム５０が、ファントム保持体５０ａを用いて寝台１７に載置され、回転撮影される状態を示す概念図である。
【図５】図４のファントムが、回転撮影される状態を示す概念図である。
【図６】本発明の実施の形態における幾何学パラメータ計算手段３００が、幾何学パラメータ計算処理を行う手順を示すフローチャートである。
【図７】回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の前半の投影データから再構成される円柱ファントムの再構成断面について示した図である。
【図８】回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の後半の投影データから再構成される円柱ファントムの再構成断面について示した図である。
【図９】回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面について示した図である。
【図１０】回転中心軸の投影パラメータが−ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面について示した図である。
【図１１】本発明の実施の形態における回転中心軸投影位置計算手段３７０が、回転中心軸投影位置計算処理を行う手順を示すフローチャートである。
【図１２】本発明の実施の形態における検出器取付け補正角計算手段３８０が、検出器取付け補正角計算処理を行う手順を示すフローチャートである。
【図１３】本発明の第二の実施の形態におけるＣ型アーム方式のコーンビームＸ線ＣＴ装置１の概略構成を示すブロック図である。
【図１４】本発明の第二の実施形態における回転中心軸投影位置及び検出器取付け補正角計算手段４７０が、回転中心軸投影位置及び検出器取付け補正角計算処理を行う手順を示すフローチャートである。
【図１５】回転中心軸投影位置及び検出器取付け補正角計算処理において、回転中心軸投影位置ＣＭと検出器取付け角ｂが、関係式：Ｃ（Ｄ）＝ＣＭ＋Ｄ・tanｂへのフィッティングにより求まる様子を示す図である。
【符号の説明】
【００９９】
１…コーンビームＸ線ＣＴ装置、１０…撮影部、１１…Ｘ線源、１１ｔ…Ｘ線管、１１ｃ…コリメータ、１２…２次元Ｘ線検出器、１２ｉ…Ｘ線イメージインテンシファイア、１２ｃ…テレビカメラ、１３…Ｃ型アーム、１４…Ｃ型アーム保持体、１５…天井支持体、１６…天井レール、１７…寝台、１８…ホールチャート、１８ｈ…ホール、２０…制御演算部、２１…ホールチャート１８に固定された横軸、２２…ホールチャート１８に固定された縦軸、２４…検出器取付け角、３０…回転軌道面（ミッドプレーン）、３１…回転中心軸、３２…回転中心軸の投影、４０…被検体、５０…円柱ファントム、５０ａ…ファントム保持体、５１…幾何学的に正しい回転中心軸で計測されるファントム端の投影位置、５２…現在使用している回転中心軸で計測されるファントム端の投影位置、５３…ファントム再構成像、５４…回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の前半の投影データから再構成される円柱ファントムの再構成断面、５５…回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、回転の後半の投影データから再構成される円柱ファントムの再構成断面、５６…回転中心軸の投影パラメータが＋ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面、５７…回転中心軸の投影パラメータが−ｕ方向に誤差を持つ場合に、全回転撮影の投影データから再構成される円柱ファントムの再構成断面、８０…表示装置、１００…撮影部制御手段、１０１…撮影系回転制御手段、１０２…撮影系位置制御手段、１０３…Ｘ線照射制御手段、１０４…寝台制御手段、１０５…検出系制御手段、１１０…画像収集手段、２００…再構成手段、２０１…前処理手段、２０２…画像歪み補正手段、２０３…フィルタリング手段、２０４…逆投影手段、２１０…画像表示手段、３００…幾何学パラメータ計算手段、３２０…画像歪み補正テーブル生成手段、３３０…画像歪み補正テーブル格納手段、３５０…回転軌道面算出手段、３７０…回転中心軸投影位置計算手段、３８０…検出器取付け補正角計算手段
【図１】
2006130060.tif 000003
【図２】
2006130060.tif 000004
【図３】
2006130060.tif 000005
【図４】
2006130060.tif 000006
【図５】
2006130060.tif 000007
【図６】
2006130060.tif 000008
【図７】
2006130060.tif 000009
【図８】
2006130060.tif 000010
【図９】
2006130060.tif 000011
【図１０】
2006130060.tif 000012
【図１１】
2006130060.tif 000013
【図１２】
2006130060.tif 000014
【図１３】
2006130060.tif 000015
【図１４】
2006130060.tif 000016
【図１５】
2006130060.tif 000017