(57)【要約】【課題】  ヘリカルスキャンにおける直線移動距離全体のうちの加・減速のための直線移動距離分も画像再構成に利用できるようにする。【解決手段】  ヘリカルスキャンを行う時に直線移動の加・減速中でも投影データを収集し、収集した投影データを画像再構成に利用する。また、直線移動の加速中には管電流を増加させながら投影データを収集し、直線移動の減速中には管電流を減少させながら投影データを収集する。【効果】  ヘリカルスキャンにおける直線移動距離全体のうちの加・減速のための直線移動距離分も画像再構成に利用できる。直線移動の加・減速中に収集した投影データの直線移動方向におけるＸ線密度を一定に出来る。【選択図】  図５
2005040582.tif 000002
【特許請求の範囲】
【請求項１】
ヘリカルスキャンを行う時にテーブルの直線移動の加・減速中でも投影データを収集し、収集した投影データを画像再構成に利用することを特徴とするＸ線ＣＴ撮像方法。
【請求項２】
ヘリカルスキャンを行う時にテーブルの直線移動の加・減速中でも投影データを収集し、スキャン中の体軸方向（以下ｚ軸という）の座標情報を各ビューデータもしくは数ビューデータごとに付加し、または別情報として保持し、収集した投影データを各ビューごとに若しくは数ビューごとに同期したｚ軸座標情報とともに画像再構成に利用することを特徴とするＸ線ＣＴ撮像方法。
【請求項３】
請求項１または請求項２に記載のＸ線ＣＴ撮像方法において、投影データの収集と並行して画像再構成を進めることを特徴とするＸ線ＣＴ撮像方法。
【請求項４】
請求項３に記載のＸ線ＣＴ撮像方法において、あるビューの投影データが収集される前に、前記投影データを画像再構成に用いるためのパラメータを予め予測して保持しておくこと又は投影データ収集時に予測しながら得ることを特徴とするＸ線ＣＴ撮像方法。
【請求項５】
請求項４に記載のＸ線ＣＴ撮像方法において、テーブルの位置の変化を表す直線移動情報を予め保持し、あるビューの投影データが収集される前に、該投影データに対応するｚ軸座標を前記直線移動情報から推定し、該推定したｚ軸座標を基に前記投影データを画像再構成に用いるためのパラメータを算出しておくことを特徴とするＸ線ＣＴ撮像方法。
【請求項６】
請求項１から請求項５のいずれかに記載のＸ線ＣＴ撮像方法において、テーブルの直線移動の加速中には管電流を増加させながら投影データを収集し、直線移動の減速中には管電流を減少させながら投影データを収集することを特徴とするＸ線ＣＴ撮像方法。
【請求項７】
請求項１から請求項６のいずれかに記載のＸ線ＣＴ撮像方法において、直線移動の加・減速を時間に対してリニアに行うことを特徴とするＸ線ＣＴ撮像方法。
【請求項８】
請求項１から請求項６のいずれかに記載のＸ線ＣＴ撮像方法において、直線移動の加・減速を時間に対してノンリニアに行うことを特徴とするＸ線ＣＴ撮像方法。
【請求項９】
請求項１から請求項８のいずれかに記載のＸ線ＣＴ撮像方法において、マルチ検出器により投影データを収集することを特徴とするＸ線ＣＴ撮像方法。
【請求項１０】
請求項９に記載のＸ線ＣＴ撮像方法において、画像再構成面をｘ軸およびｙ軸に平行なｘｙ面とし、マルチ検出器の検出器列が並ぶ方向をｚ軸方向とするとき、あるビュー位置でのマルチ検出器のｚ軸方向の中心を通るｘｙ面から画像再構成面までのｚ軸方向の距離と画像再構成面上の画素点の位置とに基づいて、前記ビューでの前記画素点の画素値を求めるのに使用する投影データを抽出することを特徴とするＸ線ＣＴ撮像方法。
【請求項１１】
請求項９または請求項１０に記載のＸ線ＣＴ撮像方法において、３次元画像再構成法により画像再構成を行うことを特徴とするＸ線ＣＴ撮像方法。
【請求項１２】
請求項１１に記載のＸ線ＣＴ撮像方法において、前記３次元画像再構成法が、収集した投影データを各ビューのｚ軸方向位置に基づいて配置し画像再構成領域上の１本のラインまたは複数画素間隔あけた複数本の平行なラインに対応する投影データを抽出し、前記投影データにコーンビーム再構成荷重を乗算して投影ライン・データを作成し、前記投影ライン・データにフィルタ処理を施して画像各位置ライン・データを作成し、前記画像各位置ライン・データを基に再構成領域上の各画素の逆投影画素データを求め、画像再構成に用いる全ビューの逆投影画素データを画素対応に加算して逆投影データを求める３次元逆投影方法であることを特徴とするＸ線ＣＴ撮像方法。
【請求項１３】
Ｘ線管と、Ｘ線検出器と、前記Ｘ線管または前記Ｘ線検出器の少なくとも一方を撮影対象の周りに回転させると共に両方を撮影対象に対して直線状に相対移動しながら且つ直線移動の加・減速中でも投影データを収集するスキャン手段と、収集した投影データを基にＣＴ画像を生成する画像再構成手段とを具備したことを特徴とするＸ線ＣＴ装置。
【請求項１４】
Ｘ線管と、Ｘ線検出器と、前記Ｘ線管または前記Ｘ線検出器の少なくとも一方を撮影対象の周りに回転させると共に両方を撮影対象に対して直線状に相対移動しながら且つ直線移動の加・減速中でも投影データを収集するとともにスキャン中の体軸方向（以下ｚ軸という）の座標情報を各ビューデータもしくは数ビューデータごとに付加するか又は別情報として保持するスキャン手段と、収集した投影データと各ビューごとに若しくは数ビューごとに同期したｚ軸座標情報とを基にＣＴ画像を生成する画像再構成手段とを具備したことを特徴とするＸ線ＣＴ装置。
【請求項１５】
請求項１３または請求項１４に記載のＸ線ＣＴ装置において、前記スキャン手段による投影データの収集と並行して前記画像再構成手段による画像再構成を進めることを特徴とするＸ線ＣＴ装置。
【請求項１６】
請求項１５に記載のＸ線ＣＴ装置において、あるビューの投影データが収集される前に、前記投影データを画像再構成に用いるためのパラメータを予め予測して保持しておくか又は投影データ収集時に予測しながら得て保持しておくパラメータ保持手段を具備したことを特徴とするＸ線ＣＴ装置。
【請求項１７】
請求項１６に記載のＸ線ＣＴ装置において、前記直線移動の位置の変化を表す直線移動情報を予め保持する直線移動情報保持手段と、あるビューの投影データが収集される前に該投影データに対応するｚ軸座標を前記直線移動情報から推定し該推定したｚ軸座標を基に前記投影データを画像再構成に用いるためのパラメータを算出しておくパラメータ推定手段を具備したことを特徴とするＸ線ＣＴ装置。
【請求項１８】
請求項１３から請求項１７のいずれかに記載のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加速中には管電流を増加させながら投影データを収集し、直線移動の減速中には管電流を減少させながら投影データを収集することを特徴とするＸ線ＣＴ装置。
【請求項１９】
請求項１３から請求項１８のいずれかに記載のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加・減速を時間に対してリニアに行うことを特徴とするＸ線ＣＴ装置。
【請求項２０】
請求項１３から請求項１８のいずれかに記載のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加・減速を時間に対してノンリニアに行うことを特徴とするＸ線ＣＴ装置。
【請求項２１】
請求項１３から請求項２０のいずれかに記載のＸ線ＣＴ撮像方法において、前記Ｘ線検出器が、マルチ検出器であることを特徴とするＸ線ＣＴ装置。
【請求項２２】
請求項２１に記載のＸ線ＣＴ装置において、前記画像再構成手段は、画像再構成面をｘ軸およびｙ軸に平行なｘｙ面とし、マルチ検出器の検出器列が並ぶ方向をｚ軸方向とするとき、あるビュー位置でのマルチ検出器のｚ軸方向の中心を通るｘｙ面から画像再構成面までのｚ軸方向の距離と画像再構成面上の画素点の位置とに基づいて、前記ビューでの前記画素点の画素値を求めるのに使用する投影データを抽出することを特徴とするＸ線ＣＴ装置。
【請求項２３】
請求項２１または請求項２２に記載のＸ線ＣＴ装置において、前記画像再構成手段は、３次元画像再構成法により画像再構成を行うことを特徴とするＸ線ＣＴ装置。
【請求項２４】
請求項２３に記載のＸ線ＣＴ装置において、前記３次元画像再構成法が、収集した投影データを各ビューのｚ軸方向位置に基づいて配置し画像再構成領域上の１本のラインまたは複数画素間隔あけた複数本の平行なラインに対応する投影データを抽出し、前記投影データにコーンビーム再構成荷重を乗算して投影ライン・データを作成し、前記投影ライン・データにフィルタ処理を施して画像各位置ライン・データを作成し、前記画像各位置ライン・データを基に再構成領域上の各画素の逆投影画素データを求め、画像再構成に用いる全ビューの逆投影画素データを画素対応に加算して逆投影データを求める３次元逆投影方法であることを特徴とするＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（Computed Tomography）撮像方法およびＸ線ＣＴ装置に関し、更に詳しくは、ヘリカルスキャン（helical scan）におけるテーブルの直線移動距離全体のうちの加・減速のための直線移動距離分も画像再構成に利用できるようにしたＸ線ＣＴ撮像方法およびＸ線ＣＴ装置に関する。
【背景技術】
【０００２】
ヘリカルスキャンを行う場合、Ｘ線管およびＸ線検出器を撮影対象の周りに回転させると共に撮影対象を載せたテーブルを直線移動させている。この直線移動では、静止しているテーブルを所定速度まで加速し、投影データを収集したい範囲でテーブルを所定速度に維持し、投影データの収集を終わるとテーブルを減速し、静止させている。なお、ある撮影部位では低い所定速度Ｖ１とし、別の撮影部位では高い所定速度Ｖ２とするように、異なる所定速度を設定する場合もある（特許文献１参照。）。
【０００３】
【特許文献１】特開平１０−３１４１６２号公報（［００４９］〜［００５１］、図５）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
従来、画像再構成に利用する投影データの収集は、直線移動速度が一定に維持されている間だけ行われ、直線移動の加・減速中には行われなかった。
つまり、従来のＸ線ＣＴ装置では、全体の直線移動距離のうちの加・減速のための直線移動距離分が画像再構成に利用されていない無駄があった。
そこで、本発明の目的は、全体の直線移動距離のうちの加・減速のための直線移動距離分も画像再構成に利用できるようにしたＸ線ＣＴ撮像方法およびＸ線ＣＴ装置を提供することにある。
【課題を解決するための手段】
【０００５】
第１の観点では、本発明は、ヘリカルスキャンを行う時にテーブルの直線移動の加・減速中でも投影データを収集し、収集した投影データを画像再構成に利用することを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第１の観点によるＸ線ＣＴ撮像方法では、直線移動速度が一定に維持されている間だけでなく、直線移動の加・減速中にも投影データを収集し、収集した投影データを画像再構成に利用する。このため、全体の直線移動距離のうちの加・減速のための直線移動距離分も画像再構成に利用できるようになる。
なお、画像再構成は、２次元画像再構成法でもよいし、３次元画像再構成法でもよい。
【０００６】
第２の観点では、本発明は、ヘリカルスキャンを行う時にテーブルの直線移動の加・減速中でも投影データを収集し、スキャン中の体軸方向（以下ｚ軸という）の座標情報を各ビューデータもしくは数ビューデータごとに付加し、または別情報として保持し、収集した投影データを各ビューごとに若しくは数ビューごとに同期したｚ軸座標情報とともに画像再構成に利用することを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第２の観点によるＸ線ＣＴ撮像方法では、直線移動速度が一定に維持されている間だけでなく、直線移動の加・減速中にも投影データを収集し、収集した投影データをｚ軸座標情報とともに画像再構成に利用する。このため、全体の直線移動距離のうちの加・減速のための直線移動距離分も画像再構成に利用できるようになる。
なお、画像再構成は、２次元画像再構成法でもよいし、３次元画像再構成法でもよい。
【０００７】
第３の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、投影データの収集と並行して画像再構成を進めることを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第３の観点によるＸ線ＣＴ撮像方法では、投影データの収集と並行して画像再構成を進めるため、画像が得られるまでのタイムラグ（時間遅れ）を小さくすることが出来る。
【０００８】
第４の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、あるビューの投影データが収集される前に、前記投影データを画像再構成に用いるためのパラメータを予め予測して保持しておくこと又は投影データ収集時に予測しながら得ることを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第４の観点によるＸ線ＣＴ撮像方法では、あるビューの投影データが収集される前に該投影データを画像再構成に用いるためのパラメータを予め保持しておくため、投影データが収集された後、迅速に画像再構成を進めることが出来る。
【０００９】
第５の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、テーブルの位置の変化を表す直線移動情報を予め保持し、あるビューの投影データが収集される前に、該投影データに対応するｚ軸座標を前記直線移動情報から推定し、該推定したｚ軸座標を基に前記投影データを画像再構成に用いるためのパラメータを算出しておくことを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第５の観点によるＸ線ＣＴ撮像方法では、テーブルの位置の変化を表す直線移動情報を保持し、あるビューの投影データが収集される前に、あるビューの投影データを画像再構成に用いるためのパラメータを算出するため、投影データが収集された後、迅速に画像再構成を進めることが出来る。
【００１０】
第６の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、テーブルの直線移動の加速中には管電流を増加させながら投影データを収集し、直線移動の減速中には管電流を減少させながら投影データを収集することを特徴とするＸ線ＣＴ撮像方法を提供する。
加速中は単位時間当たりの直線移動距離が次第に長くなっていくため、管電流が一定であると、直線移動方向におけるＸ線密度が次第に小さくなっていく。一方、減速中は単位時間当たりの直線移動距離が次第に短くなっていくため、管電流が一定であると、直線移動方向におけるＸ線密度が次第に大きくなっていく。換言すれば、収集した投影データに対応する実質的な管電流がビューによって異なることになり、前処理が煩雑になってしまう。
そこで、上記第６の観点によるＸ線ＣＴ撮像方法では、直線移動の加速中には管電流を増加させながら投影データを収集する。また、直線移動の減速中には管電流を減少させながら投影データを収集する。これにより、直線移動方向におけるＸ線密度は一定になる。換言すれば、直線移動の加・減速中に収集した投影データに対応する実質的な管電流をビューによらずに一定にでき、前処理を簡単化できる。
【００１１】
第７の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、直線移動の加・減速を時間に対してリニアに行うことを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第７の観点によるＸ線ＣＴ撮像方法では、直線移動の加速および減速を時間に対してリニアに行うため、制御が簡単になる。
【００１２】
第８の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、直線移動の加・減速を時間に対してノンリニアに行うことを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第８の観点によるＸ線ＣＴ撮像方法では、直線移動の加速および減速を時間に対してノンリニアに行うため、直線移動速度の変化を滑らかに出来る。
【００１３】
第９の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、マルチ検出器により投影データを収集することを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第９の観点によるＸ線ＣＴ撮像方法では、マルチ検出器により一度に多くの投影データを収集できる。
【００１４】
第１０の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、画像再構成面をｘ軸およびｙ軸に平行なｘｙ面とし、マルチ検出器の検出器列が並ぶ方向をｚ軸方向とするとき、あるビュー位置でのマルチ検出器のｚ軸方向の中心を通るｘｙ面から画像再構成面までのｚ軸方向の距離と画像再構成面上の画素点の位置とに基づいて、前記ビューでの前記画素点の画素値を求めるのに使用する投影データを抽出することを特徴とするＸ線ＣＴ撮像方法を提供する。
従来の画像再構成法は、直線移動速度が一定であることを前提にしている。このため、従来の画像再構成法を直線移動の加・減速中に収集した投影データにそのまま適用すると、アーチファクトを生じてしまう。
そこで、上記第１０の観点によるＸ線ＣＴ撮像方法では、あるビュー位置でのマルチ検出器のｚ軸方向の中心を通るｘｙ面から画像再構成面までのｚ軸方向の距離と画像再構成面上の画素点ｇの位置とに基づいて前記ビューでの前記画素点ｇの画素値を求めるのに使用する投影データを抽出するので、直線移動の加・減速中に収集した投影データでも適正に必要な投影データを抽出でき、アーチファクトを防止できる。
【００１５】
第１１の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、３次元画像再構成法により画像再構成を行うことを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第１１の観点によるＸ線ＣＴ撮像方法では、コーン角の大きなマルチ検出器により投影データを収集するが、３次元画像再構成法により画像再構成を行うため、大きなコーン角に起因するアーチファクトを防止できる。
なお、３次元画像再構成法としては、フェルドカンプ（Feldkamp）法もしくは重み付きフェルドカンプ法などがある。
【００１６】
第１２の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、前記３次元画像再構成法が、収集した投影データを各ビューのｚ軸方向位置に基づいて配置し画像再構成領域上の１本のラインまたは複数画素間隔あけた複数本の平行なラインに対応する投影データを抽出し、前記投影データにコーンビーム再構成荷重を乗算して投影ライン・データを作成し、前記投影ライン・データにフィルタ処理を施して画像各位置ライン・データを作成し、前記画像各位置ライン・データを基に再構成領域上の各画素の逆投影画素データを求め、画像再構成に用いる全ビューの逆投影画素データを画素対応に加算して逆投影データを求める３次元逆投影方法であることを特徴とするＸ線ＣＴ撮像方法を提供する。
上記第１２の観点によるＸ線ＣＴ撮像方法では、特願２００２−１４７２３１号および特願２００２−３３８９４７号で提案されている３次元画像再構成法を行うため、演算量を大幅に減らすことが出来る。
【００１７】
第１３の観点では、本発明は、Ｘ線管と、Ｘ線検出器と、前記Ｘ線管または前記Ｘ線検出器の少なくとも一方を撮影対象の周りに回転させると共に両方を撮影対象に対して直線状に相対移動しながら且つ直線移動の加・減速中でも投影データを収集するスキャン手段と、収集した投影データを基にＣＴ画像を生成する画像再構成手段とを具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１３の観点によるＸ線ＣＴ装置では、上記第１の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００１８】
第１４の観点では、本発明は、Ｘ線管と、Ｘ線検出器と、前記Ｘ線管または前記Ｘ線検出器の少なくとも一方を撮影対象の周りに回転させると共に両方を撮影対象に対して直線状に相対移動しながら且つ直線移動の加・減速中でも投影データを収集するとともにスキャン中の体軸方向（以下ｚ軸という）の座標情報を各ビューデータもしくは数ビューデータごとに付加するか又は別情報として保持するスキャン手段と、収集した投影データと各ビューごとに若しくは数ビューごとに同期したｚ軸座標情報とを基にＣＴ画像を生成する画像再構成手段とを具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１４の観点によるＸ線ＣＴ装置では、上記第２の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００１９】
第１５の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記スキャン手段による投影データの収集と並行して前記画像再構成手段による画像再構成を進めることを特徴とするＸ線ＣＴ装置を提供する。
上記第１５の観点によるＸ線ＣＴ装置では、上記第３の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２０】
第１６の観点では、本発明は、上記構成のＸ線ＣＴ装置において、あるビューの投影データが収集される前に、前記投影データを画像再構成に用いるためのパラメータを予め予測して保持しておくか又は投影データ収集時に予測しながら得て保持しておくパラメータ保持手段を具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１６の観点によるＸ線ＣＴ装置では、上記第４の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２１】
第１７の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記直線移動の位置の変化を表す直線移動情報を予め保持する直線移動情報保持手段と、あるビューの投影データが収集される前に該投影データに対応するｚ軸座標を前記直線移動情報から推定し該推定したｚ軸座標を基に前記投影データを画像再構成に用いるためのパラメータを算出しておくパラメータ推定手段を具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１７の観点によるＸ線ＣＴ装置では、上記第５の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２２】
第１８の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加速中には管電流を増加させながら投影データを収集し、直線移動の減速中には管電流を減少させながら投影データを収集することを特徴とするＸ線ＣＴ装置を提供する。
上記第１８の観点によるＸ線ＣＴ装置では、上記第６の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２３】
第１９の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加・減速を時間に対してリニアに行うことを特徴とするＸ線ＣＴ装置を提供する。
上記第１９の観点によるＸ線ＣＴ装置では、上記第７の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２４】
第２０の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記スキャン手段は、直線移動の加・減速を時間に対してノンリニアに行うことを特徴とするＸ線ＣＴ装置を提供する。
上記第２０の観点によるＸ線ＣＴ装置では、上記第８の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２５】
第２１の観点では、本発明は、上記構成のＸ線ＣＴ撮像方法において、前記Ｘ線検出器が、マルチ検出器であることを特徴とするＸ線ＣＴ装置を提供する。
上記第２１の観点によるＸ線ＣＴ装置では、上記第９の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２６】
第２２の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記画像再構成手段は、画像再構成面をｘ軸およびｙ軸に平行なｘｙ面とし、マルチ検出器の検出器列が並ぶ方向をｚ軸方向とするとき、あるビュー位置でのマルチ検出器のｚ軸方向の中心を通るｘｙ面から画像再構成面までのｚ軸方向の距離と画像再構成面上の画素点の位置とに基づいて、前記ビューでの前記画素点の画素値を求めるのに使用する投影データを抽出することを特徴とするＸ線ＣＴ装置を提供する。
上記第２２の観点によるＸ線ＣＴ装置では、上記第１０の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２７】
第２３の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記画像再構成手段は、３次元画像再構成法により画像再構成を行うことを特徴とするＸ線ＣＴ装置を提供する。
上記第２３の観点によるＸ線ＣＴ装置では、上記第１１の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【００２８】
第２４の観点では、本発明は、上記構成のＸ線ＣＴ装置において、前記３次元画像再構成法が、収集した投影データを各ビューのｚ軸方向位置に基づいて配置し画像再構成領域上の１本のラインまたは複数画素間隔あけた複数本の平行なラインに対応する投影データを抽出し、前記投影データにコーンビーム再構成荷重を乗算して投影ライン・データを作成し、前記投影ライン・データにフィルタ処理を施して画像各位置ライン・データを作成し、前記画像各位置ライン・データを基に再構成領域上の各画素の逆投影画素データを求め、画像再構成に用いる全ビューの逆投影画素データを画素対応に加算して逆投影データを求める３次元逆投影方法であることを特徴とするＸ線ＣＴ装置を提供する。
上記第２４の観点によるＸ線ＣＴ装置では、上記第１２の観点によるＸ線ＣＴ撮像方法を好適に実施できる。
【発明の効果】
【００２９】
本発明のＸ線ＣＴ撮像方法およびＸ線ＣＴ装置によれば、ヘリカルスキャンにおける直線移動距離全体のうちの加・減速のための直線移動距離分も画像再構成に利用できるようになる。
【発明を実施するための最良の形態】
【００３０】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
【実施例１】
【００３１】
図１は、実施例１にかかるＸ線ＣＴ装置１００を示す構成ブロック図である。
このＸ線ＣＴ装置１００は、操作コンソール１と、寝台装置１０と、走査ガントリ２０とを具備している。
【００３２】
操作コンソール１は、操作者の入力を受け付ける入力装置２と、画像再構成処理などを実行する中央処理装置３と、走査ガントリ２０で取得した投影データを収集するデータ収集バッファ５と、投影データから再構成したＣＴ画像を表示するＣＲＴ６と、プログラムやデータやＸ線ＣＴ画像を記憶する記憶装置７とを具備している。
【００３３】
テーブル装置１０は、被検体を乗せて走査ガントリ２０のボア（空洞部）に入れ出しするテーブル１２を具備している。テーブル１２は、テーブル装置１０に内蔵するモータで昇降および直線移動される。
【００３４】
走査ガントリ２０は、Ｘ線管２１と、Ｘ線コントローラ２２と、コリメータ２３と、マルチ検出器２４と、ＤＡＳ（Data Acquisition System）２５と、被検体の体軸の回りにＸ線管２１などを回転させる回転コントローラ２６と、制御信号などを前記操作コンソール１や寝台装置１０とやり取りする制御コントローラ２９と、スリップリング３０とを具備している。
【００３５】
図２および図３は、Ｘ線管２１とマルチ検出器２４の説明図である。
Ｘ線管２１とマルチ検出器２４は、回転中心ＩＣの回りを回転する。鉛直方向をｙ方向とし、水平方向をｘ方向とし、これらに垂直な方向をｚ方向とするとき、Ｘ線管２１およびマルチ検出器２４の回転平面は、ｘｙ面である。また、テーブル１２の移動方向は、ｚ方向である。
Ｘ線管２１は、コーンビームＣＢと呼ばれるＸ線ビームを発生する。コーンビームＣＢの中心軸方向がｙ方向に平行なときを、view＝０゜とする。
マルチ検出器２４は、例えば２５６列の検出器列を有する。また、各検出器列は、例えば１０２４チャネルのチャネルを有する。
【００３６】
図４は、Ｘ線ＣＴ装置１００の動作の概略を示すフロー図である。
ステップＳ１では、Ｘ線管２１とマルチ検出器２４とを撮影対象の周りに回転させ且つテーブル１２を直線移動させながら直線移動位置ｚとビュー角度viewと検出器列番号ｊとチャネル番号ｉとで表わされる投影データＤ０(z,view,j,i)を収集する。直線移動位置ｚは、エンコーダによりｚ軸方向位置パルスがカウントされ、制御コントローラ２９にてｚ軸座標に変換され、スリップリング３０を経由して、ＤＡＳ２５の投影データにｚ軸座標情報として付加される。
図５に、ｚ軸座標情報を付加した、あるビュー角度viewのデータのフォーマットを示す。
なお、このデータ収集処理については、図６〜図１４を参照して後述する。
【００３７】
ステップＳ２では、投影データＤ０(z,view,j,i)に対して、前処理（オフセット補正，対数補正，Ｘ線線量補正，感度補正）を行う。
ステップＳ３では、前処理した投影データＤ０(z,view,j,i)に対して、フィルタ処理を行う。すなわち、フーリエ変換し、フィルタ（再構成関数）を掛け、逆フーリエ変換する。
ステップＳ４では、フィルタ処理した投影データＤ０(z,view,j,i)に対して、３次元逆投影処理を行い、逆投影データＤ３(x,y)を求める。この３次元逆投影処理については、図１５を参照して後述する。
ステップＳ５では、逆投影データＤ３(x,y)に対して後処理を行い、ＣＴ画像を得る。
【００３８】
図６は、データ収集処理（図４のステップＳ１）の詳細を示すフロー図である。
ステップＡ１では、Ｘ線管２１とマルチ検出器２４とを撮影対象の周りに回転させる。
ステップＡ２では、図７や図９に示す直線移動開始位置までテーブル１２を低速で直線移動する。
ステップＡ３では、テーブル１２の直線移動をスタートさせる。
【００３９】
ステップＡ４では、テーブル１２の直線移動速度を所定関数に基づいて加速し、それに合わせて管電流を増加させる。所定関数が時間に対してリニアの場合を図７，図８に示し、所定関数が時間に対してノンリニアの場合を図９，図１０に示す。直線移動方向のＸ線密度つまり単位厚さ当たりのＸ線量は「管電流／直線移動速度」に比例する。従って、直線移動速度の増加に合わせて管電流を増加させることで「管電流／直線移動速度」を一定にすることができ、これにより加速中でも直線移動方向のＸ線密度を一定にすることが出来る。
なお、所定関数は、計算式形式で保持しておいてもよいし、表形式で保持しておいてもよい。
【００４０】
ステップＡ５では、加速中の投影データＤ０(z,view,j,i)を収集する。
ステップＡ６では、テーブル１２の直線移動速度が図７や図９に示す所定速度Ｖｃに達したらステップＡ７へ進み、所定速度Ｖｃに達してなかったらステップＡ４に戻ってさらに加速する。
【００４１】
ステップＡ７では、テーブル１２の直線移動速度を所定速度に維持した状態で定速の投影データＤ０(z,view,j,i)を収集する。
ステップＡ８では、テーブル１２が図７や図９に示す定速終了位置に達したらステップＡ９へ進み、定速終了位置に達してなかったらステップＡ７に戻って定速の投影データ収集を継続する。
【００４２】
ステップＡ９では、テーブル１２の直線移動速度を所定関数に基づいて減速し、それに合わせて管電流を減少させる。所定関数が時間に対してリニアの場合を図７，図８に示し、所定関数が時間に対してノンリニアの場合を図９，図１０に示す。直線移動方向のＸ線密度つまり単位厚さ当たりのＸ線量は「管電流／直線移動速度」に比例する。従って、直線移動速度の減少に合わせて管電流を減少させることで「管電流／直線移動速度」を一定にすることができ、これにより減速中でも直線移動方向のＸ線密度を一定にすることが出来る。
ステップＡ１０では、減速中の投影データＤ０(z,view,j,i)を収集する。
ステップＡ１１では、テーブル１２の直線移動速度が図７や図９に示す停止可能速度に達したらステップＡ１２へ進み、停止可能速度に達してなかったらステップＡ９に戻ってさらに減速する。
【００４３】
ステップＡ１２では、テーブル１２の直線移動をストップさせる。
【００４４】
なお、図１１〜図１４に示すように、定速開始位置＝定速終了位置に設定すれば、最も短い直線移動距離で投影データＤ０(z,view,j,i)を収集できる。
【００４５】
図１５は、３次元逆投影処理（図４のステップＳ４）の詳細を示すフロー図である。
【００４６】
ステップＲ１では、ＣＴ画像の再構成に必要な全ビュー（すなわち、３６０゜分のビュー又は「１８０゜分＋ファン角度分」のビュー）中の一つのビューに着目する。
【００４７】
ステップＲ２では、着目ビューの投影データＤ０(z,view,j,i)の中から再構成領域上の複数画素間隔あけた複数の平行なラインに対応する投影データＤｒを抽出する。
【００４８】
図１６に再構成領域Ｐ上の複数の平行なラインＬ０〜Ｌ８を例示する。
ライン数は、ラインに直交する方向の再構成領域の最大画素数の１／６４〜１／２とする。例えば、再構成領域Ｐの画素数が５１２×５１２であるとき、ライン数は９本とする。
また、−４５゜≦view＜４５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および１３５゜≦view＜２２５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）では、ｘ方向をライン方向とする。また、４５゜≦view＜１３５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および２２５゜≦view＜３１５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）では、ｙ方向をライン方向とする。
また、回転中心ＩＣを通り、ラインＬ０〜Ｌ８に平行な投影面ｐｐを想定する。
【００４９】
図１７は、ラインＬ０〜Ｌ８をＸ線透過方向に検出器面ｄｐに投影したラインＴ０〜Ｔ８を示している。
Ｘ線透過方向は、Ｘ線管２１とマルチ検出器２４とラインＬ０〜Ｌ８の幾何学的位置（マルチ検出器２４のｚ軸方向の中心を通るｘｙ面から画像再構成領域Ｐまでのｚ軸方向の距離や、画像再構成面Ｐ上の画素点の集合であるラインＬ０〜Ｌ８の位置を含む）によって決まるが、投影データＤ０(z,view,j,i)の直線移動方向の位置ｚが判っているため、加・減速中の投影データＤ０(z,view,j,i)でもＸ線透過方向を正確に求めることが出来る。
検出器面ｄｐに投影したラインＴ０〜Ｔ８に対応する検出器列ｊおよびチャネルｉの投影データを抽出すれば、それらがラインＬ０〜Ｌ８に対応する投影データＤｒである。
【００５０】
そして、図１８に示すように、ラインＴ０〜Ｔ８をＸ線透過方向に投影面ｐｐ上に投影したラインＬ０’〜Ｌ８’を想定し、それらラインＬ０’〜Ｌ８’に投影データＤｒを展開しておく。
【００５１】
図１５に戻り、ステップＲ３では、各ラインＬ０’〜Ｌ８’の投影データＤｒにコーンビーム再構成荷重を乗算し、図１９に示す如き投影ライン・データＤｐを作成する。
ここで、コーンビーム再構成荷重は、Ｘ線管２１の焦点から投影データＤｒに対応するマルチ検出器２４の検出器列ｊ，チャネルｉまでの距離をｒ０とし、Ｘ線管２１の焦点から投影データＤｒに対応する再構成領域上の点までの距離をｒ１とするとき、（ｒ１／ｒ０）2である。
【００５２】
ステップＲ５では、投影ライン・データＤｐに対してライン方向に補間処理を施し、図２０に示す如き高密度画像各位置ライン・データＤｈを作成する。
高密度画像各位置ライン・データＤｈのデータ密度は、ライン方向の再構成領域の最大画素数の８倍〜３２倍とする。例えば、１６倍として再構成領域Ｐの画素数が５１２×５１２であるとき、データ密度は８１９２点／ラインとする。
【００５３】
ステップＲ６では、高密度画像各位置ライン・データＤｈをサンプリングし且つ必要に応じて補間・補外処理して、図２１に示すように、ラインＬ０〜Ｌ８上の画素の逆投影データＤ２を得る。
【００５４】
ステップＲ７では、高密度画像各位置ライン・データＤｈをサンプリングし且つ補間・補外処理して、図２２に示すように、ラインＬ０〜Ｌ８間の画素の逆投影データＤ２を得る。
【００５５】
なお、図１８〜図２２は、−４５゜≦view＜４５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および１３５゜≦view＜２２５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）を想定しているが、４５゜≦view＜１３５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および２２５゜≦view＜３１５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）では、図２３〜図２７に示すようになる。
【００５６】
図１５に戻り、ステップＲ８では、図２８に示すように、図２２または図２７に示す逆投影データＤ２を画素対応に加算する。
ステップＲ９では、ＣＴ画像の再構成に必要な全ビュー（すなわち、３６０゜分のビュー又は「１８０゜分＋ファン角度分」のビュー）について、ステップＲ１〜Ｒ８を繰り返し、逆投影データＤ３(x,y)を得る。
【００５７】
実施例１のＸ線ＣＴ装置１００によれば、直線移動速度が一定に維持されている間だけでなく、直線移動の加・減速中にも投影データを収集し、収集した投影データを画像再構成に利用する。このため、全体の直線移動距離のうちの加・減速のための直線移動距離分も画像再構成に利用できるようになる。
【００５８】
なお、画像再構成法は、従来公知の２次元画像再構成法でもよいし、従来公知のフェルドカンプ法による３次元画像再構成法でもよい。さらに、特願２００２−０６６４２０号、特願２００２−１４７０６１号、特願２００２−１４７２３１号、特願２００２−２３５５６１号、特願２００２−２３５６６２号、特願２００２−２６７８３３号、特願２００２−３２２７５６号および特願２００２−３３８９４７号で提案されている３次元画像再構成法を用いてもよい。
【実施例２】
【００５９】
実施例１では、図４のステップＳ１で画像再構成に必要なビュー範囲のデータを全て取得した後、ステップＳ４で３次元逆投影処理を行っていた。この場合、データ収集と３次元逆投影処理とが完全に直列になるので、画像が得られるまでのタイムラグ（時間遅れ）が大きい。
そこで、実施例２では、データを収集するのと並行して３次元逆投影処理の一部の処理を進める。これにより、画像が得られるまでのタイムラグを短縮できる。
すなわち、実施例２に係るＸ線ＣＴ装置では、図２９のデータ収集処理と図３０のパラメータ推定処理と図３１の３次元逆投影処理とを並行して実行する。
【００６０】
図２９は、実施例２に係るデータ収集処理を示すフロー図である。
ステップＡ５’，Ａ７’，Ａ１０’以外は、図６のステップと同じなので、ステップＡ５’，Ａ７’，Ａ１０’についてのみ説明する。
ステップＡ５’では、加速中の投影データＤ０(z,view,j,i)を収集し、並行して実行中の３次元逆投影処理に渡す。
ステップＡ７’では、定速中の投影データＤ０(z,view,j,i)を収集し、並行して実行中の３次元逆投影処理に渡す。
ステップＡ５’では、減速中の投影データＤ０(z,view,j,i)を収集し、並行して実行中の３次元逆投影処理に渡す。
【００６１】
図３０は、実施例２に係るパラメータ推定処理を示すフロー図である。
ステップＢ１では、投影データＤ０を未取得の一つのビューに着目する。
ステップＢ２では、テーブル１２の直線移動速度を決める所定関数に基づいて、着目ビューの投影データＤ０が得られる時のテーブル１２のｚ座標を推定する。
ステップＢ３では、推定したテーブル１２のｚ座標に基づいて、着目ビューの投影データＤ０が得られる時のＸ線管２１，マルチ検出器２４および再構成領域Ｐの相対位置を推定する。
ステップＢ４では、推定したＸ線管２１，マルチ検出器２４および再構成領域Ｐの相対位置に基づいて、再構成領域Ｐ上の複数画素間隔あけた複数の平行なラインＬ０〜Ｌ８をＸ線透過方向に検出器面ｄｐに投影したラインＴ０〜Ｔ８を推定する。
ステップＢ５では、推定したラインＴ０〜Ｔ８をＸ線透過方向に投影面ｐｐ上に投影したラインＬ０’〜Ｌ８’の投影データＤｒに乗算すべきコーンビーム再構成荷重を算出する。
ステップＢ６では、画像再構成に必要な全ビューのコーンビーム再構成荷重を算出し終わったなら処理を終了し、未算出のビューが残っているならステップＢ１に戻る。
【００６２】
図３１は、実施例２に係る３次元逆投影処理を示すフロー図である。
ステップＣ１では、並行して実行中のデータ収集処理から、ＣＴ画像の再構成に必要な全ビュー（すなわち、３６０゜分のビュー又は「１８０゜分＋ファン角度分」のビュー）中の一つのビューの投影データＤ０(z,view,j,i)が渡されるのを待ち（図２９のステップＡ４’，Ａ７’，Ａ１０’）、渡されたらステップＣ２へ進む。
【００６３】
ステップＣ２では、データ収集処理から渡された投影データＤ０(z,view,j,i)に対して、前処理（オフセット補正，対数補正，Ｘ線線量補正，感度補正）を行う。
ステップＣ３では、前処理した投影データＤ０(z,view,j,i)に対して、フィルタ処理を行う。すなわち、フーリエ変換し、フィルタ（再構成関数）を掛け、逆フーリエ変換する。
【００６４】
ステップＣ４では、データ収集処理から渡された投影データＤ０(z,view,j,i)の中から、再構成領域Ｐ上の複数画素間隔あけた複数の平行なラインＬ１〜Ｌ８に対応する検出器面ｄｐのラインＴ０〜Ｔ８の投影データＤ０を抽出し、ラインＴ０〜Ｔ８をＸ線透過方向に投影面ｐｐ上に投影したラインＬ０’〜Ｌ８’に投影データＤ０を展開し、図１８に示す如き投影データＤｒを得る。
このとき、並行して実行しているパラメータ推定処理でラインＴ０〜Ｔ８を予め推定してあれば（図３０のステップＢ４）、迅速に投影データＤｒを得ることが出来る。
【００６５】
ステップＣ５では、各ラインＬ０’〜Ｌ８’の投影データＤｒにコーンビーム再構成荷重を乗算し、図１９に示す如き投影ライン・データＤｐを作成する。
このとき、並行して実行しているパラメータ推定処理でコーンビーム再構成荷重を予め推定してあれば（図３０のステップＢ５）、迅速に投影ライン・データＤｐを得ることが出来る。
【００６６】
ステップＣ７では、投影ライン・データＤｐに対してライン方向に補間処理を施し、図２０に示す如き高密度画像各位置ライン・データＤｈを作成する。
【００６７】
ステップＣ８では、高密度画像各位置ライン・データＤｈをサンプリングし且つ必要に応じて補間・補外処理して、図２１に示すように、ラインＬ０〜Ｌ８上の画素の逆投影データＤ２を得る。
【００６８】
ステップＣ９では、高密度画像各位置ライン・データＤｈをサンプリングし且つ補間・補外処理して、図２２に示すように、ラインＬ０〜Ｌ８間の画素の逆投影データＤ２を得る。
【００６９】
なお、図１８〜図２２は、−４５゜≦view＜４５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および１３５゜≦view＜２２５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）を想定しているが、４５゜≦view＜１３５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）および２２５゜≦view＜３１５゜（もしくはそれを主体とし周辺をも含むビュー角度範囲）では、図２３〜図２７に示すようになる。
【００７０】
図３１に戻り、ステップＣ１０では、図２８に示すように、図２２または図２７に示す逆投影データＤ２を画素対応に加算する。
【００７１】
ステップＣ１１では、ＣＴ画像の再構成に必要な全ビュー（すなわち、３６０゜分のビュー又は「１８０゜分＋ファン角度分」のビュー）について、ステップＣ１〜Ｃ１０を繰り返し、逆投影データＤ３(x,y)を得、ステップＣ１２へ進む。
【００７２】
ステップＣ１２では、逆投影データＤ３(x,y)に対して後処理を行い、ＣＴ画像を得る。
【００７３】
実施例２のＸ線ＣＴ装置によれば、直線移動速度が一定に維持されている間だけでなく、直線移動の加・減速中にも投影データを収集し、収集した投影データを画像再構成に利用する。このため、全体の直線移動距離のうちの加・減速のための直線移動距離分も画像再構成に利用できるようになる。
さらに、次の効果が得られる。
（１）パラメータ推定処理により、あるビューの投影データＤ０を取得する前にコーンビーム再構成のためのパラメータを算出してしまうので、投影データＤ０が得られたら迅速に処理できる。
（２）データ収集処理と３次元逆投影処理を並行して行うので、画像が得られるまでのタイムラグを減らすことが出来る。
【００７４】
なお、画像再構成法は、一般的な３次元画像再構成（コーンビーム再構成）であるフェルドカンプ（Feldkump）アルゴリズムや、他の３次元画像再構成アルゴリズムを用いてもよい。その場合でも、同じ効果が得られる。
【産業上の利用可能性】
【００７５】
本発明のＸ線ＣＴ撮像方法およびＸ線ＣＴ装置は、Ｘ線ＣＴ画像を撮影するのに利用できる。
【図面の簡単な説明】
【００７６】
【図１】実施例１にかかるＸ線ＣＴ装置を示すブロック図である。
【図２】Ｘ線管およびマルチ検出器の回転を示す説明図である。
【図３】コーンビームを示す説明図である。
【図４】実施例１に係るＸ線ＣＴ装置の概略動作を示すフロー図である。
【図５】データ収集処理の詳細を示すフロー図である。
【図６】データ収集処理の詳細を示すフロー図である。
【図７】リニアに加・減速する場合の直線移動速度の変化を示すグラフである。
【図８】リニアに加・減速する場合の管電流の変化を示すグラフである。
【図９】ノンリニアに加・減速する場合の直線移動速度の変化を示すグラフである。
【図１０】ノンリニアに加・減速する場合の管電流の変化を示すグラフである。
【図１１】定速区間なしでリニアに加・減速する場合の直線移動速度の変化を示すグラフである。
【図１２】定速区間なしでリニアに加・減速する場合の管電流の変化を示すグラフである。
【図１３】定速区間なしでノンリニアに加・減速する場合の直線移動速度の変化を示すグラフである。
【図１４】定速区間なしでノンリニアに加・減速する場合の管電流の変化を示すグラフである。
【図１５】３次元画像再構成処理の詳細を示すフロー図である。
【図１６】再構成領域上のラインをＸ線透過方向へ投影する状態を示す概念図である。
【図１７】検出器器面に投影したラインを示す概念図である。
【図１８】view=０゜における各ラインの投影データＤｒを投影面に投影した状態を示す概念図である。
【図１９】view=０゜における各ラインの投影ライン・データＤｐを投影面に投影した状態を示す概念図である。
【図２０】view=０゜における各ラインの高密度画像各位置ライン・データＤｈを投影面に投影した状態を示す概念図である。
【図２１】view=０゜における再構成領域上の各ラインの逆投影画素データＤ２を示す概念図である。
【図２２】view=０゜における再構成領域上の各画素の逆投影画素データＤ２を示す概念図である。
【図２３】view=９０゜における各ラインの投影データＤｒを投影面に投影した状態を示す概念図である。
【図２４】view=９０゜における各ラインの投影ライン・データＤｐを投影面に投影した状態を示す概念図である。
【図２５】view=９０゜における各ラインの高密度画像各位置ライン・データＤｈを投影面に投影した状態を示す概念図である。
【図２６】view=９０゜における再構成領域上の各ラインの逆投影画素データＤ２を示す概念図である。
【図２７】view=９０゜における再構成領域上の各画素の逆投影画素データＤ２を示す概念図である。
【図２８】逆投影画素データＤ２を画素対応に全ビュー加算して逆投影データＤ３を得る状態を示す説明図である。
【図２９】実施例２に係るデータ収集処理を示すフロー図である。
【図３０】実施例２に係るパラメータ推定処理を示すフロー図である。
【図３１】実施例２に係る３次元逆投影処理を示すフロー図である。
【符号の説明】
【００７７】
１                  操作コンソール
３                  中央処理装置
１０                寝台装置
１２                テーブル
２０                走査ガントリ
２１                Ｘ線管
２４                マルチ検出器
ｄｐ                検出器面
Ｐ                  再構成領域
ｐｐ                投影面
Ｖａ                加速時の直線移動速度
Ｖｂ                減速時の直線移動速度
Ｖｃ                定速時の直線移動速度
【図１】
2005040582.tif 000003
【図２】
2005040582.tif 000004
【図３】
2005040582.tif 000005
【図４】
2005040582.tif 000006
【図５】
2005040582.tif 000007
【図６】
2005040582.tif 000008
【図７】
2005040582.tif 000009
【図８】
2005040582.tif 000010
【図９】
2005040582.tif 000011
【図１０】
2005040582.tif 000012
【図１１】
2005040582.tif 000013
【図１２】
2005040582.tif 000014
【図１３】
2005040582.tif 000015
【図１４】
2005040582.tif 000016
【図１５】
2005040582.tif 000017
【図１６】
2005040582.tif 000018
【図１７】
2005040582.tif 000019
【図１８】
2005040582.tif 000020
【図１９】
2005040582.tif 000021
【図２０】
2005040582.tif 000022
【図２１】
2005040582.tif 000023
【図２２】
2005040582.tif 000024
【図２３】
2005040582.tif 000025
【図２４】
2005040582.tif 000026
【図２５】
2005040582.tif 000027
【図２６】
2005040582.tif 000028
【図２７】
2005040582.tif 000029
【図２８】
2005040582.tif 000030
【図２９】
2005040582.tif 000031
【図３０】
2005040582.tif 000032
【図３１】
2005040582.tif 000033