(57)【特許請求の範囲】
【請求項１】
Ｘ線を検出する検出素子を二次元配列した検出素子群を複数有し、該検出素子群を１画素に対応させて複数配列した検出部と、
前記検出素子の出力信号の加算率を決定する加算率決定部と、
前記検出素子群に属する前記検出素子の出力信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出する加算部と、
前記画素と当該画素に対応する検出素子群に属する前記検出素子との位置関係を示す画素位置情報と、前記検出素子群に含まれる欠陥素子の位置を示す欠陥素子位置情報とを記憶した位置情報記憶部と、を備え、
前記加算率決定部が、前記画素位置情報及び前記欠陥素子位置情報に基づいて、信号値を算出する画素に含まれる前記欠陥素子の出力信号の加算率と、信号値を算出する画素の中心に対して当該欠陥素子と対称に位置する対角検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定し、他の前記検出素子の加算率を略同一の値となるように決定するＸ線検出器。
【請求項２】
前記位置情報記憶部が、前記検出素子の面積または面積比から成る面積情報を具備し、
前記加算率決定部が、前記面積情報に基づいて、前記欠陥素子及び前記対角検出素子の前記加算率を決定する請求項１記載のＸ線検出器。
【請求項３】
前記加算部が、
前記検出素子の出力信号をデジタル信号に変換するデジタル変換部を備え、
前記検出素子群に属する前記検出素子のデジタル信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出する請求項１記載のＸ線検出器。
【請求項４】
前記加算部が、
前記検出素子からのアナログの出力信号を前記加算率に応じて画素毎に加算してアナログ出力値を算出するアナログ信号加算部と、
前記画素のアナログ出力値をデジタル信号に変換するデジタル変換部と、を備え、
前記加算率決定部が、スイッチであり該スイッチのオンまたオフ操作により、前記加算率を決定する請求項１記載のＸ線検出器。
【請求項５】
前記欠陥素子及び前記対角検出素子の加算率が０である請求項１記載のＸ線検出器。
【請求項６】
前記加算率決定部が、前記検出素子群における前記欠陥素子の配列位置を特定し、
前記検出素子群において前記欠陥素子が中央以外に位置するとき、
前記欠陥素子及び前記対角検出素子の加算率を夫々０に決定する請求項１記載のＸ線検出器。
【請求項７】
前記加算率決定部が、前記検出素子群における前記欠陥素子の配列位置を特定し、
前記検出素子群において前記欠陥素子が角に位置するとき、
前記欠陥素子及び前記対角検出素子の加算率を０に決定する請求項６記載のＸ線検出器。
【請求項８】
前記位置記憶部が、各前記検出素子の前記加算率を加算率情報として予め記憶し、
前記加算率決定部が、前記加算率情報を用いて各前記検出素子の加算率を決定する請求項１記載のＸ線検出器。
【請求項９】
前記投影像の画素毎の信号値に対して、前記検出素子群に属する検出素子の加算率の合計値に基づいて規格化する規格部をさらに備えた請求項１記載のＸ線検出器。
【請求項１０】
前記検出素子が、Ｘ線を検出して、検出したＸ線からＸ線フォトンのエネルギーに応じた信号を発生し、
前記デジタル変換部が、前記エネルギーを２以上のエネルギー範囲に分別してＸ線フォトン数のデジタル信号を出力し、
前記加算部が、前記エネルギー範囲毎に前記画素の信号値を算出する請求項３記載のＸ線検出器。
【請求項１１】
前記欠陥素子と前記対角検出素子の少なくとも一方の前記加算率が、前記エネルギー範囲に応じて異なる請求項１０記載のＸ線検出器。
【請求項１２】
請求項１から請求項１１の何れか１項に記載のＸ線検出器と、
Ｘ線を照射するＸ線発生部と、
前記Ｘ線検出器からの信号に対して再構成演算を行って再構成像を作成する再構成処理部と、
前記Ｘ線検出器、前記Ｘ線発生部、及び前記再構成処理部を制御する制御部と、を備えたＸ線ＣＴ装置。
【請求項１３】
Ｘ線を検出する複数の検出素子が二次元配列された検出素子群を投影像の１画素に対応させて複数配列した検出部により、電荷量に応じて出力された出力信号について、
前記画素と当該画素に対応する検出素子群に属する前記検出素子との位置関係を示す画素位置情報と、前記検出素子群に含まれる欠陥素子の位置を示す欠陥素子位置情報とを記憶するステップと、
前記画素位置情報及び前記欠陥素子位置情報に基づいて、信号値を算出する画素に含まれる前記欠陥素子の出力信号の加算率と、信号値を算出する画素の中心に対して当該欠陥素子と対称に位置する対角検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定すると共に、他の前記検出素子の加算率を略同一の値となるように決定するステップと、
前記検出素子群に属する前記検出素子の出力信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出するステップと、を備えたＸ線検出方法。
【請求項１４】
Ｘ線を検出する複数の検出素子が二次元配列された検出素子群を投影像の１画素に対応させて複数配列した検出部により、電荷量に応じて出力された出力信号について、
前記画素と当該画素に対応する検出素子群に属する前記検出素子との位置関係を示す画素位置情報と、前記検出素子群に含まれる欠陥素子の位置を示す欠陥素子位置情報とを記憶するステップと、
前記画素位置情報及び前記欠陥素子位置情報に基づいて、信号値を算出する画素に含まれる前記欠陥素子の出力信号の加算率と、信号値を算出する画素の中心に対して当該欠陥素子と対称に位置する対角検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定すると共に、他の前記検出素子の加算率を略同一の値となるように決定するステップと、
前記検出素子群に属する前記検出素子の出力信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出するステップと、をコンピュータに実行させるＸ線検出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本願発明は、Ｘ線検出器、特に、欠陥素子によって生じる出力信号の欠陥を補間するＸ線検出器に関する。
【背景技術】
【０００２】
複数の方向から撮影した被検体のＸ線透過像である投影像からＸ線吸収係数（線減弱計数）を算出し、被検体の断層像である再構成像を得るＸ線ＣＴ装置が知られている。
このようなＸ線ＣＴ装置に適用される積分型のＸ線検出器は、複数の検出素子を備え、検出素子毎に被検体を透過したＸ線のエネルギーを電気信号に変換し、所定時間積分した積分信号を出力し投影像を得ている。ここで、複数の検出素子には欠陥を有する検出素子（以下、「欠陥素子」という）が含まれることがあり、欠陥素子によって画素のサンプリング位置がずれることで出力値がずれ、投影像においてアーチファクトが生じる虞がある。
そこで、例えば、特許文献１の画像撮影装置では、欠陥素子によって投影像に画素欠陥が生じた場合には、正常な検出素子の出力信号から欠陥素子の出力信号を推定すると共に、推定値に対してあらかじめ定めた影響量パラメータを用いて当該欠陥素子の周辺素子が当該欠陥素子から受ける影響を補正している。
【０００３】
ところで、近年、Ｘ線フォトンの個数を計測するフォトンカウンティング型のＸ線検出器を搭載したＸ線ＣＴ装置が開発されている。Ｘ線ＣＴ装置にフォトンカウンティング型のＸ線検出器を搭載すると、積分型のＸ線検出器を搭載したＸ線ＣＴ装置では取得できなかったエネルギー毎の疑似単色の再構成像や、原子番号などの分布を示す吸収係数以外の再構成像（以下、これらの画像を「マルチエネルギー画像」という）を作成することができるという利点がある。
【０００４】
フォトンカウンティング型の検出器を適用したＸ線ＣＴ装置の一形態では、一画素に対して微細な複数の検出素子を割り当て、検出素子毎にＸ線フォトンの個数を計測し、その出力値を加算することにより投影像の画素毎の出力を求める。フォトンカウンティング型の検出器を適用することにより、Ｘ線ＣＴ装置のように非常に高いＸ線線量率を用いる装置においても、パイルアップを抑えることが可能となる。また、検出素子を小さくすると、同一面積範囲を撮影する際に検出素子数が増えてしまうものの、フォトンカウンティング型の検出器側で、すなわち、投影像でＸ線ＣＴ装置において必要なサイズまで加算しておくことにより、Ｘ線ＣＴ装置側で処理すべきデータ量、処理回路、処理工数、処理時間などの増加を抑えることができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１２−２３１２１０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述のように、フォトンカウンティング型の検出器では、一画素に対して微細な複数の検出素子を割り当てているため、積分型の検出器に比して検出素子数が多く、読出し回路や検出部等の異常等による欠陥素子の割合も高くなる。このため、欠陥素子によって投影像に画素欠陥が生じたときに、特許文献１の画像撮影装置のような画素欠陥を補間する技術をそのまま適用すると、より多くの処理時間、処理回路及び補間用データ等を必要とし、処理速度の低下、装置コストの増大及び作業工数の増加等を招来する。
【０００７】
本発明は上記実情に鑑みてなされたものであり、処理時間、処理回路及び補間用データ等を増大させることなく、欠陥素子に対する補間の精度を向上させ、簡易にアーチファクトを抑制することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明は以下の手段を提供する。
本発明の一態様は、Ｘ線を検出する検出素子を二次元配列した検出素子群を複数有し、該検出素子群を１画素に対応させて複数配列した検出部と、前記検出素子の出力信号の加算率を決定する加算率決定部と、前記検出素子群に属する前記検出素子の出力信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出する加算部と、前記画素と当該画素に対応する検出素子群に属する前記検出素子との位置関係を示す画素位置情報と、前記検出素子群に含まれる欠陥素子の位置を示す欠陥素子位置情報とを記憶した位置情報記憶部と、を備え、前記加算率決定部が、前記画素位置情報及び前記欠陥素子位置情報に基づいて、信号値を算出する画素に含まれる前記欠陥素子の出力信号の加算率と、信号値を算出する画素の中心に対して当該欠陥素子と対称に位置する対角検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定し、他の前記検出素子の加算率を略同一の値となるように決定するＸ線検出器を提供する。
【発明の効果】
【０００９】
本発明によれば、処理時間、処理回路及び補間用データ等を増大させることなく、欠陥素子に対する補間の精度を向上させ、簡易にアーチファクトを抑制することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１の実施形態に係るＸ線検出器を適用したＸ線ＣＴ装置の概略を示すブロック図である。
【図２】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示す参考図である。
【図３】本発明の第１の実施形態に係るＸ線検出器の概略を示すブロック図である。
【図４】本発明の第１の実施形態に係るＸ線検出器におけるサンプリングの様子の一例を示すグラフである。
【図５】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示す参考図である。
【図６】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、（ａ）は、欠陥素子の位置を示す欠陥素子位置情報の例であり、（ｂ）は、欠陥素子位置情報の配列マップの例である。
【図７】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、（ａ）は、画素と検出素子及び欠陥素子との位置関係を示す画素位置情報の例であり、（ｂ）は、画素位置情報の配列マップの例である。
【図８】本発明の第１の実施形態に係るＸ線検出器の加算部における加算処理の流れを示すフローチャートである。
【図９】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、特に、加算率の重心を説明する説明図である。
【図１０】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、特に、検出素子の加算率を説明する説明図である。
【図１１】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の加算率の例を示し、（ａ）は、図９における（−１，１）が欠陥素子であり、加算率を０とした場合の画素内の加算率を示し、（ｂ）は、欠陥素子の（−１，１）とその対角検出素子（１，−１）の加算率を０とした加算率を示し、（ｃ）〜（ｆ）は外挿により加算率を算出した例を示す説明図である。
【図１２】本発明の第１の実施形態に係るＸ線検出器の加算率決定部における加算率決定処理の流れを示すフローチャートである。
【図１３】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、特に、検出素子の加算率を説明する説明図である。
【図１４】本発明の第１の実施形態に係るＸ線検出器の加算率決定部における他の加算率決定処理の流れを示すフローチャートである。
【図１５】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を示し、特に、検出素子の加算率を説明する説明図である。
【図１６】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図であり、（ａ）は画素中に２つの欠陥素子を有する例であり、（ｂ）及び（ｃ）は、画素中に１つの欠陥素子を有する場合の例を示す。
【図１７】図１６のように検出素子が配列された場合の加算率の例を示し、（ａ）は欠陥素子周辺の検出素子の加算率を増加させた例であり、（ｂ）は、欠陥素子の対角検出素子の加算率を０とした例を示す。
【図１８】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図であり、（ａ）は画素中に２つの欠陥素子を有する例であり、（ｂ）及び（ｃ）は、画素中に１つの欠陥素子を有する場合の例を示す。
【図１９】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図である。
【図２０】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図である。
【図２１】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図である。
【図２２】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図である。
【図２３】本発明の第１の実施形態に係るＸ線検出器の検出部における検出素子の配列例を説明する説明図である。
【図２４】本発明の第２の実施形態に係るＸ線検出器の概略を示すブロック図である。
【図２５】本発明の第２の実施形態に係るＸ線検出器の信号処理部の概略を示すブロック図である。
【発明を実施するための形態】
【００１１】
以下、本発明の一実施形態について、図面を参照して説明する。
本発明に係るＸ線検出器は、Ｘ線を検出する検出素子を二次元配列した検出素子群を複数有し、該検出素子群を１画素に対応させて複数配列した検出部と、前記検出素子の出力信号の加算率を決定する加算率決定部と、前記検出素子群に属する前記検出素子の出力信号を前記加算率に応じて加算することにより投影像の画素毎の信号値を算出する加算部と、前記画素と当該画素に対応する検出素子群に属する前記検出素子との位置関係を示す画素位置情報と、前記検出素子群に含まれる欠陥素子の位置を示す欠陥素子位置情報とを記憶した位置情報記憶部と、を備え、前記加算率決定部が、前記画素位置情報及び前記欠陥素子位置情報に基づいて、信号値を算出する画素に含まれる前記欠陥素子の出力信号の加算率と、信号値を算出する画素の中心に対して当該欠陥素子と対称に位置する対角検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定し、他の前記検出素子の加算率を略同一の値となるように決定する。
【００１２】
以下、より具体的に本発明の実施形態について説明する。
＜第１の実施形態＞
以下、本発明の実施形態に係るＸ線検出器について図面を参照して説明する。本実施形態では、Ｘ線検出器をＸ線ＣＴ装置に適用した例について説明する。
【００１３】
図１に示すように、Ｘ線ＣＴ装置は、撮影系としての、Ｘ線源１００と、Ｘ線検出器１１１と、これらＸ線源１００及び検出器１１１の検出部１０４（後述）を対向配置し所定の回転軸を中心に回転するガントリー回転部１０１と、ガントリー回転部１０１の開口内に配置された寝台天板１０３と、これら撮影系の動作に伴いＸ線検出器１１１が取得した信号を処理する信号処理部１１２とを備えている。
【００１４】
Ｘ線源１００は、例えば管電圧で加速した電子ビームをタングステンやモリブデンなどのターゲット金属に衝突させ、その衝突位置（焦点）からＸ線を発生させる。
【００１５】
ガントリー回転部１０１は、Ｘ線源１００及び検出部１０４を互いに対向配置し、所定の回転軸を中心に回転する。ガントリー回転部１０１の中央には、被検体１０２が挿入される開口が設けられ、この開口内に、被検体１０２が寝かせられる寝台天板１０３が配置されている。寝台天板１０３とガントリー回転部１０１とは、所定の方向に相対的に移動可能となっている。
【００１６】
Ｘ線検出器１１１は、入射したＸ線フォトンを検出し、２つのエネルギー範囲に分別して計数を行うフォトンカウンティング方式の検出素子４００が複数配置された検出部１０４と、検出素子４００から出力される投影像を収集する信号収集部１０８とを備えている。Ｘ線検出器１１１の詳細は、後述する。
【００１７】
信号処理部１１２は、演算部１０５、表示部１０６、制御部１０７、主記憶部１０９及び入力部１１０を備えている。
演算部１０５は、収集した信号に所定の演算処理を行うため、信号収集部１０８で収集した信号に対して補正処理を行う補正処理部１０５２と、マルチエネルギー画像等の再構成像を作成する再構成処理部１０５３とを有している。
を備える。
【００１８】
表示部１０６は、演算部１０５により生成された再構成像などを表示する。制御部１０７は、Ｘ線源１００の発生駆動源の動作を制御するＸ線制御部、Ｘ線検出器１１１の信号読み出し動作を制御する読み出し制御部、ガントリー回転部１０１の回転と寝台天板１０３の移動を制御する撮影制御部、及びこれら各部全体を制御する全体制御部を備えている。主記憶部１０９は、演算部１０５のける演算処理に用いられるパラメータやデー等を記憶している。入力部１１０は、Ｘ線ＣＴ装置における撮影条件等の入力を行う。
【００１９】
演算部１０５及び制御部１０７は、その一部又は全部をＣＰＵ（中央処理装置）、メモリ及び主記憶部１０９を含むシステムとして構築することができ、演算部１０５及び制御部１０７を構成する各部の機能は、予め記憶部に格納されたプログラムをＣＰＵがメモリにロードし、実行することにより実現することができる。また機能の一部は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）や  ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）などのハードウェアで構成することも可能である。
【００２０】
以下の説明においては、特に説明しない限り、上述した撮影系、制御部１０及び信号処理部１１２を構成する要素は、公知のＸ線ＣＴ装置が備える要素と同様の構成を有し、同様の機能を持つ。
【００２１】
続いて、Ｘ線検出器１１１について説明する。
Ｘ線検出器１１１は、検出部１０４と、検出部１０４の各検出素子４００からの出力信号を投影像として収集する信号収集部１０８を備えている。
【００２２】
検出部１０４は、Ｘ線を検出する複数の検出素子４００が二次元配列された検出素子群を投影像の１画素に対応させて複数配列したものであり、検出部１０４に含まれる各検出素子４００は、所謂フォトンカウンティング方式の検出素子であり、入射したＸ線フォトンを検出し、例えば、２つのエネルギー範囲に分別して計数を行う。
【００２３】
検出部１０４は、例えば、図２に示すように、複数の同一のサイズの検出素子４００（図２中の実線）が、チャネル方向とスライス方向に２次元的に配置され、それぞれの方向で等間隔に配置された構造を成している。検出素子４００は、チャネル方向と回転方向を、スライス方向と回転軸方向を夫々一致させて配置されている。
【００２４】
図２の破線で示す矩形は、検出素子４００を二次元配列した検出素子群（画素）４１０を示しており、検出素子群４１０は、投影像における１画素に対応するようになっている（以下、検出素子群４１０を画素４１０と称する）。図２に示す例では、画素４１０は、チャネル方向に３素子、スライス方向に３素子の計９個の検出素子４００で構成される。一画素に対する検出素子４００の個数は一例であり、本発明を限定するものではない。
なお、図２では、検出部１０４に配置された検出素子４００の一部を示し、チャネル方向に７個、スライス方向に６個分を切り出して例として示したものである。
【００２５】
検出部１０４は、Ｘ線源１００を略中心とした円弧状に複数配置されており、ガントリー回転部１０１の回転に伴い、Ｘ線源１００との位置関係を保ちながら回転する。なお、図１では、説明の便宜上、検出部１０４が８個配列された例を示しているが、実際の装置では、例えば４０個程度配置されている。また検出部１０４の前面にはＸ線グリッド（図示せず）が設置されており、Ｘ線源１００から照射されたＸ線のうち、被検体１０２などで散乱されたＸ線が、検出部１０４に入射するのを防止する。
【００２６】
検出部１０４の各検出素子４００は、例えば図３に示すように、検出層４０１を挟むように正負の電極４０２、４０３が設けられ、各電極に読み出し回路４０５が接続された構造を有している。本実施形態では、正の電極４０２が、各検出素子４００間で共通する所謂共通電極である。
【００２７】
検出層４０１は、例えばＣｄＴｅ（テルル化カドミウム）、ＣｄＺｎＴｅ（カドミジンクテルル）、Ｓｉ（シリコン）などの半導体材料から成る。Ｘ線は、矢印４０４で示すように、正の電極４０２側から検出層４０１に入射し、Ｘ線フォトンを検出してそのエネルギーに応じた量の電荷を生じる。
【００２８】
図３に示すように、検出部１０４からのアナログ信号は、信号収集部１０８に入力される。収集部１０８は、読出し回路４０５、加算率決定部４０７、位置情報記憶部４０６、加算部４０８、及び処理制御部４０９を備えている。
【００２９】
読出し回路４０５は、Ｘ線の入射をトリガとして読み出し、エネルギー分別とデジタル変換を行う。読出し回路４０５では、入射した電荷で生じた電気信号を発生させ、所定の閾値により複数のエネルギー範囲に分別する。このとき、発生した電気信号の波高や発生量は、入射したＸ線フォトンのエネルギーに依存するため、Ｘ線フォトンに応じたエネルギー範囲に分別できる。次に分別結果を受けて、エネルギー範囲毎に入射したＸ線フォトン数をデジタルカウントし、デジタル信号を得る。
【００３０】
この分別方法は、例えば、２つのエネルギー範囲を、所定の閾値未満であるエネルギー範囲（以下、低エネルギー範囲と記す）か、所定の閾値以上のエネルギー範囲（以下、高エネルギー範囲と記す）かによって判別する。このような判別をサンプリング毎に行い、Ｘ線フォトンが入射したときに高エネルギー範囲と低エネルギー範囲に分別して、それぞれのＸ線フォトン数をビュー毎にデジタル信号でカウントする。
【００３１】
分別方法の一例について、図４を参照して説明する。図４は、発生電荷で生じた電圧１２０を示すグラフであり、横軸１２８が時間、縦軸１２９が電圧を示す。図示する例では、サンプリング時間１２３中にＸ線が入射してパルス出力１２１を生じ、サンプリング時間１２５中にＸ線が入射してパルス出力１２２を生じている。なお図４では、サンプリングはＸ線が入射するタイミングだけでなく、Ｘ線が入射しない場合（サンプリング時間１２４）でも周期的に行われる場合を示したが、Ｘ線フォトンが入射したタイミングでサンプリングが行われる場合も在り得る。
【００３２】
読み出し回路４０５は、サンプリング毎に、その区間における出力電圧の最大値と、エネルギー閾値１２６とエネルギー閾値１２７と比較して分別する。このエネルギー閾値１２６は、入射したＸ線フォトンが高エネルギー範囲か低エネルギー範囲かに分別するものである。エネルギー閾値１２７は、Ｘ線フォトンの入力無しかを判定するものである。ここで出力電圧１２０は、Ｘ線が入力しないときにも検出部１０４の回路ノイズによって変動しているため、これをＸ線による信号と誤検出しないためには、エネルギー閾値１２７はゼロより大きな値が必要である。
【００３３】
これらのエネルギー閾値を用いることで、例えば図４のサンプリング時間１２４では、出力電圧１２０はエネルギー閾値１２７以下のため、「Ｘ線フォトンの入力無し」と判断する。また、サンプリング時間１２５では、出力電圧１２０はエネルギー閾値１２６よりも大きいため、高エネルギー範囲のＸ線の「入射有り」と判断する。またサンプリング時間１２３では、出力電圧１２０はエネルギー閾値１２７よりも大きいがエネルギー閾値１２６以下のため、低エネルギー範囲のＸ線が入射したと判断する。以上のようにして、入射の有無とエネルギー範囲の分別を行う。
なおサンプリングでの最大値を用いて分別を行う代わりに、例えば、サンプリング中の出力電圧の積分値を用いてもよく、分別手法は上記手法に限定されない。
【００３４】
位置情報記憶部４０６は、画素と当該画素に対応する検出素子群に属する検出素子との位置関係を示す画素位置情報と、検出器１０４が備える複数の検出素子に含まれる欠陥素子（欠陥を有する検出素子）の位置を示す欠陥素子位置情報とを記憶している。すなわち、画素位置情報とは、画素４１０毎に、何れの検出素子４００が当該画素４１０のどこに位置しているかを示す情報であり、欠陥素子位置情報とは、欠陥素子が何れの検出素子群（画素）４１０に属し、当該画素４１０のどこに位置しているかを示す情報である。画素位置情報及び欠陥素子位置情報は、位置情報記憶部４０６に予め記憶させておくことができる。
【００３５】
図５〜図７は、図２と同様に、検出部１０４の一部を示し、投影像の画素４１０と検出素子４００の位置関係を示している。また、以下の説明において、図５〜図７に示す複数の検出素子のうち、検出素子４００−１−４、４００−５−２、４００−４−５が欠陥素子であるものとする。
図６は、記憶部４０６に記憶された欠陥素子位置情報の一例を示し、図６（ａ）は、図５の検出素子４００の配列において、正常な検出素子である場合は０、欠陥素子である場合は１として表現している。記憶部４０６は、欠陥素子位置情報として、例えば図６（ｂ）の配列マップを記憶する。
【００３６】
このような欠陥素子位置情報は、例えば被写体を設けずにＸ線を当てた際に、検出素子４００の出力を得て、出力値が基準値より小さい場合を欠陥として判断して作成する。先に述べたように本実施形態に係るＸ線ＣＴ装置ではエネルギー範囲毎に出力値を得るが、その一方でも基準値より小さい場合、欠陥と判断する。この作成は、装置が自動に行っても、人が判断して行っても良い。
【００３７】
また検出素子４００の出力を加算せずに得る場合は、画素中の検出素子数倍（本実施形態では９倍）のデータを出力する必要があるが、例えば、加算部にて画素中の１個の検出素子４００の出力値のみ出力して投影像を作り、この選択した検出素子４００を画素内で切り替えて、画素中の検出素子数個（本実施形態では９個）の投影像を得れば良い。ただし、ここに記す欠陥素子位置情報の作成方法や計測方法は一例であり、欠陥素子の位置を示すあらゆる情報を適用するができる。
【００３８】
図７は、画素位置情報の一例を示し、図７（ａ）は、図５の検出素子４００の配列において、画素の中心の検出素子４００の場合は３、角の場合は１、それ以外の場合は２として表現したものである。
【００３９】
記憶部４０６は、画素位置情報として、例えば図７（ｂ）の配列マップを記憶する。このような画素位置情報は一例であり、例えば、位置毎に、該当する検出素子４００の位置情報を座標情報として記憶することもできる。更に、座標に周期性がある場合は、数式などの形態で記憶することもできる。また、画素の中心位置を、例えば座標や数式などで与えるものであっても良い。
【００４０】
加算率決定部４０７は、位置情報記憶部４０６に記憶された欠陥素子位置情報及び画素位置情報に基づいて、画素４１０に属する検出素子４００の加算率を夫々決定する。つまり、画素位置情報及び欠陥素子位置情報に基づいて、信号値を算出する画素４１０に含まれる欠陥素子の出力信号の加算率と、信号値を算出する画素４１０の中心に対して当該欠陥素子と対称に位置する検出素子の出力信号の加算率とを、同一かつ他の検出素子の加算率よりも低い値となるように決定し、他の検出素子の加算率を略同一の値となるように決定する。
【００４１】
ここで加算率とは、０以上の数であり、１の場合は全ての出力信号を使用する場合、０の場合は全ての出力信号を使用しない場合である。例えば、欠陥素子の出力信号を全て使用しない場合は０とすれば良く、正常な検出素子の出力信号を全て使用する場合は１とすれば良い。加算率を決定する手順、加算率の具体例及び画素毎の出力値算出方法については後述する。
【００４２】
加算部４０８では、加算率決定部４０７において決定した加算率に従って、画素４１０に属する各検出素子４００の出力信号を重み付け加算することにより、投影像の画素毎の出力値として、画素４１０の出力値を算出する。本実施形態においては、一例として、チャネル方向に３素子、スライス方向に３素子の計９個の検出素子４００で構成される画素４１０を示した（図２参照）。従って、加算部８は、９個の検出素子４００の出力信号を夫々重み付け加算することにより画素４１０の出力値として算出する。
【００４３】
この重み付け加算は、例えば、画素４１０の出力値をＲｃｅｌｌ、その中の検出素子４００−ｉ−ｊ（ｉ，ｊは１から３の整数）の出力をｒ（ｉ，ｊ）とし、それらの重みである加算率α（ｉ，ｊ）とするとき、以下の数式（１）に従って算出することができる。
【００４４】
【数１】
0006619258.tif 000002
【００４５】
ここで、Ｍ（Ｍは自然数）は１つの画素中のチャネル方向の検出素子数、Ｎ（Ｎは自然数）は１つの画素中のチャネル方向の検出素子数であり、本実施形態では、ＭもＮも３である。
他の画素の出力値も同様に、その画素の中の検出素子４００の出力値に対して加算率で重み付け加算して算出する。これらの重み付け加算は、エネルギー範囲毎に行う。
【００４６】
処理制御部４０９は、上述の処理を行うために信号収集部１０８に含まれる各部の制御を行う。信号収集部１０８において算出した画素４１０の出力値は、投影像として、演算部１０５に出力される。
【００４７】
以上の構成を踏まえ、Ｘ線ＣＴ装置の一般的な撮影動作を説明する。なお、以下の説明において、説明の便宜上、フォトンのエネルギー範囲を２つとしているが、エネルギー範囲を３つ以上設けることができる。
まず撮影者が、入力部１１０から撮影条件を入力して実撮影の開始を入力すると、制御部１０７はＸ線源１００からのＸ線の照射と、ガントリー回転部１０１を制御し撮影を開始する。
【００４８】
このとき、例えば１２０ｋＶの管電圧で電子ビームを加速してＸ線源１００からＸ線が照射する。Ｘ線源１００の焦点から照射されたＸ線は、寝台天板１０３に載った被検体３００に向けて照射され、被検体３００を透過したＸ線は検出部１０４で検出される。検出部１０４は、検出素子４００毎に、入射Ｘ線のエネルギーに応じた電荷を発生する。信号収集部１０８は、前述したように、この電荷を高エネルギー範囲と低エネルギー範囲に分別し、エネルギー範囲毎及びビュー毎に、デジタルのカウント値を得る。更に投影像の画素毎に出力値を求め、演算部１０５に出力する。
【００４９】
次に制御部１０７は、このような撮影を、ガントリー回転部１０１を回転方向に回転することで、被検体３００に対するＸ線の照射角度を変化させる。このビューでも前ビューと同じように計測を行い、各エネルギー範囲でのカウント数として出力する。ここで、Ｘ線源１００から発生されるＸ線は、ビューに同期したパルスＸ線でも良いし、連続Ｘ線でも良い。更にこのように回転駆動させながら、ビュー毎に焦点位置を変更させて撮影を繰り返し行い、３６０度分のデジタル信号を取得する。撮影は、例えば０．４度ごとに複数ビューの間、行う。このような撮影により、３６０度分の投影像を得ることができる。
【００５０】
次いで、演算部１０５は、信号収集部１０８が収集した投影像に対し、所定の補正処理や演算処理を行い、マルチエネルギー投影像を作成する。補正処理では、例えばエア補正、演算処理では、例えば、密度画像作成、マルチエネルギー画像用投影像作成、再構成処理を行う。
【００５１】
以下に、加算部４０８における加算方法の一例について図８のフローチャートに従って説明する。
図８に示すように、まずステップＳ１１０にて、加算部４０８は、加算率決定部４０７において撮影前に決定された加算率を読出す。これは、例えば、システムの立ち上げ時や、撮影開始の入力指示が行われ、実撮影を開始する前に行う。
【００５２】
実撮影が開始されると、ステップＳ１１１では、加算部４０８が、読み出し回路４０５から１つのエネルギー範囲の１つの画素の全検出素子の出力を取得する。加算部４０８は、次のステップＳ１１２にて、加算率と検出素子の出力値から、１つのエネルギー範囲の１つの画素の出力値を求め、その結果を信号処理部１１２の演算部１０５に転送する。
【００５３】
加算部４０８は、次のステップＳ１１３において、全てのエネルギー範囲の出力値を算出したかを判定する。判定結果がＮＯである場合、ステップＳ１１１に戻り、同一画素の検出素子であって、まだステップＳ１１２を行っていない他のエネルギー範囲の出力を読み出し回路４０５から得て、同様にステップＳ１１２における処理を行う。一方、ステップＳ１１３の判定結果がＹＥＳである場合には、次のステップＳ１１４に移行し、全ての画素の出力値を算出したかを判定する。
【００５４】
ステップＳ１１４の判定結果がＮＯの場合、ステップＳ１１１に戻って同様に処理を行い、全エネルギー範囲の出力を読み出す。ステップＳ１１４の判定結果がＹＥＳの場合、処理を終了する。
【００５５】
上記の加算方法は一例であり、例えば、本方法の処理の順番を異ならせたり、読み出しや演算を複数のデータ毎にまとめて行ったりすることもできる。また、処理を行うデータの単位が異なる場合であっても同世の処理を行うことができる。さらに、加算部４０８が、出力値を画素毎に信号処理部１１２に転送するのではなく、全画素全エネルギー範囲のデータを取得した後にまとめて信号処理部１１２へ転送することもできる。
【００５６】
ここで、各検出素子の出力信号を重み付け加算して画素毎に出力値を得る場合における加算率α（ｉ，ｊ）とＳＮＲの関係は以下のようになる。
すなわち、検出素子４００−ｉ−ｊの出力信号のＳＮＲ（ｓｎｒ（ｉ，ｊ）と記す）がすべて同一とし、重み付け加算後の画素の出力値のＳＮＲ（ＳＮＲｃｅｌｌと記す）とすると、これらは式（１）から、式（２）のような関係となる。
【００５７】
【数２】
0006619258.tif 000003
【００５８】
（加算率決定方法について）
加算率決定部４０７における重み付け（加算率）の決定方法を、具体的に説明する。加算率の決定方法は、欠陥素子の有無とその位置によって異なる。
前提として、図５の画素４１０−１のように欠陥素子が無い場合、加算率は全ての検出素子４００で１とする。一方、画素４１０−２、画素４１０−３、画素４１０−４のように欠陥素子が有る場合、加算率決定部４０７は欠陥素子の加算率を、正常な場合から１だけ低減して０とする。
【００５９】
ところが、このままの加算率で重み付け合成を行うと、画素の出力値のサンプリング位置が画素中心からずれてしまう。つまり、例えば、画素４１０−４では、検出素子４００−４−５が欠陥であるが、この画素の加算率を０とすると、サンプリング位置がチャネル方向に、検出素子４００の１／９、画素の１／２７だけずれる。ただしこの量は、検出素子間のすき間が無い場合の値である。
【００６０】
サンプリング位置がずれてしまうと、以下のような影響が生じる。すなわち、例えば、被写体のエッジ付近などでは、隣接画素間で３倍程度出力が違う場合がある。この変化が指数的であると考え、サンプリング位置が検出素子４００の１／９（画素の１／２７）ずれた場合、出力は４％程度ずれてしまう。このようなずれは、特にＸ線ＣＴ装置では、アーチファクトの原因となる。
【００６１】
このため、加算率決定部４０７は、欠陥素子の加算率を０とすると共に、サンプリング位置が画素の中心となるように、画素内の加算率の重心が中心となるように他の加算率を変更する。
ここで、加算率の重心とは、始点からそれぞれの検出素子４００の中心へのベクトルとその加算率の積の総和を、その検出素子４００の個数で規格化したものである。より具体的には、画素がＮ（Ｎは２以上の整数）個の検出素子から成り、検出素子ｋ（ｋは１からＮの整数）への始点からのベクトルをｄ（ｋ）、その素子４００の加算率をα（ｋ）とする場合、加算率の重心Ｇは、式（６）のように表すことができる。
【００６２】
【数３】
0006619258.tif 000004
【００６３】
ベクトルの始点を画素の中心にとる場合つまり、始点が画素の中心と重なる場合には、加算率の重心はゼロベクトルとなる。以下、このように、始点を画素の中心にとったこととして説明する。
【００６４】
加算率の重心についての詳細を、図９を用いて説明する。図９は１つの画素内の９個の検出素子４００を表し、記載された（ｉ，ｊ）（ｉ，ｊは−１，０，１の整数）は、検出素子４００間の長さを単位長さとしたときに、画素中心を原点とした際の座標を表す。またｉはチャネル方向、ｊはスライス方向である。各検出素子４００の加算率をα（ｉ，ｊ）とし、加算率の重心のチャネル方向の座標をｘ、スライス方向の座標をｙとすると、それぞれは式（３）から、式（４−１）（４−２）（以下、併せて式（４）という）のように表すことができる。
【００６５】
【数４】
0006619258.tif 000005
【００６６】
上記の画素４１０−１の場合は、全ての検出素子４００は正常なので加算率は１となるので、加算率の重心（ｘ，ｙ）は（０，０）となり、画素の中心にあることが分かる。
【００６７】
次に、欠陥素子が有る場合、加算率決定部４０７は、欠陥素子が画素のどの位置にあるかを判別し、画素内の検出素子４００の位置に応じて加算率を決定する。欠陥素子が画素のどの位置にあるかについての判別は、記憶部４０６に記憶された欠陥素子位置情報及び画素位置情報に基づいて行う。
【００６８】
具体的には、例えば、図６の欠陥位置情報では欠陥素子は１で表されているので、画素４１０−２には欠陥素子４００−５−２があることが分かり、この検出素子４００は図１０の画素位置情報において数字が３であり、前述の定義から、画素の中心に位置することが分かる。
【００６９】
次に、以下、欠陥素子の位置に応じて加算率を決定する例について、図１０を用いて説明する。図１０の検出素子４００内に示す数字は欠陥素子が有る場合における各検出素子４００の加算率を表す。また、以下において、正常な検出素子の加算率値１から低減した量を、加算率の低減量と定義する。すなわち、例えば、当該検出素子が欠陥素子であり、加算率を０とした場合、加算率の低減量は１となる。同様に、正常な検出素子の加算率値１から増加した量を、加算率の増加量と定義する。
【００７０】
（欠陥素子が画素の角にある場合）
この場合は、例えば図６（ａ）における画素４１０−３は、検出素子４００−１−４が欠陥素子であり、欠陥素子が画素４１０−３の角に位置している。これは図９における（−１，１）が欠陥素子の場合に相当する。従って、欠陥素子４００−１−４の加算率のみを０とする、すなわち低減量と１とする場合、加算率の重心（ｘ，ｙ）は（１／９，−１／９）となり、画素の中心とならない。
【００７１】
そこで、加算率の重心を画素中心に合致させるために、加算率決定部４０７は、図１０に示すように、画素内において欠陥素子４００−１−４の画素中心（図１０では検出素子４００−２−５）に対して対称な位置にある対角検出素子である検出素子４００−３−６（図１０中、破線で示す検出素子）の加算率も、同様に低減量を１として加算率を０にする。
【００７２】
すなわち対角検出素子の加算率の低減量を、欠陥素子４００−１−４と同一にする。このように加算率を決定することにより、加算率の重心（ｘ，ｙ）は（０，０）となり、画素中心と一致させることができる。
【００７３】
このような方法は、所謂外挿を用いた方法よりも、ＳＮＲの点で優れる。図１１（ａ）は、図９における（−１，１）が欠陥素子であり、加算率を０とした場合の画素内の加算率を示し、図１１（ｂ）は、図９の欠陥素子の（−１，１）とその対角検出素子（１，−１）の加算率を０とした場合の加算率を示している。このように加算率を定めた場合、欠陥素子が無く、全ての検出素子４００の加算率が１の場合に対して、図１１（ａ）の場合は５．７％、図１１（ｂ）の場合は１１．８％、それぞれＳＮＲが低下する。
【００７４】
一方、図１１（ｃ）〜（ｆ）に示す外挿により求めた加算率の場合、ＳＮＲの低下の程度は以下のようなる。図１１（ｃ）は、斜め方向に、（０，０）と（１，−１）の検出素子４００の出力値を用いて、外挿により（−１，１）の欠陥素子の出力値を推定した場合の加算率である。このような加算率になることは、（０，０）の検出素子の出力値Ｐ（０，０）と（１，−１）の検出素子の出力値Ｐ（１，−１）から（−１，１）の検出素子の出力値Ｐ（−１，１）を推測するとき、出力値Ｐ（−１，１）の出力値は、以下の式（５）のように表すことができ、（０，０）の加算率を２だけ増加し、（１，−１）の加算率を１だけ低減する必要があることからも理解できる。また外挿を行った場合でも、対角検出素子の加算率は欠陥素子と同じ０となり、加算率の重心が画素の中心と一致する。
【００７５】
【数５】
0006619258.tif 000006
【００７６】
図１１（ｃ）のような加算率の場合、欠陥素子が無く、全ての検出素子４００の加算率が１の場合に対して２２．５％ほどＳＮＲが低下する。なお、ＳＮＲの低減率は、上述した式（２）から求めることができる。
【００７７】
同様に、図１１（ｄ）は、同一行の（０，１）と（１，１）の検出素子４００を用いて、外挿にて（−１，１）の出力値を求めた場合の加算率を示す。図１１（ｅ）は、同一列の（−１，１）と（−１，−１）の検出素子４００を用いて、外挿により（−１，１）の出力値を求めた場合の加算率を示す。これらの場合も、図１１（ｂ）と同様に、欠陥素子が無く、全ての検出素子４００の加算率が１の場合に対して２２．５％ほどＳＮＲが低下する。
【００７８】
同様に、図１１（ｆ）は、同一行、同一列、斜めのそれぞれの外挿をそれぞれ行って、それらを１／３ずつ足し合わせて欠陥素子の出力値を求めるときの加算率である。このとき、このとき（１，１）と（１，−１）と（−１，−１）の加算率は１／３だけ低減して２／３となり、（０，１）と（０，０）と（１，０）の加算率は２／３だけ増加して５／３となる。この場合、欠陥素子が無く、全ての検出素子４００の加算率が１の場合に対して１２．２％だけ、ＳＮＲが低下する。
【００７９】
以上述べた本実施形態におけるＸ線検出器の加算率決定部４０７において決定された加算率によって重み付け加算をして画素毎に出力値を算出した場合、所謂外挿よって定めた加算率を適用した場合に比して画素の出力値に対するＳＮＲの低減を抑えることができる。
【００８０】
これは、欠陥素子と対角検出素子以外の検出素子４００の加算率が、同一である場合に、最も良好に低減を抑制できる。このように、本実施形態における加算率決定部４０７では、欠陥素子及び対角検出素子の加算率を、正常な検出素子の加算率に対して同じだけ低減し、他の検出素子４００の加算率は一律に略同一の値、つまり、概ね一律とする。これにより、加算率の重心と画素の中心とを一致させ、ＳＮＲの低減を抑制することができる。なお、略同一の値、乃至は概ね一律とは、例えば加算率の違いが±１０％以内程度の場合をいう。
【００８１】
（欠陥素子が画素の角以外かつ中央以外にある場合）
例えば、図６（ａ）における画素４１０−４のように、検出素子４００−４−５が欠陥素子である場合は、図９における（−１，０）が欠陥素子である場合に相当する。
このような画素では、欠陥素子４００−４−５の加算率のみを０とする、すなわち低減量と１とする場合、加算率の重心（ｘ，ｙ）は（１／９，０）となり、加算率の重心と画素の中心とが一致しない。
【００８２】
そこで、欠陥素子４００−４−５の周辺素子の加算率を増加する。本実施形態では、図１０に示すように、検出素子４００−４−４と検出素子４００−４−６（実線の太線で記された検出素子）の加算率を０．５だけ増加する。すなわち、検出素子４００−４−４と検出素子４００−４−６の加算率を１．５（すなわち加算量０．５）、欠陥素子４００−４−５の加算率を０（すなわち低減量１）、他の検出素子４００の加算率を１とする。
【００８３】
これは検出素子４００−４−４と検出素子４００−４−６の出力信号から内挿により検出素子４００−４−５の出力信号を求めた場合と同じである。このように加算率を決定することで、加算率の重心（ｘ，ｙ）は（０，０）となり、画素中心と一致させることができる。
【００８４】
（欠陥素子が画素の中心にある場合）
例えば、図６（ａ）における画素４１０−３のように、検出素子４００−５−２が欠陥素子であるは、加算率の重心（ｘ，ｙ）は（０，０）であり、中心となる検出素子が欠陥素子そのものであるため、欠陥素子４００−４−５の加算率のみを０とし、他の検出素子４００の加算率は低減せずにそのままの値とする。
上記した例のように加算率を定めることで、投影像の画素４１０のサンプリング位置が実際の位置とずれることなく、各画素の出力値を算出することができる。
【００８５】
上述の各検出素子に対する加算率の決定処理は、概ね以下のように行われる。
まず、（１）全ての検出素子の加算率を１とし、全ての欠陥素子の加算率を１だけ低減する。次に、（２）欠陥素子の位置を欠陥素子位置情報に基づいて、画素中の角、中央、それ以外のどの位置にあるか判定する。続いて、（３）欠陥素子の位置が画素中の角の場合、対角検出素子の位置を決定し、対角検出素子が欠陥素子でない場合にはその加算率を１だけ低減する。欠陥素子の位置が画素中の角以外かつ中央以外の場合は、同一画素の周辺素子にて欠陥素子の出力を内挿するように、欠陥素子の周辺素子の加算率を増加させる。ただし、周辺素子が欠陥素子である場合などは、周辺素子にて欠陥素子の出力を内挿できない場合も有り得る。この場合は、欠陥素子の周辺素子の加算率は増加させずに、対角検出素子の加算率を１だけ低減する。
【００８６】
これらの工程は、欠陥素子が複数ある場合、上記（２）と（３）の処理を、全欠陥素子に対して繰り返すことで、加算率を決定することができる。
上記した例では、欠陥素子の低減量と、欠陥素子の位置が画素中の角の場合の対角検出素子の低減量とを１とする例について説明したが、例えば、低減量を０より大きく１未満の値とすることもできる。そして、この場合にも、加算率の低減量を、対角検出素子と欠陥素子とで同一にする。
【００８７】
以下、より詳細に、加算率決定部４０７において加算率を決定する処理について、図１２のフローチャートを用いて説明する。
加算率を決定するために、まずステップＳ２００において、加算率決定部４０７は、全検出素子の加算率を１とし、ステップＳ２０１に進み、全欠陥素子の加算率を１だけ低減する。このとき、検出素子が欠陥素子であるかの判定は、記憶部４０６に記憶された欠陥素子位置情報、すなわち、例えば、図６の配列マップに基づいて決定する。
【００８８】
次に、ステップＳ２０２において、画素中の何れの位置に欠陥素子が位置するかの判別を行う。これは、記憶部４０６に記憶された欠陥素子位置情報及び画素位置情報に基づいて、すなわち、例えば図６及び図７に示す配列マップを用いて決定する。欠陥素子の位置を判別した結果、検出素子が画素の角にある場合は次のステップＳ２０３に進み、欠陥素子が画素の中心及び角以外にある場合はステップＳ２０５に進み、欠陥素子が画素の中心である場合はステップＳ２０９に進む。
【００８９】
ステップＳ２０３では、欠陥素子の対角検出素子が欠陥素子であるかを判定し、対角にある検出素子が欠陥素子でない場合はステップＳ２０４に進み、欠陥素子の対角検出素子の加算率を１だけ低減し、次のステップＳ２０９に移行する。一方、ステップＳ２０３において、対角にある検出素子が欠陥素子であると判定された場合には、ステップＳ２０９に移行する。
【００９０】
また、ステップＳ２０２において、検出素子が画素の中心及び角以外にある場合は、ステップＳ２０５〜Ｓ２０８にて、周辺素子にて欠陥素子の出力を内挿できるかどうかを判定して、加算率を決定する。具体的には、ステップＳ２０５では、欠陥素子と同一画素中において、欠陥素子の上下に検出素子があり、且つその上下の検出素子の双方が欠陥素子ではないことを判別する。これがＹＥＳ、すなわち、欠陥素子と同一画素中であって上下に正常な検出素子がある場合は、ステップＳ２０６に進み、その欠陥素子の上下の検出素子の加算率を０．５だけ増加する。
【００９１】
反対にステップＳ２０５の判定結果がＮＯの場合、すなわち、欠陥素子と同一画素中であり、その上下の少なくとも一方に正常な検出素子が無い場合は、ステップＳ２０７に移行する。ステップＳ２０７では、欠陥素子と同一画素中において、欠陥素子の左右に検出素子があり、かつその左右の検出素子の双方が欠陥素子ではないことを判別する。
【００９２】
これがＹＥＳ、すなわち、欠陥素子と同一画素中であって左右に正常な検出素子がある場合、ステップＳ２０８にて、その欠陥素子の左右の検出素子の加算率を０．５だけ増加する。反対にＮＯの場合、すなわち、欠陥素子と同一画素中であって、その左右の少なくとも一方に、正常な検出素子が無い場合、ステップＳ２０４にて、欠陥素子の対角検出素子の加算率を１だけ低減する。
【００９３】
このような処理により、欠陥素子の上下の検出素子が正常な場合には、それらの出力から内挿にて欠陥素子の出力を求めることと同一となる。上下の検出素子の少なくとも一方が正常でない場合には、左右の検出素子を用いて、内挿にて欠陥素子の出力を求めることと同一となる。更に左右の検出素子にも、少なくとも一方が欠陥素子である場合は、欠陥素子が画素の角にある場合と同様に、対角検出素子の加算率を１だけ減らすことになる。
【００９４】
次にステップＳ２０９にて、ステップＳ２０３で行った位置判定を全ての欠陥素子に対して行ったか判定する。行っていない場合、ステップＳ２０３に戻って、位置判定を行っていない次の欠陥素子に対して、ステップＳ２０３からステップＳ２０８を同様に行う。
【００９５】
このようにして全ての欠陥素子に対して上記処理を行った場合、ステップＳ２１０にて移行して、加算率決定処理を終了する。以上の処理により、欠陥素子が複数存在する場合でも、全検出素子における加算率を決定することができる。
【００９６】
上記した欠陥素子に隣接する検出素子の加算率を高める方法、すなわち周辺素子にて欠陥素子の出力を内挿できるかどうかを判定する方法は一例である。使用する周辺素子の増加量の重心が、欠陥素子の中心に位置するように決定すればよく、例えば、ステップＳ２０５とステップＳ２０７の判定の順序を異ならせることもでき、３つ以上の周辺素子で内挿できるか判定することもできる他、上下左右以外の他の周辺素子を用いることもできる。
【００９７】
以上のように加算率決定部４０７により加算率を決定することで、故障などによって生じた欠陥素子によって、投影像の画素のサンプリング位置がずれて出力値がずれることを抑制し、延いてはアーチファクトが発生することを簡易に除去・低減することができ、正確な出力値を得ることができる。
【００９８】
ところで、上述のように欠陥素子やその画素中心の対称位置の検出素子４００の加算率を低下させると、欠陥素子を含む画素の出力値が低下する。従って、画素の加算率の合計を、画素毎に同一となるように、規格化する必要がある。このため、例えば、加算部４０８に規格化部を設け、規格化部において各画素に含まれる検出素子の出力信号に対する加算率を決定した後に、画素間の加算率の違いを規格化する。
【００９９】
具体的には、例えば、画素に含まれる検出素子の加算率の合計で規格化することができる。すなわち、例えば、図１０において、画素４１０−１の加算率の合計値は９であるので、加算部４０８が、上記の式（１）に示す加算方法により出力値を算出した後、規格化部により出力値を加算率の合計値としての９で除算する。同様に、例えば、画素４１０−２の加算率の合計値は８であるので、規格化部により出力値に対して加算率の合計である８で割る。このように規格化部は、加算部４０８で求めた画素の出力値に対して、すなわち投影像に対して規格化を行うことで、画素毎に使用する検出素子数が異なることに起因する画素毎の出力値の相違を補正することができる。
【０１００】
この他、規格化部を設けずに、例えば、加算率決定部４０７において、各画素の加算率の合計が１となるように、各検出素子の加算率を予め決定しても良い。すなわち、例えば、図１０において、画素４１０−１の加算率の合計値９を用いて、画素４１０−１に含まれる各検出素子の加算率を９で割って１／９とする、画素４１０−２の加算率の合計値８を用いて、画素４１０−２に含まれる各検出素子の加算率を８で割って１／８とすることができる。
【０１０１】
また、さらに他の方法として、演算部１０５における補正処理部１０５２において、エア補正と共に規格化に係る処理を行うことができる。
ここで、エア補正について簡単に説明する。エア補正とは、例えば、本撮影の事前に計測し作成して主記憶部１０９に保存しておいた感度・Ｘ線分布データを用いて、投影像をエネルギー範囲毎に除することで実現する。感度・Ｘ線分布データは、例えば、被検体を設けずにＸ線管１００からＸ線を照射してエネルギー毎に投影像を取得し、投影像に対して検出素子４００毎にビュー方向に加算平均を行い、検出部１０４での出力の平均値によって規格化することでエネルギー範囲毎に作成する。この補正処理は、それぞれのエネルギー範囲で取得した投影像毎に行う。
【０１０２】
また、他の例として、信号収集部１０８や演算部１０５が、欠陥素子と対角検出素子の加算率の低下分を規格化することもできる。これは、例えば正常な画素における検出素子の加算率の合計をＳ、欠陥素子のある画素における加算率の合計をＴとするとき、欠陥素子がある画素の出力は、無い画素に比べてＴ／Ｓ倍の出力になるため、加算後の出力に対して、例えばＳ／Ｔ倍して規格化すれば良い。
【０１０３】
従って、例えば、補正処理部１０５２においてエア補正を行う際に、検出素子の感度が同一であるとき、規格化を行わない場合、当該検出素子に同一の線量が入力した際には、画素４１０−２の出力値は画素４１０−１の出力値の８／９倍になる。そして、感度・Ｘ線分布データにもおいても出力値の比は８／９倍となり、同様の重みがかかることになるため、エア補正を行うことで、同時に加算率も規格化されて画素４１０−１と画素４１０−２の出力を同一とすることができる。
【０１０４】
（加算率の決定を行うタイミングの他の例）
上記した第１の実施形態においては、加算率決定部４０による加算率の決定を、撮影後、すなわち、投影像を作成する際に行う場合について説明した。加算率の決定は上記した例に限られず、例えば、撮影の事前に行うこともできる。このようにすることで、撮影後に、加算率を決定する必要が無く、撮影から投影像の完成までの時間を短縮することが可能となる。
【０１０５】
（欠陥素子の位置に応じた加算率決定方法の他の例）
上記した例では、欠陥素子が画素の角及び中央以外にある場合、欠陥素子の周辺素子の加算率を増加する例について説明したが、対角検出素子の加算率を低減しても良い。この場合、例えば、欠陥素子４００−４−５と、その対角検出素子４００−６−５の加算率を同程度低減すればよい。例えば、図１３に、低減量が１であって加算率を０とする例を示す。なお、図１３における検出素子４００内の数字は加算率を示す。このようにすることで、加算率の重心（ｘ，ｙ）と画素の中心とを一致させることができる。
【０１０６】
上述の各検出素子に対する加算率の決定処理は、概ね以下のように行われる。
まず、（１）全ての検出素子の加算率を１とし、全ての欠陥素子の加算率を１だけ低減する。次に、（２）欠陥素子の位置を欠陥素子位置情報に基づいて、画素中の角、中央、それ以外のどの位置にあるか判定する。続いて、（３）欠陥素子の位置が画素中の角の場合、対角検出素子の位置を決定し、対角検出素子が欠陥素子でない場合にはその加算率を１だけ低減する。欠陥素子の位置が画素中の角以外かつ中央以外の場合は、欠陥素子位置が画素中の中央以外の場合、対角検出素子の位置を決定し、その加算率を１だけ低減する。
【０１０７】
以下、本例の加算率決定部４０７において加算率を決定する処理について、図１４のフローチャートを用いて説明する。
図１４に示すフローチャートでは、図１２のフローチャートと同様にステップＳ３００及びステップＳ３０１において、加算率決定部４０７は、全検出素子の加算率を１とし、ステップＳ２０１に進み、全欠陥素子の加算率を１だけ低減し、ステップＳ３０２の処理を行う。ステップＳ３０２では、欠陥素子が画素の中心に位置するか否かを判定する。
【０１０８】
ステップＳ３０２の判定結果がＹＥＳ、すなわち欠陥素子が画素の中心に位置する場合には、ステップＳ３０５に進む。一方、ステップＳ３０２の判定結果がＮＯ、すなわち欠陥素子が画素の中心にない場合には、ステップＳ３０３に移行し、欠陥素子の対角検出素子が欠陥素子であるかを判定する。
ステップＳ３０３の結果がＮＯの場合、すなわち欠陥素子の対角検出素子は正常な検出素子である場合はステップＳ３０４に移行して、欠陥素子の対角検出素子の加算率を１だけ低減し、その後にステップＳ３０９に移行する。
【０１０９】
一方、ステップＳ３０３の結果がＹＥＳの場合、すなわち欠陥素子の対角検出素子も欠陥素子である場合はステップＳ３０５に移行する。
次にステップＳ３０５では、全ての欠陥素子に対して位置判定を実施するように処理を反復する。このようにして全ての欠陥素子に実施した後に加算率決定処理を終了する。
このような処理は、上述した欠陥素子が画素の角及び中央以外にある場合に、欠陥素子の周辺素子の加算率を増加する例と比較すると、簡易な処理であり、ＳＮＲはやや劣るものの高速な処理に適している。
【０１１０】
ＳＮＲについて、例えば、欠陥素子が１つで、全ての検出素子４００に一様にＸ線が入射した場合において、欠陥素子が画素の角及び中央以外にある場合に、欠陥素子の周辺素子の加算率を増加する例では、その画素のＳＮＲは、全て正常の検出素子から成る画素に比べて約７．４％低下する。一方、欠陥素子が画素の角及び中央以外にある場合に、対角検出素子の加算率を低減する例では１１．８％低下してしまう。
【０１１１】
従って、対角検出素子の加算率を低下して０にする方法は、検出素子４００が画素の角にあるときなど、補間では加算率の重心を画素の中心にすることが難しいときに適用することが望ましい。
【０１１２】
なお、この他、例えば、欠陥素子が画素の角に位置するものの、隣接画素の検出素子を用いて加算率の重心が画素の中心となるように内挿により欠陥素子の出力信号の推定を行うことが可能な欠陥素子に対して、欠陥素子の周辺の検出素子の加算率を増加しても良い。すなわち、隣接の検出素子４００から内挿して値を決定しても良い。ただし、加算率の重心を欠陥素子の中心にすることが難しい場合や、加算部４０８で隣接の画素の検出素子４００の出力値を使用するのが難しい場合などには、欠陥素子の対角検出素子の加算率を低減する方が望ましい。
【０１１３】
このような場合の例としては、（１）加算部４０８の複数のブロック毎に回路が構成され、加算前に回路間でのデータのやり取りが困難である場合、（２）処理の遅延を生じる場合、（３）検出部１０４の端部の画素であることから、隣接の検出部１０４に属する検出素子の出力値を加算前にやり取りするのが困難である場合、（４）Ｘ線照射野の端部の画素で隣接画素はＸ線照射野外であって、隣接画素の検出素子を用いて内挿が実質的に難しい場合、（５）隣接画素まで距離があり隣接画素の検出素子を用いて内挿しても精度が低い場合、（６）実質的に挟む画素が存在しない場合、などが考えられる。
【０１１４】
また、欠陥素子が画素の中央にある場合、欠陥素子の加算率を１だけ低減して０にする例について説明したが、加算率の重心が画素の中心となるように加算率を定めれば、例えば、周辺の検出素子４００から内挿して欠陥素子の値を求めるように、欠陥素子の周辺の加算率を増加することもできる。
【０１１５】
また、周辺の検出素子４００の加算率を上げる場合は、画素内の欠陥素子以外の検出素子４００の加算率が一律に近い方がＳＮＲの低下は少なくなるが、画素内の欠陥素子以外の全ての検出素子４００の出力値の平均から欠陥素子の出力値を求める場合、これは欠陥素子以外の画素内の全ての検出素子４００の加算率が同一となり、最もＳＮＲとして低減が少ない加算率の決定方法ということができる。
【０１１６】
（画素中の検出素子数の他の例）
以上の説明において、１つの画素の出力値を３×３個の検出素子から求めることとしたが、検出素子の数やその配列方法は種々考えられる。１画素を構成する検出素子が２個である場合には、一方の検出素子が欠陥素子であると、他方が対角検出素子となり、両者の加算率を０とすることにより画素の出力値が０となってしまう。従って、画素内の検出素子数は３個以上とすることが好ましい。
【０１１７】
ＳＮＲの観点からは、画素内の検出素子数が多いほど好ましい。例えば、図１５に、画素４１０が６×６個の検出素子４００から成る場合の例を示す。図１５の例では、（２，２）の検出素子４００を欠陥素子としている。
【０１１８】
この場合、欠陥素子のみの加算率を０とした場合、欠陥素子が無い場合に比べて、ＳＮＲは１．４％低下する。一方、欠陥素子と、対角検出素子（５，５）の加算率を０とした場合、欠陥素子が無い場合に比べて、ＳＮＲは２．８％低下する。
また、画素が３×３個の検出素子からなる場合には、ＳＮＲの低下が１２％であった。
つまり、画素内の検出素子４００の個数が多い場合は、欠陥素子と対角検出素子の両方の加算率を下げる方法を用いても、ＳＮＲの低下が小さく抑えられるということができる。
【０１１９】
（信号収集部１０８の搭載位置の他の例）
信号収集部１０８が回転部１０１に設けられている場合を記したが、一部が静止系に設けられていても良い。例えば、演算部１０５の一部であっても良い。
【０１２０】
（欠陥素子が画素中に複数個ある場合の例）
上述の例では、１画素を構成する複数の検出素子に、欠陥素子が１個だけ含まれる場合について説明した。１画素に、欠陥素子が複数含まれている場合には、例えば以下のように考えることができる。
【０１２１】
図１６（ａ）において、（−１，１）と（０，１）の２つの検出素子４００が欠陥素子であるとする。このとき、図１６（ａ）は、図１６（ｂ）と図１６（ｃ）との組み合わせであると考えることができる。
図１６（ｂ）は、欠陥素子（−１，１）と、それに対して画素中心に対して対称な（１，−１）の検出素子４００の加算率を１だけ低減する。
図１６（ｃ）は、欠陥素子の周辺の検出素子４００である（−１，−１）と（１，−１）の加算率を０．５だけ増加する。このとき（１，−１）は、低減量が１で増加率が０．５となって、加算率は０．５（すなわち低減量０．５）となる。
【０１２２】
従って、図１６（ａ）では、（１，−１）の検出素子４００の加算率を０．５、（−１，−１）の検出素子４００の加算率を１．５、（−１，１）の検出素子４００の加算率を０、他の検出素子４００の加算率を１とすれば良い。この時の加算率を図１７（ａ）に示す。
【０１２３】
また、上述したように、図１６（ｃ）の場合に対して、対角検出素子である（１，０）の検出素子４００の加算率を１下げて０としても良い。この時の加算率を図１７（ｂ）に示す。
【０１２４】
ＳＮＲの低下は、図１７（ａ）の場合は１４．８％であり、図１７（ｂ）の場合は２５．５％となり、欠陥素子の周辺の検出素子の加算率を増加させる方が、ＳＮＲの低減が少ないことが分かる。
【０１２５】
また他の例として、図１８（ａ）に示すように、１つの画素の、（−１，１）と（−１，０）の２つの検出素子４００が欠陥素子である場合について考える。
同様に、図１８（ａ）も、図１８（ｂ）と図１８（ｃ）との組み合わせであると考えることができる。
【０１２６】
図１８（ｂ）は、対角検出素子である（１，−１）の検出素子の加算率を１だけ低減する。
一方、図１８（ｃ）は、周辺素子である（−１，１）の検出素子４００が欠陥素子であって内挿による推定ができないことから、図１８（ｃ）も、対角検出素子である（１，０）の検出素子の加算率を１だけ低減すれば良い。
【０１２７】
このように考えると、欠陥素子である（−１，１）（−１，０）と、対角検出素子である（１，０）（１，−１）の検出素子の加算率を０（すなわち低減量１）とし、他の検出素子４００を１にすれば良い。
【０１２８】
更に、図１９に、画素中心に対して対称な検出素子４００同士が欠陥素子である例を示す。図１９では、（−１，１）と（１，１）の検出素子が欠陥素子である。この場合は、欠陥素子のそれぞれの加算率を０とすれば、加算率の重心は画像の中心となっているので、他の検出素子に対して加算率を低減する必要はない。
このように、他の位置に２つの欠陥素子が有る場合や、３つ以上ある場合も、１つ１つに分けて、それらの組み合わせとして考えことで、加算率を低減する検出素子４００を決定することができる。
【０１２９】
（欠陥素子の位置情報の他の例）
上記した例では、欠陥素子位置情報の例として図６に示す配列マップを挙げたが、この他、検出素子４００が正常か欠陥かを示す数値や、配列マップの形状は種々考えられる。また、配列マップに限られず、例えば、欠陥素子のみ、または欠陥素子のみの、座標のような位置情報として記憶することもできる。なお、位置情報に加えて、欠陥素子と対応させて、当該欠陥素子の加算率を併せて記憶させておくこともできる。
【０１３０】
（対角検出素子の決定方法の他の例）
上述の例では、加算率を低減する対角検出素子の位置を、記憶部４０６に記憶した欠陥素子位置情報と、画素と検出素子４００の位置関係を表す画素位置情報とに基づいて決定していた。画素位置情報としては、加算率を低減する対角検出素子の位置を、例えば、図２０に示すような配列マップの形態で保存することもできる。
【０１３１】
図２０は、欠陥素子が図６に示す位置に存在している場合の対角検出素子の位置を表す一例であり、太い矩形で示され検出素子内に「１」と示された検出素子４００が対角検出素子である。欠陥素子の位置情報と同様に、様々な形態の配列マップや、例えば、座標のような位置情報などとして、記憶することもできる。
【０１３２】
また、加算率を低減する対角検出素子の位置情報を、欠陥素子位置情報と共に保存することもできる。すなわち、画素位置情報が欠陥素子位置情報を兼ねても良い。図２１はその一例であり、欠陥素子の位置情報に、加算率を低減する対角検出素子を記したものと画素位置情報としている。図２１において、加算率を低減する正常の検出素子４００は２として示している。このとき加算率決定部４０７は、この欠陥素子の位置情報を読み込むことにより、同時に対角検出素子の位置情報を得ることができる。
【０１３３】
なお、図２１において、正常な検出素子は０、欠陥素子は１、対角検出素子は２で区別して示したが、これは一例に過ぎない。
更に、欠陥素子の位置情報や画素位置情報は、欠陥素子や対角検出素子の位置情報だけでなく、加算率を有していても良い。図２２はこの一例の画素位置情報であり、検出素子４００内の数字は、加算率を表す。このとき加算率決定部４０７は、この欠陥素子の位置情報を読み込んで対角検出素子の位置と加算率を得れば、即座に全ての検出素子の加算率を決定できる。
【０１３４】
（加算率の決定方法の他の例）
欠陥素子の加算率を０としたが、例えば、０より大きく１未満の様々な値を適用することができる。このとき、各検出素子間において加算率の重心が画素の中心と一致させるために、画素中の検出素子４００の大きさがすべて同一の場合、欠陥素子と対角検出素子の加算率に対する低減量を同じとする。
【０１３５】
また、上述の例では欠陥素子の加算率を、エネルギー範囲によらず同一としたが、欠陥素子の加算率の低減量をエネルギー範囲によって異ならせても良い。特に、フォトンカウンティング検出器では、例えば雑音が大きいため、低いエネルギー範囲では誤検出が多いため使用できないが、高いエネルギー範囲では正確に計数できる場合もある。
【０１３６】
このとき、低エネルギー範囲の加算率のみ０以上１未満に下げ、高エネルギー範囲の加算率は下げないように加算率を定めることもできる。また、高エネルギー範囲の加算率の低減量を、低エネルギー範囲の加算率の低減量よりも小さくなるように加算率を定めることもできる。
【０１３７】
この他、エネルギー範囲に応じて、対角検出素子の加算率の決定方法を異ならせても良い。例えば、図６（ａ）の画素４１０−４において、低エネルギー範囲では欠陥素子の周辺素子の加算率を増加させ、高エネルギー範囲では欠陥素子の対角に位置する検出素子の加算率を低減させることができる。このように、欠陥素子と対角検出素子の少なくとも一方の加算率が、エネルギー範囲に応じて異ならせることができる。
【０１３８】
（検出素子の配置及び形状の他の例）
本実施形態において、検出部１０４は、チャネル方向及びスライス方向のそれぞれに検出素子４００が等間隔に配列された例を示したが、検出素子の配列が等間隔でない場合でも、上述の例と同様に加算率を決定することができる。また、検出素子４００の形状についても必ずしも等方的な形状である必要はなく、非等方的な形状であっても、加算率の重心が画素の中心と一致するように加算率を決定することができる。そして、このようにすることで、処理時間、処理回路及び補間用データ等を増大させることなく、欠陥素子に対する補間の精度を向上させ、簡易にアーチファクトを抑制する
【０１３９】
（加算率の重心の他の計算方法）
上述の例では、加算率を用いて加算率の重心を算出、すなわち、式（３）を用いて算出した。加算率の重心は、加算率の低減量又は増加量からも算出することができる。この場合には、式（３）における検出素子４００−ｋ（ｋは１からＮの整数）の加算率の低減量をβ（ｋ）（＝１−α（ｋ））とすることにより、以下の、式（６）により加算率の重心を算出することができる。なお、低減量β（ｋ）が負となった場合は、増加率を表すと見なせば良い。
【０１４０】
【数６】
0006619258.tif 000007
【０１４１】
この式（６）は、全ての加算率が１の場合、加算率の重心は画素中心となりゼロベクトルとなることを利用して、式（３）から導出したものである。このような方法では、全ての検出素子を考える必要がなく、加算率を低減や増加する検出素子のみを考えれば良いので、計算がし易いという利点がある。
【０１４２】
（検出素子４００のサイズが異なる場合の一例）
本実施形態では、画素を構成する検出素子がすべて同一の大きさである例について説明した。１画素を構成する検出素子の大きさは、必ずしもすべて同一である必要はなく、例えば、図２３のように１画素に互いに異なる大きさの検出素子を混在させることもできる。
【０１４３】
この場合には、加算率の重心の求め方が異なる。１画素を構成する検出素子の大きさ、すなわち、面積がすべて同一である場合は、上述の式（３）に従って、始点からそれぞれの検出素子４００の中心へのベクトルとその加算率の積の総和を、その検出素子４００の個数で規格化することができる。
【０１４４】
しかしながら、検出素子の面積が一律でない場合、各検出素子の加算率に面積に応じた重み付けをする必要がある。すなわち、始点からそれぞれの検出素子の中心へのベクトルとその加算率の積に対して、面積の重み付けをする必要がある。従って、検出素子ｋ（ｋは１からＮの整数）の面積をＳ（ｋ）、画素内の検出素子の全面積をＳａｌｌとするとき、加算率の重心Ｇは、式（７）のように書ける。
【０１４５】
【数７】
0006619258.tif 000008
【０１４６】
ここで面積Ｓ（ｋ）に代えて、検出素子の面積比を用いても良い。このとき全面積Ｓａｌｌとして、面積比の合計を用いれば良い。このように考えると、図２３の例では、欠陥素子である検出素子４００−Ａの加算率を０とする際には、（０，−１）と（１，−１）の検出素子の加算率も０とすれば良い。このとき、（０，−１）と（１，−１）の検出素子の面積の合計は、（０，−１）と（１，−１）の検出素子の間の空間４１１だけ検出素子４００−Ａよりも小さく見える。
【０１４７】
これは、例えば、検出素子が半導体材料で作製されている場合、空間４１１に入射したＸ線によって発生した電荷は、（０，−１）と（１，−１）の何れか一方の検出素子で検出されるため、空間４１１はおよそ０と考えることができる。従って、（０，−１）と（１，−１）の検出素子の面積の総和と検出素子４００−Ａの面積とは同一と見做すことができる。
【０１４８】
すなわち、画素内の検出素子が同じ大きさであれば、対角検出素子の加算率の低減量は欠陥素子と同じにする必要がある。一方、画素内の検出素子が同じ大きさでなければ、対角検出素子の加算率の低減量と面積との積が、欠陥素子の加算率の低減量と面積との積と同じとなるように対角検出素子の加算率の低減量を決める必要がある。
【０１４９】
言い換えると、欠陥素子と等しい面積に相当する対角検出素子の加算率の低減量を欠陥素子の加算率の低減量と同一にする場合、欠陥素子の面積と等しい面積になるように対角検出素子を選択する必要がある。つまり、対角検出素子の加算率の低減量は、欠陥素子の加算量の低減量に、欠陥素子の面積と対角検出素子の面積の比を乗じた値にする必要がある。この場合、低減量が１以上とならないように、対角検出素子の総面積が欠陥素子の面積と同一又は大きくなるように選択することが望ましい。更に別の言い方をすると、対角検出素子の面積は、欠陥素子の面積に、欠陥素子と対角検出素子の加算率の低減量の比を乗じた値にする必要がある。
【０１５０】
そして、このように１画素に含まれる検出素子の大きさが互いに異なる場合でも、画素内で、画素の中心に対して対称の位置の検出素子の大きさが同一の場合には、対称に位置する検出素子の加算率を低減する方法を用いることができる。
【０１５１】
なお、面積や面積比等の面積情報を考慮して加算率を決定するとき、例えば、以下のような手順で加算率を決定することができる。
（１）全ての検出素子の加算率を１とし、欠陥素子の加算率を１低減する。
（２）欠陥素子の位置が画素の何れに位置するか（角、中央、それ以外）を判定する。
（３）欠陥素子の位置が画素中の角の場合、対角検出素子の位置を決定し、欠陥素子の面積を対角検出素子の面積で除し、欠陥素子の加算率の低減量を乗じた分（＝欠陥素子の面積÷対角検出素子の面積×欠陥素子の加算率の低減量）だけ、対角検出素子の加算率を低減する。欠陥素子の位置が画素中の角以外かつ中央以外の場合、周辺素子にて欠陥素子の出力を内挿するように、欠陥素子の周辺素子の加算率を増加させる。
（４）画素の加算率の合計値の違いを規格化する。
【０１５２】
上述の手順は、欠陥素子が複数ある場合、手順（２）と手順（３）とを、欠陥素子毎に繰り返すことで、加算率を決定できる。ここでは、手順（１）において全検出素子の加算率を１としたが、０より大きい値、且つ一律の値であれば良い。
【０１５３】
また、欠陥素子の位置が画素中の角以外かつ中央以外の場合に、対角検出素子の加算率を低減するとき、手順（３）に代えて以下の（３’）のようにすることが好ましい。
（３’）欠陥素子位置が画素中の中央以外の場合、対角検出素子の位置を決定し、欠陥素子の面積を対角検出素子の面積で割った値に欠陥素子の加算率の低減量を書けた分だけ、対角検出素子の加算率を低減する。
【０１５４】
（検出部の他の例）
本実施形態では、検出部１０４にＸ線を直接検出する半導体検出器を形容した例を示したが、例えば、シンチレータと半導体光検出器とから構成された検出器を適用することもできる。この場合には、Ｘ線を検出して光に変換し、半導体光検出器で光を電気信号に変換する。これらの変換をＸ線フォトン毎に行い、発生した電気信号にて、入射Ｘ線のエネルギーを分別する。
【０１５５】
＜第２の実施形態＞
以下、本発明の第２の実施形態について説明する。本実施形態が、上述した第１の実施形態と異なる点は、Ｘ線検出器が、アナログ信号を加算した後にデジタル信号に変換する検出部１０４と信号収集部１０８とを備えた点であり、その一例として、所謂積分型のＸ線検出器を適用した例について説明する。
【０１５６】
図２４に示すように、Ｘ線検出器１１１は、検出部１０４及び信号収集部１０８を備えている。信号収集部１０８は、処理制御部４０９、加算率決定部４０７、記憶部４０６、加算部４０８、アナログ−デジタル変換器４１４を備える。検出部１０４は、上述した第１の実施形態と同様の構成であるので、ここでの詳細な説明は省略する。また、加算率決定部４０７と記憶部４０６も第１の実施形態と同様であるのでここでの詳細な説明は省略する。
【０１５７】
加算部４０８は、複数のスイッチ４１２、複数の積分器１４３、及びスイッチ制御部４１５及びアナログ−デジタル変換器４１４を備えている。
【０１５８】
このような構成により、加算率決定部４０７で決定した加算率を実現するように、スイッチ制御部４１５が検出素子４００毎にスイッチ４１２のオン／オフを決定し、オンとした検出素子４００の電荷のみを積分器４１３が積分する。この積分によって、スイッチ４１２がオンである検出素子４００の電荷が加算され、積分器４１３の出力電圧をデジタル−アナログ変換器４１４がデジタル信号に変換し、デジタル信号が演算部１０５に出力される。
【０１５９】
このように本実施形態では、スイッチ制御部４１５が、スイッチ４１２をオフすることで加算率０を、オンすることで加算率１を、それぞれ実現する。また加算率決定部４０７は、例えば実施形態１の場合と同様に、記憶部４０６に記憶された欠陥素子の位置情報を用いて、欠陥の検出素子とその対角検出素子の加算率を０（すなわち低減量１）に、他の検出素子の加算率を１に、それぞれ決定する。
【０１６０】
信号収集部１０８の構成部は、図２５に示すように、処理制御部４０９によって制御される。すなわち、まず撮影の前に、記憶部４０６に記憶された欠陥素子の位置情報を用いて、加算率決定部４０７が加算率を決定するように制御する。次に、スイッチ制御部４１５を制御して加算率を実現し、後の撮影中も保持する。更に、実際の撮影の際には、積分器４１３でアナログ信号を積分し、アナログ−デジタル変換器４１４でデジタル信号に変化するように制御する。
【０１６１】
このような構成により、アナログ信号を加算した後にデジタル変換を行うＸ線検出器１１１を搭載するＸ線ＣＴ装置において、故障などによって生じた欠陥素子によって、投影像の画素のサンプリング位置がずれて出力値がずれてしまうこと、更にそのためにアーチファクトが発生することを簡単に除去・低減することができ、正しい出力値を得ることができる。
【０１６２】
（加算処理を行う検出素子数の他の例）
本実施形態では３つの検出素子毎に１つの積分器を設けたが、これに限られず、４個以上とすることもできる。また一次元方向のみの検出素子からの出力信号の加算に限らず、２次元的な方向での加算を行う場合も有り得る。
【０１６３】
（Ｘ線検出器の他の例）
本実施形態では、Ｘ線検出器１１１が積分型の場合を記したが、これは一例であり、本発明を限定するものではない。例えば、フォトンカウンティング方式の検出器であっても構わない。このときＸ線によって発生した電荷の積分とデジタル信号への変換は１つのＸ線フォトン毎に行う。更にＸ線検出器が、エネルギー分別型の検出器であっても良い。このとき、アナログ−デジタル変換器４１４にて、エネルギー分別とデジタル変換を行えば良い。
【０１６４】
（検出素子）のサイズが異なる場合の一例）
本実施形態では、画素を構成する検出素子がすべて同一の大きさである例について説明した。１画素を構成する検出素子の大きさは、必ずしもすべて同一である必要はない。上述した第１の実施形態と同様に、欠陥素子と対角検出素子のサイズが異なる場合、対角検出素子の加算率の低減量と面積との積が、欠陥素子の加算率の低減量と面積との積と同じとなるように対角検出素子の加算率の低減量を決めれば良い。
【０１６５】
本実施形態においては加算率を、スイッチ４１４のオンオフで実現するため、０または１しか実現できない。そのため本方法の適用は、欠陥素子と対角検出素子の面積比が整数比の場合が望ましい。ただし整数比でない場合でも、対角検出素子の加算率の低減量と面積との積が、欠陥素子の加算率の低減量と面積との積に近い値となるように対角検出素子の加算率の低減量を決めれば、加算率の重心は画素の中心に近づき、画素の出力値のずれは低減できる。
【０１６６】
（装置形態の他の例）
上述した第１及び第２の実施形態の例では、Ｘ線ＣＴ装置に、Ｘ線検出器を適用した例について説明したが、Ｘ線検出器１１１を単独で用いることができる他、他のＸ線撮影装置はもちろん、その他種々の装置に適用することができる。
【０１６７】
本発明は上記した各実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々に変形して実施することが可能である。更に、上記実施形態には様々な段階が含まれており、開示される複数の構成要素における適宜な組み合わせにより、多様な発明が抽出され得る。例えば、実施形態に示される全構成要素から幾つかの構成要素が削除されても良い。その一例として、画像再構成処理を行わない装置や、Ｘ線源を有しない装置などが挙げられる。具体的には、Ｘ線画像診断装置、Ｘ線画像撮影装置、Ｘ線透視装置、マンモグラフィー、デジタルサブトラクション装置、Ｘ線検出器、放射線検出器等がある。
【０１６８】
更にＸ線や放射線に限らず、可視光や赤外線など、様々な波長のフォトンを検出する検出器や、それを搭載した検出装置、イメージング装置などである場合も有り得る。
【０１６９】
本発明によれば、複数の検出素子の出力を加算して投影像の画素の出力値を求めるＸ線検出器において、処理時間、処理回路及び補間用データ等を増大させることなく、欠陥素子によって画素のサンプリング位置がずれることで出力値がずれてしまうことを簡易に低減乃至は除去することができ、欠陥素子に対する補間の精度を向上させ、簡易にアーチファクトを抑制することができる。
【符号の説明】
【０１７０】
１０４・・・検出部、１０５・・・演算部、１０６・・・表示部、１０７・・・制御部、１０８・・・信号収集部、１１０・・・入力部、１１１・・・Ｘ線検出器、４０５・・・読み出し回路、４０６・・・記憶部、４０７・・・加算率決定部、４０８・・・加算部
【図１】
0006619258.tif 000009
【図２】
0006619258.tif 000010
【図３】
0006619258.tif 000011
【図４】
0006619258.tif 000012
【図５】
0006619258.tif 000013
【図６】
0006619258.tif 000014
【図７】
0006619258.tif 000015
【図８】
0006619258.tif 000016
【図９】
0006619258.tif 000017
【図１０】
0006619258.tif 000018
【図１１】
0006619258.tif 000019
【図１２】
0006619258.tif 000020
【図１３】
0006619258.tif 000021
【図１４】
0006619258.tif 000022
【図１５】
0006619258.tif 000023
【図１６】
0006619258.tif 000024
【図１７】
0006619258.tif 000025
【図１８】
0006619258.tif 000026
【図１９】
0006619258.tif 000027
【図２０】
0006619258.tif 000028
【図２１】
0006619258.tif 000029
【図２２】
0006619258.tif 000030
【図２３】
0006619258.tif 000031
【図２４】
0006619258.tif 000032
【図２５】
0006619258.tif 000033