(57)【要約】【課題】  狭い視野の投影データを画像再構成して広い画面に表示してもノイズの粒状性テキスチャーに違和感を与えないことを課題とする。【解決手段】  X線断層像における関心領域の視野の大きさに対する画像再構成領域（又は画像表示領域）の大きさの拡大率に応じて、この関心領域の投影データを部分拡大画像再構成する際に、再構成領域の拡大率が増すほど投影データと重畳する再構成関数の周波数特性を高い側に移して粒状性を維持する。【選択図】  図２
2006068229.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、再構成領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にすることを特徴とするＸ線ＣＴ断層像の画像再構成方法。
【請求項２】
Ｘ線ＣＴ断層像を再構成する際の画像再構成領域の大きさに応じて、表示領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にすることを特徴とするＸ線ＣＴ断層像の画像再構成方法。
【請求項３】
画像再構成される画像のスライス厚が増すほど再構成関数の周波数特性を変更しないか、低域側に移すことを特徴とする請求項１又は２記載のＸ線ＣＴ断層像の画像再構成方法。
【請求項４】
再構成関数は、周波数領域で規定された基本関数群と、この基本関数の高域部分を抑制する第１の関数群と、この基本関数の高域部分を強調する第２の関数群とを備え、これらの関数の周波数特性を定めるパラメータを調整して合成することにより周波数特性を調整することを特徴とする請求項１又は２記載のＸ線ＣＴ断層像の画像再構成方法。
【請求項５】
あらかじめ異なる周波数特性を有する複数の再構成関数を備え、部分拡大画像再構成の関心領域の大きさで決められる拡大率に応じて、これらの再構成関数の１つを選択することを特徴とする請求項１又は２記載のＸ線ＣＴ断層像の画像再構成方法。
【請求項６】
あらかじめ異なる周波数特性を有する複数の再構成関数を備え、部分拡大画像再構成の関心領域の大きさで決められる拡大率に応じて、これらの再構成関数の２つを選択し、この２つを補間して再構成関数を求めることを特徴とする請求項１又は２記載のＸ線ＣＴ断層像の画像再構成方法。
【請求項７】
被検体をはさんで相対向するＸ線発生手段及びＸ線検出手段を備え、被検体のＸ線投影データに基づきＣＴ断層像を画像再構成するＸ線ＣＴ装置において、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、再構成領域の拡大率が増すほど投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にする再構成関数調整手段を備えることを特徴とするＸ線ＣＴ装置。
【請求項８】
被検体をはさんで相対向するＸ線発生手段及びＸ線検出手段を備え、被検体のＸ線投影データに基づきＣＴ断層像を再構成するＸ線ＣＴ装置において、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、表示領域の拡大率が増すほど投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にする再構成関数調整手段を備えることを特徴とするＸ線ＣＴ装置。
【請求項９】
再構成関数調整手段は、画像再構成される画像のスライス厚が増すほど再構成関数の周波数特性を低域側に移すことを特徴とする請求項７又は８記載のＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（Computed Tomography）断層像の画像再構成方法及びＸ線ＣＴ装置に関し、更に詳しくは、被検体をはさんで相対向するＸ線発生手段及びＸ線検出手段を備え、被検体のＸ線投影データに基づきＸ線ＣＴ断層像を画像再構成するＸ線ＣＴ断層像の画像再構成方法及びＸ線ＣＴ装置に関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置では、被検体を透過したＸ線量子ノイズによって、画像再構成された画素のＣＴ値にもＸ線量子ノイズなどに起因するＸ線ＣＴ断層像の画像ノイズとなって表れる。このＸ線量子ノイズはポアソン分布に従うと共に、平均Ｎ個のＸ線量子を検出するような計測では、その標準偏差が√（Ｎ）になることが知られている。このため、Ｘ線量子ノイズによる画質の粒状性はＸ線断層像では避けられない。またこの粒状性はＸ線断層像とともに認知されて来た。
【０００３】
また、例えば肺野部の診断で片側の肺のみを部分拡大再構成することがある。このように、診断の用途によっては小さな画像再構成領域を部分拡大画像再構成する場合があり、この時、Ｘ線量子ノイズなどに起因する画像の粒状性（粒子）も単純に拡大されてしまうため、同じ再構成関数を用いると、視覚的な違和感を与えていた。
【０００４】
このような状況の下、従来は、撮影範囲の全域を画質を各部位の診断目的に最適な状態でＸ線断層撮影できるように、例えば、頭部Ｍａには通常よりも骨の輪郭が明瞭な画像になるような再構成関数を使用し、また肺野部Ｍｃには通常よりもシャープな画像になるような再構成関数を使用し、さらに、腹部Ｍｄには癌などの病変部が目立つ通常よりもスムーズな画像になるような再構成関数を使用するものが知られている（特許文献１）。しかし、各撮影部位の領域の大きさに最適な再構成関数については、撮影部位の領域の大きさにかかわらず粒状性テキスチャーの視覚的な感覚・印象を維持させることについて具体的な記載はなかった。
【０００５】
また、従来は、広い視野（大きい画像再構成領域）の画像と狭い視野（小さな画像再構成領域）の画像を同時に得ることを課題とし、広い視野（大きい画像再構成領域）の画像については高周波抑制タイプの再構成関数と重畳することによりノイズの少ない滑らかな画像を得、また狭い視野（小さな画像再構成領域）の画像については高周波強調タイプの再構成関数と重畳することにより解像度が高い画像を得ると共に、必要なら、狭い視野（小さな画像再構成領域）の画像については再構成後の画像を拡大して画面に表示するものが知られている（特許文献２）。
【特許文献１】特開２００２−１９１５９２（要約，「００２８」）
【特許文献２】特開２０００−１７５９０３（要約，「００３２」，「００３３」）
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかし、上記小さな画像再構成領域の画像について部分拡大画像再構成を行うと、画質、特に粒状性テキスチャーに対する観察者の視覚的感覚・印象を一定に保つには限界があり、画像の粒状性テキスチャーが大きくぼやけた感じで違和感を感じる不自然な画像となっていた。
【０００７】
本発明は、上記従来技術の問題点に鑑みなされたもので、その目的とする所は、小さな画像再構成領域の投影データを部分拡大画像再構成して表示しても粒状性テキスチャーに違和感を与えないＸ線ＣＴ断層像の画像再構成方法及びＸ線ＣＴ装置を提供することにある。
【課題を解決するための手段】
【０００８】
第１の観点では、本発明は、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、再構成領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にすることを特徴とするＸ線ＣＴ断層像の画像再構成方法を提供する。
【０００９】
上記第１の観点における画像再構成方法では、投影データに対してその拡大率に応じた高い周波数特性（即ち、高域の延長又は強調された）の再構成関数を重畳するため、小さな部分拡大画像再構成領域の投影データを拡大断層像に画像部分拡大画像再構成しても、断層像の粒状性テキスチャーに別段の違和感を与えない。即ち、ノイズ画像を常にノイズと認識できるように保つことが可能となる。
【００１０】
第２の観点では、本発明は、Ｘ線ＣＴ断層像を再構成する際の画像再構成領域の大きさに応じて、表示領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にすることを特徴とするＸ線ＣＴ断層像の画像再構成方法を提供する。
【００１１】
上記第２の観点における画像再構成方法では、投影データに対してその拡大率に応じた高い周波数特性の再構成関数を重畳するため、小さな部分拡大画像再構成領域の投影データを拡大断層像に画像部分拡大画像再構成しても、断層像の粒状性テキスチャーに別段の違和感を与えない。即ち、ノイズ画像を常にノイズと認識できるように保つことが可能となる。
【００１２】
第３の観点では、本発明は、画像再構成される画像のスライス厚が増すほど再構成関数の周波数特性を変更しないか、低域側に移すことを特徴とする第１又は第２の観点に記載のＸ線ＣＴ断層像の画像再構成方法を提供する。
【００１３】
上記第３の観点における画像再構成方法では、再構成画像のスライス厚が増すと、再構成画像のＳ／Ｎも改善される傾向にあるため、ノイズが強調されない方向に再構成関数の周波数特性を変更する。又は変更をしなくても充分であれば変更しない。
【００１４】
第４の観点では、本発明は、再構成関数は、周波数領域で規定された基本関数群と、この基本関数の高域部分を抑制する第１の関数群と、この基本関数の高域部分を強調する第２の関数群とを備え、これらの関数の周波数特性を定めるパラメータを調整して合成することにより周波数特性を調整することを特徴とする第１又は第２の観点に記載のＸ線ＣＴ断層像の画像再構成方法を提供する。
【００１５】
上記第４の観点における画像再構成方法では、例えば図５に示すように、再構成関数Ｋ（ω）は、周波数領域で規定された基本関数Ｋ１（ω）と、この基本関数の高域部分を抑制する第１の関数Ｋ２（ω）と、この基本関数の高域部分を強調する第２の関数Ｋ３（ω）とを備え、これらの内の１又は２以上の関数の周波数特性を規定するパラメータ（ωn，ωd，ωr，ωb，ωa，Ｒ等）を調整して合成することにより周波数特性を調整するものである。このように、再構成関数の周波数特性の任意調整が行え、最適化も可能である。
【００１６】
第５の観点では、本発明は、あらかじめ異なる周波数特性を有する複数の再構成関数を
備え、部分拡大画像再構成の関心領域の大きさで決められる拡大率に応じて、これらの再構成関数の１つを選択することを特徴とする第１又は第２の観点に記載のＸ線ＣＴ断層像の画像再構成方法を提供する。
【００１７】
上記第５の観点における画像再構成方法では、あらかじめ異なる周波数特性を有する複数の再構成関数を備え、前記拡大率に応じてこれらのうちの１つを選択するものである。従って、簡単な処理で周波数特性の適正な選択が行える。
【００１８】
第６の観点では、本発明は、あらかじめ異なる周波数特性を有する複数の再構成関数を備え、部分拡大画像再構成の関心領域の大きさで決められる拡大率に応じて、これらの再構成関数の２つを選択し、この２つを補間して再構成関数を求めることを特徴とする第１又は第２の観点に記載のＸ線ＣＴ断層像の画像再構成方法を提供する。
【００１９】
上記第６の観点における画像再構成方法では、あらかじめ異なる周波数特性を有する複数の再構成関数を備え、前記拡大率に応じてこれらのうちの２つを選択し、２つを補間して１つの再構成関数を作るものである。
【００２０】
第７の観点では、本発明は、被検体をはさんで相対向するＸ線発生手段及びＸ線検出手段を備え、被検体のＸ線投影データに基づきＣＴ断層像を画像再構成するＸ線ＣＴ装置において、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、再構成領域の拡大率が増すほど投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にする再構成関数調整手段を備えることを特徴とするＸ線ＣＴ装置を提供する。
【００２１】
上記第７の観点におけるＸ線ＣＴ装置では、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、再構成領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にするため、小さな部分拡大画像再構成領域の投影データを拡大断層像に画像部分拡大画像再構成しても、断層像の粒状性テキスチャーに別段の違和感を与えない。即ち、ノイズ画像を常にノイズと認識できるように保つことが可能となる。
【００２２】
第８の観点では、本発明は、被検体をはさんで相対向するＸ線発生手段及びＸ線検出手段を備え、被検体のＸ線投影データに基づきＣＴ断層像を再構成するＸ線ＣＴ装置において、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、表示領域の拡大率が増すほど投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にする再構成関数調整手段を備えることを特徴とするＸ線ＣＴ装置を提供する。
【００２３】
上記第８の観点におけるＸ線ＣＴ装置では、Ｘ線ＣＴ断層像を部分拡大画像再構成する際の画像再構成領域上の関心領域の大きさに応じて、表示領域の拡大率が増すほど、投影データと重畳するための再構成関数の周波数特性を高域側に移して最適にするため、小さな部分拡大画像再構成領域の投影データを拡大断層像に画像部分拡大画像再構成しても、断層像の粒状性テキスチャーに別段の違和感を与えない。即ち、ノイズ画像を常にノイズと認識できるように保つことが可能となる。
【００２４】
第９の観点では、本発明は、再構成関数調整手段は、画像再構成される画像のスライス厚が増すほど再構成関数の周波数特性を低域側に移すことを特徴とする第７又は第８の観点に記載のＸ線ＣＴ装置を提供する。
【００２５】
上記第９の観点におけるＸ線ＣＴ装置では、再構成画像のスライス厚が増すと、再構成
画像のＳ／Ｎも改善される傾向にあるため、ノイズが強調されない方向に再構成関数の周波数特性を変更する。又は変更をしなくても充分であれば変更しない。
【発明の効果】
【００２６】
以上のように、本発明によれば、小さな部分拡大画像再構成領域の投影データを画像再構成して拡大断層像に表示しても、断層像の粒状性テキスチャーに違和感を与えないため、部分拡大画像再構成を用いた画像診断においても適正に行える。
【発明を実施するための最良の形態】
【００２７】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
【００２８】
図１は、本発明の一実施形態にかかるＸ線ＣＴ装置の構成ブロック図である。このＸ線ＣＴ装置１００は、被検体を載せて体軸方向（z方向）に移動させる撮影テーブル１０、Ｘ線ファンビームにより被検体のアキシャルスキャン／ヘリカルスキャンデータ収集を行う走査ガントリ２０、撮影テーブル１０および走査ガントリ２０の遠隔制御を行うとともに、操作者が各種の設定操作を行う操作コンソール１とを具備している。
【００２９】
操作コンソール１は、操作者の入力を受け付ける入力装置２と、画像再構成処理などを実行する中央処理装置３と、走査ガントリ２０で取得した投影データを収集するデータ収集バッファ５と、投影データから再構成したＣＴ画像を表示するモニタ６と、プログラムやデータやＸ線ＣＴ画像を記憶する記憶装置７とを具備している。
【００３０】
テーブル装置１０は、被検体を乗せて走査ガントリ２０のボア（空洞部）に入れ出しするクレードル１２を具備している。クレードル１２は、撮影テーブル１０に内臓するモータで昇降およびテーブル直線移動される。
【００３１】
走査ガントリ２０は、Ｘ線管２１と、Ｘ線管の管電圧・管電流を制御するＸ線コントローラ２２と、Ｘ線ファンビームのz方向の厚さ、つまりスライス厚さを制御するコリメータ２３と、多列Ｘ線検出器２４と、データ収集装置ＤＡＳ（Data Acquisition System）２５と、被検体の体軸の回りにＸ線管２１などを回転させる回転部コントローラ２６と、制御信号などを前記操作コンソール１や投影テーブル１０とやり取りする制御コントローラ２９とを具備している。
【００３２】
実施形態におけるＸ線ＣＴシステムの構成は概ね上記の通りである。この構成のＸ線ＣＴ装置において、投影データの収集は例えば次のように行われる。まず、被検体を走査ガントリ２０の回転部１５の空洞部に位置させた状態でz軸方向の位置を固定し、Ｘ線管２１からのＸ線ビームを被検体に照射し（Ｘ線の投影）、その透過Ｘ線を多列Ｘ線検出器２４で検出する。そして、この透過Ｘ線の検出を、Ｘ線管２１と多列Ｘ線検出器２４を被検体の周囲を回転させながら｛すなわち、投影角度（ビュー角度）を変化させながら｝複数Ｎ（例えば、Ｎ＝１,０００）のビュー方向で、３６０度分のデータ収集を行う。
【００３３】
検出された各透過Ｘ線は、ＤＡＳ（データ収集部）２５でディジタル値に変換されて投影データd（ch，view）（ただし、ch：チャネル、view：ビュー）としてデータ収集バッファ５を介して操作コンソール１に転送される。この動作を１スキャンとよぶ。そして、順次z軸方向にスキャン位置を所定量だけ移動して、次のスキャンを行っていく。このようなスキャン方式はコンベンショナルスキャン方式（又はアキシャルスキャン方式）とよばれる。また、このコンベンショナルスキャン（アキシャルスキャン）方式において、連続して複数回転分のスキャンを行うスキャン方式はシネスキャン方式とよばれる。また、投影角度の変化に同期して撮影テーブル１０を所定速度で移動させ、スキャン位置を移動
させながら（Ｘ線管４と多列Ｘ線検出器２４とが被検体の周囲をらせん状に周回することになる）投影データを収集する方式を、いわゆるヘリカルスキャン方式とよぶ。本発明はコンベンショナルスキャン方式、シネスキャン方式、ヘリカルスキャン方式いずれにも適用できる。
【００３４】
操作コンソール１は、走査ガントリ２０から転送されてくる投影データを中央処理装置３の記憶装置７に格納するとともに、例えば、所定の再構成関数と重畳演算を行い、逆投影処理により断層像を再構成する。ここで、操作コンソール１は、スキャン処理中に走査ガントリ２０から順次転送されてくる投影データからリアルタイムに断層像を再構成し、常に最新の断層像をモニタ６に表示させることが可能である。さらに、記憶装置７に格納されている投影データを呼び出して改めて画像再構成を行わせることも可能である。
【００３５】
以下、実施の形態によるＸ線ＣＴ撮影処理の流れを詳細に説明する。図２は実施の形態によるＸ線ＣＴ撮影処理のフローチャートで、例えば頭部中の注目する直径２５cmの領域のＣＴ断層像を５１２×５１２画素（５０cm×５０cmの表示領域に相当）の再構成領域に画像再構成する場合を示している。
【００３６】
図２において、事前に被検体のスカウトスキャンを行った後、この処理に入力する。ステップＳ１１では、続く被検体アキシャル／ヘリカルスキャンのためのスキャン計画画面をモニタ６に表示する。ステップＳ１２では、操作者が、スキャンパラメータを設定する。図３にスキャンパラメータの設定画面を示す。一例のスキャンパラメータは、以下の通りである。
【００３７】
スキャンタイプ[Scan Type]＝アキシャルスキャン
体軸上のスキャン開始位置[Start Loc]＝z１
体軸上のスキャン終了位置[End Loc]＝z２
スキャン回数[Ｎo．of Scan]＝３２
スキャン幅[Thick]＝１mm
スキャン時間[Sec]＝１秒／ガントリ１回転
Ｘ線管の管電圧[kV]＝１２０kV
Ｘ線管の管電流[mA]＝２８０mA
この状態で、操作者が表示エリア１３aの制御アイコン[Show Localizer]をクリックすると、画面上側の表示エリア１３bには被検体のスカウト像｛図は側面（９０度方向スカウト像）を示しているが、正面（０度方向スカウト像）でも良い｝が表示され、その上に各スキャン位置を示す線（スライスライン）が重ねて表示される。図の実太線はスキャンの開始位置及び終了位置をそれぞれ表し、また細線は中間の各スキャン位置を表す。操作者は、スカウト像上のスライスライン等を確認すると共に、必要なら、マウスでドラッグしてスキャン範囲・スライスライン数等を変更可能である。
【００３８】
次に、操作者が表示エリア１３a上のプロスペクティブリコンタグ[P-RecoＮ]をクリックすると、Ｘ線ＣＴ断層像の画像再構成パラメータを入力可能となる。そして、図２のステップＳ１３では、操作者がＣＴ断層像の画像再構成パラメータを設定する。図４に画像再構成パラメータの設定画面を示す。ＣＴ断層像「Ｒ」を得るための一例の画像再構成パラメータは、以下の通りである。
【００３９】
画像再構成するスライス開始位置[Start Loc]＝z３
画像再構成するスライス終了位置[End Loc]＝z４
スライス（画像）枚数[ＮO．of Images]１６枚
被検体のスライス厚[Thick]＝１mm
再構成関数＝A−１（詳細は後述）
画像フィルタ＝タイプ３（詳細は後述）
画像再構成のマトリクスサイズ＝５１２×５１２画素
更に、画面下側には被検体の注目領域（再構成領域）を指定するためのアイコン１３cが表示されており、操作者は例えば最初のスキャンで得られた全体の肺野部断層像の注目領域である片側の肺野の中心座標（ｘ，ｙ）と、この領域の直径（例えば２５cm）を画像再構成パラメータ入力又は関心領域のマウス操作で指定可能である。
【００４０】
なお、被検体の肺野部断層像は、事前に全体断層像をスキャンし画像再構成パラメータを変更して部分拡大画像再構成断層像を表示するのが通常である。こうすれば、疾患部（関心領域）を正確に指定可能である。あるいは、事前のスカウト像から被検体頭部のおおよそのサイズ、形状を推定しても良いし、又は、小児や成人の頭部の指定によって選択された典型的な頭部の断面像（模型像）であっても良い。これらのいずれにしても、疾患部を含むようなおおよその関心領域を指定可能である。
【００４１】
図２に戻り、ステップＳ１４では設定確認ボタン「COＮFIRM」の入力を待ち、このボタンが入力されると、ステップＳ１５では上記設定されたスキャンパラメータに従って被検体のスキャン制御を行う。ステップＳ１６では被検体の投影データを収集・蓄積する。ステップＳ１７では撮影条件の全スキャンを完了したか否かを判別し、完了でない場合はステップＳ１５に戻り、スキャンを繰り返す。
【００４２】
こうして、やがて、ステップＳ１７の判別で全スキャンを完了すると、ステップＳ１８では投影データに対して所定の前処理（対数変換、Ｘ線検出器の感度補正等）を行う。
図２に戻り、ステップＳ１９では、再構成関数が“調整して作成される”、又は“既に用意されているものから選択される”。再構成関数は、まず選択されている部位ごとに再構成関数は大きく変えるのが通常である。またさらに、上記“スキャン領域の１画素の大きさと、表示領域又は再構成領域の１画素の大きさとの比”、その他、スライス厚などに応じて所定の再構成関数（再構成カーネル、フィルタ関数等とも呼ばれる）の周波数特性を微調整する。ここで、所定の再構成関数とは、例えば上記図４の再構成関数で選択された再構成関数A−１である。
【００４３】
図５は再構成関数の調整処理を説明している。ここでは再構成関数の周波数領域における調整処理を示している。一般に、投影データ領域，空間領域における画像再構成は、投影データd（ch，view）と再構成関数h（ch）との重畳積分を全ビューについて行う。しかし、一般的には空間領域の重畳演算が周波数領域では乗算になることを利用して、処理の高速化の観点から投影データをフーリエ変換して周波数領域に変換し、周波数領域に変換された再構成関数を乗算し、逆フーリエ変換で空間領域に投影データを戻すことがよく行われる。
【００４４】
図５では再構成関数の空間周波数軸における調整を示し、説明する。図５（a）は周波数領域における関数Ｋ１（ω）、図５（b），（c）の関数Ｋ２（ω），Ｋ３（ω）、そして、図５（d）はこれらを組合せて構成した再構成関数Ｋ（ω）を示している。各関数Ｋ１（ω）〜Ｋ３（ω）は例えば記憶装置７に記憶されており、パラメータを変更することにより様々な周波数特性が得られる。
【００４５】
図５（a）の関数Ｋ１（ω）は次式、
【００４６】
【数１】
2006068229.tif 000003
【００４７】
で表される。この関数Ｋ１（ω）の周波数特性はゲインが０となる周波数パラメータωn（Ｘ線検出器のチャネルピッチによってきまるナイキスト周波数に相当）を指定することにより変更可能である。なお、通常はデフォルト値としてのＸ線検出器のチャネルピッチによってきまるナイキスト周波数ωnを使用する。
【００４８】
図５（b）の関数Ｋ２（ω）は、高周波成分を抑制するための関数であり、具体的な１例としては次式、
【００４９】
【数２】
2006068229.tif 000004
【００５０】
で表される。この関数Ｋ２（ω）の周波数特性はダンピングが開始される周波数パラメータωd及び又はゲインが１／２に落ちる時点の周波数パラメータωrを指定することにより変更可能である。
【００５１】
また、図５（c）の関数Ｋ３(ω)は、上記とは逆に高周波成分を強調するための修正関数であって、具体的な１例としては次式、
【００５２】
【数３】
2006068229.tif 000005
【００５３】
で表される。この関数Ｋ３(ω)は、(ω)の周波数特性はゲインが上昇開始する周波数パラメータωb及び又はゲインが{(R＋１)／２}に上がる時点の周波数（ωa＋ωb）を特定するための周波数パラメータωaを指定することにより変更可能である。
【００５４】
そして、図５（d）の再構成関数Ｋ(ω)は次式、
【００５５】
【数４】
2006068229.tif 000006
【００５６】
但し、scaleは所定のスケーリング定数
で表される。このように、再構成関数Ｋ(ω)を関数Ｋ１(ω)と、周波数特性を調整するための関数Ｋ２(ω)，Ｋ３(ω)の組合せで作ることにより、実際に使用する再構成関数Ｋ(ω)の高周波成分の抑制又は強調を、修正関数Ｋ２(ω)，Ｋ３(ω)に対する簡単なパラメ
ータ（ωd，ωr，ωb，ωa，Ｒ等）の調整で容易に行える。図５（d）に高域を強調していない再構成関数（１）と、高域を強調した再構成関数（３）の例を示す。
【００５７】
次に、“スキャン領域の１画素の大きさと、表示領域又は再構成領域の１画素の大きさとの比”（拡大率）に応じて再構成関数の周波数特性を微調整（又は選択）する関係を説明する。まず、関心投影領域ＲＯＩに対する再構成領域ＤＦＯＶの拡大率が増すほど、つまり“スキャン領域の１画素の大きさと、表示領域又は再構成領域の１画素の大きさとの比”が大きくなるほど投影データと重畳するための再構成関数の周波数特性を高域に移動させる。これにより、高い周波数特性の再構成関数と重畳することにより量子ノイズのランダム性が強調され、これにより粒状性テキスチャーのザラザラ感が強調され、全体の断層像と部分拡大断層像のテキスチャーの差がより少なくなる。
【００５８】
又は、関心投影領域に対する表示領域の拡大率が増すほど、つまり、“スキャン領域の１画素の大きさと、表示領域又は再構成領域の１画素の大きさとの比”が大きくなるほど投影データと重畳するための再構成関数の周波数特性を高域に移動させる。従って、上記同様の作用・効果が得られる。
【００５９】
また、再構成画像のスライス厚が増すほど、ＳＮが良くなるため、上記の再構成関数の周波数特性を低域に移動させるか、変えなくても充分であれば変えないことで全体の断層像と部分拡大断層像のテキスチャーの差がより少なくなる。
【００６０】
本実施の形態では、上記画像再構成パラメータの設定に従って再構成関数の周波数特性が調整（又は選択）される。従って、“スキャン領域の１画素の大きさと、表示領域又は再構成領域の１画素の大きさとの比”、画像表示領域の大きさ（即ち、部分拡大画像再構成関心領域と画像表示領域との比）やスライス厚が変わっても、ＣＴ断層像の画質、特に粒状性テキスチャーが感覚的にほぼ一定に保たれる。
【００６１】
図２に戻り、ステップＳ２０では、投影データと上記の調整（又は選択）された再構成関数との重畳積分を行う。今、投影データをＤ(ch,view)、再構成関数をＫ(ω)とすると、これらの重畳積分後の投影データＤflt(ch,view)は次式、
【００６２】
【数５】
2006068229.tif 000007
【００６３】
ここで、Ｆはフーリエ変換
Ｆ−１はフーリエ逆変換
により求められる。この時のフーリエ変換は、通常、チャンネル数が約１０００chであることから、例えば約２０００点のフーリエ変換を用いるのも可能であるが、再構成画像の粒状性を細かくしたい場合には、約４０００点のフーリエ変換を用いてもよい。ステップＳ２１では、再構成関数を重畳された投影データを逆投影処理することで断層像が再構成される。ステップＳ２２では、このように画像再構成されたＣＴ断層像をモニタ６に表示する。
【００６４】
なお、上記実施の形態における再構成関数Ｋ(ω)の部分関数Ｋ１(ω)，Ｋ２(ω)，Ｋ３(ω)を上記とは異なる部分関数で構成して、高周波領域低減、あるいは周波数特性の調整を行ってもよい。
【００６５】
また、上記実施の形態では、本発明の処理を操作コンソール１の中央処理装置３により実行したが、この操作コンソール１とは独立してスクリーニング検査等のための画像観察用操作コンソール端末を上記Ｘ線ＣＴ装置に接続し、上述した処理をこの端末に行わせることも可能である。操作コンソール１および上記画像観察用操作コンソール端末の構成自体は汎用の画像処理装置（ワークステーションやパーソナルコンピュータ等）で実現できるので、上記のようなソフトウェアを同装置にインストールし、それでもって実現することも可能である。
【００６６】
従って、本発明の機能処理をコンピュータで実現するための、コンピュータにインストールされるプログラム自体および、そのプログラムを格納した記録媒体そのものも、本発明を実現するものである。つまり、本発明の特許請求の範囲には、本発明の機能処理を実現するためのプログラム自体および、そのプログラムを格納したコンピュータ読み取り可能な記録媒体も含まれる。プログラムを供給するための記録媒体としては、例えば、フレキシブルディスク、光ディスク（CD-ROM、CD-R、CD-RW、DVD等）、光磁気ディスク（MO）、磁気テープ、メモリカード等がある。その他、プログラムの供給方法としては、インターネットを介して本発明のプログラムをファイル転送によって取得する態様も含まれる。
【００６７】
また、画像再構成法は、従来公知の２次元画像再構成法でもよいし、従来公知のフェルドカンプ法による３次元画像再構成法でもよい。さらに、特願２００２−０６６４２０号、特願２００２−１４７０６１号、特願２００２−１４７２３１号、特願２００２−２３５５６１号、特願２００２−２３５６６２号、特願２００２−２６７８３３号、特願２００２−３２２７５６号および特願２００２−３３８９４７号などで提案されている３次元画像再構成法を用いてもよい。
【００６８】
また、上記実施の形態では、スキャン直前に設定したり画像再構成条件に従ってスキャン直後に画像再構成を行う、いわゆるレトロスペクティブリコン（Retrospective Recon）時への本発明の適用例を述べたが、これに限らない。本発明はスキャン後に設定した条件に従って被検体を再びスキャンするところの、いわゆるＣＴ透視（フルオロ）モードなどに適用してもよい。
【００６９】
また、上記本発明の１つの実施の形態を述べたが、本発明思想を逸脱しない範囲内で各部の構成、制御、処理及びこれらの組み合わせの様々な変更が行えることは言うまでもない。
【図面の簡単な説明】
【００７０】
【図１】実施の形態によるＸ線ＣＴ装置の構成図である。
【図２】実施の形態によるＸ線ＣＴ撮影処理のフローチャートである。
【図３】実施の形態によるスキャンパラメータの設定画面を示す図である。
【図４】実施の形態による画像再構成パラメータの設定画面を示す図である。
【図５】実施の形態による再構成関数の調整処理を説明する図である。
【符号の説明】
【００７１】
１  操作コンソール
２  入力装置
３  中央処理装置
５  データ収集バッファ
６  モニタ
７  記憶装置
１０  撮影テーブル
１２  クレードル
１５  回転部
２０  走査ガントリ
２１  Ｘ線管
２２  Ｘ線コントローラ
２３  コリメータ
２４  多列Ｘ線検出器
２５  ＤＡＳ（データ収集装置）
２６  回転部コントローラ
２９  制御コントローラ
３０  スリップリング
【図１】
2006068229.tif 000008
【図２】
2006068229.tif 000009
【図３】
2006068229.tif 000010
【図４】
2006068229.tif 000011
【図５】
2006068229.tif 000012