(57)【要約】【課題】ヘリカルスキャンに際しての寝台移動速度の設定の自由度を維持しながら、コーンビーム・アーチファクトを低減することを可能とする。【解決手段】主制御部２７は、寝台移動速度をユーザ指示に応じて定まる設定速度としつつ第１の投影データを収集するよう各部を制御する。主制御部２７は、第１の投影データから有効素子列のうちの一部の素子列の出力に基づく投影データを抽出することで特定条件下にて収集される投影データに近似した第２の投影データを取得する。主制御部２７は、第１の投影データに基づいての第１の画像の再構成および第２の投影データに基づいての第２の画像の再構成をそれぞれ画像再構成部２３に行わせる。主制御部２７は、第１の画像と第２の画像との差分画像を生成し、この差分画像に基づいて第１の画像をそのコーンビーム・アーチファクトを低減するように補正する。【選択図】      図１
2011036473.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線を錐状に放射するＸ線源と、
前記Ｘ線源に対向して配置されるとともに、それぞれが前記Ｘ線を検出する多数の検出素子を含む素子列を第１の方向に沿って複数列並べた状態で備えるＸ線検出器と、
被検体を載置するための天板を有し、当該天板を前記Ｘ線源と前記Ｘ線検出器との間隙を通して前記第１の方向に沿って移動可能な寝台と、
前記Ｘ線検出器に備えられた複数の素子列のうちの全てまたは一部として定められ有効素子列の出力に基づいて前記被検体に関する投影データを収集する収集手段と、
前記Ｘ線源を前記第１の方向に沿った回転軸周りで回転させる回転機構と、
ヘリカルピッチおよび前記有効素子列の数に応じて定まる寝台移動速度をユーザ指示に応じて定まる設定速度としつつ、前記収集手段による前記投影データの収集が行われるように前記Ｘ線源、前記寝台、前記収集手段および前記回転機構を制御する収集制御手段と、
第１のコーンビームヘリカル再構成法を使用して前記被検体に関する第１の画像を再構成する第１の再構成手段と、
前記寝台移動速度が特定速度である特定条件下にて収集され投影データに基づくことで前記第１のコーンビームヘリカル再構成法よりもアーチファクトを低減可能な第２のコーンビームヘリカル再構成法を使用して前記被検体に関する第２の画像を再構成する第２の再構成手段と、
前記収集制御手段による制御の下に前記収集手段が収集した投影データから前記有効素子列のうちの一部の素子列の出力に基づく投影データを抽出することで前記特定条件下にて前記収集手段が収集する投影データに近似した投影データを取得する取得手段と、
前記収集制御手段による制御の下に前記収集手段が収集した投影データに基づく再構成および前記取得手段により取得された投影データに基づく再構成をそれぞれ行うように前記第１の再構成手段および第２の再構成手段を制御する再構成制御手段と、
前記第１の画像と前記第２の画像との差分画像を生成する生成手段と、
前記第１の画像を前記差分画像に基づいて、前記第１の画像に生じているコーンビーム・アーチファクトを低減するように補正する補正手段とを具備したことを特徴とするＸ線コンピュータ断層撮影装置。
【請求項２】
前記取得手段は、前記有効素子列のうちの中央に位置する一部の素子列の出力に基づく投影データを抽出することを特徴とする請求項１に記載のＸ線コンピュータ断層撮影装置。
【請求項３】
前記第２の再構成手段は、前記第１の再構成手段よりも解像度の低い画像を再構成し、
前記生成手段は、前記第２の画像に含まれる複数の画素のそれぞれを前記第１の画素に含まれる複数の画素に対応付けて、この対応付けた画素どうしの画素値の差をそれぞれ求めることで前記差分画像を生成することを特徴とする請求項１に記載のＸ線コンピュータ断層撮影装置。
【請求項４】
前記差分画像における画素値のうちの予め定められた閾値を上回るものを予め定められた上限値に置換する手段をさらに備えることを特徴とする請求項１に記載のＸ線コンピュータ断層撮影装置。
【請求項５】
前記差分画像に含まれた画像ノイズを低減する手段をさらに備えることを特徴とする請求項１に記載のＸ線コンピュータ断層撮影装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コーンビームを用いたヘリカルスキャンによって収集した投影データに基づいて被検体に関する画像を得るＸ線コンピュータ断層撮影装置（Ｘ線ＣＴ装置）に関する。
【背景技術】
【０００２】
従来よりコーンビームを用いたヘリカルスキャンによって収集した投影データに基づく画像再構成には、F.D.K再構成法を拡張したＴＣＯＴ（true cone-beam tomography）などの再構成法が用いられている。この従来の再構成法は、近似解再構成法であり、コーン角が大きくなるほどコーンビーム・アーチファクトが増大し、画質が大きく劣化する。
【０００３】
一方、University of Central FloridaのKatsevich氏により提案された再構成法も知られている。この再構成法は完全解再構成法であり、コーンビーム・アーチファクトの少ない、高画質な画像を得ることができる。なお、この再構成法に関連する技術については、例えば特許文献１に開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−０９５２９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
前述したKatsevich氏提案の再構成法（以下、Katsevich法と称する）でコーンビーム・アーチファクトの少ない画像を得るためには、投影データの収集時における寝台移動速度が制限されてしまう。具体的には、寝台移動速度が１ＰＩおよび３ＰＩと称される条件のいずれかを満たす必要がある。１ＰＩおよび３ＰＩとは、Ｘ線検出器の幅に対する寝台移動速度を表しているもので、１ＰＩとは約180度分の投影データが使用可能となる寝台移動速度であり、３ＰＩとは約540度分の投影データが使用可能となる寝台移動速度である。かくしてKatsevich法における寝台移動速度は、１ＰＩまたは３ＰＩが成立する２つの速度からしか選択できない。ちなみに寝台速度は、３ＰＩにおいて１ＰＩよりも遅い。なお、寝台移動速度とは、寝台の天板のヘリカルスキャンに際しての移動速度であり、ヘリカルピッチ、収集スライス厚およびＸ線検出器の列数から定まる。
【０００６】
一方、ＴＣＯＴのような近似解再構成法では、寝台移動速度が任意に設定可能であるが、列数が多くなる（コーン角が大きくなる）ほどにコーンビーム・アーチファクトが増大してしまう。
【０００７】
本発明はこのような事情を考慮してなされたものであり、その目的とするところは、ヘリカルスキャンに際しての寝台移動速度の設定の自由度を維持しながら、コーンビーム・アーチファクトを低減することが可能なＸ線コンピュータ断層撮影装置を提供することにある。
【課題を解決するための手段】
【０００８】
本発明の一態様によるＸ線コンピュータ断層撮影装置は、Ｘ線を錐状に放射するＸ線源と、前記Ｘ線源に対向して配置されるとともに、それぞれが前記Ｘ線を検出する多数の検出素子を含む素子列を第１の方向に沿って複数列並べた状態で備えるＸ線検出器と、被検体を載置するための天板を有し、当該天板を前記Ｘ線源と前記Ｘ線検出器との間隙を通して前記第１の方向に沿って移動可能な寝台と、前記Ｘ線検出器に備えられた複数の素子列のうちの全てまたは一部として定められ有効素子列の出力に基づいて前記被検体に関する投影データを収集する収集手段と、前記Ｘ線源を前記第１の方向に沿った回転軸周りで回転させる回転機構と、ヘリカルピッチおよび前記有効素子列の数に応じて定まる寝台移動速度をユーザ指示に応じて定まる設定速度としつつ、前記収集手段による前記投影データの収集が行われるように前記Ｘ線源、前記寝台、前記収集手段および前記回転機構を制御する収集制御手段と、第１のコーンビームヘリカル再構成法を使用して前記被検体に関する第１の画像を再構成する第１の再構成手段と、前記寝台移動速度が特定速度である特定条件下にて収集され投影データに基づくことで前記第１のコーンビームヘリカル再構成法よりもアーチファクトを低減可能な第２のコーンビームヘリカル再構成法を使用して前記被検体に関する第２の画像を再構成する第２の再構成手段と、前記収集制御手段による制御の下に前記収集手段が収集した投影データから前記有効素子列のうちの一部の素子列の出力に基づく投影データを抽出することで前記特定条件下にて前記収集手段が収集する投影データに近似した投影データを取得する取得手段と、前記収集制御手段による制御の下に前記収集手段が収集した投影データに基づく再構成および前記取得手段により取得された投影データに基づく再構成をそれぞれ行うように前記第１の再構成手段および第２の再構成手段を制御する再構成制御手段と、前記第１の画像と前記第２の画像との差分画像を生成する生成手段と、前記第１の画像を前記差分画像に基づいて補正する補正手段とを備えた。
【発明の効果】
【０００９】
本発明によれば、ヘリカルスキャンに際しての寝台移動速度の設定の自由度を維持しながら、コーンビーム・アーチファクトを低減することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】本発明の一実施形態に係るＸ線コンピュータ断層撮影装置の主要部の構成を示す図。
【図２】図１中の主制御部の処理手順を示したフローチャート。
【発明を実施するための形態】
【００１１】
以下、図面を参照して本発明の実施の形態について説明する。
【００１２】
図１は本実施形態に係るＸ線コンピュータ断層撮影装置（以下、ＣＴ装置と称する）１００の主要部の構成を示す図である。このＣＴ装置１００は、スキャンガントリ１とコンピュータ装置２とから構成される。スキャンガントリ１は、被検体２００に関する投影データを収集するための構成要素である。スキャンガントリ１で収集された投影データは、コンピュータ装置２での画像再構成等の処理に供される。
【００１３】
スキャンガントリ１は、寝台１０、Ｘ線管装置１１、Ｘ線検出器１２、架台回転駆動部１３、高電圧発生部１４、寝台駆動部１５およびデータ収集部１６を含む。
【００１４】
Ｘ線管装置１１およびＸ線検出器１２は、図示しない円環状の回転架台に対向関係で搭載されている。回転架台は、架台回転駆動部１３により駆動されて回転する。このとき、Ｘ線管装置１１とＸ線検出器１２とが、同一の回転軸ＲＡの軸周りを回転する。スキャンガントリ１は、Ｘ線管装置１１およびＸ線検出器１２の回転軌道の内側に、空洞（撮影空間）を形成している。
【００１５】
Ｘ線管装置１１は、Ｘ線管１１ａおよびＸ線フィルタ１１ｂを含む。Ｘ線管１１ａは、高電圧発生部１４から電力供給を受けて、Ｘ線検出器１２に向けてＸ線を放射する。Ｘ線フィルタ１１ｂは、被曝低減のために低エネルギー成分を除去する。高電圧発生部１４は、高電圧変圧器、フィラメント電流発生器および整流器を備える。この他に高電圧発生部１４は、管電圧およびフィラメント電流を任意にまたは段階的に調整するために、管電圧切換器およびフィラメント電流切換器等を備えている。なおＸ線管装置１１は、錐状（円錐状または角錐状）のいわゆるコーンビームとしてＸ線を放射する機能を備える。
【００１６】
Ｘ線検出器１２は、複数のＸ線検出素子をそれぞれ含んだ素子列を備える。これら複数の素子列は、回転軸ＲＡに沿った方向に配列されている。複数のＸ線検出素子はそれぞれ、入射するＸ線の強度に応じた電気信号を出力する。
【００１７】
被検体２００は、寝台１０の天板１０ａに載置される。寝台１０は、寝台駆動部１５により駆動されて、天板１０ａをその長手方向（図１中の左右方向）に移動する。通常、この長手方向が回転軸ＲＡと平行になるように寝台１０が設置される。また通常、被検体２００は、その体軸が回転軸ＲＡに沿うように天板１０ａに載置される。かくして被検体２００は、天板１０ａの移動に伴ってスキャンガントリ１の空洞内に挿入される。
【００１８】
データ収集部１６は、Ｘ線検出器１２の出力をスキャンデータとして収集し、コンピュータ装置２に供給する。なお、Ｘ線検出器１２とデータ収集部１６との間には、スリップリングや光通信などを用いたインタフェースが介挿される。これによりデータ収集部１６は、回転架台を連続回転させながらＸ線検出器１２の出力を収集できる。
【００１９】
コンピュータ装置２は、ガントリ制御部２１、前処理部２２、画像再構成部２３、記憶部２４、表示部２５、操作卓２６および主制御部２７を備える。これらのガントリ制御部２１、前処理部２２、画像再構成部２３、記憶部２４、表示部２５、操作卓２６および主制御部２７は、データ／制御バス２８を介して互いに接続されている。
【００２０】
ガントリ制御部２１は、診断に必要なスキャンデータを得るためのスキャンが行われるようにスキャンガントリ１の動作を制御する。
【００２１】
スキャンガントリ１からコンピュータ装置２に供給されたスキャンデータは、前処理部２２を介して投影データとして画像再構成部２３に供給される。画像再構成部２３は、上記の投影データに基づいて断層画像データを再構成する。なお画像再構成部２３は、近似解ヘリカル再構成法により断層画像データを再構成する機能と、完全解ヘリカル再構成法により断層画像データを再構成する機能とを備える。なお、近似解ヘリカル再構成法としては、コーンビーム・アーチファクトが完全解ヘリカル再構成法よりも大きいものの、投影データを得る際の寝台移動速度の自由度が完全解ヘリカル再構成法よりも高い周知の種々のヘリカル再構成法が適用可能である。また完全解ヘリカル再構成法としては、投影データを得る際の寝台移動速度の自由度が近似解ヘリカル再構成法よりも低いものの、コーンビーム・アーチファクトが近似解ヘリカル再構成法よりも小さい周知の種々のヘリカル再構成法を適用可能である。本実施形態では、近似解ヘリカル再構成法としてF.D.K.ヘリカル再構成法を、かつ完全解ヘリカル再構成法としてKatsevich法をそれぞれ採用することとする。
【００２２】
記憶部２４は、投影データや断層画像データを一時的に記憶する。記憶部２４としては、例えばＨＤＤ（hard disk drive）やＲＡＭ（random access memory）が利用できる。
【００２３】
表示部２５は、記憶部２４に記憶された断層画像データに基づいて断層画像を表示する。
【００２４】
操作卓２６は、操作者が例えば撮影条件などの様々な情報や各種指示を入力するために設けられている。操作卓２６は、操作画面を備える。
【００２５】
主制御部２７は、既存のＣＴ装置で実現されているような各種の動作を実現するようにＣＴ装置１００の各部の動作を総括制御する機能を備える。主制御部２７は他に、次のような各種の機能を備える。主制御部２７は、例えばプロセッサとメモリとを組み合わせて実現される。そして各種の機能は、上記プロセッサがメモリに格納されたプログラムを実行することによって実現される。
【００２６】
上記の機能の１つは、寝台移動速度をユーザ指示に応じて定まる設定速度とした状態で投影データの収集が行われるようにスキャンガントリ１の動作を制御する。上記の機能の１つは、上記の制御の下にデータ収集部１６が収集した投影データからその投影データの収集のために使用された素子列のうちの一部の素子列の出力に基づく投影データを抽出することで、Katsevichヘリカル再構成法に適合する特定条件下にてデータ収集部１６が収集する投影データに近似した投影データを取得する。上記の機能の１つは、上記の制御の下にデータ収集部１６が収集した投影データに基づくF.D.K.ヘリカル再構成法による再構成と、上記のように抽出した投影データに基づくKatsevichヘリカル再構成法による再構成をそれぞれ行うように画像再構成部を制御する。上記の機能の１つは、上記の制御の下にF.D.K.ヘリカル再構成法により再構成された画像（以下、近似解再構成画像と称する）とKatsevichヘリカル再構成法により再構成された画像（以下、完全解再構成画像と称する）とのとの差分画像を生成する。上記機能のもう１つは、近似解再構成画像を差分画像に基づいて補正する。
【００２７】
次に以上のように構成されたＣＴ装置１００の動作について説明する。
【００２８】
このＣＴ装置１００は、既存のＣＴ装置で実現されている各種の撮影を行うことが可能であるが、これについての説明は省略する。そしてここでは、既存のＣＴ装置では実現されていないＣＴ装置１００に特有の撮影動作について説明することとする。
【００２９】
図２は主制御部２７の処理手順を示したフローチャートである。
【００３０】
ステップＳａ１において主制御部２７は、撮影条件をユーザの指示に応じて設定する。ここで設定する撮影条件には、寝台移動速度を含む。なお、寝台移動速度は本実施形態においては、次式のように定義することとする。なお次式において、ＣＳは寝台移動速度［mm/rot］を、ＨＰはヘリカルピッチを、inNsegは投影データの収集のために使用する素子列（以下、有効素子列と称する）の列数を、ＳＷは収集スライス厚［mm］をそれぞれ示す。
【００３１】
ＣＳ＝（ＨＰ／inNseg）・ＳＷ・inNseg
主制御部２７は、ヘリカルピッチＨＰ、有効素子列数inNsegおよび収集スライス厚ＳＷをそれぞれユーザに指定させても良いし、これらの各値の一部のみをユーザに指定させて、他は固定値を適用しても良い。例えば、Ｘ線検出器１２が備える素子列の全てを常に有効素子列とするのであれば、有効素子列数inNsegはＸ線検出器１２が備える素子列の列数を固定的に適用し、ヘリカルピッチＨＰおよび収集スライス厚ＳＷをそれぞれユーザに指定させれば良い。
【００３２】
なお主制御部２７は、１ＰＩ条件および３ＰＩ条件を考慮せずに撮影条件を設定する。そこでユーザは、F.D.K.ヘリカル再構成法の特性を考慮して希望の撮影条件が設定されるように主制御部２７に指示する。ちなみに、F.D.K.ヘリカル再構成法の場合、64列0.5mm収集時に推奨されるヘリカルピッチは一般に５１であり、寝台移動速度は26.5［mm／rot］となる。これに対して64列0.5mm収集時には、１ＰＩ条件および３ＰＩ条件が満たされる寝台移動速度はそれぞれ45.7［mm／rot］および17.8［mm／rot］である。従って、寝台移動速度が上記の推奨値に設定されたならば、そのままでは１ＰＩ条件および３ＰＩ条件のいずれも満たされない。
【００３３】
ステップＳａ２において主制御部２７は、ステップＳａ１で設定した撮影条件でコーンビームヘリカルスキャンを行うようにスキャンガントリ１の動作を制御する。このときにデータ収集部１６が収集した投影データは、前処理部２２で周知の前処理が施された上で、記憶部２４に記憶される。
【００３４】
ステップＳａ３において主制御部２７は、上記のように収集されて記憶部２４に記憶された投影データの全てに基づくF.D.K.ヘリカル再構成法による再構成を行うように画像再構成部２３を制御する。この制御の下に画像再構成部２３がF.D.K.ヘリカル再構成法により再構成した近似解再構成画像は、記憶部２４に記憶される。
【００３５】
一方、ステップＳａ４において主制御部２７は、抽出列数outNsegを次式によって算出する。
【００３６】
outNseg＝｛（π＋2γmax）／（2π・cosγmax）｝・（ＣＳ／ＳＷ）
ただし、γmaxは次式により定義される。次式において、ＦＡはファン角度［deg］である。
【００３７】
γmax＝（ＦＡ／2）／（π・180）
ステップＳａ５において主制御部２７は、ステップＳａ２でのスキャンにおける有効素子列のうちのoutNseg列の素子列（以下、抽出素子列と称する）で取得された投影データをステップＳａ２にて取得されて記憶部２４に記憶されている投影データから抽出する。なお、抽出素子列としては、互いに連続していれば有効素子列のうちのどの素子列を割り当てても良い。ただし、抽出素子列は、有効素子列の中央付近に存在する素子列を割り当てることが望ましい。さらに望ましくは、抽出素子列の中の中央に位置する素子列と有効素子列の中の中央に位置する素子列とを一致させる。
【００３８】
具体的には、64列、収集スライス厚0.5mm、ヘリカルピッチが５３の場合、抽出列数は37列となる。
【００３９】
ステップＳａ６において主制御部２７は、上記の抽出した投影データに基づいてKatsevichヘリカル再構成法による再構成を行うように画像再構成部２３を制御する。この制御の下に画像再構成部２３が、Katsevichヘリカル再構成法により再構成した完全解再構成画像は、記憶部２４に記憶される。
【００４０】
ところで、37列、収集スライス厚0.5mm、寝台移動速度26.5［mm/rot］の条件下では１ＰＩ条件が満たされる。従って、ここで得られる完全解ヘリカル再構画像は、コーンビーム・アーチファクトがほとんど生じていない画像となる。
【００４１】
ステップＳａ７において主制御部２７は、ステップＳａ３で再構成されて記憶部２４に記憶されている近似解再構成画像とステップＳａ６で再構成されて記憶部２４に記憶されている完全解再構成画像との画素毎の差分を求める。さて、近似解再構成画像および完全解再構成画像はいずれも、同一の投影データに基づいて再構成されているため、両者の差分はコーンビーム・アーチファクト成分が支配的となる。そして完全解再構成画像にはコーンビーム・アーチファクトがほとんど生じていないので、近似解再構成画像に生じているコーンビーム・アーチファクト成分が抽出されたこととなる。またこのステップＳａ７において主制御部２７は、上記のように求めた差分値のうちで予め定められた閾値を上回るものを、予め定められた上限値に置き換える閾値処理を行う。このように閾値処理を行ったのちの差分値を各画素の画素値とした画像を、以下においてはアーチファクト成分画像と称する。
【００４２】
ステップＳａ８において主制御部２７は、アーチファクト成分画像に対し、スムージングを施す。このスムージングは、アーチファクト成分画像に含まれる画像ノイズ成分を低減するための処理である。画像ノイズ成分は、近似解再構成画像および完全解再構成画像のそれぞれに生じている画像ノイズの差により生じたものである。スムージングは例えば、ＸＹ面でのスムージングフィルタ（例えば、ガウシアンフィルタ）を適用して、画像ノイズ成分を低減させることによって実現できる。
【００４３】
ステップＳａ９において主制御部２７は、近似解再構成画像をステップＳａ８でスムージングが施された後のアーチファクト成分画像を用いて補正する。すなわち、スムージングが施された後のアーチファクト成分画像は、近似解再構成画像に含まれたコーンビーム・アーチファクト成分とほぼ同一の成分の分布を示すから、このアーチファクト画像を近似解再構成画像から減じることによって、コーンビーム・アーチファクト成分を低減するように近似解再構成画像を補正することができる。なお、ここでアーチファクト画像を近似解再構成画像から減じるに当たって、アーチファクト画像に異常に大きな値や画像ノイズ成分が含まれていると、適正な補正が行えない恐れがある。そこでステップＳａ７における閾値処理やステップＳａ８におけるスムージングを施しておくことによって、補正を適正に行えるようにしている。このため、ステップＳａ７における閾値処理で使用する閾値や、ステップＳａ８におけるスムージングに使用するスムージングフィルタのゲインは、ステップＳａ２におけるスキャンの特性、あるいはステップＳａ３における再構成で生じ得るコーンビーム・アーチファクト成分の大きさなどを考慮して、ステップＳａ９での補正が適正に行われるような値に設定しておくことが望ましい。この値の設定は、例えばシミュレーションや実験などに基づいて行うことが可能であるが、設計者の経験則に基づいて行われても良い。
【００４４】
このようにＣＴ装置１００によれば、画像再構成法は実質的にF.D.Kヘリカル再構成法であるので、寝台移動速度の設定の自由度が高い。これにより、検査目的に応じた柔軟なＣＴ装置１００の運用が可能となる。
【００４５】
また、スムージング処理や閾値処理を施していることにより、F.D.Kヘリカル再構成法に特有の画像ノイズなどは除去されず、これまで広く使用されてきたF.D.Kヘリカル再構成法の画質特性を持った再構成画像を得ることが可能である。
【００４６】
さらにＣＴ装置１００では、１ＰＩ条件が成り立たない寝台移動速度で収集された投影データのうちから１ＰＩ条件が成り立つように投影データを抽出して、この抽出した投影データに基づいてKatsevichヘリカル再構成法によりコーンビーム・アーチファクトの少ない完全解再構成画像を得る。そしてＣＴ装置１００では、F.D.Kヘリカル再構成法により得られた近似解再構成画像に生じているコーンビーム・アーチファクトを上記の完全解再構成画像を用いて補正する。このためＣＴ装置１００によれば、Katsevichヘリカル再構成法と同程度にコーンビーム・アーチファクトの少ない再構成画像を得ることが可能である。
【００４７】
この実施形態は、次のような種々の変形実施が可能である。
【００４８】
近似解再構成画像と完全解再構成画像とでマトリクス数を同じにしておけば、近似解再構成画像の全てを個別の補正値で補正することができるために最適な補正を行うことができる。しかしながら、完全解再構成画像は補正量の算出用として用いられるだけであり、そのマトリクス数は補正後の再構成画像の解像度には寄与しないので、完全解再構成画像のマトリクス数を近似解再構成画像のマトリクス数よりも小さくしても良い。このようにすれば、完全解再構成画像の再構成に要する時間を短縮することが可能となる。
【００４９】
抽出列数は、３ＰＩ条件が成立するように求めても良い。
【００５０】
ステップＳａ７における閾値処理やステップＳａ８におけるスムージング処理は、その一方または両方を省略しても良い。
【００５１】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。
【符号の説明】
【００５２】
１…スキャンガントリ、２…コンピュータ装置、１０…寝台、１０ａ…天板、１１…Ｘ線管装置、１１ｂ…Ｘ線フィルタ、１１ａ…Ｘ線管、１２…Ｘ線検出器、１３…架台回転駆動部、１４…高電圧発生部、１５…寝台駆動部、１６…データ収集部、２１…ガントリ制御部、２２…前処理部、２３…画像再構成部、２４…記憶部、２５…表示部、２６…操作卓、２７…主制御部、２８…制御バス、１００…Ｘ線コンピュータ断層撮影装置（ＣＴ装置）。
【図１】
2011036473.tif 000003
【図２】
2011036473.tif 000004