(57)【特許請求の範囲】
【請求項１】
計算機式断層写真法イメージング・システム（１０）の各部を制御手段が制御して再構成される画像の雑音の低減を促進する方法であって、
前記システムが生成した投影データの雑音分布の特性を決定する工程（６４）と、
前記システムが前記雑音分布の特性決定を用いて前記投影データに適応型雑音低減演算を施す工程（６６）と、
前記システムが画像を再構成する工程（６８）とを備え、
前記投影データの雑音分布の特性を決定する工程（６４）は、
空気較正ベクトルを生成する工程と、
ボウタイ形ベクトルを生成するために前記空気較正ベクトルの多項式フィットを実行する工程と、
前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似する工程とを含んでおり、
前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似する工程は、
σ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）
に従って前記雑音分布の標準偏差を近似する工程を含んでおり、ここで、
γは特定のビュー及び特定の投影角の投影標本についての検出器角度であり、
σ（γ）は特定の検出器行及び投影角についての検出器角度γでの投影データの標準偏差であり、
Ｂ（γ）は前記ボウタイ形ベクトルであり、
ｐ（γ）は特定の検出器行及び投影角についての検出器角度γでの較正後の投影データであり、
ｐ（γ0）は参照チャネル・データである、方法（６０）。
【請求項２】
再構成される画像の雑音の低減を促進する計算機式断層写真法（ＣＴ）イメージング・システム（１０）であって、
検出器アレイ（１８）と、
１以上の放射線源（１４）と、
前記検出器アレイ及び放射線源に結合されているコンピュータ（３６）とを備えており、
該コンピュータ（３６）は、投影データを生成し（６２）、該投影データの雑音分布の特性を決定し（６４）、該雑音分布の特性決定を用いて前記投影データに適応型雑音低減演算を施して（６６）、画像を再構成するように構成されており、
前記投影データの雑音分布の特性を決定するために、前記コンピュータ（３６）は、空気較正ベクトルを生成し、ボウタイ形ベクトルを生成するために前記空気較正ベクトルの多項式フィットを実行して、前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似するようにさらに構成されており、
前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似するために、前記コンピュータ（３６）は、
σ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）に従って前記雑音分布の標準偏差を近似するようにさらに構成されており、ここで、
γは特定のビュー及び特定の投影角の投影標本についての検出器角度であり、
σ（γ）は特定の検出器行及び投影角についての検出器角度γでの投影データの標準偏差であり、
Ｂ（γ）は前記ボウタイ形ベクトルであり、
ｐ（γ）は特定の検出器行及び投影角についての検出器角度γでの較正後の投影データであり、
ｐ（γ0）は参照チャネル・データである、
計算機式断層写真法（ＣＴ）イメージング・システム（１０）。
【請求項３】
計算機式断層写真法イメージング・システム（１０）において再構成される画像の雑音の低減を促進するシステムにより実行可能なプログラムで符号化されているコンピュータ読み取り可能な媒体であって、前記プログラムは、
投影データを生成し（６２）、
該投影データの雑音分布の特性を決定し（６４）、
該雑音分布の特性決定を用いて前記投影データに適応型雑音低減演算を施して（６６）、画像を再構成す（６８）べく前記コンピュータに指示するように構成されており、
前記投影データの雑音分布の特性を決定するために、前記プログラムは、空気較正ベクトルを生成し、ボウタイ形ベクトルを生成するために前記空気較正ベクトルの多項式フィットを実行して、前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似するようにさらに構成されており、
前記ボウタイ形ベクトルを用いて前記雑音分布の標準偏差を近似するために、前記プログラムは、
σ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）に従って前記雑音分布の標準偏差を近似するようにさらに構成されており、ここで、
γは特定のビュー及び特定の投影角の投影標本についての検出器角度であり、
σ（γ）は特定の検出器行及び投影角についての検出器角度γでの投影データの標準偏差であり、
Ｂ（γ）は前記ボウタイ形ベクトルであり、
ｐ（γ）は特定の検出器行及び投影角についての検出器角度γでの較正後の投影データであり、
ｐ（γ0）は参照チャネル・データである、
コンピュータ読み取り可能な媒体。
【発明の詳細な説明】
【０００１】
【発明の背景】
本発明は、断層写真法撮像に関し、さらに具体的には、計算機式断層写真法（ＣＴ）イメージング・システムにおいて再構成される画像の雑音を低減させる方法及び装置に関する。
【０００２】
少なくとも一つの公知のＣＴイメージング・システム構成においては、Ｘ線源がファン（扇形）形状のビームを投射し、このビームは、デカルト座標系のＸＹ平面であって、一般に「イメージング（撮像）平面」と呼ばれる平面内に位置するようにコリメートされる。Ｘ線ビームは患者等の被撮像物体を透過する。ビームは物体によって減弱された後に放射線検出器のアレイに入射する。検出器アレイで受光される減弱したビーム放射線の強度は、物体によるＸ線ビームの減弱量に依存している。アレイ内の各々の検出器素子が、検出器の位置でのビーム減弱の測定値である別個の電気信号を発生する。すべての検出器からの減弱測定値を別個に取得して透過プロファイル（断面）を形成する。
【０００３】
公知の第３世代ＣＴシステムでは、Ｘ線源及び検出器アレイは、Ｘ線ビームが物体と交差する角度が定常的に変化するように撮像平面内で被撮像物体の周りをガントリと共に回転する。一つのガントリ角度での検出器アレイからの一群のＸ線減弱測定値すなわち投影データを「ビュー」と呼ぶ。物体の「走査（スキャン）」は、Ｘ線源及び検出器が一回転する間に様々なガントリ角度すなわちビュー角度において形成される一組のビューで構成される。
【０００４】
アキシャル・スキャン（軸方向走査）では、投影データを処理して、物体を通して得られる２次元スライスに対応する画像を構成する。一組の投影データから画像を再構成する一つの方法は、当業界でフィルタ補正逆投影法と呼ばれている。この手法は、走査からの減弱測定値を「ＣＴ数」又は「ハンスフィールド(Hounsfield)単位」と呼ばれる整数へ変換し、これらの整数を用いて陰極線管表示器上の対応するピクセルの輝度を制御するものである。
【０００５】
【発明が解決しようとする課題】
幾つかの公知のイメージング・システムでは、ミリメートル未満の単位での走査速度を用いると、Ｘ線源のピーク電力及び持続電力に対する要件が大幅に強まって、結果的に検出器が信号損失を生ずる。加えて、信号損失は、個々の検出器セルの間の間隙によって、検出されるＸ線フォトンの数が二分の一未満に減少することにより生じ、画質の低下を助長する場合がある。加えて、Ｘ線量が増大すると人体に悪影響を及ぼす虞があり、従って、患者が受けるＸ線量は最小限に抑えなければならないので、Ｘ線フォトン束を減少させることにより患者が受けるＸ線量の減少を促す。このことは一般的には、Ｘ線管電流を減少させることにより行なわれており、これにより典型的にはＸ線フォトンの数が減少する。
【特許文献１】
米国特許第６０７３０４０号
【０００６】
【課題を解決するための手段】
一つの側面では、計算機式断層写真法イメージング・システムにおいて再構成される画像の雑音の低減を促進する方法を提供する。この方法は、投影データを生成する工程と、投影データの雑音分布の特性を決定する工程と、雑音分布の特性決定を用いて投影データに適応型（アダプティブ）雑音低減演算を施す工程と、画像を再構成する工程とを含んでいる。
【０００７】
もう一つの側面では、計算機式断層写真法イメージング・システムにおいて再構成される画像の雑音の低減を促進する方法を提供する。この方法は、投影データを生成する工程と、オフセット補正、一次速度補正、参照チャネル補正、及び空気較正の１以上を用いて投影データを前処理する工程とを含んでいる。この方法はまた、投影データの雑音分布の特性を決定する工程を含んでおり、雑音分布の特性を決定する工程は、空気較正ベクトルを生成する工程と、空気較正ベクトルの多項式フィットを実行してボウタイ形ベクトルを生成する工程と、ボウタイ形ベクトルを用いてσ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）に従って雑音分布の標準偏差を近似する工程とを含んでいる。この方法はまた、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施す工程と、ｐr（γ）＝ｐ（γ）−α（γ）［ｐ（γ）−ｐf（γ）］に従って投影データをフィルタ処理する工程と、画像を再構成する工程とを含んでいる。
【０００８】
さらなる側面では、再構成される画像の雑音の低減を促進する計算機式断層写真法（ＣＴ）イメージング・システムを提供する。この計算機式断層写真法（ＣＴ）イメージング・システムは、検出器アレイと、１以上の放射線源と、検出器アレイ及び放射線源に結合されており投影データを生成するように構成されているコンピュータとを含んでいる。コンピュータはまた、投影データの雑音分布の特性を決定し、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施して、画像を再構成するように構成されている。
【０００９】
さらにもう一つの側面では、再構成される画像の雑音の低減を促進する計算機式断層写真法（ＣＴ）イメージング・システムを提供する。この計算機式断層写真法（ＣＴ）イメージング・システムは、検出器アレイと、１以上の放射線源と、検出器アレイ及び放射線源に結合されており投影データを生成すると共にオフセット補正、一次速度補正、参照チャネル補正、及び空気較正の１以上を用いて投影データを前処理するように構成されているコンピュータとを含んでいる。コンピュータはまた、投影データの雑音分布の特性を決定するように構成されており、雑音分布の特性を決定するために、空気較正ベクトルを生成し、空気較正ベクトルの多項式フィットを実行してボウタイ形ベクトルを生成し、ボウタイ形ベクトルを用いてσ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）に従って雑音分布の標準偏差を近似することを含んでいる。コンピュータはまた、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施し、ｐr（γ）＝ｐ（γ）−α（γ）［ｐ（γ）−ｐf（γ）］に従って投影データをフィルタ処理して、画像を再構成するように構成されている。
【００１０】
さらにもう一つの側面では、計算機式断層写真法イメージング・システムにおいて再構成される画像の雑音の低減を促進するコンピュータにより実行可能なプログラムで符号化されているコンピュータ読み取り可能な媒体を提供する。このプログラムは、投影データを生成し、投影データの雑音分布の特性を決定し、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施して、画像を再構成すべくコンピュータに指示するように構成されている。
【００１１】
さらにもう一つの側面では、計算機式断層写真法イメージング・システムにおいて再構成される画像の雑音の低減を促進するコンピュータにより実行可能なプログラムで符号化されているコンピュータ読み取り可能な媒体を提供する。このプログラムは、投影データを生成すると共にオフセット補正、一次速度補正、参照チャネル補正、及び空気較正の１以上を用いて投影データを前処理すべくコンピュータに指示するように構成されている。プログラムはまた、投影データの雑音分布の特性を決定するように構成されており、雑音分布の特性を決定するために、空気較正ベクトルを生成し、空気較正ベクトルの多項式フィットを実行してボウタイ形ベクトルを生成し、ボウタイ形ベクトルを用いてσ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）に従って雑音分布の標準偏差を近似することを含んでいる。プログラムはまた、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施し、ｐr（γ）＝ｐ（γ）−α（γ）［ｐ（γ）−ｐf（γ）］に従って投影データをフィルタ処理して、画像を再構成すべくコンピュータに指示するように構成されている。
【００１２】
【発明の実施の形態】
本書で用いられる用語として、単数形で記載されており単数不定冠詞を冠した要素又は工程という用語は、排除すると明記していない限りかかる要素又は工程を複数備えることを排除しないものと理解されたい。さらに、本発明の「一実施形態」に対する参照は、所載の特徴を同様に組み入れている他の実施形態の存在を排除しないものと解釈されたい。
【００１３】
また、本書で用いられる「画像を再構成する」という句は、画像を表わすデータが生成されるが可視画像は形成されないような本発明の実施形態を排除するものではない。但し、多くの実施形態は１以上の可視可能な画像を形成する（か又は形成するように構成されている）。
【００１４】
図１及び図２には、計算機式断層写真（ＣＴ）イメージング・システム１０が、「第３世代」ツイン・ビームＣＴスキャナに典型的なガントリ１２を含むものとして示されている。ガントリ１２はＸ線放射線源１４を有しており、Ｘ線放射線源１４は、Ｘ線放射線ビーム１６をガントリ１２の対向する側に設けられている検出器アレイ１８に向かって投射する。検出器アレイ１８は検出器素子２０によって形成されており、検出器素子２０は一括で、物体２２、例えば患者を透過する投射Ｘ線を感知する。各々の検出器素子２０は、入射Ｘ線ビームの強度を表わし従って患者２２を透過する際のビームの減弱を表わす電気信号を発生する。Ｘ線投影データを取得するための一回の走査の間に、ガントリ１２及びガントリ１２に装着されている構成部品は、回転中心２４の周りを回転する。検出器アレイ１８は、検出器アレイ１８が検出器素子又は検出器セル２０の行を複数で有するマルチスライス構成で作製されており、その一行のみを図２に示す。ツイン・ビーム螺旋走査時には、二行の検出器行からデータが同時に取得される。かかる構成を有する検出器素子２０の１以上の追加行が図示の行に平行に配列されており、各々の行が患者２２の並進方向（すなわちｚ軸又は患者の体軸）を横断している。
【００１５】
ガントリ１２の回転及びＸ線源１４の動作は、ＣＴシステム１０の制御機構２６によって制御されている。制御機構２６はＸ線制御器２８とガントリ・モータ制御器３０とを含んでおり、Ｘ線制御器２８はＸ線源１４に電力信号及びタイミング信号を供給し、ガントリ・モータ制御器３０はガントリ１２の回転速度及び位置を制御する。制御機構２６内に設けられているデータ取得システム（ＤＡＳ）３２が検出器素子又は検出器セル２０からのアナログ・データをサンプリングして、後続の処理のためにこのデータをディジタル信号へ変換する。画像再構成器３４が、サンプリングされてディジタル化されたＸ線データをＤＡＳ３２から受け取って高速画像再構成を実行する。再構成された画像はコンピュータ３６への入力として印加され、コンピュータ３６は記憶装置３８に画像を記憶させる。コンピュータ３６はまた、キーボードを有するコンソール４０を介して操作者から指令及び走査用パラメータを受け取る。陰極線管表示器又は液晶表示器のような付設されている表示器４２によって、操作者は、再構成された画像及びコンピュータ３６からのその他のデータを観測することができる。操作者が供給した指令及びパラメータはコンピュータ３６によって用いられて、ＤＡＳ３２、Ｘ線制御器２８及びガントリ・モータ制御器３０に制御信号及び情報を供給する。加えて、コンピュータ３６は、モータ式テーブル４６を制御するテーブル・モータ制御器４４を動作させて、患者２２をガントリ１２内で配置する。具体的には、テーブル４６は患者２２の各部分をガントリ開口４８を通して移動させる。
【００１６】
一実施形態では、コンピュータ３６は着脱式媒体５２に読み書きを行なう装置５０を含んでいる。例えば、装置５０はフロッピー（商標）・ディスク・ドライブ、ＣＤ−Ｒ／Ｗドライブ又はＤＶＤドライブである。これに対応して、媒体５２は、フロッピー・ディスク、コンパクト・ディスク又はＤＶＤのいずれかとなる。一実施形態では、装置５０及び媒体５２を用いて、イメージング・システム１０から他のコンピュータへ取得された投影データを転送してさらなる処理を行ない、或いは他の実施形態では、コンピュータ３６によって処理される機械読み取り可能な命令を入力する。
【００１７】
イメージング・システム１０のコンピュータ３６及び／又は画像再構成器３４は、単独で又は組み合わせて、本発明の１以上の実施形態において本書に記載する計算工程を実行するのに必要な処理能力を提供する。計算工程を実行するための命令は、記憶装置３８、読み取り専用メモリ若しくは読み書きメモリ（図１では別個には示されていない）のような付設メモリ、又は媒体５２に記憶されている。
【００１８】
投影データに平滑化演算を施す際の一つの問題点は、統計的ゆらぎから実際の信号変動を分離するのが困難であることである。しかしながら、統計的雑音のみによって生ずる信号の特定の変動を予め決定できる場合には、空間分解能に影響を与えずに投影データに低域通過フィルタ演算を施すことができる。例えば、測定された信号すなわち投影データの雑音特性によって、統計的ゆらぎに起因する変動と物体２２の実際の構造に起因する変動との間の弁別を容易にする。
【００１９】
図３は、計算機式断層写真法イメージング・システムの統計的雑音の低減を促進する方法６０の流れ図である。一実施形態では、雑音はＸ線フォトン雑音である。もう一つの実施形態では、雑音は電子雑音とフォトン雑音との組み合わせである。方法６０は、投影データを生成する工程６２、投影データの雑音分布の特性を決定する工程６４、雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施す工程６６、及び画像を再構成する工程６８を含んでいる。
【００２０】
一実施形態では、投影データを生成する工程６２は、計算機式断層写真法（ＣＴ）イメージング・システム１０（図１に示す）を用いて物体２２（図１に示す）を走査する工程を含んでいる。投影データは前処理される。一実施形態では、前処理は、オフセット補正、一次速度補正、参照チャネル補正、及び空気較正を含むがこれらに限定されない。一実施形態では、投影データはマイナス対数演算の工程まで前処理される。投影データを前処理してイメージング・システム１０を較正することにより、画像のアーティファクトの減少が促進される。
【００２１】
投影データの雑音分布の特性を決定する工程６４は、空気較正ベクトルを生成する工程、空気較正ベクトルの多項式フィットを実行してボウタイ形ベクトルを生成する工程、及びボウタイ形ベクトルを用いて雑音分布の標準偏差を近似する工程を含んでいる。
【００２２】
一実施形態では、空気較正ベクトルはボウタイ・フィルタを用いて生成される。従って、空気較正後の信号は、実際の測定されたＸ線フォトン束を表わすものでは最早なくなる。次いで、ボウタイ効果を推定するボウタイ形ベクトルを生成する。ボウタイ形ベクトルは、式１に従って空気較正ベクトルの多項式フィットを実行することにより生成される。
【００２３】
Ｂ（γ）＝ｐｏｌｙ＿ｆｉｔ（ｘ2，Ａｉｒｃａｌ（γ））          式１
ここで、Ｂ（γ）はボウタイ形ベクトルであり、ｘは検出器チャネル番号であり、イソチャネルをゼロに等しくする。一実施形態では、五次多項式フィットを用いる。もう一つの実施形態では、ｎ次多項式を用い、すなわちｎ＞５とする。代替的には、ボウタイ・ベクトルＢ（γ）をＸ線フォトン束測定値から直接決定することもできる。利用時には、空気走査すなわち走査平面内に物体を含まない走査の第一の集合をボウタイを用いて収集する。空気走査の第二の集合はボウタイを用いずに収集する。空気走査の第一の集合と空気走査の第二の集合との比がボウタイ・ベクトルＢ（γ）に相当する。
【００２４】
投影データを前処理した後には、投影データの雑音分布はポアソン分布に従うものと仮定することができる。従って、投影データの標準偏差σ（γ）を下記の式によって近似することができる。
【００２５】
σ（γ）＝［Ｂ（γ）・（ｐ（γ））1/2］／ｐ（γ0）              式２
ここで、ｐ（γ）はマイナス対数処理前の較正後の投影データであり、ｐ（γ0）は参照チャネル信号である。この方法はすべての行で同一であるので検出器行番号は用いない。同様に、簡単化のために投影角番号も省略されている。実施形態の一例では、データ集合のあらゆる投影標本、すなわち各々のチャネル、各々の検出器行及び各々の投影ビューに式２を適用する。
【００２６】
図４は、品質保証（ＱＡ）ファントムの較正後の投影のグラフ図である。図５は、式（２）に基づいて近似された標準偏差のグラフ図であり、同図においては、近似標準偏差を実線で表わし、対応する測定標準偏差を点線で表わしている。
【００２７】
雑音分布の特性決定を用いて投影データに適応型雑音低減演算を施す工程６６は、上限閾値ｔhigh及び下限閾値ｔlowを画定する工程を含んでいる。投影標本の近傍に跨がる投影標本の測定変動がσ（γ）ｔlowよりも小さい場合には、完全平滑化演算を適用する。測定変動がσ（γ）ｔhighよりも大きい場合には、平滑化演算を行なわない。変動が近似的にσ（γ）ｔhighとσ（γ）ｔlowとの間にある場合には、変動の水準に基づいて平滑化の量を調節する。一実施形態では、平滑化はチャネル毎の方式で行なわれる。変動は、小さい関心領域における標準偏差を計算することにより測定することができる。一実施形態では、変動測定のために複数の異なる配向を用いる。一実施形態では、チャネル、ビュー及び検出器行に跨がる分散を算出することができる。もう一つの実施形態では、複数の中間的な方向、例えば両チャネル及びビューに跨がる斜方向等であるがこれに限定されない方向を用いることができる。さらにもう一つの実施形態では、複数の多重分解能形式のアプローチを用いることができる。例えば、２、３又は４の隣接する検出器チャネルを加算して単一のチャネルを形成することができ、圧縮された投影に変動解析及びフィルタ処理を施すことができる。
【００２８】
投影データをフィルタ処理する工程は、変動が前段で画定された閾値の範囲内すなわちｔlowとｔhighとの間になるような方向に沿ってフィルタ処理を行なうことを含んでいる。フィルタ処理された投影を導き出すためには複数のフィルタ処理方法を用いることができる。例えば、任意の形式の低域通過フィルタ又は統計方式フィルタを用いてよい。導き出されたフィルタ処理後の投影ｐf（γ）を用いて、修正後の投影ｐr（γ）を下記の式に従って得ることができる。
【００２９】
ｐr（γ）＝ｐ（γ）−α（γ）［ｐ（γ）−ｐf（γ）］            式３
ここで、α（γ）は前段で算出した変動及び閾値に基づいて決定される。
【００３０】
図６は、フィルタ処理方法の実効性を説明するのに用いられる再構成したファントムの複数の画像である。画像１００は、走査して従来の再構成アルゴリズムで再構成したＱＡファントムの例を示す。画像１０２は、本書に記載した一つの方法に従って処理した同じＱＡファントムを示す。適応型処理によって画像雑音が大幅に減少している。画像１０４は、画像１００と画像１０２との差画像を示す。一実施形態では、差画像には構造情報が存在しておらず、空間分解能を保存する本アルゴリズムの性能を示している。画像１０６は、１．８倍の管電流で走査して従来の再構成アルゴリズムで再構成したＱＡファントムを示す。画像雑音レベルは画像１０２と同等であることに留意されたい。
【００３１】
様々な特定の実施形態によって本発明を説明したが、当業者であれば、特許請求の範囲の要旨及び範囲に含まれる改変を施して本発明を実施し得ることを理解されよう。
【図面の簡単な説明】
【図１】ＣＴイメージング・システムの見取り図である。
【図２】図１に示すシステムのブロック概略図である。
【図３】再構成される画像の雑音の低減を促進する方法の流れ図である。
【図４】品質保証（ＱＡ）ファントムの較正後の投影のグラフ図である。
【図５】近似した標準偏差のグラフ図である。
【図６】フィルタ処理方法の実効性を示すのに用いられる再構成後のファントムの複数の画像の図である。
【符号の説明】
１０  ＣＴシステム
１２  ガントリ
１４  Ｘ線源
１６  Ｘ線ビーム
１８  検出器アレイ
２０  検出器素子
２２  患者
２４  回転中心
２６  制御機構
４２  表示器
４６  モータ式テーブル
４８  ガントリ開口
５０  読み書き装置
５２  着脱式媒体
６０  ＣＴイメージング・システムの統計的雑音の低減を促進する方法
１００  従来のアルゴリズムで再構成したＱＡファントム画像
１０２  本発明の方法で再構成したＱＡファントム画像
１０４  画像１００と画像１０２との差画像
１０６  １．８倍の管電流を用いて従来のアルゴリズムで再構成したＱＡファントム画像
【図１】
0004150900.tif 000002
【図２】
0004150900.tif 000003
【図３】
0004150900.tif 000004
【図４】
0004150900.tif 000005
【図５】
0004150900.tif 000006
【図６】
0004150900.tif 000007