JPH11511375A

JPH11511375A - コンピュータ断層撮影用リング抑圧フィルタ

Info

Publication number: JPH11511375A
Application number: JP9532630A
Authority: JP
Inventors: ゴードン，バーナード，エム．; チン―ミン，ライ
Original assignee: アナロジックコーポレーション
Priority date: 1996-03-13
Filing date: 1997-02-25
Publication date: 1999-10-05
Also published as: NL1005511C2; KR19990087775A; BR9708163A; WO1997033517A1; US5745542A; NL1005511A1; CN1218377A; DE19781653T1

Abstract

(57)【要約】コンピュータ断層撮影システムに用いられるリング抑圧フィルタ（７４）は、高域通過フィルタ（３１２）、ヒストグラム生成器（３１８）、および結合器（３２２）を含む。高域通過フィルタ（３１２）は、高周波成分を表す各投影データ信号（ＰＤＳ）に関するエラー信号（ＥＳ）を生成する。ヒストグラム生成器（３１８）は、エラー信号（ＥＳ）と投影データ信号（ＰＤＳ）とを対応させるヒストグラム（HIST）を生成する。結合器（３２２）は、エラー信号（ＥＳ）とヒストグラム（HIST）とを結合し、リングを修正された信号（ＲＣＳ）を生成する。

Description

【発明の詳細な説明】コンピュータ断層撮影用リング抑圧フィルタ関連する特許出願本発明は、１９９６年１月１７日付で出願された本発明の譲受人に譲渡された米国特許出願第０８/５８７,４６８号“STREAK SUPPRESSION FILTER FOR USE IN COMPUTED TOMOGRAPHY SYSTEMS（コンピュータ断層撮影システムに用いるしま抑圧フィルタ）”、１９９６年３月１３日付で本発明の譲受人に譲渡された米国特許出願第０８/６１４,６２３号“MOTION ARTIFACT SUPPRESSION FILTER FOR USE IN COMPUTED TOMOGRAPHY SYSTEMS（コンピュータ断層撮影システムに用いる移動ファーティファクト抑圧フィルタ）”、および１９９６年３月１３日付で出願された本発明の譲受人に譲渡された米国特許出願第０８/６１４,６６０号“SELF -CALIBRATING RING SUPPRESSION FILTER FOR USE IN COMPUTED TOMOGRAPHY SYST EMS（コンピュータ断層撮影システムに用いる自己校正リング抑圧フィルタ）” と関連するために、これらについても参照されたい。発明の分野本発明は、一般的には、人の患者のなどのＣＴ画像を生成するための医療技術において用いられるコンピュータ断層撮影（ＣＴ）システムに関する。より詳細には、本発明は、ＣＴ画像内のリングを低減するための改良されたリング抑圧フィルタに関する。発明の背景第三世代のタイプのＣＴシステムは、Ｘ−線源、およびＸ−線検出器システムを含み、これからがおのおの環状のディスクの直径的に反対側に固定される。ディスクは、橋形サポート内に回転できるように搭載され、走査の際、ディスクは、回転軸のまわりを連続的に回転し、このとき、Ｘ−線が線源からディスクの開口部内に配置された対象を貫通して検出器システムに至る。検出器システムは、典型的には、単一の行として円弧状に配列された検出器アレイを含む。この円は、曲率中心を“焦点スポット（focal spot）”と呼ばれる点に持ち、Ｘ−線源からの放射線は、ここを通過する。Ｘ−線源と、検出器アレイとは、線源と各検出器との間のＸ−線の経路がディスクの回転軸に垂直な共通平面（以降、“スライス平面”あるいは“走査平面”と呼ばれる）内にくるように配置される。Ｘ−線の経路は実質的に点源であるＸ−線源から出発し、異なる角度を通って、検出器に到達するために、Ｘ−線の経路は扇に類似し、このために、しばしば、任意の時間における全てのＸ−線の経路を記述するために“扇形光線（fan beam）”という用語が用いられる。ある走査のある測定瞬間において単一の検出器上に当たるＸ−線は、通常、“光線（ray）”と呼ばれ、各検出器は、対応する光線の強度を表す出力信号を生成する。各光線は、その経路内の全ての物質によって部分的に減衰される。各検出器によって生成される出力信号は、実際には、検出器とＸ−線源との間に位置する全ての物質の密度（つまり、検出器の対応する光線路内に位置する物質の密度）を表すとともに、この関数である。Ｘ−線検出器によって生成された出力信号は、通常は、ＣＴシステムの信号処理部によって処理される。この信号処理部は、通常は、データ取得システム（ＤＡＳ）を含み、これは、Ｘ−線検出器によって生成された出力信号をフィルタリングすることで、これらの信号対雑音比を改善する。ＤＡＳによって生成されたフィルタリングされた出力信号は、通常は、“生のデータ信号（raw data signa l）”と呼ばれる。信号処理部は、通常、投影フィルタを含む。これは、この生のデータ信号を対数的に処理することで、セットの投影データ信号を生成するが、各投影データ信号は、対応する光線路内に位置する物質の密度を表す。ある測定瞬間における全ての投影データ信号の集合は通常“投影（projection）” あるいは“ビュー（view）”と呼ばれる。単一の走査において、ディスクが回転するにことに伴い複数の投影が生成されるが、各投影はディスクの異なる角度位置において生成される。ディスク上の線源と検出器の特定の投影に対応する角度方位は、“投影角度（projection angle）”と呼ばれる。周知のアルゴリズム、例えば、Rodon アルゴリズムを用いることで、ＣＴ画像が各投影角度にて収集された全ての投影データ信号から生成される。ＣＴ画像は、走査された対象の走査平面に沿っての二次元“スライス(slice)”の密度を表す。投影データ信号からＣＴ画像を生成するプロセスは、ＣＴ画像は投影データから再生されるために、通常、“フィルタ後の逆投影（filtered back projecti on）”あるいは“再生（reconstruction）”と呼ばれる。ＣＴシステムの信号処理部は、通常、投影データ信号から再生されたＣＴ画像を生成するための逆投影器を含む。ＣＴシステムと関連する一つの問題は、様々なノイズおよびエラー源が潜在的に再生されたＣＴ画像へのノイズあるいはファーティファクトとして寄与することである。ＣＴシステムは、従って、典型的には、信号対雑音比を改善し、再生されたＣＴ画像内に存在するファーティファクトを低減するための様々な信号処理技術を用いる。再生されたＣＴ画像内に望ましくないファーティファクトが出現する原因となる一つの重要な要因は、Ｘ−線検出器の一様性と安定性に関係する。ある単一の検出器の調節が、検出器アレイ内の他の検出器に対してずれている場合は、その単一の検出器が原因で再生されたＣＴ画像内にファーティファクトが出現する。このファーティファクトは、再生されたＣＴ画像の“中心（center）”を中心とする円形リングとして、あるいは、一つあるいはそれ以上の円弧として出現する（ここで、再生されたＣＴ画像の“中心”は、ディスクの回転軸の位置に対応する）。さらに、一つ以上の検出器の調節がずれている場合は、これらが集合的に働いて、再生されたＣＴ画像内に、一群の同心円状のリングあるいは円弧が出現する。これらファーティファクトは、典型的には、“リング（rings）”と呼ばれ、“リング除去（deringing）”あるいは“リング抑圧（ring suppressing） ”という用語は、再生されたＣＴ画像内に出現するこれらリングを低減あるいは除去するための方法および装置との関連で用いられる。理想的には、Ｘ−線検出器は、それらの伝達関数、より簡単には、それらの“ 応答（responses）”が全て等しくなるように構成される。ただし、実際には、これを達成することは困難である。このため、多くのＣＴシステムにおいては、投影フィルタは、“応答校正テーブル（response calibration table）”を含み、これが投影データ信号を調節するために用いられ、こうして検出器の応答の差が補償される。この応答校正テーブルは、典型的には、通常、“幻影（phantoms ）”と呼ばれる既知の密度の対象を走査することによって生成され、この応答校正テーブルは定期的に更新される。この応答校正テーブルは、通常は、テーブルが更新された直後は、リングを抑圧するのに有効である。ただし、検出器の応答は、温度の変動、放射線の損傷、その他の要因により、時間と共にドリフトする傾向があり、通常は、任意の単一の検出器の応答が他の検出器のある与えられた公差内に留まることを保証することは困難であり、しばしば、一つあるいは複数の検出器の応答が応答校正テーブルが更新される前に大きくドリフトし、再生されたＣＴ画像内にリングを出現させる。図１は、人の頭部の再生されたＣＴ画像であり、リングと関連する問題を示す。当業者においては理解できるように、図１に示すＣＴ画像は画像の解釈を妨げるリングを含む。リングを抑圧するための従来の技術による方法は、典型的には、投影データ信号を高域通過フィルタリングすることによりセットの高周波信号を生成し、次に、この高周波信号を複数の隣接する投影角度に渡って平均化することで各検出器に対するオフセットを決定する。各高周波信号は、データ成分上に重なったエラー成分を含む（ここで、データ成分は、患者の密度プロフィルの高周波成分を表し、エラー成分はある検出器の他の検出器に対する応答のオフセットを表す）。データ成分は、どちらかというと投影角度に渡って見た場合ある程度ランダムであり、エラー成分はどちらかというと一定にとどまる傾向があるために、高周波成分を平均化した場合、データ成分は相殺され、エラー成分のみが残る傾向がある。こうして平均化された信号が、次に、投影データ信号を調節するためのオフセット値として用いられ、これによって検出器の応答の偏差（ずれ）が補償される。このような従来の技術による方法は、幾つかの短所を持つ。第一に、殆どのＣＴスキャナにおいて、リングエラーは、投影角度には関係なく、従って、高周波信号を複数の隣接する投影角度に渡って平均化しても満足できるオフセット値は生成されない。第二に、走査対象内の高いコントラストを持つ形状は、しばしば、大きな振幅を持つ高周波信号を生成し、このような大きな振幅の高周波信号が平均化プロセスに含まれると不正確なオフセット値が生成される原因となる。従って、ＣＴ画像内のリングを低減するための改善された方法および装置に対する必要性が存在する。発明の目的本発明の一つの目的は、従来の技術の上述の問題を大きく低減あるいは克服することにある。本発明のもう一つの目的は、各チャネル内のエラーを、投影の振幅の関数として、記述する各チャネルに対するヒストグラム信号を生成するリング抑圧フィルタを提供することにある。本発明のさらにもう一つの目的は、リングを表す高品質のデータが得られない局所領域においてさえも、大局的なデータを用いることで、リングを抑圧することができるリング抑圧フィルタを提供することにある。本発明のさらにもう一つの目的は、現在の走査におけるリングを抑圧するために、一つあるいは複数の過去の走査からのデータを用いることができるコンピュータ断層撮影システムを提供することにある。発明の要約これらおよびその他の目的がＣＴシステムにおいて用いるための改善されたリング抑圧フィルタによって達成される。本発明によるシステムは、複数の投影角度に対する投影を生成するが、各投影は、複数の投影データ信号の一つの測定値を含む。リング抑圧フィルタは、高域通過フィルタ、ヒストグラム生成器、および結合器（減算器）を含む。高域通過フィルタは、複数の投影データ信号の一つの測定値に対応するエラー信号の測定値を生成するが、各エラー信号の測定値は、複数の投影の一つのある部分の高周波成分を表す。ヒストグラム生成器は、エラー信号の測定値と、一つの投影データ信号の測定値との間の関係を表すヒストグラム信号を生成する。結合器は、一つの投影データ信号の測定値と、こうして生成されたヒストグラム信号とを結合することにより、リングを修正された信号の複数の測定値を生成する。リング抑圧フィルタは、複数のリングを修正された信号の測定値を生成するが、各リングを修正された信号は、投影データ信号の一つと対応する。投影データ信号ではなく、リングを修正された信号を用いることで、ＣＴシステムは、リングの低減されたＣＴ画像を生成することができる。当業者においては以下の詳細な説明から本発明のその他の目的および長所が明きからになると思われる。以下では、複数の実施例について説明されるが、これらは単に本発明の最良の実施形態を説明するためのものである。理解できるように本発明は他の異なる実施形態を持つことも可能であり、後に説明される本発明の実施例の細部は、本発明から逸脱することなく、様々な点で修正することが可能である。従って、図面および説明は単に解説を目的とするもので、制約あるいは限定を意図するものではない。本発明の範囲は請求の範囲によって規定される。図面の簡単な説明本発明の特徴および目的は以下の詳細な説明を、付録の図面を参照しながら読むことによって一層理解できるものである。図面中、同一の参照符号は同一あるいは類似する要素を示す。図１は、ＣＴ画像であり、リングと関連する問題を示し；図２は、本発明に従って構成された一つの好ましいＣＴシステムを軸方向に示す図面であり；図３は、図２に示すリング抑圧フィルタをより詳細に示すブロック図であり；図４Ａは、本発明に従って構成されたリング抑圧フィルタによって生成される一例としてのヒストグラム信号を示すグラフであり；図４Ｂ〜４Ｃは、本発明に従って構成されたリング抑圧フィルタによってヒストグラム信号を生成するために用いられる中間変数のグラフであり；図５は、本発明に従って構成されたリング抑圧フィルタによってヒストグラム信号を生成するために用いられるプロセスを示す流れ図であり；図６は、図１に示す画像を生成するのに用いたのと同一のデータを用いて本発明に従って構成されたリング抑圧フィルタを用いるＣＴシステムによって生成されたＣＴ画像の一例を示し；図７は、本発明に従って構成された以前の走査からのデータを現在の走査におけるリングを抑圧するために用いるＣＴシステムを示すブロック図であり；図８は、本発明に従って構成された一つの好ましいＣＴシステムの信号処理部のブロック図であり；図９Ａ〜９Ｂは、それぞれ、扇形の光線の投影と、平行な光線の投影の生成を示し；図１０Ａ〜１０Ｂは、本発明に従って構成されたＣＴシステムに用いることができる平行な光線の投影を生成するための一つの方法を示し；図１１Ａ〜１１Ｂは、それぞれ、０度と１８０度の投影角度に対するＸ−線源と検出器との間の空間的な関係を示し；図１２は、０度と１８０度の投影角度における検出器アレイの間の空間的な関係を示し；図１３は、本発明に従って構成されたもう一つの好ましいＣＴシステムの信号処理部をブロック図にて示す。図面の詳細な説明図２は、本発明の原理を組み込む一例としてのＣＴシステム、すなわち、スキャナ４０を示す。スキャナ４０は、ディスク４６に搭載されたＸ−線源４２と、検出器アセンブリ４４を含み、検出器アセンブリは、検出器のアレイから構成される。Ｘ−線源４２と検出器アセンブリ４４は、回転軸４８（図２に示す面に対して垂直に延びる）のまわりを回転し、ＣＴ走査の際に、ディスク４６の中央開口部を貫通して延びる対象５０のまわりを回転する。対象５０は、生きた患者の一部分、例えば、頭あるいは胴であり得る。Ｘ−線源４２は、走査平面（回転軸４８に対して垂直）内で、Ｘ−線の連続的な扇形の光線５２を放出し、これら光線が対象５０を貫通して通過した後に、検出器アセンブリ４４によって検出される。好ましくは、散乱防止プレートのアレイ５４が、散乱した光線が検出器によって検出されることを防止するために、対象５０と検出器アセンブリ４４との間に配置される。一つの好ましい実施例においては、検出器の数は、３８４とされ、４８°の弧をカバーする。ただし、検出器の数および弧の角度はこれとは異なってもかまわない。ディスク４６は、好ましくは、軽量の材料、例えば、アルミニウムとされ、軸４８のまわりを高速に滑らかに回転するようにされる。ディスク４６は、対象５０をディスクの開口部を貫通して配置できるように開放型のフレーム構造とされる。対象５０は、例えば、好ましくは、Ｘ−線に対して可能な限り透明なテーブル５６てに支持される。検出器アセンブリ４４によって生成される出力信号は、ＤＡＳ７０（ブロック図の形式にて示す）に加えられ、ＤＡＳ７０は、この出力信号からセットの生のデータ信号を生成する。この生のデータ信号は投影フィルタ７２に加えられ、投影フィルタ７２は、セットの投影データ信号を生成する。ディスク４６が回転することで、多数の投影角度からの投影データ信号が得られる。次に、これら投影データ信号は、リング抑圧フィルタ７４に加えられ、リング抑圧フィルタ７４は、再生されたＣＴ画像内のリングを低減させるために、本発明による方法で、投影データ信号をフィルタリングする。リング抑圧フィルタ７４によって生成された出力信号は、“リングを修正された投影データ信号（ring corrected project ion data signal）”、あるいは、単に、“リングを修正された信号（ring corr ected signal）”と呼ばれる。この信号は、次に、逆投影器７６に加えられ、逆投影器７６は、リングを修正された信号からＣＴ画像を生成する。逆投影器７６は逆投影のためにデータの畳み込を遂行する畳み込みフィルタを含む入力段を持つ。スキャナ４０は、マルチチャネルデバイスであり、チャネル毎に、データとして、一つの生のデータ信号、一つの投影データ信号、および一つのリングを修正された信号を送信する。各チャネルにおいて、これら信号のおのおのが、好ましくは、複数の投影角度の各角度に対して生成される。後に詳しく説明するように、投影フィルタ７２によって生成された投影データ信号を用いるのではなく、リング抑圧フィルタ７４によって生成されたリングを修正された信号を用いることで、逆投影器７６によって生成されるＣＴ画像は、目立つリングがより少ない、より鮮明となる。上述のように、再生されたＣＴ画像のリングは、一つの原因としては、ある検出器の調節が検出器アレイ内の隣接する検出器に対してずれている場合に発生する。同様に、あるチャネル内の任意の要素が原因で、そのチャネル調節が、スキャナ内の他のチャネルに対してずれる場合にもリングが発生する。さらに、リングは、勿論、単一のチャネルでなく、複数のチャネルの調節が他のチャネルに対してずれている場合にも発生する。加えて、一つあるいは複数の共通の要素が検出器アセンブリを構成する複数のサブセットの検出器からの出力を処理するために用いられる場合、例えば、あるマルチプレクサを用いて、これら一つあるいは複数の共通の要素を通じて処理された検出器出力が処理される場合、これら共通の要素の一つの調節が他の要素に対してずれている場合にも、最終的な画像内に複数のリングが発生する。このように、ここでは“チャネル（channel）”という用語は、ある検出器から出発し、その後、説明の様々な要素を通じて処理されるある信号の経路を指すために用いられ、ここでは、必ずしも必須ではないが、複数のチャネルが一つの共通の信号経路の少なくとも一部分を、例えば、複数のチャネルが一つの共通の要素を共有するものと理解される。上述の事象のいずれかが存在する場合、チャネルの応答には上述のように偏差（ずれ）が発生し、このために、各チャネル内の投影データ信号は、患者のある部分の密度を表すデータ成分に加えて、この上に重なった“リングエラー（ring error）”と呼ばれるエラー成分を含むこととなる。リング抑止フィルタ７４の動作は、リングを修正された信号がデーダ成分のみを含むように各チャネル内のリングを除去することにある。こうして、リングを修正された信号から再生されたＣＴ画像は、より少数のリングを含み、品質が改善される。図３は、本発明に従って構成されたリング抑圧フィルタ７４の一つの実施例をブロック図にて示す。リング抑圧フィルタ７４は、投影フィルタ７２によって生成された投影データ信号を受信し、これから、リングを修正された信号を生成し、これを逆投影器７６に加える。リング抑圧フィルタ７４は、低域通過フィルタ３１０、高域通過フィルタ３１２、ヒストグラム生成器３１８、リミッタ３２０、および、減算器（結合器）３２２を含む。投影データ信号は、低域通過フィルタ３１０と、減算器３２２の正の入力端子とに加えられる。低域通過フィルタ３１０は、投影データ信号をフィルタリングすることで、セットの低周波投影信号を生成し、これが、高域通過フィルタ３１２と、ヒストグラム生成器３１８の２つの入力の一つとに加えられる。高域通過フィルタ３１２は、低周波投影信号をフィルタリングすることで、セットのエラー信号を生成し、これらエラー信号がヒストグラム生成器３１８の他の入力に加えられる。ヒストグラム生成器３１８は、低周波投影信号とエラー信号からセットのヒストグラム信号を生成する。これらヒストグラム信号は、リミッタ３２０に加えられ、リミッタ３２０は、これからセットのオフセット信号を生成し、このオフセット信号が、減算器３２２の負の入力端子に加えられる。減算器３２２は、投影データ信号からオフセット信号を減ずることで、リングを修正された信号を生成する。単一の走査（つまり、ディスクの一回転）の際に収集される投影データ信号の測定値は、以下の式（１）によって示されるようにマトリックスＰＤＳに編成される：このＰＤＳマトリックス内の各要素ＰＤＳ(ｉ、θ) は、ｉ番目のチャネルにおけるθに等しい投影角度に対する投影データ信号の測定値を表す。式（１）において、Ｎは、スキャナ４０内のチャネルの数を表す。上述のように、この好ましい実施例においては、検出器アレイ４４内には、３８４個の検出器が存在し、このため、この好ましい実施例においては、スキャナ４０内には３８４個のチャネルが存在し、このため、Ｎは、３８４となる。上述のように、チャネルの０からＮ −１への通し番号は、検出器の配列順に対応し、検出器アセンブリ４４の検出器の一端から他端に番号が与えられたとき、チャネル０は、端の検出器１に対応し、チャネル１は、検出器１と検出器３との間に位置する検出器２に対応し、以下同様となる。Δθは、一連の投影の間のディスク４６の回転の量（つまり、一連の投影の間の投影角度の角度増分）を表す。この好ましい実施例においては、ディスク４６は、各投影間で、１°の八分の一だけ回転し、スキャナ４０は、単一の走査の際に、２８８０個（つまり、１°当たり８個の投影×３６０°＝２８８０個）を生成し、このため、この好ましい実施例においては、Δθは、０．１２５ °（１°÷８＝０．１２５°）に等しい。ＰＤＳマトリックスの行は、単一の投影角度において収集される投影データ信号の全ての測定値を表す。従って、この好ましい実施例においては、ＤＰＳマトリックスは、２８８０個の行を含む。ＤＰＳマトリックスの列は、単一の走査の際に収集される一つのチャネルを通じての全ての測定値（投影データ信号の値）を表し、この好ましい実施例においては、ＤＰＳマトリックスは、３８４個の列を含む。空間的には、ＰＤＳマトリックスは、巡回的な性質を持ち、第一の行は、最後の行に続く。つまり、ＰＤＳ(ｉ、０)＝ＰＤＳ(ｉ、３６０)となる。これは、これら２個の値が同一の投影角度にて取られるためである。ただし、これらの値はこれらが異なる時間において測定されるために異なることもある。低域通過フィルタ３１０は、低周波投影信号を、好ましくは、投影データ信号を複数の投影角度に渡って（つまり、ＰＤＳマトリックスの複数の行に渡って）低域通過フィルタリングすることによって生成する。低域通過フィルタ３１０は、好ましくは、低周波投影信号の測定値のマトリックスＰＬＦを生成するが、ＰＬＦマトリックスの各要素ＰＬＦ(ｉ、θ)は、ｉ番目のチャネルのθに等しい投影角度に対する低周波投影信号の測定値を表す。当業者においては理解できるように、リング抑止フィルタ７４内に用いられる低域通過フィルタ３１０は、この特定の形式に限定されるものではない。例えば、低域通過フィルタ３１０は、有限インパルス応答（ＦＩＲ）フィルタ、無限インパルス応答（ＩＩＲ）フィルタ、反復あるいは非反復フィルタ、さらには、例えば、フーリエ変換技術を用いる周波数領域フィルタを用いて実現することもできる。一つの実施例においては、低域通過フィルタ３１０は、ＰＬＦマトリックスの各要素ＰＬＦ(ｉ、θ)を、これがＰＤＳマトリックスの要素の加重平均となるように以下の式（２）に従って生成する：ここで、Ｍは、好ましくは、５０（Δθが０．１２５度の場合）とされ、重み係数Ａ_k は、好ましくは、全て１とされる。こうして、この好ましい実施例においては、チャネルｉ、特定の投影角度θに対して、ＰＬＦ(ｉ、θ)が、投影角度θ と、前の５０および次の５０の投影角度、すなわち、全体で約１２．５度の角度に対するＰＤＳ値から決定される。ただし、勿論、重み係数Ａ_kに対して、他の値（の配分）を用いることも、Ｍに対して他の値を用いることもできる。式（２）においては、好ましくは、ＰＤＳマトリックスの巡回的な性質が利用され、変数θはモジュロ３６０として解釈される（例えば、θ＝（３６０＋Δθ）は、θ ＝Δθとして解釈される）。こうして、０度の投影角度に対して収集されたデータは、３６０度に対して収集されたデータに“ラップアラウンド（wraps around ）”される。高域通過フィルタ３１２は、エラー信号を、好ましくは、低周波投影信号を複数のチャネルに渡って（つまり、ＰＬＦマトリックスの複数の列に渡って）高域通過フィルタリングすることによって生成する。各エラー信号は、従って、ある一つのチャネルの投影データ信号が、そのチャネルと隣接するチャネルの投影データ信号とどの程度異なるかの測定値となる。高域通過フィルタ３１２を実現するために用いられるフィルタは、各エラー信号が、そのチャネルのリングエラーの良好な測定値を与えるように選択される（これらエラー信号は、そのチャネルのデータ成分の測定値上に重なり、データ成分は患者の密度プロフィルの高周波成分を表す）。高域通過フィルタ３１２は、好ましくは、エラー信号の測定値のマトリックスＥＳを生成し、マトリックスＥＳの各要素ＥＳ(ｉ、θ)は、ｉ番目のチャネルの θに等しい投影角度に対するエラー信号の測定値を表す。一つの好ましい実施例においては、高域通過フィルタ３１２は、２パスフィルタとして実現され、式（３）に従ってマトリックスＥＳの各要素ＥＳ(ｉ、θ)を生成する。第一のパスの結果ＰＨＦ(ｉ、θ)が、さらに、高域通過フィルタリングされ、第二のパスの結果がＥＳ(ｉ、θ)として得られる。当業者においては理解できるように、フィルタリングのこれら２つのパスは、単一のパスに結合することもでき、あるいは、他の形式の高域通過フィルタを用いることもできる。例えば、２つより多数のパスを含むフィルタ、各パス内により少数あるいは多数の隣接する検出器を含むフィルタ、あるいは、隣接する検出器の出力に異なる重みを与えるフィルタ、その他を用いることも考えられる。 PHF(i,θ)＝PLF(i,θ)− 1/6[PLF(i-2,θ)+2PLF(i-1,θ)+2PLF(i+1,θ)+PLF(i+2,θ)] （３） ES(i,θ)＝PHF(i,θ)− 1/5[PHF(i-2,θ)+PHF(i-1,θ)+PHF(i,θ)+PHF(i+1,θ)+PHF(i+2,θ)] 式（３）から分かるように、各高域通過フィルタリングされた信号ＰＨＦ(ｉ、 θ)は、そのチャネルに対する低域通過フィルタリングされた信号ＰＬＦ(ｉ、θ) から、チャネル（ｉ）の両側の次の２個の隣接するチャネル（検出器）の低域通過フィルタリングされた信号（ここで、この２個の隣接するチャネル（ｉ−１）および（ｉ＋１）の低域通過フィルタリングされた信号にはおのおの２の係数が掛けられる）と、これら２つのチャネル（検出器）の外側に隣接するチャネル（ｉ−２）および（ｉ＋２）の低域通過フィルタリングされた信号（ここで、この２つの低域通過フィルタリングされた信号にはおのおの係数１が掛けられる）との加重総和を６で割った値を引くことによって得られる。容易に理解できるように、５個の低域通過フィルタリングされた信号のおのおのの値が同一の場合は、ＰＨＦ(ｉ、θ)の値は、零となり、差は存在しないことを示す。ＥＳ(ｉ、θ)の値は、各チャネルに対するＰＨＦ(ｉ、θ)の値から、ＰＨＦ（i- 2、θ）、ＰＨＦ(i-1、θ)、ＰＨＦ(ｉ、θ)、ＰＨＦ(i+1、θ)、およびＰＨＦ(i+ 2、θ)の値（ここで、これらには等しい重みが与えられる）の総和を５で割った値を引くことによって決定される。これらエラー信号と低周波投影信号が、ヒストグラム生成器３１８に加えられ、ヒストグラム生成器３１８は、これからヒストグラム信号を生成する。この好ましい実施例においては、ヒストグラム生成器３１８は、量子化／蓄積デバイス３３０、２つの低域通過フィルタ３３２、３３４、およびディバイダ３３６を含む。量子化／蓄積デバイス３３０は、各チャネルに対して、２つの中間信号HIST₁ および HIST₂を生成する。HIST₁信号は、低域通過フィルタ３３２に加えられ、低域通過フィルタ３３２は、これから、セットの低域通過フィルタリングされた信号 HISTLF₁を生成する。HIST₂信号は、低域通過フィルタ３３４に加えられ、低域通過フィルタ３３４は、これから、セットの低域通過フィルタリングされた信号 HISTLF₂を生成する。ディバイダ３３６は、これらのフィルタリングされた信号を受信し、各チャネルのこうしてフィルタリングされた信号 HISTLF₁をこうしてフィルタリングされた信号 HISTLF₂によって割ることによってヒストグラム信号を生成する。図４Ａ、４Ｂ、４Ｃは、それぞれ、ヒストグラム生成器３１８によって生成されるヒストグラム信号 HIST、中間信号 HIST₁、および HIST₂の例を示すグラフである。より具体的には、図４Ａ、図４Ｂ、および図４Ｃは、それぞれ、全て第一のチャネルに対するヒストグラム信号 HIST（1、j）、中間信号 HIST₁（1、j）、および HIST₂（1、j）の例を示すグラフである。このグラフは、第一のチャネルに対して生成された信号を示すが、当業者においては理解できるように、他のチャネルのおのおのに対して生成される信号も、同様に、これらグラフのように表すことができる。図４Ａ〜４Ｃにおいて、各グラフの横軸は、第一のチャネルの投影振幅（つまり、第一のチャネルの低周波投影信号の振幅）を表すが、これらは、好ましくは、幾つかのビンに量子化される。変数ｊは、量子化された投影振幅を表し、変数ｊの各値は、投影振幅の量子化間隔、すなわち、レンジを表す。ｊの最大値は、任意の投影振幅の期待される最大の値を含むレンジに対応し、ｊの最小値は、任意の投影振幅の期待される最小の値を含むレンジに対応し、この２つの間の差によって動的レンジが定義される。ヒストグラム生成器３１８は、ヒストグラム信号を、量子化変数ｊの関数となるように生成する。図４Ａ〜４Ｃにおいては、投影振幅は、１００個の間隔に対応する１００個のビンに量子化されるように示されるが、ただし、他の任意の個数のビンを用いることもできる。図４Ａにおいて、縦軸は、第一のチャネルのリングエラーを表す。こうして、図４Ａに示される信号 HIST（1、j）は第一のチャネルのリングエラーを投影振幅の関数として記述する。図５は、量子化／蓄積デバイス３３０内で、全てのチャネルに対して中間ヒストグラム信号HIST₁、HIST₂を生成するために、用いられる一つの好ましいプロセス５００をブロック図にて示す。以下では、簡単のために、プロセス５００は、第一のチャネルに対する中間ヒストグラム信号 HIST₁（1、j）、HIST₂（1、j）の生成との関連で説明されるが、当業者においては理解できるように、他のチャネルに対する中間ヒストグラムも類似の方法にて計算される。プロセス５００の第一のステップは、初期化ステップ５１０であり、ここで、このプロセスにおいて用いられる変数が０にセットされる。より詳細には、投影角度を表す変数が０度にセットされ、HIST₁（1、j）と HIST₂（1、j）の値がｊの全ての値に対して０にセットされる。ステップ５１０に続いて、比較ステップ、あるいは、閾値ステップ５１２が実行される。ステップ５１２において、エラー信号マトリックスＥＳの要素ＥＳ(1 、θ)が、所定の閾値の値ＴＨＲと比較される。要素ＥＳ(1、θ)の絶対値がＴＨＲの値より小さな場合は、要素ＥＳ(1、θ)は、“閾値以下（under-threshold） ”であると呼ばれ、要素ＥＳ(1、θ)の絶対値がＴＨＲの値より小さくない場合は、要素ＥＳ(1、θ)は、“閾値以上（over-threshold）”であると呼ばれる。要素ＥＳ(1、θ)が閾値以下である場合は、その後、量子化ステップ５１４、更新ステップ５１６、および増分ステップ５１８が順次に実行され、プロセス５００は、要素ＥＳ(1、θ)とＰＬＦマトリックスの対応する要素ＰＬＦ(1、θ)を用いて、中間ヒストグラム信号を生成する。反対に、要素ＥＳ(1、θ)が閾値以上である場合は、量子化ステップ５１４および更新ステップ５１６はスキップされ、プロセス５００は、要素ＥＳ(1、θ)と対応する要素ＰＬＦ(1、θ)を用いることなく、中間ヒストグラムを生成する。閾値ステップ５１２は、隣接する検出器間で画像データ（例えば、骨と隣接する柔らかな組織とのエッジ）に起因して生成されるコントラスト（この場合は、ＥＳ(i、θ)の絶対値は閾値以上となる）と、リングアーティフクトに起因するコントラスト（この場合は、ＥＳ(i、θ)の絶対値は閾値以となる）とを区別するために遂行される。ＥＳ(i、θ)の絶対値は経験的に決定される。量子化ステップ５１４において、ビン番号ｊは、好ましくは、要素ＰＬＦ(1、θ)を以下の式（４）に示す量子化関数ＱＵＡＮＴを用いて量子化することによって計算される：ｊ＝ＱＵＡＮＴ[ＰＬＦ(1,θ)] （４）当業者においては理解できるように、この量子化関数QUANTは、多くの異なる方法にて実現することができるが、一つの好ましい量子化関数は、以下の式（５）によって与えられる：式（５）において、ｍａｘと、ｊ_maxの値は、定数であり、ｍａｘは、低周波投影信号に対して期待される最大値に等しく、ｊ_maxは、ｊの最大値に等しい（こうして、図４Ａ〜４Ｃに示すようなヒストグラム信号を生成する実施例においてはｊは１００となる）。式（５）において、ＩＮＴは、周知の整数関数であり、一般的には、ＩＮＴ(ｘ)は、引数ｘの整数部分に等しい（例えば、ＩＮＴ(３.６ )は３に等しい）。量子化ステップ５１４においてビンの番号を決定した後に、更新ステップ５１６が実行される。更新ステップ５１６において、HIST₁（1、j）および HIST₂（1、 j）の値が更新される。HIST₁（1、j）の新たな値は HIST₁（1、j）の古い値に要素ＥＳ(１、θ)の値を加えた値にセットされ、HIST₂（1、j）の新たな値は HIST₂（ 1、j）の古い値に１を加えた値にセットされる。更新ステップ５１６に続いて、増分ステップ５１８が実行される。増分ステップ５１８において、変数θの現在の値が量Δθだけ増分される。増分ステップ５１８に続いて、テストステップ５２０が実行される。テストステップ５２０において、変数θが３６０度と比較される。θが３６０より小さな場合は、ステップ５１２およびその後のステップが再び実行され、θが３６０より小さくない場合は、プロセス５００は完了し、１回の走査の全体に対する HIST₁（1、j）およびH IST₂（1、j）の値の計算が完結する。図４Ｂおよび４Ｃは、それぞれ、プロセス５００によって計算される HIST₁（ 1、j）および HIST₂（1、j）を示すグラフである。ｊの各値に対して、HIST₁（1、j ）の値は、ＥＳ（１、θ）が閾値以下で、かつ、QUANT［PLF(1、θ)］がｊに等しいときのθの全ての値に対する要素ＥＳ(１、θ)の総和に等しく、HIST₂（1、j）の値は、HIST₁（1、j）を形成するために加算された要素ＥＳ(１、θ)の数に等しい。中間ヒストグラム信号 HIST₁（1、j）およびHIST₂（1、j）は、それぞれ、図３に示す低域通過フィルタ３３２および３３４に加えられ、これらフィルタは、これら中間信号を複数の投影振幅に渡って平滑化する。一つの実施例においては、低域通過フィルタ３３２、３３４は、それぞれ、フィルタリングされたヒストグラム信号 HISTLF₁（1、j）、HISTLF₂（1、j）を、それぞれ、以下の式（６）、式（７）に従って生成する：ここでは、フィルタリングの幅Ｍは、８とされ、重み係数Ｂ_kは、全て１とされる。ただし、重み係数Ｂ_kとフィルタリングの幅Ｍに対して他の値を用いることもでき、さらに、別の実施例として、別の低域通過フィルタリング関数を用いることもできる。フィルタリングされたヒストグラム信号 HISTLF₁（1、j）、HISTLF₂（1、j）は、次に、図３に示すディバイダ３３６に加えられる。ディバイダ３３６は、全てのｉおよびｊに対して、最終的なヒストグラム信号HIST（i、j）を、以下の式（８）に従って生成する：ｊの各値に対して、HIST（i、j）は、閾値以下で、かつ、ビンｊ内に入る対応する投影振幅を持つ全ての要素ＥＳ(１、θ)の平均となる。各要素ＥＳ(１、θ)は、データ成分上に重なったリングエラーを表し、また、データ成分は、要素ＥＳ( １、θ)が複数の投影角度θに渡って、かつ、ｊに隣接する複数の振幅ビンに渡って平均化されたとき互いに相殺される傾向を持つため、HIST（i、j）はｉ番目のチャネルの投影振幅の関数としてのリングエラーの測定値となる。ディバイダ３３６は、好ましくは、HISTLF₂（i、j）が０となったときを検出し、HIST（i、j）の対応する値を０に設定するようにされ、これによって、０による割算を試みることによってエラーあるいはオーバーフロー状態が発生することが回避される。さらに好ましくは、ディバイダ３３６は、式（８）の割算を遂行する前に、HISTLF₂（i、j）を１あるいは他の小さな値だけ増分するようにされ、これによって、HIST（i、j）に対して誤った大きな値が生成されるのが回避される。ヒストグラム信号は、データが得られない投影振幅に対応するギャップを含むことがある。フィルタ３３２、３３４によって用いられる式（６）、（７）によって定義されるフィルタリングの幅Ｍが比較的小さな場合は、特にこれが顕著となる。例えば、図４Ａ〜４Ｃに示す中間ヒストグラム信号HIST₁（1、j）およびHI ST₂（1、j）は、ビン番号４（つまり、ｊ＝４）の所にギャップを含む。このようなギャップは、ある特定のビンに量子化された要素ＰＬＦ(ｉ、θ)が存在しない場合、あるいは、ある特定のビンに量子化された要素ＰＬＦ(ｉ、θ)の全てが閾値以上である要素ＥＳ(ｉ、θ)に対応する場合に、ｉ番目のチャネルに対する中間ヒストグラム信号内に発生する。ギャップは、図４Ｂ、４Ｃに示すように、単一の振幅ビンに相当する狭いものも、複数の振幅ビンに相当する広いものもある。低域通過フィルタ３３２、３３４は、このギャップを満たす傾向がある。例えば、図４Ｂ〜４Ｃに示すギャップは満たされ、このため、図４Ａに示すヒストグラム信号 HIST（i、j）にはギャップは存在しない。フィルタ３３２、３３４はギャップが比較的小さな場合はこれを満たす傾向があるが、ギャップの幅が比較的大きな場合はこれを完全に満たすことはない。ただし、ヒストグラム信号内のギャップは、リング抑止フィルタ７４の性能の妨害となることはない。これは、ヒストグラム信号内のギャップは特定の振幅における投影データの不在に対応し、従って、ギャップは、単に存在しないデータの修正に寄与するためである。好ましくは、低域通過フィルタ３３２、３３４は、投影振幅の両端の付近（例えば、ｊが１に等しいビンの付近）においては、フィルタリングされたヒストグラム信号の計算に異なる方法を用いる。式（６）、（７）において、パラメータＭは、フィルタリングされたヒストグラム信号を生成するのに用いる平均化ウインドウのサイズを定義する。好ましくは、投影信号の両端の付近のフィルタリングされたヒストグラム信号を計算するときには、平均化ウインドウが両端を超えないように、平均化ウインドウのサイズが（例えば、Ｍの値を小さくすることによって）小さくされる。こうして、例えば、HISTLF₁（i、3）、HISTLF₂（i、3）を計算するこきは、好ましくは、Ｍは２にセットされる。ただし、別の方法として、平均化ウインドウのサイズは一定に維持し、重みＢ_kを調節することで、両端に来るヒストグラム信号の値を無視することもできる（例えば、式（６）、（７）において、Ｂ_kを、ｊ＋ｋが１より小さな場合、０にセットすることもできる）。さらに、好ましくは、低域通過フィルタ３３２、３３４は、最小の投影振幅（つまり、ｊ＝１）に対応するビンを他のビンとは異なる方法にて扱う。一般的には、このビンは、低域通過フィルタ３３２、３３４によって遂行される平均化には含まない方が好ましい。このために、HISTLF₁（i、1）、HISTLF₂（i、1）は、全てのチャネルｉに対して、それぞれ、HIST₁（i、1）、HIST₂（i、1）と等しくなるようにセットされる。さらに、このビン内の値は、他のビンに対して遂行される平均化にも含まないことが望ましい。一般的に、リングエラーは、投影振幅が０に接近すると著しく異なる振る舞いを示し、このため、このビンに対して測定されたリングエラーは、他のビンに対して測定されたリングエラーと平均化すべきではなく、このビンは、異なる方法で扱うことが望ましい。このビンは、Ｘ− 線光子の吸収が最小であった光線、つまり、（患者の部分ではなく）空気のみを通過した光線に対応する。殆どの走査において、検出器アレイの端付近の検出器に対応するチャネルの投影データ信号は、殆ど常に、最小振幅のビンに量子化され（つまり、患者は、殆どの場合、Ｘ−線源と検出器アレイの端付近の検出器との間には位置せず、このため検出器アレイの端付近の検出器は患者を殆ど“見る（see）”ことはななく）、反対に、検出器アレイの中央付近のチャネルの投影データ信号は、殆ど決して、最小の振幅ビンに量子化されることはない（つまり、検出器アレイの中央付近の検出器は殆ど常に患者を“見る（see）”）。このため、典型的には、あるチャネル内に、最小の振幅ビンに対する多数のエラー信号の測定値が存在する場合と、それらが全く存在しない場合とがあるが、最小の振幅ビンに対する多数のエラー信号の測定値を持つチャネルについては、これら測定値を他のビンからのエラー信号と平均化する必要はなく；単に最小の振幅ビン内のの多数のエラー信号の測定値を平均化することで、エラー信号の高品質な測定値が得られる。最小の振幅ビンに対するエラー信号の測定値を持たないチャネルについては、最小の振幅ビンに対するリングエラーを推定する必要はない。上述のように、従来の技術によるシステムは、リングエラーは、投影データ信号の高周波フィルタリングされたバージョン（つまり、高周波信号）を複数の隣接する投影角度に渡って平均化することによって測定される。例えば、従来の技術によるシステムにおいては、一つのチャネルの２２．５度の投影角度に対するリングエラーは、そのチャネルの高周波信号を０から４５度の間の全ての投影角度に対して平均化することによって測定された。このような従来の技術によるシステムは、リングエラーの測定値を生成するために、“局所的なデータ（local data）”を用いるために制約がある。このような従来の技術によるシステムとは対照的に、図３に示すリング抑圧フィルタ７４は、リングエラーの測定値を生成するために、“大域的なデータ（global data）”を用いる。つまり、各投影振幅に対して、ヒストグラム信号が全ての投影角度からの閾値以下のエラー信号を平均化することによって生成される。このために、リング抑圧フィルタ７４は、良好なリング修正データを得ることができない画像の部分のリングを、画像の他の部分からのデータを用いて修正することが可能となる。このことは、従来の技術によるシステムでは不可能であった。さらに、従来の技術によるシステムでは、単一のオフセット測定値が、各チャネルに対して高周波信号を隣接する投影角度に渡って平均化することによって生成されるが、これとは対照的に、本発明によるリング抑圧フィルタフィルタ７４では、この平均化プロセスが投影振幅の異なる複数のビンに分離され、こうして、リングエラーの測定値が各ビンに対して得られる。これは検出器の応答がその検出器に当たるＸ−線の強度の関数として変動するために有益である。図３に示すリミッタ３２０は、ディバイダ３３６からヒストグラム信号を受信し、これからオフセット信号を生成する。リミッタ３２０は、ヒストグラム信号の規模が閾値より小さな場合は、オフセット信号を対応するヒストグラム信号と等しくし、規模が閾値を超える場合は、オフセット信号を閾値に制限する。これは図２に示すような走査対象５０がリング状の構造を含み、このため、これに対応する検出器および振幅に対して、偽のHIST（i、j）値が生成される場合があるために有益である。リミッタ３２０は、好ましくは、ヒストグラム信号をＭＡＸ＿ＥＲＲのオフセット限界に飽和させ、以下の式（９）に従ってオフセット信号 OFFSET（i、J）を生成する：＝ＭＡＸ_ＥＲＲ， ifHIST(i,j)＞ＭＡＸ_ＥＲＲＯＦＦＳＥＴ(i,j)＝−ＭＡＸ_ＥＲＲ，ifHIST(i,j)＜−ＭＡＸ_ＥＲＲ (９) ＝HIST(i,j)， otherwise 当業者においては理解できるように、オフセット信号 OFFSET（i、j）は、HIST（ i、j）に類似するヒストグラム信号の形式を持つ。また、用いられるＭＡＸ＿ＥＲＲの実際の値は経験的に決定される。次に、オフセット信号と投影データ信号は、減算器３２２に加えられ、減算器は、投影データ信号から適当なオフセット信号を引くことにより、リングを修正された信号を生成する。減算器３２２は、リングを修正された信号のマトリックスＲＣＳを生成するが、個々の要素ＲＣＳ(ｉ、θ)は、以下の式（１０）に従って生成される：ＲＣＳ(i,θ)＝ＰＤＳ(i,θ)−ＯＦＦＳＥＴ(i,ＱＵＡＮＴ[ＰＤＳ(i,θ)]) （１０）次に、リングを修正された信号は、逆投影器７６に加えられ、逆投影器は、これから再生されたＣＴ画像を生成する。各投影データ信号は、データ成分とリングエラーを含み、オフセット信号はリングエラーの測定値であるために、リング修正された信号は、データ成分のみを含む。従って、投影データ信号ではなく、リング修正された信号を用いることで、逆投影器７６は、改善された再生ＣＴ画像を生成する。当業者においては理解できるように、リング抑止フィルタ７４は、多数の異なる方法を用いて実現することができる。リング抑止フィルタ７４の各要素は、離散要素を用いて実現することも、あるいは、リング抑止フィルタ７４の一つあるいは複数の要素を、デジタルコンピュータ、例えば、アレイプロセッサ上に実現することもできる。例えば、減算器３２２のある一つの実現においては、メモリ、例えば、検索テーブル（ＬＵＴ）が含まれ、このメモリにオフセット信号のテーブルが格納される。そして、減算器は、投影データ信号が受信されたときに、この検索テーブル（ＬＵＴ）を用いて、リングを修正された信号を生成する。さらに、リング抑止フィルタ７４は、デジタルシステムであり、各信号がサンプリングされる（例えば、投影データ信号がマトリックスＰＤＳを形成するためにサンプリングされる）ものとして説明されたが、当業者においては理解できるように、リング抑止フィルタ７４は、別の方法とて、信号をサンプリングするのではなく、これらを連続した信号として扱うアナログフィルタとして実現することもできる。図６は、図１に示すＣＴ画像を再生するのに用いたのと同一のデータを用いて再生されたＣＴ画像を示す。ただし、図６においては、投影データ信号は、最初に、リング抑止フィルタ７４を用いて処理された。この画像の生成においては、リング抑止フィルタ７４は、０．００４の閾値ＴＨＲ、および０．００２のオフセット限界ＭＡＸ＿ＥＲＲを用いたが、患者の最大投影振幅は、典型的には、５〜１０のレンジであった。図６に示す画像は、図１に示す画像と比較して、目立つリングは遙かに少なく、鮮明度が大幅に改善されている。上では、リング抑止フィルタ７４は、一回の走査において収集される投影データ信号からセットのオフセット信号を生成し、これらオフセット信号を、その同一走査における投影データ信号を修正するために（つまり、セットのリングを修正された信号を生成するために）用いるものと説明されたが、リング抑止フィルタ７４は、他のモードにて用いることもできる。例えば、一つあるいは複数の前の走査において生成されたセットのオフセット信号を、現在の走査における投影データ信号を修正するために用いることの方が望ましい場合がある。このモードでは、オフセット信号をＣＴシステムがアイドルのときに“フオライン”にて生成することができ、この場合は、リング抑止フィルタ７４の速度は、それほど重要でなくなる。さらに別のモードとして、リング抑止フィルタ７４は、オフセット信号を生成するために、一つあるいは複数の前の走査からの投影データ信号を用いることもできる。これは、ある走査においてオフセット信号を計算するための良好なデータが不足するような場合（例えば、エラー信号の異常に大きな割合が閾値以上となるような場合）に、他の走査の結果をこの不足を補うために用いることができ有益である。さらに、複数の走査の結果を平均化することもできるが、この場合は、遷移現象などに対して強いより正確なオフセット信号を生成することが可能となる。一つの好ましい実施例においては、ヒストグラム生成器３１８は、前の走査からの前に生成されたヒストグラム信号の加重平均である信号を生成し、リング抑止フィルタ７４は、オフセット信号を生成するためにこの加重平均を用いる。このような実施例においては、図３のヒストグラム生成器３１８は、ディバイダ３３６と、低域通過フィルタ３３４の出力の所に、それぞれ、信号 H_AVE₁（i、j）と、H_AVE₂（i、j）を提供するが、これら２つの信号は、それぞれ、ヒストグラム信号 HIST（i、j）と、低域通過フィルタリングされた信号 HISTLF₂（i、j）の加重平均である。最初、信号 H_AVE₁（i、j）と、H_AVE₂（i、j）は、ｉとｊの全ての値に対して、０にセットされ、各走査の後に、ヒストグラム生成器３１８は、セットの重みＷ_old（i、j）と、Ｗ_new（i、j）を以下のセットの式（１１）、（１２）、（１３）に従って更新する：ＳＵＭ(i,j)＝Ｈ_ＡＶＥ₂(i,j)＋HISTLF₂(i,j) (１３) 重みＷ_old（i、j）と、Ｗ_new（i、j）がいったん計算されると、信号 H_AVE₁（i、 j）と、H_AVE₂（i、j）が以下の式（１４）、（１５）に従って更新される：Ｈ_ＡＶＥ₁(i,j)＝Ｈ_ＡＶＥ₁(i,j)*Ｗ_old(i,j)+HIST(i,j)*Ｗ_new(i,j)(１４) ここで、Ｔは、一定の閾値を表す。式（１４）、（１５）は、信号 H_AVE₂(i、j) がＴより大きくなるのを阻止する。当業者においては理解できるように、これは、前の走査からの古いヒストグラム信号の信号 H_AVE₁（i、j）への寄与を制限する効果を持つ。ただし、古いヒストグラム信号の信号 H_AVE₁（i、j）への寄与を制限するための他の機能を用いることも可能である。各走査の後に、信号 HIST（i、j）が、式（８）と同一の以下の式（１６）に従って生成され、その後、H_AVE₁（i、j）信号が、その走査の際に生成された信号を用いて生成される：当業者においては理解できるように、信号 H_AVE₁（i、j）は、過去のヒストグラムの加重平均を含む。ヒストグラム生成器３１８によって、いったん信号H_AVE₁ （i、j）が生成されると、次に、オフセット信号が、H_AVE₁（i、j）信号をリミッタ３２０に加えることによって生成され、次に、リングを修正された信号が式（１０）に従って生成される。リング抑止フィルタ７４が、複数の走査からの過去の投影データ信号を用いて、リングを修正された信号を生成する場合も、リング抑止フィルタ７４は、複数の動作モードにて用いることができる。一つの動作モードにおいては、リング抑止フィルタ７４は、過去のリングエラー情報を保持し、オフセット信号を、この過去のエラー情報を用いて生成する。このモードの一例においては、リング抑止フィルタ７４は、H_AVE₁（i、j）信号を生成し、この信号を用いてオフセット信号を生成する。図７は、リング抑圧フィルタ７４に対するもう一つの動作モードを示す。図７は、本発明に従って構成されたＣＴシステム７００の信号処理部のブロック図である。システム７００においては、投影フィルタ７２は、応答校正テーブル７１０とリングオフセットテーブル７１２を含む。上述のように、応答校正テーブルは、検出器の応答、より厳密には検出器の伝達関数の変動を補償するために用いられ、通常は、稀にしか更新されず、周知の方法にて“幻像（phantoms）”を走査することによって生成される。システム７００においては、応答校正テーブル７１０によって生成された出力信号は、リングオフセットテーブル７１２に加えられ、リングオフセットテーブル７１２は、これから投影データ信号を生成し、投影データ信号はリング抑圧フィルタ７４に加えられる。リング抑圧フィルタ７４は、リングを修正された信号を生成し、リングを修正された信号は逆投影器７６に加えられる。リング抑圧フィルタ７４は、また、リングオフセットテーブル７１２をフィードバック経路を介して更新する。リングオフセットテーブル７１２と応答校正テーブル７１０は、両方とも、メモリ、例えば、検索テーブル（ＬＵＴ）を含むように実現することも、あるいは、当業者においては理解できるように、リングオフセットテーブル７１２と応答校正テーブル７１０は互いに結合して単一のＬＵＴとして実現することもできる。最初は、応答校正テーブル７１０の更新の後に、リングオフセットテーブル７１２は、好ましくは、投影データ信号を、応答校正テーブル７１０によって生成された出力信号と等しくなるように生成する。ただし、各走査の後には、リング抑圧フィルタ７４は、リングオフセットテーブル７１２の内容を更新し、投影データ信号は、前の走査から得られた古いリングエラー情報を用いて生成される。リングオフセットテーブル７１２は、リングエラーのセットの推定値を格納するが、各推定値は、対応するチャネルの応答の応答校正テーブル７１０が最後の更新された後の校正値からの偏差（ずれ）の測定値を表す。こうして、システム７００においては、投影フィルタ７２は、前の走査の際にリング抑圧フィルタ７４によって観測された検出器応答のエラーを補償し、リング抑圧フィルタ７４は、現在の走査に際に観測された検出器応答のエラーを補償する。システム７００を実現する一つの方法においては、リング抑圧フィルタ７４は、低域通過フィルタリングされたヒストグラム信号 HISTLF₁、HISTLF₂を投影フィルタ７２に加え、投影フィルタフィルタ７２は、これから H_AVE₁（i、j）信号を生成し、これら信号をリングオフセットテーブル７１２内に格納する。リングオフセットテーブル７１２も、H_AVE₁（i、j）信号の値をある閾値、例えば、ＭＡＸ＿ＥＲＲに制限する図３に示すリミッタ３２０に類似するリミッタを含む。投影フィルタ７２は、次に、投影データ信号を以下の式（１７）に従って生成する：ＰＤＳ(i,θ)＝ＰＤＳ_pre(i,θ)−Ｈ_ＡＶＥ₁(i,ＱＵＡＮＴ[ＰＤＳ(i,θ)]) （１７）ここで、ＰＤＳ_pre（i、θ）は、リングオフセットテーブル７１２によって修正される前の投影データ信号の測定値を表す（つまり、ＰＤＳ_pre（i、θ）は応答校正テーブル７１０によって生成される測定値を表す）。こうして、システム７００においては、投影フィルタ７２によって生成される投影データ信号は、リング抑圧フィルタ７４によって前の走査の際に観測されたリングエラーを補償すように修正される（つまり、リングオフセットテーブル７１２は、応答校正テーブル７１０の最後の更新と、最も最近の走査との間に観測されたチャネル応答の偏差を補償する）。リング抑圧フィルタ７４は、次に、投影データ信号を、式（１０）に従ってリングを修正された信号を生成することで、現在の走査の際に観測されたリングエラーについて補償する。リング抑圧フィルタ７４は、リングエラーの新たなセットの推定値も生成するが、これらの各推定値は、１つのチャネルの応答の現在の走査と最も最近の前の走査との間で発生した偏差の測定値を表す。投影フィルタ７２は、次に、リング抑圧フィルタ７４からのこの新たな推定値を用いて、リングオフセットテーブル７１２内に格納されている推定値を更新する。複数の走査の際に収集された投影データ信号の過去の情報をリングを修正された信号を生成するために用いる任意の動作モードにおいては、ＣＴシステムは、リング抑圧フィルタの入力に接続された信号処理チャネルを形成する検出器アセンブリ４４を含む複数の要素の伝達関数を自動的に校正するデバイスであるとみなすこともできる。ＣＴシステムは、勿論、この場合でも、ＣＴ画像の生成を行なうが、ただし、これに加えて、ＣＴシステムは、患者の走査の際に収集されたデータを用いて、これら検出器アセンブリの校正を維持する。上では、リング抑圧フィルタ７４が、ＣＴシステムと共に、投影フィルタ７２によって供給される投影データ信号からリングを修正された信号を生成するために用いられるものとしての説明された。ただし、本発明に従って構成されるＣＴシステムは、投影データ信号に、ＣＴ画像の再生の前に、他の多くの異なるタイプのフィルタリングを適用することもできる。図８は、本発明に従って構成された一つの好ましいＣＴシステム９００の信号処理部をブロック図にて示す。システム９００は、リング抑圧フィルタ７４に加えて、平行光線コンバータ９１０、モーション補償（抑止）フィルタ９１２、しま抑圧（補償）フィルタ９１４、および挿間フィルタ９２０を含む。平行光線コンバータ９１０、モーション補償フィルタ９１２、およびしま抑圧フィルタ９１４に関しては、上述の米国特許出願第０８/５８７,４６８号“STREAK SUPPRESSION FILTER FOR USE IN COMPUTED TO MOGRAPHY SYSTEM”および米国特許出願第０８/６１４,６２３号“MOTION ARTIFA CT SUPPRESSION FILTER FOR USE IN COMPUTED TOMOGRAPHY SYSTEM”において詳細に説明されているため、これらを参照されたい。これらは、ＣＴ画像から、他のタイプのアーティファクトを抑圧するために用いられる。システム９００の一つの好ましい実施例においては、投影フィルタ７２によって生成された投影データ信号は、平行光線コンバータ９１０に加えられ、平行光線コンバータ９１０は、これから、セットの平行な光線信号を生成する。この平行な光線信号は、リング抑圧フィルタ７４に加えられ、リング抑圧フィルタ７４はリングを修正された信号を生成する。このリングを修正された信号は、モーション補償フィルタ９１２に加えられ、モーション補償フィルタ９１２は、セットのモーションを補償された信号を生成する。このモーションを補償された信号は、次に、しま抑圧フィルタ９１４に加えられ、しま抑圧フィルタ９１４は、これからセットのしまを抑圧された信号を生成する。このしまを抑圧された信号は、次に、挿間フィルタ９２０に加えられ、挿間フィルタ９２０の出力は、逆投影器７６に加えられ、逆投影器７６は、この出力信号から、ＣＴ画像を再生する。システム９００の他の幾つかの実施例においては、並列光線コンバータ、リング抑圧フィルタ、モーション補償フィルタ、しま抑圧フィルタのどれか１つ、２つ、あるいは３つが除去される。ただし、システム９００の好ましい実施例は、これら４つのフィルタ９１０、７４、９１２、９１４の全てを含む。他の幾つかの実施例においては、リング抑圧フィルタ７４は、図７との関連で上で説明したように、フィードバック経路を介して、投影フィルタ７２に結合される。さらに別の幾つかの実施例においては、リング抑圧フィルタ７４および／あるいはしま抑圧フィルタ９１４は、平行光線コンバータ９１０と逆投影器７６との間ではなく、投影フィルタ７２と平行光線コンバータ９１０との間に結合される。平行光線コンバータ９１０は、リオーダコンバータ９１６とインターリーブコンバータ９１８を含む。リオーダコンバータ９１６は、投影フィルタ７２から投影データ信号を受信し、これからセットのリオーダされた（再編成された）信号を生成する。こうしてリオーダされた信号は、インターリーブコンバータ９１８に加えられ、インターリーブコンバータ９１８は、これから並列な光線信号を生成する。投影フィルタ７２によって生成された投影は、全ての投影が図２に示すような扇形の光線５２を用いて生成されるために、“扇形の光線（fan beam）” である。平行光線コンバータ９１０は、これら扇形の投影を再編成して、平行な光線の投影を形成する。図９Ａは、患者５０の断面の単一の投影ビューの一部分を形成するセットの光線１０１０を示す。各光線は、Ｘ−線源４２から放射されるが、これは、本質的にて点源であるために、光線１０１０のどれも平行ではなく、こうして得られた投影は、扇形光線の投影となり、ＰＤＳマトリックスの各行は、単一の扇形光線の投影に対応する。リオーダコンバータ９１６は、この扇形の投影データ信号を再編成するが、こうして、再編成された投影は、セットの平行な光線、例えば、図９Ｂに示す光線１０２０から形成される。図１０Ａ〜１０Ｂは、リオーダコンバータ９１６によって再編成された投影を生成するために用いるための一つの好ましい方法を示す。図１０Ａ〜１０は、２つの連続する投影ビューを生成する際のＸ−線源42と検出器アレイ44の位置を示す。ある一つの走査の際に、Ｘ−線源４２と検出器アレイ４４は、円１１２０の中心１１１０のまわりを反時計まわりに回転する。図１０Ａに示す第一の投影の際は、光線１１３０は、検出器４４：４（つまり、アレイ４４の第四のチャネルの部分を形成する検出器）に当たる。そして、図１０Ｂに示す次の投影の際は、光線１１３２は、検出器４４：３（つまり、アレイ４４の第三のチャネルの部分を形成する検出器）に当たる。この好ましい実施例においては、検出器間の間隔は、一連の投影を生成する間の回転の量と一致し、光線１１３０は、光線１１３２と平行であるが、少しオフセットする。この好ましい実施例においては、この基本的な関係は、全ての検出器に対して適用する。こうして、一連の投影の際に隣接する検出器に当たる任意の２つの光線は、平行で、互いに少しオフセットする。上述のように、この好ましい実施例においては、Δθは、０．１２５度とされ、このため、この好ましい実施例においては、アレイ４４内の各検出器は、隣接する検出器と、０．１２５度だけ離れる。リオーダコンバータ９１６は、この基本的な関係を用いてデータを再編成し、再編成された投影を生成する。リオーダコンバータ９１６は、好ましくは、ＰＤＳマトリックスを再編成することで、ＲＥマトリックスの各行が平行な光線によって形成される投影と等価となるようなリオーダされた信号のマトリックスＲＥを形成する。リオーダコンバータ９１６は、好ましくは、ＲＥマトリックスを、ＲＥマトリックスの各要素ＲＥ(i、θ)が、以下の式（１８）に従って選択されるように生成する：ＲＥ(i,θ)＝ＰＤＳ(i,[i-k][Δθ]＋θ) （１８）ここで、ｋ番目のチャネルは、検出器アレイの幾何中心に最も近いチャネルを表し、ＲＥマトリックスの各要素ＲＥ(ｉ、θ)は、ｉ番目のチャネルの平行な光線の投影角度θに対するリオーダされた信号の測定値を表す。リオーダコンバータ９１６は、さらに、各チャネルに対して、隣接する投影角度の投影を平均化するための低域通過フィルタを含む。平均化された、すなわちデシメートされた平行な光線のマトリックスＲＥ(ｉ、θ)は、より大きな角度間隔Δθだけ離された、より少数の行を含む。このようにＲＥマトリックスをデシメートすることは、これによってその後の動作に対する計算が低減されるために有益である。（図８に示す）インターリーブコンバータ９１８は、こうしてリオーダされた信号を受信し、これから平行な光線の信号を生成する。インターリーブコンバータ９１８は、好ましくは、１８０度の間隔だけ離れたペアの平行な光線の投影を結合し、より高密度の投影を形成する。図１１Ａと図１１Ｂは、Ｘ−線源４２、患者の断面５０、および検出器アレイ４４の間の関係を、それぞれ、０度と１８０度の投影角度に対して示す。図１１Ａ〜１１Ｂにおいては、検出器アレイ４４は７個の検出器を含み、第四のチャネルの部分を形成する検出器４４：４がアレイ４４の中央の検出器であるように示される。上述のように、この好ましい実施例においては、実際には、検出器アレイ４４は、３８４個の検出器を含むが、以下では、説明を簡単にするために、７個の検出器を含む実施例について説明する。この好ましい実施例においては、検出器アレイ４４は、ディスク４６の中央１２１０から少しオフセットされ、Ｘ−線源４２の焦点とディスクの中央１２１０とを交差する線１２３４は、中央の検出器４４：４の中央を交差しない。このような検出器システムの構成が１９９４年２月３日付けで出願され(弁理士ドケット番号ANA-044)、本発明の譲受人に譲渡された米国特許出願第０８/１９１,４２８号“X-RAY TOMOGRAPHY SYSTEM FOR AND METHOD OF IMPROVING THE QUALITY OF A SCANNED IMAGE”において開示されているために、これを参照されたい。図１２は、０度および１８０度の投影角度における、検出器アレイ４４と３つの検出器に当たる光線１３１０、１３１２、１３１４との間の空間的な関係を示す。検出器アレイ４４とディスク４６の中央１２１０との間のオフセットのために、０度の投影角度における検出器アレイ４４は、１８０度における検出器アレイ４４から少しオフセットする。このため、１８０度の投影角度において６番目のチャネルの検出器４４：６に当たる光線１３１０は、それぞれ、０度の投影角度において検出器４４：２と４４：３に当たる光線１３１２と１３１４の間に来る。この例では、検出器４４：６は、“中央（central）”検出器であると呼ばれ、検出器４４：２および４４：３は、“反対の隣接する（opposite-adjacent ）”検出器であると呼ばれる。各投影角度において、各検出器は、患者の一部分の密度を測定するが、一般的には、反対の隣接する検出器によって測定される部分は、任意の他の検出器によって測定される部分より、中央検出器によって測定される部分に近い（例えば、０度の投影角度において検出器４４：２、４４：３によって測定される部分は、１８０度の投影角度において検出器４４：５、４４：７によって測定される部分より、１８０度の投影角度において検出器４４：６によって測定される部分に近い）。１８０度だけ離れた任意の投影が中央検出器と反対の隣接する検出器との間のこの関係を用いてインターリーブされ、これによって単一のより高密度の投影が形成される。例えば、図１２に示される構成では、一つのこのようなインターリーブされた投影は、量［ＲＥ(1、0)、ＲＥ(7、18 0)、ＲＥ(2、0)、ＲＥ(6、180)、ＲＥ(3、0)、ＲＥ(5、180)、ＲＥ(4、0)、ＲＥ(4、18 0)、ＲＥ(5、0)、ＲＥ(3、180)、ＲＥ(6、0)、ＲＥ(2、180)、ＲＥ(7、0)、ＲＥ(1、18 0)］から構成され、ここで、ＲＥ(ｉ、θ)は、ｉ番目のチャネルの検出器によって、投影角度θに対して、生成されるリオーダされた信号を表す。インターリーブコンバータ９１８は、リオーダされた信号をこのようにインターリーブすることで、より高密度の投影を形成する。インターリーブコンバータ９１８は、好ましくは、平行な光線の信号の測定値のマトリックスＰＡＲを生成するが、ここで、ＰＡＲマトリックスの各要素ＰＡＲ(ｉ、θ)は、ｉ番目のチャネルのθに等しい平行な光線の投影角度に対する測定値を表す。ＰＡＲマトリックスの構造は、以下の式（１９）によって表される：式（１９）に示されるように、ＰＡＲマトリックスは、ＰＤＳマトリックスの二倍の列と、半分の行を持つ。このため、ＰＡＲマトリックスの各行は、ＰＤＳマトリックスの行の二倍のデータを含む平行な光線の投影を表す。このため、各平行な光線の投影は、扇形の光線の投影のチャネルの二倍の数のチャネルからの画像データを持つ。ＰＡＲマトリックスは、ＤＰＳマトリックスとは少し異なる巡回的な特性を持ち、ここでは、最後の行は、逆の順番にて最初に行に続き、ＰＡＲ（0、180）＝ＰＡＲ（2n-1、0）；ＰＡＲ（1、180）＝ＰＡＲ（2n-2、0）、となる。この好ましい実施例においては、インターリーブコンバータ９１８は、ＰＡＲマトリックスの要素を以下の式（２０）に従って生成する：ＰＡＲ(2i,θ)＝ＲＥ(i,θ) ＰＡＲ(2i+1,θ)＝ＲＥ(Ｎ-1-i,θ+180) （２０） for0≦ｉ＜Ｎ周知のように、扇形の光線のデータを平行な光線のデータに変換するための平行光線コンバータ、例えば、図８のコンバータ９１０は、通常は、挿間フィルタ、例えば、（図８の）フィルタ９２０を含む。ただし、挿間フィルタは、通常は、インターリーブコンバータ９１８の直後に配置される。検出器は、通常は、隣接する検出器間のＸ−線源に対する角度オフセットが全ての検出器に対して同一となるような間隔に配置されるが、これら検出器は線形的な意味では等距離の配置ではない。従って、ＰＡＲマトリックスの各行（つまり、各平行な光線の投影）は等間隔ではないデータ点を含み、むしろ、各投影の中央付近、つまり、等角点を通る中央光線付近の要素は各投影の両端付近の要素より広い間隔を持つ。挿間フィルタはデータを挿間することで、各投影の全ての要素が等間隔のとなるような平行光線データの新たなマトリックスを生成する。この好ましい実施例においては、挿間フィルタ９２０は、データを挿間することで等間隔の要素を含む投影を生成するために周知の技術を用いる。好ましくは、挿間フィルタ９２０は、インターリーブコンバータ９１８の直後に配置するのではなく、図８に示すように、しま補償フィルタ９１４の後に配置される。ただし、本発明は、従来の方法において通常行なわれるように、挿間フィルタ９２０をインターリーブコンバータ９１８の直後に配置しても問題なく機能する。さらに、周知のように、扇形の光線のデータを、平行な光線のデータに変換した場合、通常は、少しの回転が導入され、このため、０度の平行な光線の投影角度は、０度の扇形の光線の投影角度とは完全には一致しなくなり、修正しない場合は、この回転のために、水平から少し回転した再生画像が生成される。この回転は、一般的には、式（８）において用いられる中央の検出器“ｋ”が光源から等角点を経て中央の検出器ｋ’へと通る線の中央に完全にはこないために導入される。ただし、回転の量は、一般的には、Δθ／２より小さく、挿間フィルタ９２０あるいは逆投影器７６によって周知の技術を用いて修正することも、あるいは、単に無視することもできる。こうして、図８に示す好ましい実施例においては、リング抑圧フィルタ７４は、平行光線コンバータ９１０によって生成される平行な光線のデータに関して動作する。ただし、このデータには、モーションによるアーティファクトが存在する可能性があり、リング抑圧フィルタ７４は、実際には、インターリーブされてない信号の関数として（つまり、リオーダコンバータ９１６によって生成されるリオーダされた信号を用いて）平行な光線のヒストグラムおよびオフセット信号を生成する方が好ましい場合がある。図１３は本発明に従って構成されたもう一つの好ましいＣＴシステム１３００の信号処理部をブロック図にて示す。システム１３００は、（図８に示す）システム９００と類似するが、システム１３００では、リング抑圧フィルタ７４が、インターリーブコンバータ９１８とモーション抑圧フィルタ９１２との間ではなく、リオーダコンバータ９１６とインターリーブコンバータ９１８との間に結合される点が異なる。この実施例においてもリング抑圧フィルタ７４は、実質的に上述と同様に動作する。ただし、この実施例では、リング抑圧フィルタ７４は投影フィルタ７２によって生成された投影データ信号ではなく、リオーダフィルタ９１６によって生成されたリオーダされた信号に関して動作する点が異なる。上述の装置に対して幾らかの変更を本発明の範囲から逸脱することなく加えることが可能ある。このため上述の説明あるいは付録の図面に含まれる全ての事項は解説を意図し、制限を意図するものではないと解釈されるべきである。

Claims

【特許請求の範囲】１．患者のコンピュータ断層画像を生成するためのコンピュータ断層撮影システムであって、このシステムが、（ａ）Ｘ−線源、および（ｂ）対応する伝達関数を持つ複数の信号処理チャネルを形成（定義）するチャネル形成（定義）手段を含み、この複数の信号処理チャネルが、前記Ｘ−線源から放出されるＸ−線を検出するための検出器アレイの対応する複数の検出器を含み、これら複数の検出器が、断層走査の複数の各投影角度において対応する検出器によって検出されるＸ−線の関数として、複数の画像データ信号を生成し、これら複数の画像データ信号が、患者の断層画像を表し、このシステムがさらに：患者の一つの完全な走査の際に生成される画像データ信号に応答してエラー修正データを生成する手段；および前記チャネルの伝達関数を、前記画像データ信号から生成されるコンピュータ断層画像内のリングファーティファクトを最小に押さえるために、前記エラー修正データの関数として自動的に校正するための手段を含むことを特徴とするコンピュータ断層撮影システム。２．患者のコンピュータ断層画像を生成するためのコンピュータ断層撮影システムであって、このシステムが、（ａ）Ｘ−線源、および（ｂ）対応する伝達関数を持つ複数の信号処理チャネルを形成するチャネル形成手段を含み、この複数の信号処理チャネルが、前記Ｘ−線源から放出されるＸ−線を検出するための検出器アレイの対応する複数の検出器を含み、これら複数の検出器が、断層走査の複数の各投影角度において対応する検出器によって検出されるＸ−線の関数として、複数の画像データ信号を生成し、これら複数の画像データ信号が、患者の断層画像を表し、このシステムがさらに：前記チャネルの伝達関数を、前記画像データ信号から生成されるコンピュータ断層画像内のリングファーティファクトを最小に押さえるために、患者の走査の際に生成される画像データ信号から得られるエラー修正データの関数として、自動的に校正する手段を含み、前記チャネルの伝達関数を自動的に校正する手段が、前記エラー修正データを、各チャネルを通じて処理される各画像データ信号から、そのチャネルを通じて処理されるデータ信号の振幅の関数として、決定するための手段を含むことを特徴とする請求の範囲１に記載のコンピュータ断層撮影システム。３．前記エラー修正データを各チャネルを通じて処理される各画像データ信号からそのチャネルを通じて処理されるデータ信号の振幅の関数として決定する手段が、前記各チャネルの前記各投影角度における伝達関数の推定エラーを表す推定信号を生成するための手段、前記各推定信号の振幅の値を所定の閾値の値と比較するための手段、および前記推定信号の値が閾値以下である場合に限って、前記各推定信号に対して、エラー修正データを生成するための手段を含むことを特徴とする請求の範囲２に記載のコンピュータ断層撮影システム。４．前記チャネルの伝達関数を自動的に校正する手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、前記対応する検出器に隣接する前記検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅との関数として決定するための手段を含むことを特徴とする請求の範囲１に記載のコンピュータ断層撮影システム。５．前記チャネルの伝達関数を自動的に校正する手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、前記対応する検出器に隣接する前記検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅の関数として、前記投影角度の選択された一つに対して、決定するための手段を含むことを特徴とする請求の範囲１に記載のコンピュータ断層撮影システム。６．前記チャネルの伝達関数を自動的に校正する手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅の関数として、所定の数の各投影角度に対して、決定するための手段を含むことを特徴とする請求の範囲１に記載のコンピュータ断層撮影システム。７．前記チャネルの伝達関数を自動的に校正する手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、（ａ）前記各投影角度において対応する検出器によって提供され、そのチャネルを通じて処理されるデータ信号の振幅と、（ｂ）所定の選択された数の隣接する投影角度におけるその検出器のデータ信号の振幅の関数として決定するための手段を含むことを特徴とする請求の範囲１に記載のコンピュータ断層撮影システム。８．前記チャネルの伝達関数を自動的に校正する手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、前記対応する検出器に隣接する前記検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅との関数として決定するための手段を含むことを特徴とする請求の範囲７に記載のコンピュータ断層撮影システム。９．対象を走査するためのコンピュータ断層撮影システムであって、このシステムが：（ａ）複数の信号処理手段を含み、これら信号処理手段が、複数の各投影データ信号を生成するための対応する複数の検出器を含み、各信号処理手段が、複数の伝達関数の対応する一つによって特性化され、各投影データが、データ成分とエラー成分を含み、エラー成分は、信号処理手段の伝達関数間の相対的な差を表し、データ成分は対象の部分の密度を表し；このシステムがさらに（ｂ）フィルタを含み、このフィルタがエラー成分の推定値を生成するための手段を含み、これら各推定値が、投影データ信号の対応する一つの振幅の関数として生成されることを特徴とするコンピュータ断層撮影システム。１０．対象のコンピュータ断層画像を生成するコンピュータ断層撮影システムであって、このシステムが、複数の投影を生成するための手段を含み、各投影が対応する投影角度にて生成され、各投影が、複数の各投影データ信号の測定値を含み、投影データ信号の各測定値が、対象の部分の密度を表し、このシステムがさらに、断層画像内のリングファーティファクトを抑圧するためのリング抑圧フィルタを含み、このフィルタが：（Ａ）複数の測定エラー信号を、投影データ信号の対応する測定値の関数として、各投影に対して、生成するための高域通過フィルタ手段を含み、各測定エラー信号が、投影の一つの少なくとも一部分と関連する高周波成分を表し；このシステムがさらに（Ｂ）投影データ信号の各測定値をある閾値と比較し、閾値より小さな投影データ信号の各測定値を、閾値以下として指定する閾値手段；および（Ｃ）閾値以下の測定値を持つ投影データ信号に対応する測定エラー信号の全ての測定値を平均化するためのヒストグラム手段を含むことを特徴とするコンピュータ断層撮影システム。１１．対象の画像を生成するコンピュータ断層撮影システム内で用いるためのリング抑圧フィルタであって、このシステムが、複数の投影（データ信号）を生成するための手段を含み、各投影が、対応する投影角度にて生成され、各投影が、複数の各投影データ信号の一つの測定値を含み、投影データ信号の各測定値が、対象の部分の密度を表し、このリング抑圧フィルタが：（Ａ）複数のエラー信号の測定値を生成するための高域通過フィルタ手段を含み、このエラー信号の各測定値が、投影の一つの一部分と関連する高周波成分を表し；このシステムがさらに（Ｂ）投影データ信号の一つの測定値と、エラー信号の測定値とに応答して、ヒストグラム信号を生成するヒストグラム手段を含み、このヒストグラム信号が、エラー信号の選択された測定値と、各対応する投影データ信号との間の関係を表し；このシステムがさらに：（Ｃ）対応する投影データ信号の測定値と、ヒストグラム信号とに応答して、これから、リングを修正された信号の複数の測定値を生成するための結合手段を含むことを特徴とするコンピュータ断層撮影システム。１２．前記ヒストグラム手段が、一つの投影データ信号の各測定値とそのエラー信号の対応する測定値とに同一のビン番号が割り当てられるように、複数のビン番号の一つを、一つの投影データ信号の各測定値とそのエラー信号とに割り当てるための量子化手段を含むことを特徴とする請求の範囲１１に記載のリング抑圧フィルタ。１３．前記量子化手段が、同一の振幅を持つ一つの投影データ信号の測定値に、同一のビン番号を割り当てることを特徴とする請求の範囲１２に記載のリング抑圧フィルタ。１４．前記ヒストグラム手段が、ヒストグラム信号を、ヒストグラム信号が、ビン番号の各値を持つように、生成することを特徴とする請求の範囲１３に記載のリング抑圧フィルタ。１５．前記ヒストグラム手段が、１つのビン番号に割り当てられたエラー信号の全ての測定値を平均化することで、その１つのビン番号に対するヒストグラム信号の値を生成することを特徴とする請求の範囲１４に記載のリング抑圧フィルタ。１６．前記ヒストグラム手段が、１つのビン番号に割り当てられ、かつ、ある閾値より小さな規模を持つエラー信号の全ての測定値を平均化することで、その１つのビン番号に対するヒストグラム信号の値を生成することを特徴とする請求の範囲１４に記載のリング抑圧フィルタ。１７．前記ヒストグラム手段が、エラー信号の各測定値をある閾値と比較するための手段、および、閾値より小さな規模を持つエラー信号の各測定値を閾値以下の測定値として指定するための手段を含むことを特徴とする請求の範囲１４に記載のリング抑圧フィルタ。１８．前記ヒストグラム手段が、１つのビン番号に割り当てられた全ての閾値以下の測定値を平均化することで、その１つのビン番号に対するヒストグラム信号の値を生成することを特徴とする請求の範囲１７に記載のリング抑圧フィルタ。１９．前記結合手段が、リングを修正された信号の各測定値を、どのビン番号が一つの投影データ信号の一つの測定値に割り当てられているかを決定し、次に、その一つの投影データ信号のその１つの測定値からそのビン番号に対するヒストグラム信号の値を減ずることによって生成することを特徴とする請求の範囲１８に記載のリング抑圧フィルタ。２０．前記ヒストグラム手段が、第一の中間信号と第二の中間信号を生成するための手段を含み、前記第一の中間信号が、各ビン番号に対して、そのビン番号に割り当てられたエラー信号の全ての閾値以下の測定値の総和と実質的に等しい値を持ち、前記第二の中間信号が、各ビン番号に対して、実質的にそのビン番号に割り当てられたエラー信号の全ての閾値以下の測定値の総数に等しい値を持つことを特徴とする請求の範囲１７に記載のリング抑圧フィルタ。２１．前記ヒストグラム手段が、前記第一の中間信号を受信し、これを低域通過フィルタリングすることで、これから第一のフィルタリングされた信号を生成するための第一の低域通過フィルタ手段、および前記第二の中間信号を受信し、これを低域通過フィルタリングすることで、これから第二のフィルタリングされた信号を生成するための第二の低域通過フィルタ手段を含むことを特徴とする請求の範囲１７に記載のリング抑圧フィルタ。２２．前記ヒストグラム手段が、前記第一の中間信号を前記第二の中間信号にて割ることで、これからヒストグラム信号を生成するためのディバイダ手段を含むことを特徴とする請求の範囲２１に記載のリング抑圧フィルタ。２３．さらに、ヒストグラム信号を受信し、これからオフセット信号を生成するためのオフセット手段を含み、このオフセット信号が、ヒストグラム信号の以下の規模を持つことを特徴とする請求の範囲２２に記載のリング抑圧フィルタ。２４．さらに、各投影を低域通過フィルタリングし、低域通過フィルタリングされた投影を高域通過フィルタ手段に加える低域通過フィルタ手段を含むことを特徴とする請求の範囲１１に記載のリング抑圧フィルタ。２５．前記各投影が扇形の光線の投影であることを特徴とする請求の範囲１１に記載のリング抑圧フィルタ。２６．前記各投影が平行な光線の投影であることを特徴とする請求の範囲１１に記載のリング抑圧フィルタ。２７．前記ヒストグラム手段が、投影データ信号の測定値および２回あるいはそれ以上の走査の際に生成されたエラー信号に応答することを特徴とする請求の範囲１１に記載のリング抑圧フィルタ。２８．走査対象の断層画像を生成するためのコンピュータ断層撮影システムであって、このシステムが、Ｘ−線源と、対応する複数の伝達関数にて特性化される複数の信号処理チャネルを形成するチャネル形成手段を含み、前記複数の信号処理チャネルが、前記Ｘ−線源から放出されるＸ−線を検出するための検出器アレイの対応する複数の検出器を含む手段を含み、これら複数の検出器が、複数のデータ信号を、対応する検出器によって、断層走査の複数の各投影角度に対して、対応する信号処理チャネルを通じて、検出されるＸ−線の関数として生成し、このシステムがさらに：前記チャネルの伝達関数を、患者の少なくとも一回の前の走査の際に生成されたデータ信号から得られたエラー修正データの関数として、自動的に校正するための手段を含み、これによって前記データ信号から生成される断層画像内のリングファーティファクトが最小に押えられることを特徴とするコンピュータ断層撮影システム。２９．前記チャネルの伝達関数を自動的に校正するための手段が、前記エラー修正データを、患者の少なくとも一回の前の走査の際に生成されたデータ信号と、患者の現在の走査の際に生成されたデータ信号の両方から、決定するための手段を含み、これによって、前記データ信号から生成される断層画像内のリングファーティファクトが最小に押さえられることを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３０．前記チャネルの伝達関数を自動的に校正するための手段が、前記エラー修正データを、一人あるいは複数の患者の複数の前の走査の際に生成されたデータ信号から、決定するための手段を含み、これによって、前記データ信号から生成される断層画像内のリングファーティファクトが最小に押さえられることを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３１．前記チャネルの伝達関数を自動的に校正するための手段が、前記エラー修正データを、各チャネルを通じて処理された各画像データ信号から、そのチャネルを通じて処理されたデータ信号の振幅の関数として、決定するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３２．前記エラー修正データを各チャネルを通じて処理された各画像データ信号からそのチャネルを通じて処理されたデータ信号の振幅の関数として決定するための手段が、各チャネルの各投影角度における伝達関数の推定エラーを表す推定信号を生成するための手段、各推定信号の振幅の値を所定の閾値の値と比較するための手段、および各推定信号に対して、推定信号の値が閾値以下である場合に限って、エラー修正データを生成するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３３．前記チャネルの伝達関数を自動的に校正するための手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、前記対応する検出器に隣接する前記検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅との関数として、決定するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３４．前記チャネルの伝達関数を自動的に校正するための手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、前記対応する検出器に隣接する検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅との関数として、前記投影角度の選択された一つに対して、決定するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３５．前記チャネルの伝達関数を自動的に校正するための手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅の関数として、所定の個数の投影角度に対して、決定するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３６．前記チャネルの伝達関数を自動的に校正するための手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、（ａ）前記各投影角度に対して、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、（ｂ）所定の数の隣接する投影角度におけるその検出器のデータ信号の振幅との関数として決定するための手段を含むことを特徴とする請求の範囲２８に記載のコンピュータ断層撮影システム。３７．前記チャネルの伝達関数を自動的に校正するための手段が、エラー修正データを、各チャネルを通じて処理された前記各画像データ信号から、対応する検出器によって提供され、そのチャネルを通じて処理されたデータ信号の振幅と、対応する検出器に隣接する前記検出器アレイの検出器と関連するチャネルを通じて処理されたデータ信号の振幅との関数として、決定するための手段を含むことを特徴とする請求の範囲３７に記載のコンピュータ断層撮影システム。３８．複数の投影角度を通じて対象を走査するためのコンピュータ断層撮影システムであって、このシステムが複数の信号処理手段を含み、これら各信号処理手段が検出器アレイの複数の検出器を含み、前記各処理手段が複数の投影データ信号の対応する一つを生成し、前記各処理手段が複数の信号伝達関数の対応する一つによって特性化され、前記各投影データ信号がデータ成分とエラー成分を含み、エラー成分は、前記処理手段の信号伝達関数間の相対的な差を表し、前記複数の処理手段が、各走査に対して、おのおのの投影データ信号を生成し、これらデータ成分が、走査される対象の部分の密度を表わし、このシステムがさらに：（Ａ）少なくとも一回の前の走査において、前記各検出器によって、前記各投影角度に対して収集されたエラー成分の前のセットの推定値を格納するための手段を含み、この格納手段が、現在の走査からのエラー成分の新たなセットの推定値を受信して、前のセットの推定値と新たなセットの推定値の関数として、更新されたセットの推定値を生成するための手段を含み、このシステムがさらに（Ｂ）投影データ信号を受信し、これから、エラー成分の新たなセットの推定値を生成し、この新たなセットの推定値を、前記格納手段に供給するリング抑圧フィルタ手段を含むことを特徴とするコンピュータ断層撮影システム。３９．対象を走査するためのコンピュータ断層撮影システムであって、このシステムが、複数の信号処理手段を含み、これら信号処理手段が複数の投影データ信号の対応する一つを生成し、前記各処理手段が複数の信号伝達関数の対応する一つによって特性化され、前記各投影データ信号がデータ成分とエラー成分を含み、エラー成分は、前記複数の処理手段の信号伝達関数間の相対的な差を表し、前記複数の処理手段が、各走査に対して、おのおのの投影データ信号を生成し、データ成分は、走査される対象の部分の密度を表し、このシステムがさらに：（Ａ）少なくとも一回の前の走査の際に得られたエラー成分の前のセットの推定値を格納するための手段を含み、この格納手段が、現在の走査からのエラー成分の新たなセットの推定値を受信し、古いセットの推定値と新たなセットの推定値の関数として、セットの推定値を更新するための手段を含み、このシステムがさらに（Ｂ）ある走査の投影データ信号と古いセットのエラー推定値とを受信し、前記投影データ信号と古いセットのエラー推定値との関数として、複数の修正された投影データ信号を生成するための結合手段；および（Ｃ）前記修正された投影データ信号を受信し、これから、エラー成分の新たなセットの推定値を生成し、この新たなセットの推定値を、前記格納手段に供給するリング抑圧フィルタ手段を含むことを特徴とするコンピュータ断層撮影システム。