(57)【要約】【課題】放射線画像検出器の放射線検出素子が異常画素や欠陥画素であるか否かを的確に判定することが可能な異常画素判定方法を提供する。【解決手段】放射線画像撮影では実写画像データＦ（ｘ，ｙ）を出力し、ダーク読取ではダーク読取値Ｄ（ｘ，ｙ）を出力する複数の放射線検出素子（ｘ，ｙ）を備える放射線画像検出器１における異常画素判定方法であって、一の放射線検出素子（ｘ，ｙ）から、放射線画像撮影で出力された実写画像データＦ（ｘ，ｙ）又は放射線画像撮影の前又は後に少なくとも１回行われたダーク読取で出力されたダーク読取値Ｄ（ｘ，ｙ）と、その放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）又はダーク読取値Ｄ（ｘ´，ｙ´）とに基づいて、一の放射線検出素子（ｘ，ｙ）が異常画素であるか否かを判定する。【選択図】図６
2010074644.tif 000002
【特許請求の範囲】
【請求項１】
放射線が照射されて行われる放射線画像撮影では実写画像データを出力し、放射線が照射されない状態で行われるダーク読取ではダーク読取値を出力する、２次元状に配置された複数の放射線検出素子を備える放射線画像検出器における異常画素判定方法であって、
前記放射線画像検出器の一の前記放射線検出素子から、放射線画像撮影で出力された前記実写画像データ又は放射線画像撮影の前又は後に少なくとも１回行われたダーク読取で出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データ又は前記ダーク読取値とに基づいて、当該一の放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定することを特徴とする異常画素判定方法。
【請求項２】
放射線画像撮影の前又は後に少なくとも１回行われた前記ダーク読取において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項１に記載の異常画素判定方法。
【請求項３】
放射線画像撮影の前又は後に少なくとも１回行われた前記ダーク読取において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値との第１差分を算出し、
過去のキャリブレーション時に複数回行われた前記ダーク読取において当該一の放射線検出素子から出力された前記ダーク読取値の時間的平均値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値の時間的平均値との第２差分を算出し、
前記第１差分と前記第２差分との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項１に記載の異常画素判定方法。
【請求項４】
被写体が存在する状態で行われた前記放射線画像撮影において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項１に記載の異常画素判定方法。
【請求項５】
被写体が存在する状態で行われた前記放射線画像撮影において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項１に記載の異常画素判定方法。
【請求項６】
被写体が存在しない状態で行われた前記放射線画像撮影において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項１に記載の異常画素判定方法。
【請求項７】
放射線が照射されて行われる放射線画像撮影では実写画像データを出力し、放射線が照射されない状態で行われるダーク読取ではダーク読取値を出力する、２次元状に配置された複数の放射線検出素子を備える放射線画像検出器における欠陥画素判定方法であって、
請求項１から請求項６のいずれか一項に記載の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に、前記信号値の異常の度合いに応じて当該一の放射線検出素子が欠陥画素であるか否かを判定することを特徴とする異常画素判定に基づく欠陥画素判定方法。
【請求項８】
前記一の放射線検出素子が前記異常画素であると判定した場合に、前記異常な信号値の大きさに応じて当該一の放射線検出素子が欠陥画素であるか否かを判定することを特徴とする請求項７に記載の異常画素判定に基づく欠陥画素判定方法。
【請求項９】
前記一の放射線検出素子が前記異常画素であると判定した場合、当該一の放射線検出素子が前記異常な信号値を出力する確率が予め設定された発生確率の閾値を越えた場合に当該一の放射線検出素子が欠陥画素であると判定することを特徴とする請求項７または請求項８に記載の異常画素判定に基づく欠陥画素判定方法。
【請求項１０】
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを記憶する記憶手段と、
請求項１から請求項６のいずれか一項に記載の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に当該一の放射線検出素子の情報を前記記憶手段に記憶させる判定手段と、
を備えることを特徴とする放射線画像検出器。
【請求項１１】
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと、前記ダーク読取値と、欠陥画素である放射線検出素子の情報を登録する欠陥画素マップと、を記憶する記憶手段と、
請求項７から請求項９のいずれか一項に記載の異常画素判定に基づく欠陥画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が欠陥画素であるか否かを判定し、欠陥画素であると判定した場合に当該一の放射線検出素子の情報を前記欠陥画素マップに登録する判定手段と、
を備えることを特徴とする放射線画像検出器。
【請求項１２】
バッテリが内蔵されており、前記バッテリから各部材に給電されることを特徴とする請求項１０又は請求項１１に記載の放射線画像検出器。
【請求項１３】
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを送信する通信手段と、
を備える放射線画像検出器と、
請求項１から請求項６のいずれか一項に記載の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に当該一の放射線検出素子の情報を記憶手段に記憶させる異常画素判定装置と、
を備えることを特徴とする異常画素判定システム。
【請求項１４】
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを送信する通信手段と、
を備える放射線画像検出器と、
請求項７から請求項９のいずれか一項に記載の異常画素判定に基づく欠陥画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が欠陥画素であるか否かを判定し、欠陥画素であると判定した場合に当該一の放射線検出素子の情報を、欠陥画素である放射線検出素子の情報を登録する欠陥画素マップに登録する欠陥画素判定装置と、
を備えることを特徴とする欠陥画素判定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、異常画素判定方法、異常画素判定に基づく欠陥画素判定方法、放射線画像検出器、異常画素判定システム及び欠陥画素判定システムに関するものである。
【背景技術】
【０００２】
通称フラットパネルディテクタ（Flat Panel Detector：ＦＰＤ）と呼ばれる固体撮像素子を２次元的に配置した放射線画像検出器には、検出素子として、ａ−Ｓｅ（アモルファスセレン）のような光導電物質を用いて放射線エネルギーを直接電荷に変換し、この電荷を２次元的に配置されたＴＦＴ（Thin Film Transistor：薄膜トランジスタ）等の信号読み出し用のスイッチ素子によって画素単位に電気信号として読み出す直接方式や、放射線エネルギーをシンチレータ等で光に変換し、この光を２次元的に配置されたフォトダイオード等の光電変換素子で電荷に変換してＴＦＴ等によって電気信号として読み出す間接方式等が良く知られている。
【０００３】
そして、いずれの方式においても、被写体を透過してきた放射線を放射線画像検出器で検出して得られた実写画像データに対してゲイン補正やオフセット補正等を行い、実写画像データを補正する必要があることが知られている。
【０００４】
一般的に、実写画像データの補正では、下記（１）式に示される通り、放射線画像検出器の各放射線検出素子（センサパネル部における座標は（ｘ，ｙ））から出力された実写画像データＦ（ｘ，ｙ）からオフセット補正値Ｏ（ｘ，ｙ）を差し引き、その差分にゲイン補正値Ｇ（ｘ，ｙ）を乗算することにより、最終的な画像データＦＯ（ｘ，ｙ）を得るようにして補正が行われる。なお、放射線画像検出器の１つの放射線検出素子から出力される実写画像データＦ（ｘ，ｙ）に基づいて生成される画像データＦＯ（ｘ，ｙ）は、撮影画像（すなわち放射線画像）における１画素分の画像データに相当するため、以下、放射線検出素子のことを画素と呼ぶ場合がある。
ＦＯ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）−Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）    …（１）
【０００５】
このように、実写画像データの補正では、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）を得ることが必要となるため、放射線画像検出器に対してキャリブレーションを定期的に行い、経時的に特性が変動し得るオフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）を更新するのが一般的である。
【０００６】
しかしながら、この方法では、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が一致していることが前提となるため、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が異なる場合は、特に温度依存性の高いオフセット補正値Ｏ（ｘ，ｙ）が適正な値からずれてしまい、最終的な画像データＦＯ（ｘ，ｙ）のＳＮ比を劣化させてしまうという問題がある。
【０００７】
この問題を解決するために、放射線画像撮影ごとに、撮影の直前や直後に放射線を照射しない状態で各放射線検出素子からの出力値（以下、ダーク読取値Ｄ（ｘ，ｙ）という。また、ダーク読取値をダーク画像と呼ぶ場合がある。）を検出して、当該放射線画像撮影におけるオフセット補正値Ｏ（ｘ，ｙ）が算出される場合もある。これは、放射線画像撮影で実写画像データＦ（ｘ，ｙ）が得られた時点における放射線検出素子の温度特性とできるだけ同じ温度条件下でオフセット補正値Ｏ（ｘ，ｙ）を得るための処理である。
【０００８】
しかし、実写画像データＦ（ｘ，ｙ）を取得する場合と同様に、ダーク読取値Ｄ（ｘ，ｙ）を取得する際にも、各種電気ノイズ、すなわち、フォトダイオードの暗電流ノイズ、ＴＦＴ過渡ノイズ、ＴＦＴサーマルノイズ、ＴＦＴリークノイズ、ＴＦＴから電荷を読み出すデータラインの寄生容量によって生じるサーマルノイズ、読み出し回路内部のアンプノイズ、Ａ／Ｄ変換によって生じる量子化ノイズなどが影響するため、仮に同一の温度条件下でダーク読取値Ｄ（ｘ，ｙ）を読み取ったとしても、ダーク読取値Ｄ（ｘ，ｙ）には、これら電気ノイズに起因する信号値のゆらぎ（ばらつき）が生じる。
【０００９】
そのため、放射線画像撮影の直前や直後にダーク読取値Ｄ（ｘ，ｙ）を読み取るとしても、読み取られたダーク読取値Ｄ（ｘ，ｙ）の値がすなわちその温度条件等の撮影条件におけるオフセット補正値Ｏ（ｘ，ｙ）の真値であるとは必ずしも言えない。そこで、放射線画像撮影の直前や直後にダーク読取を複数回行って各ダーク読取値Ｄ（ｘ，ｙ）の平均値を算出し、その平均値をオフセット補正値Ｏ（ｘ，ｙ）として採用することがしばしば行われる（例えば特許文献１〜３等参照）。
【００１０】
これは、放射線画像撮影の直前や直後に、放射線画像撮影における放射線検出素子等の温度特性と同じ温度条件下で複数回読み出されたダーク読取値Ｄ（ｘ，ｙ）の平均値を算出すれば、ダーク読取値Ｄ（ｘ，ｙ）のゆらぎが緩和、若しくは相殺されるため、その平均値は、その撮影条件下におけるオフセット補正値Ｏ（ｘ，ｙ）の真値に等しく、或いは少なくともそれに近い値となるという考えに基づくものである。そして、その平均値であるオフセット補正値Ｏ（ｘ，ｙ）を用いて実写画像データＦ（ｘ，ｙ）を補正すれば、補正後の最終的な画像データＦＯ（ｘ，ｙ）のＳＮ比を良好なものとすることができる。
【００１１】
一方、このようなオフセット補正をはじめとする補正処理等の画像処理は、従来から、放射線画像検出器等の撮影装置とは異なる画像処理プロセッサやコンソール等の処理装置で行われることも多い（例えば、特許文献４等参照）。また、近年、バッテリを内蔵し、ケーブルを介さずに無線方式により外部の処理装置等との間で実写画像データＦ（ｘ，ｙ）等の送受信を行うポータブルの放射線画像検出器が開発されている（例えば、特許文献５等参照）。
【００１２】
また、放射線画像検出器では、製造過程において、正常な画素値を出力できない欠陥画素が発生することがある。例えば、欠陥画素の例としては、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）などがある。
【００１３】
これら欠陥画素は、ＴＦＴ製造過程で発生する場合もあるし、ＴＦＴ上に光電変換素子を形成する際に発生する場合もある。放射線画像検出器では、照射した放射線量に比例して画素値が適正に変化することが求められるが、これら欠陥画素では、照射した放射線量に応じた正常な画素値の変化は期待できない。
【００１４】
診断画像に対してこのような欠陥画素が発生すると、誤診の原因になったり、診断行為の妨げになるため、欠陥画素が存在する場合は、この欠陥画素を周囲の画素値を用いて補間処理により補正を行い、最終的な画像データＦＯ（ｘ，ｙ）中に欠陥画素が視認できないように処理することが求められる。
【００１５】
このような補正処理を行うためには、予め欠陥画素の画素位置が分かっていなければならない。このような欠陥画素は、製造後の出荷検査や、定期的なキャリブレーション時に放射線照射を伴う検査を実施することで検出され、検出された欠陥画素の画素位置が欠陥画素マップに登録される（例えば、特許文献６参照）。この欠陥画素マップは放射線画像検出器毎に構築され、以後、個々の放射線画像検出器と対応付けて運用、管理される。
【００１６】
前述したオフセット補正、ゲイン補正を行う際にも、この欠陥画素マップが参照され、欠陥画素の画素値をオフセット補正やゲイン補正に誤って使用しないように管理される。
【特許文献１】米国特許第５４５２３３８号明細書
【特許文献２】米国特許第６２２２９０１号明細書
【特許文献３】米国特許第７０４１９５５号明細書
【特許文献４】特開平１１−１１３８８９号公報
【特許文献５】特開平７−１４０２５５号公報
【特許文献６】特開２００１−８１９８号公報
【発明の開示】
【発明が解決しようとする課題】
【００１７】
欠陥画素の例として、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）を上げたが、これらの欠陥画素は非常に容易に検出することができる。
【００１８】
しかしながら、上記のいずれにも属さないが、ある確率で異常な信号値を出力する画素（以下、異常画素と呼ぶ。）が存在する。これらの画素は欠陥画素マップに登録されていないため、ある確率で最終的な画像データＦＯ（ｘ，ｙ）中に異常な値を持つ画素を生成させ、診断の妨げとなる危険性があった。
【００１９】
また、これまでの欠陥画素判別方法では、欠陥画素を判別する際、放射線の照射を行った画像情報を用いて欠陥画素を検出していたため、画像生成のために使用される時間の関係で、欠陥画素を判別するために使用する画像数が数枚からせいぜい十数枚に限定されてしまうという問題があった。従って、定常的に異常値を出力する欠陥画素は検出できてもある確率で発生する異常画素を捕まえられなかったり、標準偏差などの統計値を用いて欠陥画素の判別を行う際に、精度の良い統計値を求めることが困難であった。そのため、放射線の照射回数を増やそうとすると生産性が低下し、生産コストが増加するという問題点があった。
【００２０】
また、このような欠陥画素や異常画素は、経時的に増加する傾向にある、そのため、工場出荷時に欠陥画素を登録しておけば安全ということではなく、ユーザー環境に設置後も、定期的にキャリブレーションを行い、新たな欠陥画素や異常画素が発見された場合には、その画素位置を欠陥画素マップに登録したり、異常画素として適切に処理することが好ましい。しかしながら、ユーザー先に設置後に行うキャリブレーション作業においても、定期的にユーザーに複数回の放射線照射を義務づけることになり、ユーザーの作業負荷が増大するという課題があった。
【００２１】
また、キャリブレーション処理により欠陥画素や異常画素を判別する際に、放射線画像検出器内の各素子の温度変動を考慮した欠陥画素や異常画素を判別する方法が提案されていないため、温度変化の影響を受けずに精度良く欠陥画素や異常画素を判定することが困難であった。
【００２２】
本発明は、以上のような課題を解決するためになされたものであり、放射線画像検出器の放射線検出素子が異常画素や欠陥画素であるか否かを的確に判定することが可能な異常画素判定方法、異常画素判定に基づく欠陥画素判定方法、放射線画像検出器、異常画素判定システム及び欠陥画素判定システムを提供することを目的とする。
【課題を解決するための手段】
【００２３】
前記の問題を解決するために、本発明の異常画素判定方法は、
放射線が照射されて行われる放射線画像撮影では実写画像データを出力し、放射線が照射されない状態で行われるダーク読取ではダーク読取値を出力する、２次元状に配置された複数の放射線検出素子を備える放射線画像検出器における異常画素判定方法であって、
前記放射線画像検出器の一の前記放射線検出素子から、放射線画像撮影で出力された前記実写画像データ又は放射線画像撮影の前又は後に少なくとも１回行われたダーク読取で出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データ又は前記ダーク読取値とに基づいて、当該一の放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定することを特徴とする。
【００２４】
また、本発明の異常画素判定に基づく欠陥画素判定方法は、
放射線が照射されて行われる放射線画像撮影では実写画像データを出力し、放射線が照射されない状態で行われるダーク読取ではダーク読取値を出力する、２次元状に配置された複数の放射線検出素子を備える放射線画像検出器における欠陥画素判定方法であって、
上記の本発明の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に、前記信号値の異常の度合いに応じて当該一の放射線検出素子が欠陥画素であるか否かを判定することを特徴とする。
【００２５】
また、本発明の放射線画像検出器は、
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを記憶する記憶手段と、
上記の本発明の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に当該一の放射線検出素子の情報を前記記憶手段に記憶させる判定手段と、
を備えることを特徴とする。
【００２６】
また、本発明の放射線画像検出器は、
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと、前記ダーク読取値と、欠陥画素である放射線検出素子の情報を登録する欠陥画素マップと、を記憶する記憶手段と、
上記の本発明の異常画素判定に基づく欠陥画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が欠陥画素であるか否かを判定し、欠陥画素であると判定した場合に当該一の放射線検出素子の情報を前記欠陥画素マップに登録する判定手段と、
を備えることを特徴とする。
【００２７】
また、本発明の異常画素判定システムは、
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを送信する通信手段と、
を備える放射線画像検出器と、
上記の本発明の異常画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定し、異常画素であると判定した場合に当該一の放射線検出素子の情報を記憶手段に記憶させる異常画素判定装置と、
を備えることを特徴とする。
【００２８】
また、本発明の欠陥画素判定システムは、
２次元状に配置された複数の放射線検出素子と、
放射線が照射されて行われる放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
前記実写画像データと前記ダーク読取値とを送信する通信手段と、
を備える放射線画像検出器と、
上記の本発明の異常画素判定に基づく欠陥画素判定方法により前記放射線画像検出器の一の前記放射線検出素子が欠陥画素であるか否かを判定し、欠陥画素であると判定した場合に当該一の放射線検出素子の情報を、欠陥画素である放射線検出素子の情報を登録する欠陥画素マップに登録する欠陥画素判定装置と、
を備えることを特徴とする。
【発明の効果】
【００２９】
本発明のような方式の異常画素判定方法、異常画素判定に基づく欠陥画素判定方法、放射線画像検出器、異常画素判定システム及び欠陥画素判定システムによれば、放射線画像検出器の放射線検出素子から出力されたダーク読取値や実写画像データの信号値が、それと対応付けられた複数の放射線検出素子から出力された信号値の空間的平均値等とかけ離れた異常な値である場合に、的確に当該放射線検出素子を異常画素であると判定することが可能となる。そして、異常画素であると判定した当該一の放射線検出素子の情報を記憶手段に一時的に記憶させておけば、その画素が欠陥画素として欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データの中に異常画素が発生することを有効に防止することが可能となる。
【００３０】
また、異常と判断された画素が、次の撮影に係る判定時に、正常状態へ戻ったと判断された時には、周辺画素による補間値では無く、当該画素の出力値を使用して診断用の画像を生成するので、高精細の診断用画像を提供可能となる。
【００３１】
また、上記のようにして異常画素であると判定された放射線検出素子から出力された異常な信号値の大きさに応じて欠陥画素であると判定したり、異常画素であると判定した放射線検出素子から異常な信号値を出力する確率に応じて欠陥画素であると判定するように構成することで、欠陥画素であると判定した放射線検出素子の情報を的確に欠陥画素マップに登録して、前述した定常的な欠陥画素と同程度に異常な信号値を出力する異常画素を欠陥画素として的確に判定して管理することが可能となるとともに、そのような異常画素（欠陥画素）について、周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データの中に異常画素に起因するアーチファクト（画像欠陥）が発生することを有効に防止することが可能となる。
【発明を実施するための最良の形態】
【００３２】
以下、本発明に係る異常画素判定方法、異常画素判定に基づく欠陥画素判定方法、放射線画像検出器、異常画素判定システム及び欠陥画素判定システムの実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。
【００３３】
［放射線画像検出器の基本的な構成］
まず、以下の本発明に用いられる放射線画像検出器の基本的な構成について説明する。
【００３４】
放射線画像検出器（フラットパネルディテクタ）１は、図１に示すように、内部を保護する筐体２を備えており、筐体２の放射線入射面Ｘの内側には、照射された放射線を光に変換する図示しないシンチレータ層が形成されている。シンチレータ層は、例えばＣｓＩ：ＴｌやＧｄ２Ｏ２Ｓ：Ｔｂ、ＺｎＳ：Ａｇ等の母体内に発光中心物質が付活された蛍光体を用いて形成されたものを用いることができる。
【００３５】
シンチレータ層の放射線が入射する側の面とは反対側の面側には、図２の等価回路図に示すように、放射線検出素子として、シンチレータ層から出力された光を電気信号に変換する複数のフォトダイオード１４が２次元状に配置されたセンサパネル部４が設けられている。１つのフォトダイオード１４から出力される電荷（信号値）は１つの画素を形成する。また、後で詳しく説明するように、各フォトダイオード１４にはそれぞれ信号読み出し用のスイッチ素子であるＴＦＴ１５が接続されている。
【００３６】
なお、以下では、上記のようにシンチレータ層で放射線を光に変換してフォトダイオード等の光電変換素子で検出する、いわゆる間接方式の放射線画像検出器１を用いる場合について説明するが、放射線画像検出器は、この他にも、前述したシンチレータ層を介さず検出素子で入射した放射線を直接電気信号に変換する、いわゆる直接方式の放射線画像検出器を用いることも可能であり、その場合にも本発明を適用することが可能である。
【００３７】
また、以下、これらの各方式の放射線画像検出器に用いられる検出素子を、あわせて放射線検出素子という。すなわち、放射線検出素子は、例えば本実施形態のような間接方式の放射線画像検出器１では１個のフォトダイオード１４、それに接続されたＴＦＴ１５及びシンチレータ層の当該フォトダイオード１４に対応する部分で構成され、例えば直接方式の放射線画像検出器では検出素子とそれに接続されたＴＦＴ等のスイッチ素子とで構成される。
【００３８】
放射線画像検出器１には、バッテリ２１（図２参照）が内蔵されている。また、図１に示すように、本実施形態では、放射線画像検出器１の筐体２の側面部分に設けられたバッテリ交換用の蓋部材１０には、無線通信手段であるアンテナ装置３が埋め込まれて設けられている。さらに、筐体２の側面部分には、放射線画像検出器１の電源スイッチ１１や各種の操作状況等を表示するインジケータ１２、ユーザーからの指示を入力するための入力ボタン２３等が設けられている。
【００３９】
図２に示すように、センサパネル部４の近傍には、センサパネル部４の各放射線検出素子の出力値を読み取る読取部５が設けられている。読取部５は、マイクロコンピュータ等からなる制御手段６や、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）、フラッシュメモリ等からなる記憶手段７、走査駆動回路８、読み出し回路９等で構成されている。
【００４０】
センサパネル部４に２次元状に配置された複数の放射線検出素子には、図３に示すように、それぞれセンサパネル部４における放射線検出素子の行方向の位置ｘと列方向の位置ｙとを各成分とする座標（ｘ，ｙ）が当該放射線検出素子の番号（ｘ，ｙ）として予め割り当てられている。以下、個々の放射線検出素子を特定する場合には、放射線検出素子（ｘ，ｙ）という。
【００４１】
なお、図３では、８×１６個の放射線検出素子（ｘ，ｙ）が記載されているが、これは簡略化して表現したものであり、実際にはさらに多くの放射線検出素子（ｘ，ｙ）が２次元状に配置されていてそれぞれ番号が割り当てられている。また、座標（ｘ，ｙ）（すなわち放射線検出素子の番号（ｘ，ｙ））を画素の番号（ｘ，ｙ）もしくは画素位置（ｘ．ｙ）と標記する場合もある。
【００４２】
センサパネル部４及び読取部５の構成についてさらに説明すると、図２の等価回路図に示すように、センサパネル部４の各放射線検出素子（ｘ，ｙ）の一方の電極にはそれぞれ信号読み出し用のスイッチ素子であるＴＦＴ１５のソース電極が接続されている。また、各放射線検出素子（ｘ，ｙ）の他方の電極にはバイアス線Ｌｂが接続されており、バイアス線Ｌｂはバイアス電源１６に接続されていて、バイアス電源１６から各放射線検出素子（ｘ，ｙ）にバイアス電圧が印加されるようになっている。
【００４３】
各ＴＦＴ１５のゲート電極はそれぞれ走査駆動回路８から延びる走査線Ｌｌに接続されており、各ＴＦＴ１５のドレイン電極はそれぞれ信号線Ｌｒに接続されている。各信号線Ｌｒは、それぞれ読み出し回路９内の増幅回路１７に接続されており、各増幅回路１７の出力線はそれぞれサンプルホールド回路１８を経てアナログマルチプレクサ１９に接続されている。また、アナログマルチプレクサ１９にはＡ／Ｄ変換器２０が接続されており、読み出し回路９はＡ／Ｄ変換器２０を介して制御手段６に接続されている。制御手段６には、記憶手段７が接続されている。
【００４４】
放射線画像検出器１では、図示しない被写体を撮影する放射線画像撮影において、被写体を透過した放射線がシンチレータ層に入射すると、シンチレータ層からセンサパネル部４に光が照射され、光の照射を受けた量に応じて、放射線検出素子（ｘ，ｙ）内に電荷が蓄積される。
【００４５】
そして、放射線画像撮影を終了し、放射線画像検出器１から実写画像データを電気信号として読み出す際には、走査線ＬｌからＴＦＴ１５のゲート電極に読み出し電圧を印加して各ＴＦＴ１５のゲートを開き、放射線検出素子（ｘ，ｙ）からＴＦＴ１５を介して蓄積された電荷を電気信号として信号線Ｌｒに取り出す。そして、電気信号を増幅回路１７で増幅する等して、アナログマルチプレクサ１９から順次Ａ／Ｄ変換器２０を介して制御手段６に出力する。
【００４６】
制御手段６は、放射線検出素子（ｘ，ｙ）から出力され増幅された電気信号を、前述した放射線検出素子（すなわち画素）の番号（ｘ，ｙ）と対応付け、実写画像データＦ（ｘ，ｙ）として記憶手段７に保存するようになっている。
【００４７】
ＴＦＴ１５に読み出し電圧を印加する走査線Ｌｌを順次走査して上記の読み出し処理を走査線Ｌｌごとに行うことで、センサパネル部４の全放射線検出素子（ｘ，ｙ）から電気信号をそれぞれ読み出し、各電気信号にそれぞれ画素の番号（ｘ，ｙ）を対応付け、各実写画像データＦ（ｘ，ｙ）として記憶手段７に保存するようになっている。
【００４８】
このように、放射線画像検出器１では、複数の放射線検出素子（ｘ，ｙ）が２次元状に配置されたセンサパネル部４や、制御手段６や走査駆動回路８、読み出し回路９等で構成された読取部５等で、実写画像データＦ（ｘ，ｙ）を取得する信号値取得手段が形成されている。
【００４９】
また、放射線画像検出器１の信号値取得手段では、上記のような実写画像データＦ（ｘ，ｙ）の取得のみならず、ダーク読取も行われるようになっている。
【００５０】
ダーク読取では、放射線画像検出器１の複数の放射線検出素子（ｘ，ｙ）をリセットして電荷を放出させた後、各ＴＦＴ１５のゲートを閉じて、放射線画像検出器１を放射線が照射されない状態に保つ。そして、所定時間経過後（通常は、放射線照射時の放射線検出素子（ｘ，ｙ）への電荷蓄積時間と同じ時間経過後）、走査線ＬｌからＴＦＴ１５のゲート電極に読み出し電圧を印加して各ＴＦＴ１５のゲートを開き、各放射線検出素子（ｘ，ｙ）に溜まった電荷（暗電荷）を信号線Ｌｒに取り出し、上記と同様に、出力値を増幅回路１７で増幅する等してアナログマルチプレクサ１９から順次Ａ／Ｄ変換器２０を介して制御手段６に出力する。
【００５１】
このようにして、放射線が曝射されない各放射線検出素子（ｘ，ｙ）から出力される電気信号がダーク読取値である。制御手段６は、各放射線検出素子（ｘ，ｙ）から出力された各電気信号を各画素の番号（ｘ，ｙ）と対応付けてダーク読取値Ｄ（ｘ，ｙ）として記憶手段７に保存するようになっている。なお、ＴＦＴ１５に読み出し電圧を印加する走査線Ｌｌを順次走査して、全放射線検出素子（ｘ，ｙ）からダーク読取値Ｄ（ｘ，ｙ）が読み出される。
【００５２】
なお、前述したように、放射線画像検出器１では、各放射線検出素子（ｘ，ｙ）からの出力値特性変動把握のために各放射線検出素子（ｘ，ｙ）のキャリブレーションが行われる。通常の場合、キャリブレーションでは、被写体が介在しない状態で放射線画像検出器１に放射線を照射し、各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）ごとにゲイン補正値が算出される。
【００５３】
また、キャリブレーションでは複数回ダーク読取も行われ、ダーク読取では、放射線画像検出器１に放射線を照射せずに各放射線検出素子（ｘ，ｙ）からダーク読取値Ｄ（ｘ，ｙ）が出力され、出力されたダーク読取値Ｄ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）ごとにオフセット補正値が算出される。
【００５４】
さらに、放射線画像撮影で実写画像データＦ（ｘ，ｙ）が得られた時点における放射線検出素子（ｘ，ｙ）の温度特性等の撮影条件と同じ撮影条件下でオフセット補正値を得るために、放射線画像撮影ごとに、撮影の直前又は直後に放射線を照射しない状態で各放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）を複数回検出して、その平均値を当該放射線画像撮影におけるオフセット補正値として算出する場合もあることは前述した通りである。
【００５５】
［放射線画像検出器の各放射線検出素子のオフセット補正値取得の原理］
さて、本発明に係る異常画素判定方法、異常画素判定に基づく欠陥画素判定方法、放射線画像検出器、異常画素判定システム及び欠陥画素判定システムの説明の前に、放射線画像検出器を用いた画像生成における放射線画像検出器の各放射線検出素子（各画素）のオフセット補正値取得の原理について説明する。
【００５６】
各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）を上記（１）式に従って補正するために、オフセット補正値（以下、Ｏ（ｘ，ｙ）と表す。）の真値が必要となる。しかし、前述したように、ダーク読取を行って実際に得られるデータは、前述した電気ノイズ等の影響のためにゆらぐ（変動する、もしくは誤差を持つ）ダーク読取値Ｄ（ｘ，ｙ）であり、オフセット補正値Ｏ（ｘ，ｙ）の真値を直接的には得ることはできない。
【００５７】
そこで、従来は、前述したように、放射線画像撮影の直前や直後に、時間的に連続して複数回のダーク読取を行って一の放射線検出素子（ｘ，ｙ）から出力される複数のダーク読取値Ｄ（ｘ，ｙ）の平均値を算出してゆらぎを緩和、若しくは、相殺し、その平均値をオフセット補正値Ｏ（ｘ，ｙ）とする手法を採用した。
【００５８】
具体的には、ダーク読取をＫ回行うとし、図４に示すように、ｋ回目（ｋ＝１〜Ｋ）のダーク読取で当該放射線検出素子（ｘ，ｙ）から出力されるダーク読取値をＤｋ（ｘ，ｙ）とすると、例えば図５に示すように、各ダーク読取で当該放射線検出素子（ｘ，ｙ）から出力される各ダーク読取値Ｄｋ（ｘ，ｙ）は平均値Ｄｋave（ｘ，ｙ）、標準偏差σＤｋ（ｘ，ｙ）の正規分布を形成するようにゆらぐことが知られている。これは、当該一の放射線検出素子（ｘ，ｙ）に対して繰り返しダーク読取を行うと、取得されるダーク読取値の値が、同じ画素位置であるにも関わらず読取り毎にゆらぐ（変動する、もしくは誤差を持つ）ことから、このゆらぎを「時間的ゆらぎ」と定義する。ゆらぎの大きさを表す統計的指標には、通常、ゆらぎの分布の標準偏差（ダーク読取値Ｄｋ（ｘ，ｙ）の場合は標準偏差σＤｋ（ｘ，ｙ））が使用される。
【００５９】
また、平均値Ｄｋave（ｘ，ｙ）は各ダーク読取値Ｄｋ（ｘ，ｙ）の時間軸上の平均値であるから、この平均値を「時間的平均値」と定義する。この平均値Ｄｋave（ｘ，ｙ）をこの放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）として使用すれば、各ダーク読取値Ｄｋ（ｘ，ｙ）が持つ値の時間的ゆらぎを緩和、若しくは相殺することができるので、一般的には、各ダーク読取値Ｄｋ（ｘ，ｙ）の時間的平均値をオフセット補正値Ｏ（ｘ，ｙ）として使用する場合が多い。オフセット補正値Ｏ（ｘ，ｙ）を下記（２）式に示す。
【数１】
2010074644.tif 000003
【００６０】
ダーク読取値Ｄｋ（ｘ，ｙ）は、標準偏差σＤｋ（ｘ，ｙ）でゆらぐが、これをＫ回で平均化すると、ゆらぎの大きさは、一般に１／√Ｋになる。すなわち、オフセット補正値Ｏ（ｘ，ｙ）は標準偏差 {１／√Ｋ・σＤｋ（ｘ，ｙ）} のゆらぎを持つ分布となる。従って、上記（２）式は、ダーク読取値Ｄｋ（ｘ，ｙ）そのものをオフセット補正値として使用するよりも、ゆらぎの大きさが１／√Ｋの値をオフセット補正値として使用することを示している。
【００６１】
しかしながら、このようにしてオフセット補正値Ｏ（ｘ，ｙ）を算出すると、ダーク読取やダーク読取値の送信を複数回行わなければならないため、電力を消費してしまう等の問題があることは前述した通りである。
【００６２】
また、この問題を回避するために、キャリブレーション時に複数のダーク読取を行っておき、上記（２）式を用いてあらかじめ算出されたオフセット補正値Ｏ（ｘ，ｙ）を使用すると、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が異なる場合は、オフセット補正値Ｏ（ｘ，ｙ）が適正な値からずれてしまうという問題があることも前述した通りである。
【００６３】
そこで、オフセット補正方法の実施例を通して、キャリブレーションの手法や温度補償の手法を交えて説明する。
【００６４】
この説明の中で取り上げられるキャリブレーションの手法や温度補償の手法は、本発明である異常画素判定方法や異常画素判定に基づく欠陥画素判定方法でも使用するため、その基本的考え方については、以下の説明の中で論ずるものとする。また、オフセット補正方法の説明には、放射線画像撮影の直前又は直後に行うダーク読取を少なくとも１回だけ行って、たとえダーク読み取りが１回だけであっても、精度の良いオフセット補正値Ｏ（ｘ，ｙ）を算出する手法を用いる。ただし、本発明に係る異常画素判定方法や異常画素判定に基づく欠陥画素判定方法は、以下に説明するオフセット補正方法に限定されるものではない。
【００６５】
本手法（放射線画像撮影の直前又は直後に行うダーク読取を少なくとも１回だけ行って、たとえダーク読み取りが１回だけであっても、精度の良いオフセット補正値Ｏ（ｘ，ｙ）を算出する手法）の基本的な考え方は以下の通りである。
【００６６】
まず、キャリブレーション時に、各放射線検出素子（ｘ，ｙ）に対応する各画素位置（ｘ，ｙ）に対して、温度変化のない所定の温度条件下で複数のダーク画像ｄｋ（ｘ，ｙ）（以下、上記のように放射線画像撮影の直前や直後に行われる１回のダーク読取で得られるダーク読取値Ｄ（ｘ，ｙ）と区別するためにダーク読取値ｄｋ（ｘ，ｙ）と表す。）を取得し、取得された複数のダーク画像ｄｋ（ｘ，ｙ）を用いて、オフセット補正値δ（ｘ，ｙ）（以下、最終的に得られるオフセット補正値Ｏ（ｘ，ｙ）と区別するためにオフセット補正値δ（ｘ，ｙ）と表す。）を求めておく。なお、ダーク画像ｄｋ（ｘ，ｙ）を取得している間は、各画素位置（ｘ，ｙ）の温度特性は同一であると仮定する。
【００６７】
この時、そのオフセット補正値δ（ｘ，ｙ）を求めた時の各画素位置（ｘ，ｙ）の信号値の温度変化に伴う変化量を代弁する温度補償変数を定義し、これを各画素位置（ｘ，ｙ）毎に計算して記憶しておく。
【００６８】
ここで、温度補償変数の位置づけについて説明する。温度補償変数は、各画素位置（ｘ，ｙ）における放射線検出素子（ｘ，ｙ）の温度そのものを測定することを目的としていない。知りたいのは、各放射線検出素子（ｘ，ｙ）の温度そのものではなく、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したかであるため（例えば、各放射線検出素子（ｘ，ｙ）の温度そのものが変化しても、その放射線検出素子（ｘ，ｙ）に対応する画素の信号値が変化していなければ、その放射線検出素子（ｘ，ｙ）の温度変化はなかったものと見なしても問題ない。）、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したか、その変化量を定量的に示唆する変数であることが望ましい。
【００６９】
次に、放射線画像撮影の直前や直後に、ダーク読取を１回だけ行って、ダーク読取値Ｄ（ｘ，ｙ）を取得する。この時、ダーク読取値Ｄ（ｘ，ｙ）に対しても、上述と同様の温度補償変数を各画素位置（ｘ，ｙ）毎に計算して記憶する。
【００７０】
次に、オフセット補正値δ（ｘ，ｙ）を求めた時の各画素位置（ｘ，ｙ）における温度補償変数とダーク読取値Ｄ（ｘ，ｙ）の各画素位置（ｘ，ｙ）における温度補償変数を比較し、その結果に基づいて、あらかじめ求めてあったオフセット補正値δ（ｘ，ｙ）を実際に放射線撮影された時の温度環境下におけるオフセット補正値Ｏ（ｘ，ｙ）に変換して、最終的な画像データＦＯ（ｘ，ｙ）を求める演算に使用する。
【００７１】
通常は、過去のキャリブレーション時に取得したオフセット補正値Ｏ（ｘ，ｙ）を使用すると、温度変動の影響を受けるため、放射線画像撮影の直前や直後に取得したダーク読取値Ｄ（ｘ，ｙ）そのものをオフセット補正値Ｏ（ｘ，ｙ）として使用することが多い（すなわち、Ｏ（ｘ，ｙ）＝Ｄ（ｘ，ｙ）とする）。一方、本手法では、放射線画像撮影毎に取得したダーク読取値Ｄ（ｘ，ｙ）はオフセット補正値そのものとして直接的に使用せず、放射線撮影された時の温度環境下における各画素位置の温度補償変数の算出に使用する点が従来手法と異なっている。そして、あらかじめ複数のダーク画像を用いて算出されたオフセット補正値δ（ｘ，ｙ）をそのまま使用するのではなく、上述の温度補償変数を用いて、実際に放射線画像撮影が実施された時点への温度補正を実施した後にオフセット補正値Ｏ（ｘ，ｙ）として使用するという点が、従来手法にはない点である。
【００７２】
この方法の利点は以下の通りである。
【００７３】
まず第１に、最終的な画像データＦＯ（ｘ，ｙ）を求めるのに使用されるオフセット補正値Ｏ（ｘ，ｙ）を算出するための基準となるオフセット補正値δ（ｘ，ｙ）が、温度変化がないと見なせる所定の温度条件下の複数のダーク画像ｄｋ（ｘ，ｙ）から算出されているため、電気ノイズ等の誤差成分（信号値の時間的なゆらぎ）が相殺された、真のオフセット補正値に近い値になっていることである。従って、放射線画像撮影時に取得された１回のダーク読取値Ｄ（ｘ，ｙ）そのものをオフセット補正値として使用する場合に比べて、補正誤差が少なく、ＳＮ比が良好な最終的な画像データＦＯ（ｘ，ｙ）を求めることができる。
【００７４】
第２に、各画素位置（ｘ，ｙ）毎に温度補償変数を計算することで、温度変化による画素値の変動を各画素位置（ｘ，ｙ）毎に個別に補正しているため、各画素位置（ｘ，ｙ）の信号値がどのように温度変化の影響を受けても、これを補正し、良好なＳＮ比の最終的な画像データＦＯ（ｘ，ｙ）が得られることである。そして、この画素単位の温度補正が、たった１回のダーク画像の取得で行えることである。
【００７５】
ここで重要になるのが、各画素位置（ｘ，ｙ）毎の温度補償変数に何を使用するかである。各画素位置（ｘ，ｙ）にあらかじめ温度センサを内蔵できれば良いが、この方法は現実的ではない。また、現実的な個数の温度センサを放射線画像検出器１に内蔵しておくことも考えられるが、これでは、各画素位置（ｘ，ｙ）毎に正確な温度補正を実現することができない。また、求めたいのは、放射線検出素子（ｘ，ｙ）の温度そのものではなく、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したかである。
【００７６】
上記より、温度補償変数の要件は、以下の４項目であると考えられる。
要件１）各画素単位（もしくは、画素に近い小領域単位毎に）に求められる変数であること。
要件２）温度変化に伴う信号値の変化を把握できる変数であること。
要件３）ダーク画像Ｄ（ｘ，ｙ）から求められる変数であること。
要件４）時間的なゆらぎが少ない変数であること。
【００７７】
この４項目の要件を満足する温度補償変数として、本件発明者は、ダーク読取値Ｄ（ｘ，ｙ）の注目画素（ｘ，ｙ）に対する周辺画素の画素値の平均値（もしくはメディアン値等）に注目した。
【００７８】
放射線検出素子（ｘ，ｙ）や電気回路などの各素子が温度変化の影響を受けた結果は、ダーク読取値Ｄ（ｘ，ｙ）に反映される。我々が知りたいのは、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって各画素の信号値がどのように変化したかであるので、温度補償変数にダーク読取値Ｄ（ｘ，ｙ）を用いるのは理にかなっており、かつ、要件１〜３を満足している。
【００７９】
次に選択した温度補償変数が要件４を満足するか否かについて、図６を用いて説明する。
【００８０】
図６に示すように、例えば複数の放射線検出素子（ｘ，ｙ）として当該一の放射線検出素子（ｘ，ｙ）すなわち、画素位置（ｘ，ｙ）を中心とする５×５の正方領域内に存在する複数の画素位置（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）の信号値を温度補償変数算出に使用する。なお、当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ，ｙ）の選択の仕方等については後で考察する。
【００８１】
しかし、当該一の放射線検出素子（ｘ，ｙ）を含む２５個の各放射線検出素子（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）から出力されるダーク読取値Ｄ（ｘ−２，ｙ−２）〜Ｄ（ｘ＋２，ｙ＋２）の信号値の時間的なゆらぎの分布は、全てが図５に示したような分布であるとは限らず、実際には、例えば図７に示すように、平均値μＤや標準偏差σＤが異なる分布となるのが一般的である。
【００８２】
本手法では、下記（３）式に示すように、当該１回のダーク読取でそれらの放射線検出素子（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）から出力された全ダーク読取値Ｄ（ｘ−２，ｙ−２）〜Ｄ（ｘ＋２，ｙ＋２）の平均値Ｗ（ｘ，ｙ）を算出して、これを温度補償変数として採用する。なお、この場合の平均値（温度補償変数）Ｗ（ｘ，ｙ）はダーク読取値Ｄ（ｘ，ｙ）の画素位置（ｘ，ｙ）に対する空間的平均値である。（３）式中のＮは、この場合２５（＝５×５）である。
【数２】
2010074644.tif 000004
【００８３】
ここで、この温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）が要件４）を満足するか否かについて考察する。
【００８４】
ダーク読取値Ｄ（ｘ−２，ｙ−２）〜Ｄ（ｘ＋２，ｙ＋２）は、それぞれが独立変数であり、かつそのゆらぎは時間的にも空間的にもランダムであるから、ダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの時間的平均値をμＤ（ｘ，ｙ）、時間的ゆらぎの分布の標準偏差をσＤ（ｘ，ｙ）と表し、温度補償変数Ｗ（ｘ，ｙ）の時間的ゆらぎの時間的平均値をμＷ（ｘ，ｙ）、時間的ゆらぎの分布の標準偏差をσＷ（ｘ，ｙ）と表すと、
【数３】
2010074644.tif 000005
【数４】
2010074644.tif 000006
となる（図８参照）。
【００８５】
ここで、σＤ（ｘ，ｙ）の中に極めて大きな値（異常値）を持つ欠陥画素や多くの異常画素が含まれていなければ、通常、
σＷ（ｘ，ｙ）＜σＤ（ｘ，ｙ）
の関係式が成立する。
【００８６】
例えば、σＤ（ｘ−２，ｙ−２）〜σＤ（ｘ＋２，ｙ＋２）のそれぞれの値は異なるが、平均を取るとσＤ（ｘ，ｙ）にほぼ等しいと仮定すると、上記（５）式より、
【数５】
2010074644.tif 000007
となり、σＷ（ｘ，ｙ）は、σＤ（ｘ，ｙ）の値の１／√Ｎとなる。Ｎ＝２５の場合は、１／√Ｎ＝１／５であるから、温度補償変数Ｗ（ｘ，ｙ）は、時間的ゆらぎの分布の広がり（標準偏差）が、ダーク読取値Ｄ（ｘ，ｙ）にくらべて１／５の分布となる。すなわち、図８に示すように、広がりが狭い（標準偏差が小さい）分布となり、要件４）を満足する。
【００８７】
極端な例として、画素位置（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）の内、注目画素である（ｘ，ｙ）以外の（Ｎ−１）の画素位置におけるダーク読取値時間的ゆらぎの標準偏差値が、注目画素位置におけるダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの標準偏差値σＤ（ｘ，ｙ）に対して、Ｎ倍の大きさを持ったと仮定しても、Ｎ＞２であれば、
【数６】
2010074644.tif 000008
となり、
σＷ（ｘ，ｙ）＜σＤ（ｘ，ｙ）
が成立する。
【００８８】
Ｎ＝２５の場合、
ＳＱＲＴ［｛１＋（Ｎ−１）・Ｎ｝／Ｎ２］
＝ＳＱＲＴ｛（１＋２４×２５）／（２５×２５）｝
≒０．９８１
となるから、
σＷ（ｘ，ｙ）≒０．９８１・σＤ（ｘ，ｙ）＜σＤ（ｘ，ｙ）
通常、σＤ（ｘ−２，ｙ−２）〜σＤ（ｘ＋２，ｙ＋２）のそれぞれの値は異なるが、平均を取るとσＤ（ｘ，ｙ）にほぼ等しい値を取ることから、上記方法は要件４）を満足すると言って良い。
【００８９】
温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布は、当該一の放射線検出素子（ｘ，ｙ）を中心とする正方領域内の複数の放射線検出素子１４（以下、所定の領域Ｒに属する放射線検出素子１４という意味で放射線検出素子１４Ｒと表記する。）の数が多くなればなるほど、すなわちＮの値が大きい程、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの標準偏差σＷ（ｘ，ｙ）が小さくなる。これは、すなわち、σＷ（ｘ，ｙ）が十分に小さい値を取るとき、
Ｗ（ｘ，ｙ）≒μＷ（ｘ，ｙ）    …（６）
と近似できることを意味している。
【００９０】
すなわち、放射線画像撮影の直前又は直後に行われる１回のダーク読取で放射線検出素子１４Ｒから出力されるダーク読取値Ｄ（ｘ−２，ｙ−２）〜Ｄ（ｘ＋２，ｙ＋２）（以下、所定の領域Ｒに属する放射線検出素子１４Ｒから出力されるダーク読取値という意味でダーク読取値ＤＲと表記する。）は、図７に例示したように各時間的ゆらぎの分布内でその時間的ゆらぎの平均値μＤ（ｘ−２，ｙ−２）〜μＤ（ｘ＋２，ｙ＋２）を中心としてゆらぐが、上記（３）式に従って算出される空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）は、空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の分布の時間的平均値μＷ（ｘ，ｙ）にほぼ等しい値が算出される。
【００９１】
従って、当該一の放射線検出素子（ｘ，ｙ）を含む各放射線検出素子１４Ｒから出力されるダーク読取値ＤＲの時間的ゆらぎの分布の平均値μＤや標準偏差σＤが異なるとしても（図７参照）、それらの空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出すれば、各放射線検出素子１４Ｒから出力された各ダーク読取値ＤＲの時間的ゆらぎが有効に相殺されて数値の変動が非常に小さくなり、時間的ゆらぎの分布の平均値μＷ（ｘ，ｙ）にほぼ等しい空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出することができる。
【００９２】
言い換えれば、複数のダーク読取を実施しなければ求められない時間的ゆらぎの平均値を、１回のダーク読取で得られたダーク読取値Ｄ（ｘ，ｙ）の空間的平均値Ｗ（ｘ，ｙ）で、少ない誤差をもって代用できるということである。
【００９３】
ここで、一の放射線検出素子（ｘ，ｙ）について温度補償変数Ｗ（ｘ，ｙ）を算出する際の当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）の選択の仕方について説明する。
【００９４】
放射線検出素子（ｘ，ｙ）では、通常、放射線検出素子（ｘ，ｙ）や、この信号値を読み出すための読み出し回路の素子の温度が低いと当該放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）は小さな値となり、放射線検出素子（ｘ，ｙ）の温度Ｔが上昇するにつれて、出力されるダーク読取値Ｄ（ｘ，ｙ）の値も大きくなっていくことが知られている。そのため、例えば、図７に示したような各放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μＤ（ｘ，ｙ）も、温度が高くなると、同図のグラフ上で右方向にシフトしていく。
【００９５】
また、例えば放射線画像検出器１の電源がＯＮされるとセンサパネル部４の温度が次第に上昇していくが、センサパネル部４上の各放射線検出素子（ｘ，ｙ）や読み出し回路の素子の温度変動は一様ではなく、例えば、周囲に多数の放射線検出素子（ｘ´，ｙ´）が存在するセンサパネル部４の中央部分の放射線検出素子（ｘ，ｙ）と、センサパネル部４の周縁部分の放射線検出素子（ｘ，ｙ）とでは、温度の変動の仕方が異なるし、読み出し回路を構成する読み出しＩＣが異なれば、温度の変動の仕方が異なる。また、同じ読み出しＩＣ内であっても、増幅回路１７が異なれば、温度の変動の仕方が異なる。
【００９６】
このような状況において、一の放射線検出素子（ｘ，ｙ）について空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出する際の当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）として、例えばセンサパネル部４上で当該一の放射線検出素子（ｘ，ｙ）から遠く離れた放射線検出素子（ｘ´，ｙ´）を選ぶと、その温度変動が当該一の放射線検出素子（ｘ，ｙ）の温度変動と異なることが予測される。
【００９７】
そのため、例えば、図７に示したような各放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の時間的ゆらぎの分布の温度による左右方向へのシフトの度合が、当該一の放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布のシフトの度合と異なるものとなり、グラフ上の相対的な位置がそれぞれの画素位置の温度に依存して変動する。
【００９８】
このように、温度によって放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布と他の放射線検出素子（ｘ´，ｙ´）の時間的ゆらぎの分布とのグラフ上の位置が、それぞれの画素位置の温度に依存して相対的にずれてしまうと、図９に示すように、それらの空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の分布はブロードになり、図８に示したような標準偏差σＷ（ｘ，ｙ）が小さい分布が得られなくなる。そのため、上記（３）式に従って算出される空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）が、空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の時間的平均値μＷ（ｘ，ｙ）にほぼ等しい値となるという有益な効果を活用することができなくなってしまう。
【００９９】
そこで、本手法では、放射線画像撮影の直前又は直後に行われる１回のダーク読取の結果から一の放射線検出素子（ｘ，ｙ）について空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出するための当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）として、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動する放射線検出素子（ｘ´，ｙ´）が予め選択されて、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けるようになっている。
【０１００】
その際、センサパネル部４上の全ての放射線検出素子（ｘ，ｙ）について同じ選択の仕方とされる必要はなく、放射線検出素子ごとに選択の仕方が異なってもよい。また、選択の仕方としては、例えば上記のように当該一の放射線検出素子（ｘ，ｙ）を中心とするｎ×ｎの正方領域内に存在する複数の放射線検出素子（ｘ，ｙ）Ｒを用いることが可能であるが、後述する図１６（Ａ）〜（Ｃ）や図１９（Ｃ）、（Ｄ）に示すように、正方領域は必ずしも当該一の放射線検出素子（ｘ，ｙ）が中心に存在するように設定されなくてもよい。
【０１０１】
また、正方領域のサイズは、本実施例では５×５や７×７の例を用いたが（７×７の例は後述する図２４〜図２６の説明文を参照）、このサイズに限定するものではない。１×１よりも大きいサイズであれば、２×２や３×３であっても良いし、９×９もしくはそれ以上のサイズであっても良い。また、正方領域の代わりに長方形状の領域を設定したり、不定形状の領域とすることも可能である。
【０１０２】
さらに、当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）は当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動する放射線検出素子であればよく、上記のような当該一の放射線検出素子（ｘ，ｙ）を含む空間的に連続する領域内に存在する放射線検出素子（ｘ，ｙ）Ｒでなく、センサパネル部４上に点在する放射線検出素子（ｘ，ｙ）であってもよい。なお、これらの同一の温度特性を示す検出素子群は、工場検査等において、予め把握することができる。また、施設設置後にも再設定することが可能である。
【０１０３】
次に、上記の温度補償変数として使用する空間的平均値Ｗ（ｘ，ｙ）と当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値との関係について説明する。
【０１０４】
放射線画像撮影の直前又は直後に１回だけ行われるダーク読取では、図１０に示したように、一の放射線検出素子（ｘ，ｙ）からは、時間的なゆらぎを含む１個のダーク読取値Ｄ（ｘ，ｙ）しか得られず、オフセット補正値Ｏ（ｘ，ｙ）の真値は分からない。また、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）から出力される１個ずつの各ダーク読取値Ｄ（ｘ−２，ｙ−２）〜Ｄ（ｘ＋２，ｙ＋２）も時間的なゆらぎを含むものであり、それらに基づいて直接的に当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値を得ることはできない。
【０１０５】
そこで、本手法では、以下で詳しく説明するように、過去の（当回の被写体の放射線画像撮影以前に行われた）キャリブレーション時にダーク読取を複数回行って当該一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄ（ｘ，ｙ）を得ておき、まず、一の放射線検出素子（ｘ，ｙ）から出力された複数回分のダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）を算出する。
【０１０６】
また、温度補償変数として、過去のキャリブレーション時の各回のダーク読取で、当該一の放射線検出素子（ｘ，ｙ）を中心とする前記正方領域内の複数の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄの空間的平均値（以下、上記の放射線画像撮影の直前や直後に行われる１回のダーク読取で得られる温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）と区別するためにキャリブレーション時に得られる温度補償変数としての空間的平均値をｗ（ｘ，ｙ）と表す。）を算出し、さらに各回の温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の複数回分の時間的平均値ω（ｘ，ｙ）を算出する。
【０１０７】
そして、上記の当該一の放射線検出素子（ｘ，ｙ）自身のダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）と、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）とを用いて、当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するようになっている。
【０１０８】
なお、本発明でキャリブレーションという場合、前述したような定期的に（例えば毎月）行われるゲイン補正値及びオフセット補正値の両方を更新するためのキャリブレーションだけでなく、ゲイン補正値に比べ、比較的温度変動が大きいオフセット補正値のみを対象とし、必要に応じて実行されるオフセットキャリブレーションも含まれる。
【０１０９】
すなわち、本発明では、キャリブレーションとは、放射線画像検出器１の工場出荷時、施設への導入時、或いは放射線が曝射されていない待機時に、必要に応じて、ダーク読取値Ｄ（ｘ，ｙ）を含む実写画像データＦ（ｘ，ｙ）を補正するためのデータの全部又は一部を予め取得して作成するための作業をいい、上記のように定期的に行われるキャリブレーションよりも広い概念である。
【０１１０】
従って、施設への導入後のキャリブレーションも、定期的なキャリブレーションよりも短い周期で（例えば毎日）行われる場合もある。また、１回のキャリブレーションやオフセットキャリブレーションでは、通常、複数回のダーク読取が行われる。
【０１１１】
なお、あまりに古いデータを用いて時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗ（ｘ，ｙ）、その時間的平均値ω（ｘ，ｙ）を算出するとそれらの値の信頼性が低下するため、それらの値の算出の元となるダーク読取値ｄ（ｘ，ｙ）は、現在から所定回数前までのダーク読取で得られた値や、最新の（すなわち前回の）キャリブレーションやオフセットキャリブレーションで行われた複数回のダーク読取で得られた値を用いることが好ましい。
【０１１２】
また、ダーク読取値ｄ（ｘ，ｙ）から時間的平均値δ（ｘ，ｙ）等を算出する処理は、放射線画像検出器１で行ってもよく、或いは放射線画像検出器１からコンソール３１（図２８参照）等の外部の装置にダーク読取値ｄ（ｘ，ｙ）を送信して外部装置で行うように構成することも可能である。
【０１１３】
さて、上記の過去の複数回のダーク読取で得られた当該一の放射線検出素子（ｘ，ｙ）自身の各ダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）は、下記（７）式に従って算出される。なお、下記（７）式で、Ｍは時間的平均値δ（ｘ，ｙ）等の算出に用いられるダーク読取値ｄｍ（ｘ，ｙ）の数、すなわち前記所定回数を表す。
【数７】
2010074644.tif 000009
【０１１４】
図１１に示すように、ｍ回目（ｍ＝１〜Ｍ）のダーク読取で当該一の放射線検出素子（ｘ，ｙ）からダーク読取値ｄｍ（ｘ，ｙ）が出力され、それらの値をヒストグラムにまとめて表すと、図１２に示すように、時間的平均値δ（ｘ，ｙ）を中心とし標準偏差σｄ（ｘ，ｙ）とする正規分布状に分布する。この分布は、ダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎを表している。
【０１１５】
また、過去のキャリブレーション時のｍ回目（ｍ＝１〜Ｍ）のダーク読取で、当該一の放射線検出素子（ｘ，ｙ）を中心とする前記正方領域内の複数の放射線検出素子（ｘ，ｙ）Ｒから出力される各ダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）は、上記（３）式と同様の下記（８）式に従って算出される。なお、この場合、（８）式中のＮも２５（＝５×５）である。
【数８】
2010074644.tif 000010
【０１１６】
そして、各回の空間的平均値ｗｍ（ｘ，ｙ）のＭ回分の時間的平均値ω（ｘ，ｙ）は、下記（９）式に従って算出される。
【数９】
2010074644.tif 000011
【０１１７】
図１１に示すように、ｍ回目（ｍ＝１〜Ｍ）のダーク読取でそれぞれ空間的平均値ｗｍ（ｘ，ｙ）が算出され、それらの値をヒストグラムにまとめて表すと、図１２に示すように、時間的平均値ω（ｘ，ｙ）を中心とし、標準偏差σｗ（ｘ，ｙ）の正規分布状に分布する。
【０１１８】
なお、上記（９）式に（８）式と（７）式を代入して変形すると、
【数１０】
2010074644.tif 000012
となることから分かるように、過去のキャリブレーション時のＭ回のダーク読取での各放射線検出素子（ｘ，ｙ）についての時間的平均値δ（ｘ，ｙ）が算出されれば、各回ごとにそれぞれ空間的平均値ｗｍ（ｘ，ｙ）を算出しなくても、各放射線検出素子（ｘ，ｙ）ごとの時間的平均値δ（ｘ，ｙ）を用いてダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）のＭ回分の時間的平均値ω（ｘ，ｙ）を算出することができる。
【０１１９】
ここで、過去のキャリブレーション時の複数回のダーク読取で当該一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）と、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値との関係について考察する。
【０１２０】
キャリブレーションは、一般的に、患者撮影を行わない始業準備中、又は、終業準備中に行われることが多く、キャリブレーション時のダーク読取は、通常、予め設定された同一の温度Ｔ０の条件下で行われると見做すことができる。従って、過去のキャリブレーション時のダーク読取で得られる当該一の放射線検出素子（ｘ，ｙ）に関する時間的平均値δ（ｘ，ｙ）及び温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）は、上記温度Ｔ０における値である。
【０１２１】
一方、放射線画像撮影が行われる際の放射線検出素子（ｘ，ｙ）やセンサパネル部４上の当該放射線検出素子（ｘ，ｙ）の近傍の温度Ｔは、図１３に示すように、必ずしもキャリブレーション時の温度Ｔ０ではない。
【０１２２】
特に、放射線画像検出器１では、電源ＯＮの状態でも、放射線画像撮影に使用されない場合には、電力の消費を抑えるために自動的又は手動の操作で電源消費状態がスリープモードとされるようになっているものも多く、図１３に示すように、撮影可能モードでは放射線検出素子（ｘ，ｙ）等の温度Ｔが上昇し、スリープモード（図中のＳの部分参照）ではセンサパネル部４への通電がされなくなるため、放射線検出素子（ｘ，ｙ）等の温度Ｔが低下する。
【０１２３】
また、前述したように、放射線検出素子（ｘ，ｙ）では、通常、放射線検出素子（ｘ，ｙ）等の温度Ｔが低いと出力されるダーク読取値Ｄ（ｘ，ｙ）は小さな値となり、放射線検出素子（ｘ，ｙ）等の温度Ｔが上昇するにつれて出力されるダーク読取値Ｄ（ｘ，ｙ）の値も大きくなる。
【０１２４】
そのため、図１４に示すように、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の中心値（平均値）や、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ，ｙ）Ｒから出力されるダーク読取値ＤＲの平均値である温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布も、温度Ｔによって大きくなったり小さくなったりする。
【０１２５】
すなわち、図１４中のダーク読取値Ｄ（ｘ，ｙ）の分布や、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布が、温度Ｔによって、図１４中の横軸上を左右に変動する。一般には、温度Ｔが上昇すると、分布は横軸上を右側にシフトする傾向にあり、温度Ｔが下降すると、分布は横軸上を左側にシフトする傾向にある。なお、今回の放射線画像撮影の直前や直後に行われた１回のダーク読取で一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）等のゆらぎの分布はあくまで推定されるものであるため、図１４では、破線で表示されている。
【０１２６】
しかし、前述したように、その場合でも、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子（ｘ´，ｙ´）Ｒは同じように温度変動する。そのため、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）について推定されるゆらぎの分布の図１４のグラフ上の位置と、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ，ｙ）Ｒから出力されるダーク読取値ＤＲの平均値である温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布のグラフ上の位置との相対的な位置関係（すなわち、それぞれの分布の平均値の距離）は温度が変動しても変わらないはずである。
【０１２７】
ところで、もし今回の放射線画像撮影が温度Ｔ０の環境下で行われた場合、図１４に示したダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μＤ（ｘ，ｙ）や、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μＷ（ｘ，ｙ）は、図１２に示した過去のキャリブレーション時の複数回のダーク読取で得られた複数のダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布の平均値δ（ｘ，ｙ）や、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子１４Ｒから出力される各ダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）の時間的ゆらぎの分布の平均値ω（ｘ，ｙ）にほぼ等しくなる。
【０１２８】
また、上記のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μＤ（ｘ，ｙ）と温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μＷ（ｘ，ｙ）との相対的な位置関係は上記のように、温度Ｔが変化しても変わらないため、各分布の平均値の差分、すなわちダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μＤ（ｘ，ｙ）と温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μＷ（ｘ，ｙ）との差分も、温度Ｔが変化しても変わらない。
【０１２９】
そのため、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる当該一の放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の分布の時間的平均値μＤ（ｘ，ｙ）と、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布の時間的平均値μＷ（ｘ，ｙ）との差分をε（ｘ，ｙ）と置くと（図１４参照）、
ε（ｘ，ｙ）＝μＤ（ｘ，ｙ）−μＷ（ｘ，ｙ）    …（１１）
は、温度Ｔにかかわらず、図１２に示した過去のキャリブレーション時の複数回のダーク読取で得られた複数のダーク読取値ｄｍ（ｘ，ｙ）の分布の時間的平均値δ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子１４Ｒから出力される各ダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）の分布の時間的平均値ω（ｘ，ｙ）との差分、
δ（ｘ，ｙ）−ω（ｘ，ｙ）
に等しくなる。
【０１３０】
すなわち、
ε（ｘ，ｙ）＝μＤ（ｘ，ｙ）−μＷ（ｘ，ｙ）
＝δ（ｘ，ｙ）−ω（ｘ，ｙ）    …（１２）
が成立する。
【０１３１】
（１２）式は、
μＤ（ｘ，ｙ）＝ε（ｘ，ｙ）＋μＷ（ｘ，ｙ）    …（１３）
これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ）の値が分かっていれば、後は、μＷ（ｘ，ｙ）さえ分かれば、μＤ（ｘ，ｙ）を推定することができることを表している。
【０１３２】
ここで、前述したように、（６）式から
Ｗ（ｘ，ｙ）≒μＷ（ｘ，ｙ）    …（６）
これを（１３）式に代入すると、
μＤ（ｘ，ｙ）≒ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）    …（１４）
の関係が得られる。
【０１３３】
これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ）の値が分かっていれば、後は、放射線画像撮影時に取得される１つのダーク画像から算出される温度補償変数Ｗ（ｘ，ｙ）さえ計算できれば、μＤ（ｘ，ｙ）を推定することができることを表している。
【０１３４】
そして、前述した従来の場合と同様に、ダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μＤ（ｘ，ｙ）は、まさに当該放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値とみなすことができるから、結局、上記（１４）式は、
Ｏ（ｘ，ｙ）≒ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）    …（１５）
と書き換えることができる。
【０１３５】
これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ）の値が分かっていれば、後は、放射線画像撮影時に取得される１つのダーク画像から算出される温度補償変数Ｗ（ｘ，ｙ）さえ計算できれば、Ｏ（ｘ，ｙ）、すなわちオフセット補正値の真値を推定することができることを表している。言い換えれば、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ）の値が分かっていれば、放射線画像撮影時の放射線検出素子や電気回路などの各素子がその後どのように温度変化していても、たった１回のダーク画像における各画素位置毎の温度補償変数Ｗ（ｘ，ｙ）を計算すれば、放射線画像撮影時の放射線検出素子や電気回路などの各素子の温度状態における理想的なオフセット補正値Ｏ（ｘ，ｙ）が求められることを表している。
【０１３６】
そこで、本手法では、
Ｏ（ｘ，ｙ）＝ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）    …（１６）
（ただし、ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ））
に従って放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するようになっている。
【０１３７】
これにより、前述した課題を解決するための骨格となる手法が提案できた。
【０１３８】
ところで、これまでの説明では、キャリブレーションを行った時の放射線画像検出器１内の各素子の温度は一定と見なして来たが、キヤリブレーション実施時、ダーク画像の取得数（上述のＭの値）が大きくなると、キャリブレーション中にも放射線画像検出器１内の各素子の温度が変動してしまう。
【０１３９】
上記（７）式〜（１０）式からも明白なように、Ｍの値が大きいほど、ｄｍ（ｘ，ｙ）、ｗｍ（ｘ，ｙ）の時間的なゆらぎが相殺され、（１６）式の精度を左右する、ω（ｘ，ｙ）やδ（ｘ，ｙ）の精度が向上するが、Ｍの値が大きいほど、ω（ｘ，ｙ）やδ（ｘ，ｙ）の値は温度変動の影響を受けやすくなることを意味している。
【０１４０】
そこで、この問題に対しても、これまで説明してきた温度補正の考え方を適応する。
【０１４１】
（１２）式に（７）式、（９）式を代入すると、
【数１１】
2010074644.tif 000013
【０１４２】
ここで、
λｍ（ｘ，ｙ）＝ｄｍ（ｘ，ｙ）−ｗｍ（ｘ，ｙ）    …（１８）
とおくと、
【数１２】
2010074644.tif 000014
【０１４３】
ここで、λｍ（ｘ，ｙ）は、キャリブレーション時のダーク読取値ｄｍ（ｘ，ｙ）の温度変化を温度補償変数ｗｍ（ｘ，ｙ）によって温度補正された、温度補正済みダーク読取値λｍ（ｘ，ｙ）と考えることができる。
【０１４４】
キャリブレーション時に放射線画像検出器１内の各素子の温度が変動してしまうと図１５に示すように、ｄｍ（ｘ，ｙ）の分布は非常にブロードな分布となり、従って、上記（７）式で求められるｄｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）は、正しい値を示せなくなる。
【０１４５】
そこで、一般的にキャリブレーション時にダーク読取を複数回実施し、その平均値を求める場合（従来は、この方法でオフセット補正値を直接求めることは前述した通りである）、直接ｄｍ（ｘ，ｙ）の平均値を計算するのではなく、上述のように、温度補正済みダーク読取値λｍ（ｘ，ｙ）の平均値（時間的平均値）であるε（ｘ，ｙ）を求めるようにすれば、温度変動が補正された正確なダーク読取の時間的平均値を求めることができるはずである。すなわち、ε（ｘ，ｙ）は温度補正済みダーク読取の時間的平均値と呼ぶことが可能である。
【０１４６】
温度補正済みダーク読取値λｍ（ｘ，ｙ）や、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）は、キャリブレーション時、さらには、放射線画像撮影時においても、ダーク読取を複数回実施し、その平均（時間的平均）を取る作業全てに適応することができる。また、時間的平均回数であるＭの値が大きくとも、温度変動の影響を受けにくくすることができる。
【０１４７】
また、キャリブレーション時にゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、同様な温度補償を適応することが可能である。ゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、あらかじめ定めておいた所定の放射線量で、被写体を置かない状態での放射線画像データを取得し、これにオフセット補正を施してゲイン補正値Ｇ（ｘ，ｙ）を計算するため、この際に取得する放射線画像データを、本手法の実写画像データＦ（ｘ，ｙ）と見なせば、本手法と同様の温度補正のかかったゲイン補正データＧ（ｘ，ｙ）を得ることができる。また、ゲイン補正値Ｇ（ｘ，ｙ）を求める際に使用されるオフセット補正値に対しても、同様に前述の温度補正を行うことが好ましい。これにより、精度の高いゲイン補正値Ｇ（ｘ，ｙ）を求めることができ、結果として、ＳＮ比の高い良好な画像データを得ることが可能となる。
【０１４８】
本手法において、キャリブレーション時のＭ回のダーク読取で放射線検出素子（ｘ，ｙ）から出力された各ダーク読取値ｄｍ（ｘ，ｙ）の分布（図１２参照）においては、時間的平均値δ（ｘ，ｙ）が算出されるが、さらに、分布の標準偏差σｄ（ｘ，ｙ）を算出し、
ｄｍ（ｘ，ｙ）−δ（ｘ，ｙ）    …（２０）
ｄｍ（ｘ，ｙ）／σｄ（ｘ，ｙ）    …（２１）
或いは、
（ｄｍ（ｘ，ｙ）−δ（ｘ，ｙ））／σｄ（ｘ，ｙ）    …（２２）
等の演算を行って各ダーク読取値ｄｍ（ｘ，ｙ）を正規化し、正規化された各ダーク読取値に基づいて上記の処理を構成することも可能である。
【０１４９】
［放射線検出素子に予め対応付ける他の放射線検出素子の選択方法の例］
前述したように、一の放射線検出素子（ｘ，ｙ）について温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）を算出するために当該一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子（ｘ´，ｙ´）については、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動するものであればよく、それらの選択の仕方としては、上記のように種々の手法があり得る。ここでは、放射線画像検出器１の実際的な構成に即した選択の仕方の一例について説明する。
【０１５０】
例えば、図６では、放射線画像検出器１のセンサパネル部４上で一の放射線検出素子（ｘ，ｙ）の周囲の行方向や列方向に十分な数の他の放射線検出素子（ｘ´，ｙ´）がある場合に、当該一の放射線検出素子（ｘ，ｙ）を中心として例えば５×５の正方領域内に存在する他の放射線検出素子（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）を選択する場合を示した。しかし、当該一の放射線検出素子（ｘ，ｙ）がセンサパネル部４の周縁部分にある場合には、このように例えば５×５の正方領域の中心に当該一の放射線検出素子（ｘ，ｙ）を位置させることができない場合もある。
【０１５１】
従って、そのような場合には、例えば図１６（Ａ）〜（Ｃ）に示すように、図中で濃く着色されて表される当該一の放射線検出素子（ｘ，ｙ）の位置などを選択すれば良い。
【０１５２】
なお、この場合、キャリブレーション時のダーク読取における空間的平均値ｗｍ（ｘ，ｙ）の算出やそれらの時間的平均値ω（ｘ，ｙ）の算出においても、当該放射線検出素子（ｘ，ｙ）についての領域Ｒは同じ要領で設定される。また、図１６（Ａ）〜（Ｃ）において、格子状に表される各放射線検出素子の上方及び左方に記載されている数字は、それぞれ各放射線検出素子（ｘ，ｙ）の行方向及び列方向の位置を表す数字であり、左方の座標（ｘ，ｙ）は当該放射線検出素子（ｘ，ｙ）の座標を表すものである。
【０１５３】
一方、各放射線検出素子（ｘ，ｙ）から電気信号を読み出し増幅するための読み出し回路９（図２参照）は、図１７に示すように、通常、行方向に１２８画素や２５６画素ずつの放射線検出素子（ｘ，ｙ）ごとに１個ずつの読み出しＩＣ９１、９２、…が接続されており、それが必要な数だけ並設されて構成されている。そして、読み出し回路９側の読み出しＩＣ９１、９２、…ごとの温度特性やノイズ特性等が必ずしも同一ではないため、各放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）が読み出しＩＣ９１、９２、…ごとに異なる場合がある。
【０１５４】
そのため、図１８に示すように、放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するために設定される領域Ｒに属する各放射線検出素子から電気信号を読み出す読み出しＩＣが、隣接する読み出しＩＣ９１、９２で別々に読み出される場合、その領域Ｒに属する各放射線検出素子（ｘ，ｙ）については温度補償変数Ｗ（ｘ，ｙ）やω（ｘ，ｙ）の値が誤差を持ってしまい、オフセット補正値Ｏ（ｘ，ｙ）を良好に算出できない場合がある。
【０１５５】
そこで、このような場合には、例えば図１９（Ｃ）、（Ｄ）に示すように、読み出しＩＣ９１に接続されている放射線検出素子（ｘ，ｙ）と読み出しＩＣ９２に接続されている放射線検出素子（ｘ，ｙ）とを分けて考える。そして、読み出しＩＣ９１、９２の境界Ｌの近傍に存在する放射線検出素子（ｘ，ｙ）を、上記のセンサパネル部４の周縁部分にある場合と同様に扱うことで、一の放射線検出素子（ｘ，ｙ）に予め対応付けられる複数の放射線検出素子を同一の読み出しＩＣ９１に接続される放射線検出素子の中から選択するように構成することが可能となる。
【０１５６】
なお、図１９（Ａ）、（Ｂ）には、図６等に示した通常の仕方で一の放射線検出素子（ｘ，ｙ）に予め対応付けられる複数の放射線検出素子が選択される場合が示されている。
【０１５７】
［欠陥画素判定方法や異常画素判定方法等の特徴］
以下、本発明に係る欠陥画素判定方法や異常画素判定方法等の特徴について説明する。
【０１５８】
まず、欠陥画素は、製造後の出荷検査や、定期的なキャリブレーション時に放射線照射を伴う検査を実施することで検出され、検出された欠陥画素の画素位置が欠陥画素マップに登録し、個々の放射線画像検出器と対応付けて運用、管理されることは前述した通りである。
【０１５９】
欠陥画素の例として、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）については、容易に判別することが可能であるため、このような欠陥画素については、既に判別されて欠陥画素マップに登録されているものとする。
【０１６０】
ここで解決しようとしているのは、ある確率で異常な信号値を出力する異常な画素についてであり、これらの異常な画素を放射線を照射せずとも精度良く検出し、欠陥画素であるか否かを判定する方法、もしくは、異常な画素を欠陥画素マップに登録する方法、もしくは異常な画素を欠陥画素マップに登録せずに異常画素として適切に処理する方法、もしくは、放射線を照射した画像データ中の異常画素について対処する方法についてである。
【０１６１】
そのため、本発明の主たる方法では、キャリブレーション時や放射線画像撮影の前後に取得するダーク読取値（放射線照射を必要としない読取信号）等を使用する。また、欠陥画素や異常画素を判別する際に、放射線画像検出器内の各素子の温度変動を考慮した欠陥画素を判別方法についても説明を行う。また、放射線を照射して得られた画像データの中から異常画素を検出し適切に対処する方法についても説明を行う。
【０１６２】
まず、本発明で解決したい欠陥画素と異常画素についての説明を行う。（１）式に示したように、最終的な画像データＦＯ（ｘ，ｙ）を得るには、実写画像データＦ（ｘ，ｙ）からオフセット補正値Ｏ（ｘ，ｙ）を差し引き、その差分にゲイン補正値Ｇ（ｘ，ｙ）を乗算する。
ＦＯ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）−Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）    …（１）
【０１６３】
従って、実写画像データＦ（ｘ，ｙ）が正常な画素値を有していても、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）に異常があれば、その演算結果である最終的な画像データＦＯ（ｘ，ｙ）は異常値を示すため、これは異常な画素（本来補正すべき画素）ということになる。
【０１６４】
また、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）が正常であっても、実写画像データＦ（ｘ，ｙ）の値が異常であれば、その演算結果である最終的な画像データＦＯ（ｘ，ｙ）は異常値を示すため、これも異常な画素（本来補正すべき画素）ということになる。
【０１６５】
ある確率でオフセット補正値Ｏ（ｘ，ｙ）、ゲイン補正値Ｇ（ｘ，ｙ）、実写画像データＦ（ｘ，ｙ）で異常画素値を出力する画素位置（ｘ，ｙ）は、その異常値の強度（どれくらい異常かを示す指針）や異常値を出力する頻度に応じて分類され、許容できない、もしくは許容すべきでないと判断された異常な画素は、欠陥画素として欠陥画素マップに登録を行う。
【０１６６】
また、突発的に異常値を出力する画素であっても、その異常値の強度（どれくらい異常かを示す指針）や異常値を出力する頻度が許容されるレベルと判断される場合は、異常画素として欠陥画素と同様の処理（補正処理など）を実施するが、欠陥画素として欠陥画素マップに登録しないことも許容する。もちろん、異常画素を全て欠陥画素として欠陥画素マップに登録しても良い。
【０１６７】
ここで、ゲイン補正値Ｇ（ｘ，ｙ）については、キャリブレーション時に放射線を照射して算出するものであるが、オフセット補正値Ｏ（ｘ，ｙ）については、放射線を照射せずに算出できるパラメータである。
【０１６８】
本発明者らは、このオフセット補正値の算出に使用するダーク読取値Ｄ（ｘ，ｙ）もしくはダーク読取値ｄｍ（ｘ，ｙ）に着目し、ある確率で発生する異常画素をダーク読取値Ｄ（ｘ，ｙ）もしくはダーク読取値ｄｍ（ｘ，ｙ）から見出し、判別し、欠陥画素や異常画素として処理する（場合によっては異常画素を欠陥画素として欠陥画素マップに登録する）方法を見出した。また、放射線画像データ中の異常画素についても検出し、適切に処理する方法についても見出した。
【０１６９】
［欠陥画素判定方法］
以下、まず、欠陥画素判定方法の例について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。
【０１７０】
前述した図５等に示したように、同じ温度Ｔの環境下にあったとしても、放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）にはゆらぎ（ばらつき）が生じ、ゆらぎの分布は平均値μ（ｘ，ｙ）を中心に標準偏差σ（ｘ，ｙ）を有する正規分布状の分布となる。
【０１７１】
しかし、放射線検出素子（ｘ，ｙ）によっては、分布の平均値μ（ｘ，ｙ）が他の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の分布の平均値μ（ｘ´，ｙ´）から大きくかけ離れていたり、或いはゆらぎが非常に大きい場合、すなわち分布の標準偏差σ（ｘ，ｙ）が他の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の分布の標準偏差σ（ｘ´，ｙ´）に比べて非常に大きいような場合がある。
【０１７２】
そのため、そのような放射線検出素子（ｘ，ｙ）では、最終的に得られる画像データＦＯ（ｘ，ｙ）の安定性も悪くなり、欠陥画素と見なすことができる程の異常な値となってしまう場合がある。以下、このように、欠陥画素と見なすことができる放射線検出素子（ｘ，ｙ）を欠陥画素（ｘｓ，ｙｓ）と表す。
【０１７３】
また、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）は正常な範囲に入るものの、時折異常な値を出力する画素が存在する。このような画素については、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）のみに注目していても検出することはできない。このように、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）は正常な範囲に入るものの、時折異常な値を出力する画素も、許容範囲を越えるものは異常画素と見なす。
【０１７４】
そこで、欠陥画素判定方法では、放射線検出素子（ｘ，ｙ）がある確率で異常画素となる場合、その発生確率や異常度合いから、欠陥画素（ｘｓ，ｙｓ）に登録するか否かを判定するようになっている。
【０１７５】
前述したように、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）の複数回分（Ｍ回分）のデータが、各放射線検出素子（ｘ，ｙ）ごとに蓄積される。そして、それらのダーク読取値ｄｍ（ｘ，ｙ）をヒストグラムにまとめて表すと、図１２に示したように、各放射線検出素子（ｘ，ｙ）ごとにダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布の時間的平均値δ（ｘ，ｙ）と標準偏差σｄ（ｘ，ｙ）等の時間的統計値を算出することができる。
【０１７６】
そこで、欠陥画素判定方法では、このように、過去のキャリブレーション時に複数回行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の値そのもの、もしくはダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）や標準偏差σｄ（ｘ，ｙ）等の時間的統計値に基づいて、下記の各判定手法に従って、当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定するようになっている。
【０１７７】
［判定手法１］
まず、時間的統計値として、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｄ（ｘ，ｙ）に着目し、標準偏差σｄ（ｘ，ｙ）が予め設定された閾値σｄthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）を欠陥画素（ｘｓ，ｙｓ）として判定し、登録するように構成することが可能である。
【０１７８】
ダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｄ（ｘ，ｙ）が大きいということは、ダーク読取値ｄｍ（ｘ，ｙ）がある確率で異常な値を出力する可能性が大きいことを示しているので、予め設定された閾値σｄthよりも大きい場合に欠陥画素（ｘｓ，ｙｓ）として判定し、欠陥画素マップに登録することで本発明が課題として認識している問題を解決することができる。
【０１７９】
すなわち、過去のキャリブレーション時に複数回行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）のゆらぎの分布における標準偏差σｄ（ｘ，ｙ）が閾値σｄth以下、すなわち、
σｄ（ｘ，ｙ）≦σｄth    …（２３）
であれば、当該一の放射線検出素子（ｘ，ｙ）のダーク読取値ｄｍ（ｘ，ｙ）のゆらぎ（ばらつき）は許容範囲内であり、正常画素であると判定される。この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素としては登録されない。
【０１８０】
しかし、図２０の右側の正規分布に示すように、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｄ（ｘ，ｙ）が、図２０の左側の正規分布に示す通常の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における平均的な標準偏差σｄ（ｘ，ｙ）よりも大きく、かつ閾値σｄthを越える場合、すなわち、
σｄ（ｘ，ｙ）＞σｄth    …（２４）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）のダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎ（ばらつき）は許容範囲を越えるものであり、欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０１８１】
上記の閾値σｄthは、例えば、放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における平均的な標準偏差σｄ（ｘ，ｙ）の３倍〜１０倍程度の値に設定されることが好ましいが、この値には限定されない。
【０１８２】
なお、標準偏差σｄ（ｘ，ｙ）の代わりに分散σｄ２（ｘ，ｙ）を用い、分散σｄ２（ｘ，ｙ）に閾値を設定して、上記と同様に処理するように構成することも可能である。
【０１８３】
［判定手法２］
また、時間的統計値として、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）に着目し、時間的平均値δ（ｘ，ｙ）が予め設定された閾値δth１（以下、第１閾値δth１という。）よりも大きいか、或いは前記閾値δth１よりも小さい値に設定された別の閾値δth２（以下、第２閾値δth２という。）よりも小さい場合に、当該一の放射線検出素子（ｘ，ｙ）は何からの欠陥を有しているとのと見なせるため、これを欠陥画素（ｘｓ，ｙｓ）として判定し、登録するように構成することが可能である。
【０１８４】
ダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）が大きすぎる場合、もしくは小さすぎる場合は、ダーク読取値ｄｍ（ｘ，ｙ）が今後もある確率で異常な値を出力する可能性が大きいことを示しているので、予め設定された第１閾値δth１よりも大きいか、或いは第２閾値δthよりも小さい場合に欠陥画素（ｘｓ，ｙｓ）として判定し、欠陥画素マップに登録することで本発明が課題として認識している問題を解決することができる。
【０１８５】
すなわち、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）が、他の放射線検出素子（ｘ＊，ｙ＊）と同様に、第２閾値δth２以上第１閾値δth１以下、すなわち、
δth２≦δ（ｘ，ｙ）≦δth１    …（２５）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず許容範囲内であり、正常画素であると判定される。この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素としては登録されない。
【０１８６】
しかし、図２１に示すように一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）が大きく、閾値δth１より大きい場合、すなわち、
δ（ｘ，ｙ）＞δth１    …（２６）
が成り立つ場合、或いは、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）が小さく、閾値δth２より小さい場合、すなわち、
δ（ｘ，ｙ）＜δth２    …（２７）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０１８７】
［判定手法３］
また、時間的統計値として、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の値の大きさそのものに着目し、ダーク読取値ｄｍ（ｘ，ｙ）の値が予め設定された閾値ｄｍthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）は何からの欠陥を有しているとのと見なせるため、これを欠陥画素（ｘｓ，ｙｓ）として判定し、登録するように構成することが可能である。
【０１８８】
ダーク読取値ｄｍ（ｘ，ｙ）そのものの値が大きすぎるということは、今後もある確率でダーク読取値ｄｍ（ｘ，ｙ）が異常な値を出力する可能性が大きいことを示しているので、予め設定された閾値ｄｍthよりも大きい場合に欠陥画素（ｘｓ，ｙｓ）として判定し、欠陥画素マップに登録することで本発明が課題として認識している問題を解決することができる。
【０１８９】
すなわち、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の値が、他の放射線検出素子（ｘ＊，ｙ＊）と同様に、閾値ｄｍth以下、すなわち、
ｄｍ（ｘ，ｙ）≦ｄｍth    …（２８）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず許容範囲内であり、正常画素であると判定される。この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素としては登録されない。
【０１９０】
しかし、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の値が閾値ｄｍthより大きい場合、すなわち、
ｄｍ（ｘ，ｙ）＞ｄｍth    …（２９）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０１９１】
さて、上記［判定手法１］〜［判定手法３］では、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）のみを用いて、その時間的統計値から当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定するように構成する場合を説明した。すなわち、複数回（Ｍ回）行われたダーク読取期間中、放射線検出素子（ｘ，ｙ）から出力される信号値は温度変化の影響を受けていない、もしくは、信号値の温度変化による影響は無視できるほど小さいと見なされた。
【０１９２】
しかしながら、複数回（Ｍ回）のダーク読取期間中に、放射線検出素子（ｘ，ｙ）から出力される信号値は温度変化の影響を受けている場合、この温度変化を補償する方が、精度の良い判定を行うことができる。
【０１９３】
そこで、前述した放射線画像検出器の各放射線検出素子のオフセット補正値手法におけるキャリブレーション時の温度補償の考え方を用いて、一の放射線検出素子（ｘ，ｙ）自身から出力されるダーク読取値ｄｍ（ｘ，ｙ）や、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動し、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線画像検出素子（ｘ´，ｙ´）から出力される各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）等を用いて、当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定し、登録するように構成する方法について［判定手法４］〜［判定手法６］で説明する。
【０１９４】
この温度補償の考え方を取り入れることによって、欠陥画素の判定を行う際に、温度起因の誤差を除去することができるので、欠陥画素の判定を温度の影響を受けることなく精度良く実施することが可能となる。
【０１９５】
［判定手法４］
ここで説明を行う判定手法は、上記［判定手法１］に温度補償の考え方を取り入れた手法である。まず、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において、各回のダーク読取ごとに、一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動し当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）との、下記（３０）式で表される差分ｅｍ（ｘ，ｙ）を新たに定義する。すなわち、空間的平均値ｗｍ（ｘ，ｙ）をダーク読取値ｄｍ（ｘ，ｙ）の温度補償変数として使用する。
ｅｍ（ｘ，ｙ）＝ｄｍ（ｘ，ｙ）−ｗｍ（ｘ，ｙ）    …（３０）
【０１９６】
なお、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）は、当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定するために当該一の放射線検出素子（ｘ，ｙ）自身から出力されるダーク読取値ｄｍ（ｘ，ｙ）と対比して用いられるものであるから、空間的平均値ｗｍ（ｘ，ｙ）を算出する対象である複数の放射線検出素子（ｘ´，ｙ´）からの各ダーク読取値ｄｍ（ｘ´，ｙ´）の中に、当該一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）は含めない方が好ましいという考え方と、空間的平均値ｗｍ（ｘ，ｙ）を算出する対象である複数の放射線検出素子（ｘ´，ｙ´）からの各ダーク読取値ｄｍ（ｘ´，ｙ´）の中に、当該一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）を含めても良いという考え方がある。
【０１９７】
本発明では、このいずれも許容するものであり、いずれか一方に限定するものではない。この考え方については、前述のオフセット補正値取得手法についても同様である。
【０１９８】
また、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられ、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動する他の複数の放射線検出素子（ｘ´，ｙ´）の選択の仕方は、前述した通りである。
【０１９９】
さて、このようにして、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取の各回のダーク読取ごとに算出される差分ｅｍ（ｘ，ｙ）の値をヒストグラムにまとめて表すと、当該一の放射線検出素子（ｘ，ｙ）が正常画素であれば、例えば図２２に示すように、差分ｅｍ（ｘ，ｙ）の時間的ゆらぎ（ばらつき）の分布は、平均的な標準偏差σｅ（ｘ，ｙ）を有する正規分布状の分布となる。
【０２００】
しかし、当該一の放射線検出素子（ｘ，ｙ）が欠陥画素（ｘｓ，ｙｓ）であれば、当該一の放射線検出素子（ｘｓ，ｙｓ）から出力されるダーク読取値ｄｍ（ｘｓ，ｙｓ）自体が正常画素の場合よりも格段に大きくゆらぐため、差分ｅｍ（ｘｓ，ｙｓ）の時間的ゆらぎ（ばらつき）の分布は、例えば図２３に示すように、大きな標準偏差σｅ（ｘｓ，ｙｓ）を有する正規分布状の分布となる。
【０２０１】
そこで、上記の差分ｅｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｅ（ｘ，ｙ）に着目し、標準偏差σｅ（ｘ，ｙ）が予め設定された閾値σｅthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）を欠陥画素（ｘｓ，ｙｓ）として判定し、登録するように構成することが可能である。
【０２０２】
すなわち、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において、各回のダーク読取ごとに、一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）との差分ｅｍ（ｘ，ｙ）を算出して、差分ｅｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｅ（ｘ，ｙ）を時間的統計値として算出し、標準偏差σｅ（ｘ，ｙ）が予め設定された閾値σｅthよりも大きい場合、すなわち、
σｅ（ｘ，ｙ）＞σｅth    …（３１）
が成り立つ場合に、当該一の放射線検出素子（ｘ，ｙ）は欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０２０３】
上記の閾値σｅthは、前述した閾値σｄthの場合と同様に、例えば、差分ｅｍ（ｘ，ｙ）の時間的ゆらぎの分布における平均的な標準偏差σｅ（ｘ，ｙ）の３倍〜１０倍程度の値に設定されることが好ましいが、この値には限定されない。また、標準偏差σｅ（ｘ，ｙ）の代わりに分散σｅ２（ｘ，ｙ）を用い、分散σｅ２（ｘ，ｙ）に閾値を設定して、上記と同様に処理するように構成することも可能である。
【０２０４】
なお、差分ｅｍ（ｘ，ｙ）は、（３０）式に示したように、ダーク読取値ｄｍ（ｘ，ｙ）に対してｗｍ（ｘ，ｙ）で温度補償を行った変数であるため、ダーク読取値ｄｍ（ｘ，ｙ）の温度変化による誤差が相殺されている。従って、差分ｅｍ（ｘ，ｙ）の時間的ゆらぎの標準偏差σｅ（ｘ，ｙ）は、ダーク読取値ｄｍ（ｘ，ｙ）の時間的ゆらぎの標準偏差σｄ（ｘ，ｙ）に対して、
σｅ（ｘ，ｙ）≦σｄ（ｘ，ｙ）    …（３２）
なる関係が成立している。すなわち、差分ｅｍ（ｘ，ｙ）のヒストグラムの分布の広がりは、ダーク読取値ｄｍ（ｘ，ｙ）のヒストグラムの分布の広がりに比べて同等以下になっている。これより、［判定手法４］は［判定手法１］に比べて、温度変化に対する誤差が少ない手法と言える。
【０２０５】
［判定手法５］
ここで説明を行う判定手法は、上記［判定手法２］に温度補償の考え方を取り入れた手法である。
【０２０６】
時間的統計値として、上記の差分εｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値に着目する。なお、上記の図２２及び図２３に示したいずれの場合も差分ｅｍ（ｘ，ｙ）の時間的ゆらぎの分布の時間的平均値は、前述の（１２）式
ε（ｘ，ｙ）＝δ（ｘ，ｙ）−ω（ｘ，ｙ）    …（１２）
で算出される温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）に計算上等しくなる。
【０２０７】
そのため、この場合、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）が、予め設定された閾値εth１（以下、第１閾値εth１という。）よりも大きいか、或いは前記閾値εth１よりも小さい値に設定された別の閾値εth２（以下、第２閾値εth２という。）よりも小さい場合に、当該一の放射線検出素子（ｘ，ｙ）を欠陥画素（ｘｓ，ｙｓ）として判定し、欠陥画素マップに登録するように構成する。
【０２０８】
すなわち、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において得られた差分ε（ｘ，ｙ）が、第２閾値εth２以上第１閾値εth１以下、すなわち、
εth２≦ε（ｘ，ｙ）≦εth１    …（３３）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず、許容範囲内であり正常画素であると判定される。この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素としては登録されない。
【０２０９】
しかし、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）が閾値εth１より大きい場合、すなわち、
ε（ｘ，ｙ）＞εth１    …（３４）
又は、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）が閾値εth２より小さい場合、すなわち、
ε（ｘ，ｙ）＜εth２    …（３５）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０２１０】
なお、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）は、（１２）式に示すように、ダーク読取値ｄｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）に対してω（ｘ，ｙ）で温度補償を行った変数であるため、ダーク読取値ｄｍ（ｘ，ｙ）の温度変化による誤差が相殺されている。これより、［判定手法５］は［判定手法２］に比べて、温度変化に対する誤差が少ない手法と言える。
【０２１１】
［判定手法６］
ここで説明を行う判定手法は、上記［判定手法３］に温度補償の考え方を取り入れた手法である。すなわち、判定に用いる時間的統計値として、各ダーク読取値ｄｍ（ｘ，ｙ）の変わりに、（３０）式で定義された、差分ｅｍ（ｘ，ｙ）の絶対値である｜ｅｍ（ｘ，ｙ）｜を使用する。
【０２１２】
差分ｅｍ（ｘ，ｙ）の絶対値｜ｅｍ（ｘ，ｙ）｜の値が予め設定された閾値ｅｍthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）は何からの欠陥を有しているとのと見なせるため、これを欠陥画素（ｘｓ，ｙｓ）として判定し、登録するように構成することが可能である。
【０２１３】
すなわち、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）と空間的平均値ｗｍ（ｘ，ｙ）の差分ｅｍ（ｘ，ｙ）の絶対値が、他の放射線検出素子（ｘ＊，ｙ＊）と同様に、閾値ｅｍth以下、すなわち、
｜ｅｍ（ｘ，ｙ）｜≦ｅｍth    …（３６）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず許容範囲内であり、正常画素であると判定される。この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素としては登録されない。
【０２１４】
しかし、一の放射線検出素子（ｘ，ｙ）について計算される差分ｅｍ（ｘ，ｙ）の絶対値が閾値ｅｍthより大きい場合、すなわち、
｜ｅｍ（ｘ，ｙ）｜＞ｅｍth    …（３７）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値ｄｍ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、欠陥画素であると判定される。そして、この場合、当該一の放射線検出素子（ｘ，ｙ）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録される。
【０２１５】
なお、差分ｅｍ（ｘ，ｙ）は、（３０）式に示したように、ダーク読取値ｄｍ（ｘ，ｙ）に対してｗｍ（ｘ，ｙ）で温度補償を行った変数であるため、ダーク読取値ｄｍ（ｘ，ｙ）の温度変化による誤差が相殺されている。これより、［判定手法６］は［判定手法１］に比べて、温度変化に対する誤差が少ない手法と言える。
【０２１６】
また、上記［判定手法１］〜［判定手法６］においては、予め定められた閾値を越えた画素を直ちに欠陥画素として登録するのではなく、例えば、その回数（頻度）によって、欠陥画素として登録するか否かを決定するようにしても良い。さらに、所定の間隔で閾値を複数設定しておき、その間隔ごとのヒストグラム（所定の区間に入る異常画素の個数）を算出して、その区間位置と、その区間での異常画素の発生頻度に応じて欠陥画素をとして登録するか否かを決定するようにしても良い。
【０２１７】
以上のように、本実施形態に係る欠陥画素判定方法によれば、ダーク読取値ｄｍ（ｘ，ｙ）のゆらぎの分布において、例えば、標準偏差σｄ（ｘ，ｙ）が格段に大きかったり、ゆらぎの分布における時間的平均値δ（ｘ，ｙ）が異常な値であったり、ダーク読取値ｄｍ（ｘ，ｙ）そのものが異常な値となるような放射線検出素子（ｘ，ｙ）を的確に見出して、その放射線検出素子（ｘ，ｙ）を的確に欠陥画素であると判定して欠陥画素マップに登録することが可能となる。
【０２１８】
そのため、そのような欠陥画素を予め登録しておくことで、オフセット補正値Ｏ（ｘ，ｙ）を算出したり、実写画像データＦ（ｘ、ｙ）を補正して最終的な画像データＦＯ（ｘ、ｙ）を生成する際に、そのような欠陥画素の悪影響を的確に排除することが可能となり、最終的な画像データＦＯ（ｘ，ｙ）のＳＮ比を良好なものとすることが可能となる。
【０２１９】
また、欠陥画素を予め登録しておくことで、その他の要因で欠陥画素として登録された欠陥画素群と同様に、例えば欠陥画素マップの形で一括して統一的に登録、管理することが可能となる。
【０２２０】
さらに、前述したように、放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）は温度変動によりその値が上下する。しかし、上記の［判定手法４］、［判定手法５］、［判定手法６］によれば、過去のキャリブレーション時の各回のダーク読取ごとに、一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）との差分εｍ（ｘ，ｙ）を算出し、それに基づいて当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定する。
【０２２１】
そのため、差分εｍ（ｘ，ｙ）に着目することで、ダーク読取値ｄｍ（ｘ，ｙ）の温度変動による値の変動の影響が除去された状態で、放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定することが可能となり、温度に影響されない状態で、欠陥画素をより安定的に見出すことが可能となる。
【０２２２】
なお、［判定手法１］〜［判定手法３］もしくは［判定手法４］〜［判定手法６］については、それぞれ単独で使用して欠陥画素を登録することができるが、［判定手法１］〜［判定手法３］もしくは［判定手法４］〜［判定手法６］の手法の組み合わせ、例えば、［判定手法１］と［判定手法３］もしくは［判定手法４］と［判定手法６］の組み合わせや、［判定手法１］、［判定手法２］、［判定手法３］もしくは［判定手法４］、［判定手法５］、［判定手法６］の組み合わせを用いることにより、さらに精度良く欠陥画素の登録を行うことができる。
【０２２３】
また、放射線画像撮影ごとに、撮影の直前や直後のダーク読取値Ｄ（ｘ，ｙ）（ダーク読取値が複数ある場合は複数のダーク読取値Ｄｍ（ｘ，ｙ）の平均値）をオフセット補正値Ｏ（ｘ，ｙ）として使用する場合においては、放射線画像撮影の直前や直後のダーク読取値Ｄ（ｘ，ｙ）、もしくはダーク読取値Ｄｍ（ｘ，ｙ）、もしくはダーク読取値Ｄｍ（ｘ，ｙ）の平均値に対して、上記の［判定手法３］もしくは［判定手法６］を使用して、ダーク読取値Ｄ（ｘ，ｙ）、もしくはダーク読取値Ｄｍ（ｘ，ｙ）、もしくはダーク読取値Ｄｍ（ｘ，ｙ）の平均値が異常値であるか否かを判定して、欠陥画素の判定、登録を行うことができる。
【０２２４】
［異常画素判定方法及び異常画素判定に基づく欠陥画素判定方法］
上記の［判定手法１］〜［判定手法６］においては、キャリブレーション時に複数回（Ｍ回）実施されるダーク読取や放射線画像撮影の直前又は直後に行われるダーク読取のダーク読取値を使用するため、Ｍの値を十分大きくすることで発生頻度を含めた統計的判断が可能である。そのため、欠陥画素として判定して良か否か或いは登録して良いか否かの判断ができた。
【０２２５】
しかしながら、放射線画像撮影の直前や直後のダーク読取値Ｄ（ｘ，ｙ）を使用する場合、キャリブレーション時のように必ずしも十分な数のダーク読取値Ｄ（ｘ，ｙ）を使用することができない場合があり（特に１回つもしくは２回のダーク読取で対応する場合など）、認識された異常画素値が極めて稀に発生した突発的現象なのか、欠陥画素として登録して良い程に異常値が頻度高く出てくるものであるのかの判定が難しい場合がある。
【０２２６】
そこで、本説明に係る異常画素判定方法では、このような異常画素は直ちに欠陥画素として登録を行わずに、異常画素として一時的な対応を行う場合についても説明する。これら異常画素を欠陥画素として登録するか否かについては、後で説明を行う。
【０２２７】
以下に、放射線画像撮影の直前や直後のダーク読取値Ｄ（ｘ，ｙ）に対して［判定手法６］を応用した場合の異常画素判定方法について説明を行う。
【０２２８】
［判定手法７］
まず、ダーク読取値Ｄ（ｘ，ｙ）と空間的平均値Ｗ（ｘ，ｙ）の差分Ｅ（ｘ，ｙ）を以下のように定義する。
Ｅ（ｘ，ｙ）＝Ｄ（ｘ，ｙ）−Ｗ（ｘ，ｙ）    …（３８）
【０２２９】
差分Ｅ（ｘ，ｙ）の絶対値を｜Ｅ（ｘ，ｙ）｜と表した時、差分Ｅ（ｘ，ｙ）の絶対値｜Ｅ（ｘ，ｙ）｜が予め設定された閾値Ｅthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）は何からの異常を有しているとのと見なせるため、これを異常画素として判定し、一時的に記憶しておく。
【０２３０】
すなわち、一の放射線検出素子（ｘ，ｙ）から出力され放射線画像撮影の直前や直後のダーク読取値Ｄ（ｘ，ｙ）と空間的平均値Ｗ（ｘ，ｙ）の差分Ｅ（ｘ，ｙ）の絶対値｜Ｅ（ｘ，ｙ）｜が、他の放射線検出素子（ｘ＊，ｙ＊）と同様に、閾値Ｅth以下、すなわち、
｜Ｅ（ｘ，ｙ）｜≦Ｅth    …（３９）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値Ｄ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず許容範囲内であり、正常画素であると判定される。
【０２３１】
しかし、一の放射線検出素子（ｘ，ｙ）について計算される差分Ｅ（ｘ，ｙ）の絶対値が閾値Ｅthより大きい場合、すなわち、
｜Ｅ（ｘ，ｙ）｜＞Ｅth    …（４０）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値Ｄ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、異常画素であると判定される。
【０２３２】
なお、差分Ｅ（ｘ，ｙ）は、（３８）式に示したように、ダーク読取値Ｄ（ｘ，ｙ）に対してＷ（ｘ，ｙ）で温度補償を行った変数であるため、ダーク読取値Ｄ（ｘ，ｙ）の温度変化による誤差が相殺されている。
【０２３３】
もしも異常画素と判定された場合は、その画素位置が欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なし、欠陥画素と同様に、周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことによって、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素が発生するのを防ぐことができる。
【０２３４】
上記［判定手法７］で示した方法は、さらに、以下に示す［判定手法８］に発展させることができる。
【０２３５】
［判定手法８］
まず、［判定手法７］と同様に、ダーク読取値Ｄ（ｘ，ｙ）と空間的平均値Ｗ（ｘ，ｙ）の差分Ｅ（ｘ，ｙ）を（３８）式に従って計算する。
【０２３６】
次に、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において、（１２）式に従って計算された、ダーク読取値ｄｍ（ｘ，ｙ）の分布の時間的平均値δ（ｘ，ｙ）とダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）の分布の時間的平均値ω（ｘ，ｙ）との差分ε（ｘ，ｙ）と、前記（３８）式によって計算された差分Ｅ（ｘ，ｙ）との差分Λ（ｘ，ｙ）を（４１）式に従って求める。
Λ（ｘ，ｙ）＝ε（ｘ，ｙ）−Ｅ（ｘ，ｙ）    …（４１）
【０２３７】
ここで、差分Λ（ｘ，ｙ）の絶対値を｜Λ（ｘ，ｙ）｜と表すと、この絶対値が予め設定された閾値Λthよりも大きい場合に、当該一の放射線検出素子（ｘ，ｙ）は何からの異常を有していると見なせるため、これを異常画素として判定し、一時的に記憶しておく。
【０２３８】
すなわち、差分Λ（ｘ，ｙ）の絶対値が、他の放射線検出素子（ｘ＊，ｙ＊）と同様に、閾値Λth以下、すなわち、
｜Λ（ｘ，ｙ）｜≦Λth    …（４２）
であれば、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値Ｄ（ｘ＊，ｙ＊）の大きさとさほど大きくは変わらず許容範囲内であり、正常画素であると判定される。
【０２３９】
しかし、一の放射線検出素子（ｘ，ｙ）について計算される差分Λ（ｘ，ｙ）の絶対値が閾値Λthより大きい場合、すなわち、
｜Λ（ｘ，ｙ）｜＞Λth    …（４３）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ＊，ｙ＊）から出力されるダーク読取値Ｄ（ｘ＊，ｙ＊）の大きさとの違いが大きく許容範囲を越えるものであり、異常画素であると判定される。
【０２４０】
なお、差分Λ（ｘ，ｙ）は、（４１）式に示したように、共に温度補償を行った変数である差分ε（ｘ，ｙ）と差分Ｅ（ｘ，ｙ）どうしの演算であるため、ダーク読取値Ｄ（ｘ，ｙ）の温度変化による誤差が相殺されている。
【０２４１】
もしも異常画素と判定された場合は、その画素位置が欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なし、欠陥画素と同様に、周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことによって、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素が発生するのを防ぐことができる。
【０２４２】
また、本発明に係る異常画素判定に基づく欠陥画素判定方法では、上記の［判定手法７］、［判定手法８］や後述する［判定手法９］、［判定手法１０］、［判定手法１１］で異常画素であると判定された画素（放射線検出素子）について、以下の選択肢１）〜５）の方法のいずれかの方法で取り扱われるようになっている。
【０２４３】
選択肢１）偶発的に発生した異常画素と見なし、欠陥画素としては登録しない。
選択肢２）欠陥画素（ｘｓ，ｙｓ）として、直ちに欠陥画素マップに登録する。
選択肢３）異常画素の値の大きさに応じて、欠陥画素として登録するか否かを決める。すなわち、異常画素の値が非常に大きい場合（別途定められた欠陥画素への登録を判定するための閾値を越えた場合）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録する。
選択肢４）異常画素として、欠陥画素とは別に画素位置を記憶しておき、その後、同じ画素位置の信号値がある頻度で異常画素となれば、その時点で発生確率の判定を行い、予め定められた発生確率の閾値を越えた時点で欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録する。
選択肢５）上記、選択肢３）と選択肢４）を組み合わせて実施する。
【０２４４】
上記［判定手法４］〜［判定手法８］については、温度補償変数である空間的平均値ｗｍ（ｘ，ｙ）や空間的平均値Ｗ（ｘ，ｙ）、及びこれらの温度補償変数を使用して求められた新たな変数である差分ε（ｘ，ｙ）、差分ｅｍ（ｘ，ｙ）、差分Ｅ（ｘ，ｙ）、絶対値｜Λ（ｘ，ｙ）｜などが使われている。ここで、空間的平均値ｗｍ（ｘ，ｙ）や空間的平均値Ｗ（ｘ，ｙ）の中に、欠陥画素として登録されていない突発的な異常画素が存在したならば、温度補償係数の値の精度が低下してしまう懸念がある。
【０２４５】
しかしながら、温度補償係数の場合は、ダーク読取値（Ｄ（ｘ，ｙ）やＤｍ（ｘ，ｙ）やｄｍ（ｘ，ｙ））の注目画素位置（ｘ，ｙ）に対してその周辺のＮ個の画素値の空間的平均値であるため、仮にＮ個の周辺画素の中に欠陥画素として登録されていない突発的な異常画素が存在したとしても、その影響は１／Ｎに緩和されているので大きな問題とはならないが、別途手段を設けて突発的な異常画素を検出し、温度補償変数の計算に使用する画素から除外するようにしても良い。
【０２４６】
さて、ここまでは、欠陥画素や異常画素をダーク読取値を用いて判別する方法について述べてきた。しかしながら、前述したように、良好な最終的な画像データＦＯ（ｘ，ｙ）を得るには、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）のみならず、実写画像データＦ（ｘ，ｙ）についても良好であること求められる。
【０２４７】
すなわち、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）に対して欠陥画素や異常画素の対処が適切になされていても、実写画像データＦ（ｘ，ｙ）の画素値の中に突発的な異常画素が含まれていれば、最終的な画像データＦＯ（ｘ，ｙ）にも異常画素が認識されることになるため好ましくない。従って、実写画像データＦ（ｘ，ｙ）についても、欠陥画素補正だけでなく、突発的に発生する異常画素に対して適切な処置を施すことが好ましい。
【０２４８】
しかしながら、実写画像データＦ（ｘ，ｙ）については、下記に示すような特徴があるため、これまで説明して来たダーク読取値ほど扱いが簡単ではない。
【０２４９】
１）実写画像データＦ（ｘ，ｙ）は殆どの場合、１回しか放射線撮影されない実写画像データＦ（ｘ，ｙ）として扱う必要がある。同じ被写体に対して複数回撮影される場合もあるが、その場合、複数回の放射線画像撮影の間に被写体を動かしたり被写体の向きを変えたりするため、おのおのが独立した実写画像データＦ（ｘ，ｙ）であると考えなければならない。
【０２５０】
２）実写画像データＦ（ｘ，ｙ）を取得するための放射線画像撮影とは、ゲイン補正データを取得するためのベタ画像撮影（ベタ画像撮影とは被写体を置かずに、一様な放射線を照射する撮影のこと）をのぞき、必ず被写体を置いての撮影となるため、実写画像データＦ（ｘ，ｙ）には被写体の画像が写り込んでいる。
【０２５１】
しかも、前述したように、実写画像データＦ（ｘ，ｙ）は１回しか放射線撮影されない独立したデータであるため、実写画像データＦ（ｘ，ｙ）の値が被写体情報を正確に反映しているのか否かを客観的に判断することが難しい。例えば、同じ被写体に対して同条件での放射線画像撮影が複数回行われれば（複数回の撮影の間、被写体は動いていないと仮定する）、所定の画素位置（ｘ，ｙ）について複数回の画像データのばらつきを比較することにより、統計的な見地から被写体情報をほぼ正確に反映しているか否かを判断できる。
【０２５２】
例えば、複数回の放射線画像撮影における画素位置（ｘ，ｙ）の信号値の平均値を計算し、その平均値からの各信号値のばらつきを測定した時、ばらつきが少なければ、それらの信号値は被写体情報をほぼ正確に反映していると判断できる。一方、複数回の放射線画像撮影における画素位置（ｘ，ｙ）の信号値の平均値に対して、ある信号値の値が平均値から大きく逸脱した値を有していれば、その信号値は被写体情報を正確に反映しておらず、ノイズ等の影響を受けた異常画素であると判断できる。
【０２５３】
このように、同じ被写体に対して同条件での放射線画像撮影が複数回行われれば、統計的な見地から被写体情報をほぼ正確に反映しているか否かを判断できるが、実写画像データＦ（ｘ，ｙ）は１回しか放射線撮影されない独立したデータであるため、このような統計的判断をすることができない。
【０２５４】
従って、各画素値に対して、本来どのような値であるべきか（各画素値のあるべき値、すなわち各画素値の真値）が分からない。例えば、ダーク読取値であれば、キャリブレーション時に取得した複数回のデータに対して今回のデータが正常であるか異常であるかは、上記のような統計的処理によって判断することができるが、１回しか放射線撮影されない実写画像データＦ（ｘ，ｙ）に対しては同様の判断を行う術がない。
【０２５５】
このような課題の中で、実写画像データＦ（ｘ，ｙ）中に発生した突発的な異常画素を判別し、適切な処置を施す手法について以下に説明を行う。
【０２５６】
［判定手法９］
まず、実写画像データＦ（ｘ，ｙ）の各画素値は、被写体の放射線吸収率によって様々に変化するが、それは画素サイズレベルのミクロの視点で見れば、非常になだらかに変化する信号値である（通常の画素サイズは１００〜２００ミクロンの範疇で選択されることが多く、１００ミクロン以下の選択はあっても２００ミクロンを越える画素サイズの選択は行われない場合が多い）。特に特定画素だけが大きく変化するような画素は存在しないと言って良い。また、最も画素レベルに近いサイズを有する病変としては、乳房撮影における微小石灰化があるが、これは放射線を吸収する方向に作用するため、信号値が減衰する方向で信号値が変化する。
【０２５７】
ここで問題にしたい異常画素は、殆どが信号値が大きくなるように作用するため（信号値が小さくなるように作用する異常は既に欠陥画素として登録されている）、放射線を吸収する方向に作用する病変が異常画素と混同されることはない。
【０２５８】
一方、信号値が大きくなるように作用する病変で、かつ画素レベルのサイズを持つものは存在しないため、独立した画素で周囲の画素に対して大きな値を持つものは突発的な異常画素と見なして問題ない。本発明ではこの特性を利用して、実写画像データＦ（ｘ，ｙ）の突発的な異常画素の判別を行う。
【０２５９】
さらに突発的に発生する画素は常に孤立した１画素である。そもそも突発的に発生する異常画素は、発生頻度が低く（発生頻度が高いものは、既に欠陥画素として登録されている）、かつ、１枚のフラットパネルディテクタは通常５メガ画素前後の画素数を有するため、今回問題にしている突発的な異常画素が並んで、もしくはクラスター状の塊として発生する確率は天文学的に小さい。
【０２６０】
従って、上述したように、今回課題として取り上げている突発的な異常画素は常に孤立した１画素として発生すると見なして良い。もちろん、隣り合って発生する場合や、クラスター状の塊として発生する場合を考慮しても良いことは言うまでもないが、以下の説明では、孤立した１画素として発生すると見なして説明を行う。
【０２６１】
上記とまとめると、今回課題としている実写画像データＦ（ｘ，ｙ）の突発的な異常画素は以下の性質をもつものとして考えられる。
【０２６２】
１）異常画素は必ず信号値が大きくなる（放射線が多く照射された画素のような挙動を示す）
２）異常画素は孤立した１点である。異常画素の周囲の画素は、欠陥画素として登録されている画素以外は必ず正常な画素である。
３）正常画素は、信号値が大きくなる方向にはなめらかに変化し、１画素単位で大きな画素変化を示さない。
【０２６３】
上記突発的異常画素を判別する方法の１つの実施例を、前述した図６を使って以下に示す。
【０２６４】
例えば複数の放射線検出素子（ｘ，ｙ）として当該一の放射線検出素子（ｘ，ｙ）すなわち、画素位置（ｘ，ｙ）を、実写画像データの注目画素位置とする。そして、画素位置（ｘ，ｙ）を中心とする５×５の正方領域内に存在する２５の画素位置（ｘ−２，ｙ−２）〜（ｘ＋２，ｙ＋２）を比較領域として定義し、この比較領域の信号値から、注目画素である画素位置（ｘ，ｙ）の信号値を除く２４個の信号値の平均値をＡ（ｘ，ｙ）と表すと、
【数１３】
2010074644.tif 000015
【０２６５】
ここでＬの値は、この場合２４（＝５×５−１）である。本例では、Ｌ＝２４とし、画素位置（ｘ，ｙ）中心とする５×５の正方領域を比較領域として説明を行うが、比較領域の選定方法はこれに限定するものでない。Ｌ＝８（＝３×３−１）として、画素位置（ｘ，ｙ）中心とする３×３の正方領域を比較領域としても良いし、例えば７×７もしくはそれ以上のサイズの正方領域を比較領域として採用しても本発明の効果を損なうものではない。また、比較領域として画素位置（ｘ，ｙ）囲む領域が選定されていれば、比較領域は必ずしも正方領域でなくとも良い。また、演算回数を減らすために、比較領域を画素位置（ｘ，ｙ）の上下、左右の４画素としたり、比較領域を画素位置（ｘ，ｙ）の上下もしくは左右の２画素としても良い。
【０２６６】
次に、（４５）式に従って、差分Ｖ（ｘ，ｙ）を計算する。
Ｖ（ｘ，ｙ）＝Ｆ（ｘ，ｙ）−Ａ（ｘ，ｙ）    …（４５）
【０２６７】
この差分Ｖ（ｘ，ｙ）が、予め設定された閾値Ｖth以下であれば、すなわち、
Ｖ（ｘ，ｙ）≦Ｖth    …（４６）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。
【０２６８】
一方、差分Ｖ（ｘ，ｙ）が、予め設定された閾値Ｖthよりも大きい場合、すなわち、
Ｖ（ｘ，ｙ）＞Ｖth    …（４７）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）を異常画素と判定して、（４９）式に定義される平均値Ｂ（ｘ，ｙ）で実写画像データＦ（ｘ，ｙ）の値を置き換える（（４８）式参照）。
Ｆ（ｘ，ｙ）←Ｂ（ｘ，ｙ）    …（４８）
【０２６９】
ここで、画素位置（ｘ，ｙ）を中心とする３×３の正方領域内に存在する９つの画素位置（ｘ−１，ｙ−１）〜（ｘ＋１，ｙ＋１）を計算領域として定義し、この計算領域の信号値から、注目画素である画素位置（ｘ，ｙ）の信号値を除く８個の信号値の平均値を平均値Ｂ（ｘ，ｙ）とする。
【数１４】
2010074644.tif 000016
【０２７０】
（４９）式において、Ｐの値は、この場合８（＝３×３−１）である。本例では、Ｐ＝８とし、画素位置（ｘ，ｙ）中心とする３×３の正方領域を計算領域として説明を行うが、計算較領域の選定方法はこれに限定するものでない。Ｐ＝２４（＝５×５−１）として、画素位置（ｘ，ｙ）中心とする５×５の正方領域を比較領域としても良いし、例えば７×７もしくはそれ以上のサイズの正方領域を比較領域として採用しても良い。また、計算領域として画素位置（ｘ，ｙ）囲む領域が選定されていれば、計算領域は必ずしも正方領域でなくとも良い。また、演算回数を減らすために、計算領域を画素位置（ｘ，ｙ）の上下、左右の４画素としたり、計算領域を画素位置（ｘ，ｙ）の上下もしくは左右の２画素としても良い。
【０２７１】
また、ここでは、比較領域Ａ（ｘ，ｙ）と計算領域Ｂ（ｘ，ｙ）を個別に定義したが（比較領域と計算領域は必ずしも一致している必要はない）、比較領域Ａ（ｘ，ｙ）と計算領域Ｂ（ｘ，ｙ）を同一のものとして扱っても良い。なお、比較領域や計算領域の中に欠陥画素が存在した場合は、欠陥画素に隣接する正常画素を使用して平均値Ａ（ｘ，ｙ）や平均値Ｂ（ｘ，ｙ）を計算すれば良い。
【０２７２】
このように、実写画像データＦ（ｘ，ｙ）の中に突発的に異常画素に対する対応を行えば、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素は発生することがなくなり、良好な診断画像を安定的に供給することが可能となる。なお、ここで検出された異常画素の取扱（欠陥画素として登録するか否か）については、上述の［判定手法８］の後に記述した選択肢１）〜５）の方法のいずれかを採用することが可能である。
【０２７３】
［判定手法１０］
なお、上記では、実写画像データＦ（ｘ，ｙ）中の異常画素は必ず信号値が大きくなる（放射線が多く照射された画素のような挙動を示す）として説明を行ったが、異常画素が、信号値が小さくなる方向に変化する場合も考慮しても良い。すなわち、
【０２７４】
１）異常画素は必ず周囲の画素値に対して大きくなる方向、もしくは小さくなる方向に異常値を持つ。
２）異常画素は孤立した１点である。異常画素の周囲の画素は、欠陥画素として登録されている画素以外は必ず正常な画素である。
３）正常画素は、はなめらかに変化し、１画素単位で大きな画素変化を示さない。
として判断を行っても良い。
【０２７５】
すなわち、（４５）式のＶ（ｘ、ｙ）に対する絶対値｜Ｖ（ｘ、ｙ）｜を定義し、この絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´th以下であれば、すなわち、
｜Ｖ（ｘ、ｙ）｜≦Ｖ´th    …（５０）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。
【０２７６】
一方、差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´thよりも大きい場合、すなわち、
｜Ｖ（ｘ、ｙ）｜＞Ｖ´th    …（５１）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）を異常画素と判定して、（４９）式に定義される平均値Ｂ（ｘ，ｙ）で実写画像データＦ（ｘ，ｙ）の値を置き換えるようにしても良い（（４８）式参照）。以下、［判定手法９］での説明と同様である。
【０２７７】
［判定手法１１］
また、前述したように、乳房撮影における微小石灰化の病変については病変サイズが画素サイズレベルに近いため、乳房撮影を行って周囲の画素から出力された信号値と大きく異なる実写画像データＦ（ｘ，ｙ）を出力した放射線検出素子があったとしても、それが病変部に存在する現実の微小石灰化を撮影したものなのか、単に異常値が出力されたものなのかを判別することは難しく、乳房撮影を行って得られた実写画像データＦ（ｘ，ｙ）に基づいて放射線検出素子を異常画素と見なしてよいか否かの判定を行うことは困難である。
【０２７８】
しかし、その一方で、特に乳房撮影に用いられる放射線画像検出器１については、放射線が照射された際にどの画素（放射線検出素子）がどの程度異常な信号値を出力するかを予め把握しておかないと、微小石灰化が存在するのに異常値と誤判定したり、異常値を微小石灰化と誤診してしまうという問題が生じる。
【０２７９】
そこで、放射線画像検出器１に対して、被写体が存在しない状態で、放射線画像検出器１に一様に放射線を照射して、各放射線検出素子（ｘ，ｙ）から実写画像データＦ（ｘ，ｙ）を出力させ、その被写体が存在しない状態で取得された各実写画像データＦ（ｘ，ｙ）に基づいて、上記［判定手法１０］と同様にして異常画素の判定を行うように構成することができる。
【０２８０】
すなわち、被写体が存在しない状態で放射線画像検出器１の放射線入射面Ｘ（図１参照）に対してその全面にわたってほぼ一定の線量となるように放射線を照射して、各放射線検出素子（ｘ，ｙ）からほぼ同じ信号値の実写画像データＦ（ｘ，ｙ）が出力されるべき状況とする。
【０２８１】
そして、その状態で実際に各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）について、上記（４４）式及び（４５）式に従ってＶ（ｘ，ｙ）を算出してＶ（ｘ、ｙ）に対する絶対値｜Ｖ（ｘ、ｙ）｜を算出し、この絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´th以下であれば、すなわち、
｜Ｖ（ｘ、ｙ）｜≦Ｖ´th    …（５０）
であれば、当該一の放射線検出素子は異常画素ではないと判定する。そして、実際の乳房撮影において、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。また、更に、ダーク読取値に基づく異常画素か否かの判定も行い、両方で異常画素では無いと判定された場合にのみ実写画像データをそのまま使用するとより好ましい。
【０２８２】
しかし、各放射線検出素子（ｘ，ｙ）からほぼ同じ信号値の実写画像データＦ（ｘ，ｙ）が出力されるべき状況であるにもかかわらず、差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´thよりも大きい場合、すなわち、
｜Ｖ（ｘ、ｙ）｜＞Ｖ´th    …（５１）
であれば、当該一の放射線検出素子を異常画素と判定する。このような場合には、周辺画素は正常組織を検出した出力値となっているので、当該異常画素に対し、周辺画素の値を用いた補間処理した値を使用することはできず、(1)異常画素が存在することを告知し他の放射線画像検出器を使用するよう促す、(2)当該異常画素のアドレスを表示し、異常画素の存在する領域を照射野領域外（乳房領域外）として撮影を行うよう被写体に対する放射線画像検出器の位置を変更するよう促す、等を技師等の操作者に対してコンソールを介して告知する。
【０２８３】
さて、これまでは、異常画素の判定方法、補正方法をオフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に対して個別に実施する説明を行って来たが、異常画素の判定方法、補正方法をオフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に対して個別に実施せずに、それらの演算結果である最終的な画像データＦＯ（ｘ，ｙ）に対して実施する方法もある。
【０２８４】
最終的な画像データＦＯ（ｘ，ｙ）は（１）式
ＦＯ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）−Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）    …（１）
によって計算されるため、オフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に異常画素が存在すれば、最終的な画像データＦＯ（ｘ，ｙ）も異常画素となる。
【０２８５】
異常画素は殆どの場合、オフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）の画素値が大きくなる方向に変化するため、上記（１）式の演算結果では、オフセット補正値Ｏ（ｘ，ｙ）に異常画素があれば最終的な画像データＦＯ（ｘ，ｙ）の信号値が小さくなる方向に変化し、実写画像データＦ（ｘ，ｙ）に異常画素があれば最終的な画像データＦＯ（ｘ，ｙ）の信号値が大きくなる方向に信号値が変化することになる。
【０２８６】
このため、判定方法としては、上記［判定手法１０］の実写画像データＦ（ｘ，ｙ）を最終的な画像データＦＯ（ｘ，ｙ）と置き換えて考えれば、［判定手法１０］で説明した方法と同様の方法で異常画素を判定することが可能である。
【０２８７】
［欠陥画素がある場合の空間的平均値の算出手法］
次に、一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子の中に、欠陥を有する放射線検出素子（以下、欠陥画素という。）が含まれる場合の当該一の放射線検出素子（ｘ，ｙ）の温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の算出手法について説明する。
【０２８８】
以下の説明を分かり易くするために、ここでは、図２４に示すように、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）を算出する当該一の放射線検出素子（注目画素）が放射線検出素子（４，４）であり、放射線検出素子（４，４）にそれを中心とする７×７個の放射線検出素子（１，１）〜（７，７）が対応付けられており、放射線検出素子（６，６）が欠陥画素である場合について説明する。
【０２８９】
［算出手法１］
最も単純な算出手法としては、欠陥画素（６，６）を含めて、当該一の放射線検出素子（４，４）に予め対応付けられた複数の放射線検出素子（１，１）〜（７，７）から出力されるダーク読取値Ｄ（１，１）〜Ｄ（７，７）を用いて、
【数１５】
2010074644.tif 000017
により放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出することが考えられる。しかし、この算出手法では、Ｗ（４，４）が欠陥画素（６，６）のダーク読取値Ｄ（６，６）の影響を受けてしまう可能性が残る。
【０２９０】
［算出手法２］
当該一の放射線検出素子（４，４）に予め対応付けられた複数の放射線検出素子（１，１）〜（７，７）のうち、欠陥画素（６，６）を除く放射線検出素子から出力されたダーク読取値Ｄ（１，１）〜Ｄ（７，７）（ただしＤ（６，６）を除く。）のみを用いて、
【数１６】
2010074644.tif 000018
により放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出することが可能である。
【０２９１】
［算出手法３］
欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子から出力されたダーク読取値に置き換えて空間的平均値Ｗ（４，４）の算出を行うように構成することも可能である。
【０２９２】
例えば、欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子（５，６）から出力されたダーク読取値Ｄ（５，６）を用いて、
Ｄ（６，６）←Ｄ（５，６）    …（５４）
と置き換える。そして、上記（５２）式に従って放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出する。
【０２９３】
なお、この場合、欠陥画素（６，６）から出力されるダーク読取値Ｄ（６，６）を、領域Ｒの内部で欠陥画素（６，６）に近接する放射線検出素子（５，６）等から出力されるダーク読取値Ｄで置換する代わりに、例えば図２５に示す放射線検出素子（８，６）のように、領域Ｒ外の位置に存在する放射線検出素子であって欠陥画素（６，６）に近接する放射線検出素子から出力されるダーク読取値Ｄで欠陥画素（６，６）から出力されるダーク読取値Ｄ（６，６）を置換するように構成することも可能である。
【０２９４】
［算出手法４］
欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）を、欠陥画素（６，６）の近傍の複数の放射線検出素子から出力された各ダーク読取値Ｄで補間して空間的平均値Ｗ（４，４）の算出を行うように構成することも可能である。
【０２９５】
例えば、欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子（５，５）、（６，５）、（７，５）、（５，６）、（７，６）、（５，７）、（６，７）、（７，７）から出力された各ダーク読取値Ｄ（５，５）、Ｄ（６，５）、Ｄ（７，５）、Ｄ（５，６）、Ｄ（７，６）、Ｄ（５，７）、Ｄ（６，７）、Ｄ（７，７）を用いて、例えばそれらの値の平均値を算出してＤ（６，６）の値をその平均値で補間する。そして、上記（５２）式に従って放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出する。
【０２９６】
なお、この場合、例えば欠陥画素が図２６に示すように、領域Ｒの端部に存在する放射線検出素子（７，６）であるような場合には、欠陥画素（７，６）から出力されたダーク読取値Ｄ（７，６）の値を、欠陥画素（７，６）の近傍の放射線検出素子（５，５）、（６，５）、（７，５）、（５，６）、（６，６）、（５，７）、（６，７）、（７，７）から出力された各ダーク読取値Ｄ（５，５）、Ｄ（６，５）、Ｄ（７，５）、Ｄ（５，６）、Ｄ（６，６）、Ｄ（５，７）、Ｄ（６，７）、Ｄ（７，７）を用いて補間するように構成することが可能である。
【０２９７】
また、図２６に示すように、領域Ｒ外の位置に存在する放射線検出素子であって欠陥画素（７，６）に近接する放射線検出素子（８，５）、（８，６）、（８，７）から出力されるダーク読取値Ｄ（８，５）、Ｄ（８，６）、Ｄ（８，７）をも用いて、欠陥画素（７，６）から出力されたダーク読取値Ｄ（７，６）の値を、放射線検出素子（６，５）、（７，５）、（８，５）、（６，６）、（８，６）、（６，７）、（７，７）、（８，７）から出力された各ダーク読取値Ｄ（６，５）、Ｄ（７，５）、Ｄ（８，５）、Ｄ（６，６）、Ｄ（８，６）、Ｄ（６，７）、Ｄ（７，７）、Ｄ（８，７）を用いて補間するように構成することも可能である。
【０２９８】
なお、上記の各算出手法では、領域Ｒ内に１つの欠陥画素、すなわち１つの点欠陥がある場合について説明したが、複数の点欠陥が存在する場合にも同様の算出手法を用いることができる。また、欠陥画素が線状に存在する線欠陥がある場合においても、同様の算出手法を用いることができる。
【０２９９】
［放射線画像検出器における異常画素判定や異常画素判定に基づく欠陥画素判定］
上述の［判定手法７］、［判定手法８］は、いずれもダーク読取値を使用するため、放射線の照射を伴わずに異常画素であるか否かを判定することができる。
【０３００】
すなわち、図２等に示した放射線画像検出器１において、マイクロコンピュータ等からなる制御手段６を判定手段とし、例えば上記の［判定手法７］に従って、放射線画像撮影の前又は後に少なくとも１回行われた、オフセット補正値生成のためのダーク読取で放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力されたダーク読取値Ｄ（ｘ´，ｙ´）の空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）との差分Ｅ（ｘ，ｙ）を算出し（上記（３８）式参照）、その絶対値｜Ｅ（ｘ，ｙ）｜が予め設定された閾値Ｅthを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定するように構成することが可能である。
【０３０１】
また、判定手段が、例えば上記の［判定手法８］に従って、［判定手法７］と同様に一の放射線検出素子（ｘ，ｙ）についてダーク読取値Ｄ（ｘ，ｙ）と空間的平均値Ｗ（ｘ，ｙ）の差分Ｅ（ｘ，ｙ）を（３８）式に従って第１差分として計算し、また、過去のキャリブレーション時に複数回行われたダーク読取において（１２）式に従って計算されたダーク読取値ｄｍ（ｘ，ｙ）の分布の時間的平均値δ（ｘ，ｙ）とダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）の分布の時間的平均値ω（ｘ，ｙ）との差分ε（ｘ，ｙ）を第２差分として計算し、これらの第１差分Ｅ（ｘ，ｙ）と第２差分ε（ｘ，ｙ）との差分Λ（ｘ，ｙ）（上記（４１）式参照）の絶対値｜Λ（ｘ，ｙ）｜が予め設定された閾値Λthを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定するように構成することも可能である。
【０３０２】
このように構成すれば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）が、それと対応付けられた複数の放射線検出素子から出力されたダーク読取値Ｄ（ｘ´，ｙ´）の空間的平均値Ｗ（ｘ，ｙ）等とかけ離れた異常な値である場合に、的確に当該放射線検出素子（ｘ，ｙ）を異常画素であると判定することが可能となる。
【０３０３】
そして、異常画素であると判定した当該一の放射線検出素子（ｘ，ｙ）の情報を記憶手段７に一時的に記憶させておけば、その画素が欠陥画素として欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素に起因するアーチファクト（画像欠陥）が発生することを有効に防止することが可能となる。
【０３０４】
尚、次回の異常判定時、前回異常と判定された画素が正常にもどっていれば、当該画素の出力値を使用して画像を生成するので、高精細の画像を提供可能となる。
【０３０５】
また、上述の［判定手法９］〜［判定手法１１］では、いずれも放射線を照射して得られた実写画像データを使用して異常画素であるか否かを判定することができる。
【０３０６】
すなわち、図２等に示した放射線画像検出器１において、マイクロコンピュータ等からなる制御手段６を判定手段とし、例えば上記の［判定手法９］や［判定手法１０］に従って、被写体が存在する状態で行われた放射線が照射された放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）との差分Ｖ（ｘ，ｙ）（上記（４５）式参照）が、予め設定された閾値Ｖthを越えた場合、或いはその絶対値｜Ｖ（ｘ，ｙ）｜が予め設定された閾値Ｖ´thを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定するように構成することが可能である。
【０３０７】
また、判定手段が、例えば上記の［判定手法１１］に従って、被写体が存在しない状態で行われた放射線が照射された放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）との差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ，ｙ）｜が予め設定された閾値Ｖ´thを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定するように構成することも可能である。
【０３０８】
このように構成すれば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）が、それと対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）とかけ離れた突発的に異常な値である場合に、的確に当該放射線検出素子（ｘ，ｙ）を異常画素であると判定することが可能となる。
【０３０９】
そして、異常画素であると判定した当該一の放射線検出素子（ｘ，ｙ）の情報を記憶手段７に一時的に記憶させておけば、その画素が欠陥画素として欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素が発生することを有効に防止することが可能となる。
【０３１０】
さらに、上記のようにして異常画素であると判定された放射線検出素子（ｘ，ｙ）について上記の選択肢１）〜５）の方法のいずれかの方法で取り扱うことで、放射線画像検出器１において異常画素判定に基づく欠陥画素判定方法を実現することができる。
【０３１１】
具体的には、図２等に示した放射線画像検出器１において、マイクロコンピュータ等からなる制御手段６を判定手段とし、判定手段が、異常画素であると判定した放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）である異常な信号値の大きさに応じて、その信号値の大きさが予め設定された閾値を越える場合に欠陥画素であると判定したり、異常画素であると判定した放射線検出素子（ｘ，ｙ）から異常な信号値を出力する確率が予め設定された発生確率の閾値を越えた場合に欠陥画素であると判定するように構成することができる。
【０３１２】
そして、欠陥画素であると判定した放射線検出素子の情報を欠陥画素マップに登録するように構成することで、前述した定常的な欠陥画素と同程度に異常な信号値を出力する異常画素を欠陥画素として的確に判定して管理することが可能となるとともに、そのような異常画素（欠陥画素）について、周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素に起因するアーチファクト（画像欠陥）が発生することを有効に防止することが可能となる。
【０３１３】
［放射線画像検出器における画像生成］
ここで、放射線画像検出器の各放射線検出素子のオフセット補正値取得の原理を用いて放射線画像検出器で最終的な画像データを生成する場合について説明する。本実施形態では、前述したように、放射線撮影された画像データのオフセット補正に用いるオフセット補正値を得るためのダーク読取は、当該放射線画像撮影の前又は後に少なくとも１回行われるようになっており、そのダーク読取によって取得されたダーク読取値に基づいて各放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するようになっている。
【０３１４】
なお、以下では、放射線画像撮影の後にダーク読取を１回だけ行う場合について説明するが、ダーク読取は当該放射線画像撮影の前に行ってもよく、また、ダーク読取を２回程度行うように構成することも可能である。なお、ダーク読取を２回以上行っても、消費電力量は増加するが、本実施形態のその他の効果が得られることは言うまでもない。
【０３１５】
放射線画像検出器１を用いた画像生成は、図２７に示すフローチャートに従って行われるようになっており、欠陥画素判定ステップ（ステップＳ１）と、実写画像データ取得ステップ（ステップＳ２）と、ダーク読取ステップ（ステップＳ３）と、温度補償変数算出ステップ（ステップＳ４）と、時間的平均値等算出ステップ（ステップＳ５）と、オフセット補正値算出ステップ（ステップＳ６）と、画像補正ステップ（ステップＳ７）とを経て、最終的な画像データが得られるようになっている。なお、本実施形態で、統計値とは前述した空間的平均値や時間的平均値等のことである。しかし、他の統計値のとり方を採用することも可能であり、この点については後で述べる。
【０３１６】
欠陥画素判定ステップ（ステップＳ１）では、前述した欠陥画素判定方法に基づいて、過去のキャリブレーション時に複数回行われたダーク読取の結果や放射線画像撮影の直前又は直後に１回以上行われたダーク読取の結果から、放射線画像検出器１のセンサパネル部４の２次元状に配置された複数の放射線検出素子（ｘ，ｙ）の全てについて欠陥画素であるか否かが判定され、欠陥画素であると判定された放射線検出素子の番号（ｘ，ｙ）等の情報が放射線画像検出器１の記憶手段７に保存された欠陥画素マップに登録される。
【０３１７】
また、前述した本発明に係る異常画素判定方法に基づいて異常画素と判定され、その画素位置を欠陥画素マップに登録しない場合でも、一時的に欠陥画素同等と見なして欠陥画素と同様に補正処理等を施す場合には、その番号（ｘ，ｙ）等の情報が放射線画像検出器１の記憶手段７に一時的に保存され、以下の処理では欠陥画素と同等に扱われる。
【０３１８】
なお、前述した本発明に係る異常画素判定に基づく欠陥画素判定方法で欠陥画素と判定された異常画素が、欠陥画素として欠陥画素マップに登録されることは言うまでもない。また、欠陥画素判定ステップ（ステップＳ１）をコンソール３２等の外部装置で行う場合には、有線方式や無線方式によって、キャリブレーション時のダーク読取で各放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）が外部装置に送信される。そして、欠陥画素であると判定された放射線検出素子の番号（ｘ，ｙ）等の情報がコンソール３２の図示しない記憶手段やサーバ手段３９の記憶手段３８（図２８参照）に保存されて登録される。
【０３１９】
また、欠陥画素判定ステップ（ステップＳ１）を、実写画像データ取得ステップ（ステップＳ２）やダーク読取ステップ（ステップＳ３）等の後に行うように構成することも可能である。
【０３２０】
実写画像データ取得ステップ（ステップＳ２）では、放射線画像検出器１の放射線入射面Ｘ（図１参照）側に図示しない患者の手や胸等の被写体部位を載置又は接近させた状態で、図示しない放射線発生装置から放射線を照射する通常の放射線照射が行われる。
【０３２１】
前述したように、放射線画像検出器１には、検出手段にａ−Ｓｅのような光導電物質を用いて放射線エネルギーを直接電荷に変換し、この電荷を２次元的に配置されたＴＦＴ等の信号読み出し用のスイッチ素子によって画素単位に電気信号として読み出す直接方式や、放射線エネルギーをシンチレータ等で光に変換し、この光を２次元的に配置されたフォトダイオード等の光電変換素子で電荷に変換してＴＦＴ等によって電気信号として読み出す間接方式などが良く知られている。
【０３２２】
その際、放射線画像検出器１が間接方式であれば、被写体を透過した放射線がシンチレータ層３で光に変換されて、画素単位に配置された複数の放射線検出素子（ｘ，ｙ）に入射し、放射線検出素子（ｘ，ｙ）で光電変換された後に、画素単位の電気信号として読み出され、実写画像データＦ（ｘ，ｙ）が生成される。放射線画像検出器１が直接方式であれば、被写体を透過した放射線がａ−Ｓｅ等の光導電物質によって直接電荷に変換されて、複数の放射線検出素子１４によって画素単位に電気信号として読み出され、実写画像データＦ（ｘ，ｙ）が生成される。
【０３２３】
そして、信号読み取り操作により、各放射線検出素子（ｘ，ｙ）（各画素）から出力値として実写画像データＦ（ｘ、ｙ）が順次読み取られ、放射線画像検出器１の制御手段６により各放射線検出素子の番号（ｘ，ｙ）と各実写画像データＦ（ｘ、ｙ）とがそれぞれ対応付けられて記憶手段７に保存される。
【０３２４】
ダーク読取ステップ（ステップＳ３）では、続いて、放射線画像撮影の後に少なくとも１回ダーク読取が行われる。すなわち、放射線を照射しない状態で、放射線画像検出器１の各走査線Ｌｌに順次読み出し電圧を印加して各放射線検出素子（ｘ，ｙ）内に溜まった電荷をダーク読取値Ｄ（ｘ，ｙ）として取り出し、上記と同様に、放射線画像検出器１の制御手段６により各ダーク読取値Ｄ（ｘ，ｙ）がそれぞれ各放射線検出素子の番号（ｘ，ｙ）と対応付けられて記憶手段７に保存される。
【０３２５】
また、記憶手段７に記憶されている異常画素判定プログラムが読み出され、それに従って、前述した本発明に係る異常画素判定方法（［判定手法７］〜［判定手法１１］参照）が実行され、得られた実写画像データＦ（ｘ，ｙ）やダーク読取値Ｄ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）が異常画素であるか否かが判定される。なお、異常画素判定方法が実行される際、あわせて本発明に係る異常画素判定に基づく欠陥画素判定方法を実行するように構成することも可能である。
【０３２６】
そして、異常画素と判定できる放射線検出素子（ｘ，ｙ）があれば、その番号（ｘ，ｙ）が記憶手段７に一時的に記憶されるとともに、上記（４８）式等に従って異常画素と判定された放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）が正常値に置き換えられて記憶手段に記憶される。
【０３２７】
後述する温度補償変数算出ステップ（ステップＳ４）や時間的平均値等算出ステップ（ステップＳ５）、オフセット補正値算出ステップ（ステップＳ６）、画像補正ステップ（ステップＳ７）をコンソール３１等の外部装置で行う場合には、有線方式や無線方式によって、各放射線検出素子の番号（ｘ，ｙ）と対応付けられた各実写画像データＦ（ｘ，ｙ）や各ダーク読取値Ｄ（ｘ，ｙ）がそれぞれ放射線画像検出器１から外部装置に送信される。
【０３２８】
温度補償変数算出ステップ（ステップＳ４）では、センサパネル部４上の一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動し当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値Ｄ（ｘ´，ｙ´）とを用いて、上記（３）式に従って、今回の放射線画像撮影時に取得されたダーク読取における当該一の放射線検出素子（ｘ，ｙ）に関する空間的平均値である温度補償変数Ｗ（ｘ，ｙ）を算出する。そして、この温度補償変数Ｗ（ｘ，ｙ）の算出を、センサパネル部４上に２次元状に配置された全ての放射線検出素子（ｘ，ｙ）についてそれぞれ行う。
【０３２９】
その際、一の放射線検出素子（ｘ，ｙ）について温度補償変数Ｗ（ｘ，ｙ）を算出するにあたって、当該一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子の中に欠陥画素が含まれる場合には、前述した［算出手法１］〜［算出手法４］のいずれかの手法を用いて欠陥画素から出力されたダーク読取値Ｄ（ｘｓ，ｙｓ）がその近傍の放射線検出素子から出力されたダーク読取値で置換されたり補間されるなどして、当該一の放射線検出素子（ｘ，ｙ）の温度補償変数Ｗ（ｘ，ｙ）が算出される。
【０３３０】
時間的平均値等算出ステップ（ステップＳ５）では、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において得られた一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）のＭ回分の時間的平均値δ（ｘ，ｙ）が、上記（７）式に従って算出される。また、時間的平均値等算出ステップ（ステップＳ５）では、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）のＭ回分の時間的平均値ω（ｘ，ｙ）が、上記（９）式或いは（１０）式に従って算出される。
【０３３１】
なお、この時間的平均値等算出ステップ（ステップＳ５）は、必ずしも温度補償変数算出ステップ（ステップＳ４）の後や、実写画像データ取得ステップ（ステップＳ２）、ダーク読取ステップ（ステップＳ３）の後に行われる必要はなく、適宜のタイミングで事前に行っておくように構成することも可能である。
【０３３２】
オフセット補正値算出ステップ（ステップＳ６）では、時間的平均値等算出ステップ（ステップＳ５）で算出された時間的平均値δ（ｘ，ｙ）と時間的平均値ω（ｘ，ｙ）とから、それらの差分、すなわち温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）（＝δ（ｘ，ｙ）−ω（ｘ，ｙ））を算出する。そして、算出したε（ｘ，ｙ）と、温度補償変数算出ステップ（ステップＳ４）で算出した当該一の放射線検出素子（ｘ，ｙ）に関する温度補償変数Ｗ（ｘ，ｙ）とを加算して、当該一の放射線検出素子（ｘ，ｙ）に対するオフセット補正値Ｏ（ｘ，ｙ）を算出する。また、このオフセット補正値Ｏ（ｘ，ｙ）の算出を、センサパネル部４上に２次元状に配置された全ての放射線検出素子（ｘ，ｙ）についてそれぞれ行う。
【０３３３】
画像補正ステップ（ステップＳ７）では、前述した実写画像データ取得ステップ（ステップＳ２）で各放射線検出素子（ｘ，ｙ）から読み取られた実写画像データＦ（ｘ，ｙ）から、上記のオフセット補正値算出ステップ（ステップＳ６）で各放射線検出素子（ｘ，ｙ）について算出された各オフセット補正値Ｏ（ｘ，ｙ）をそれぞれ差し引いて実写画像データＦ（ｘ，ｙ）が補正される。そして、実写画像データＦ（ｘ，ｙ）と、オフセット補正値Ｏ（ｘ，ｙ）と、各放射線検出素子（ｘ，ｙ）について予め算出されているゲイン補正値Ｇ（ｘ，ｙ）とを上記（１）式に代入して、各放射線検出素子（ｘ，ｙ）について最終的な画像データＦＯ（ｘ，ｙ）が生成される。
【０３３４】
このようにして、実写画像データ取得ステップ（ステップＳ２）で生成された実写画像データＦ（ｘ、ｙ）に対して画素（放射線検出素子（ｘ，ｙ））ごとの特性ばらつきが補正されて最終的な画像データＦＯ（ｘ，ｙ）が生成される。
【０３３５】
なお、オフセット補正値算出ステップ（ステップＳ６）で処理される、過去のキャリブレーション時のダーク読取で得られた一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄｍ（ｘ，ｙ）や、画像補正ステップ（ステップＳ７）で処理される実写画像データＦ（ｘ，ｙ）等についても、欠陥画素であると判定されて登録されている放射線検出素子から出力されたデータである場合には、上記の欠陥画素がある場合の空間的平均値の算出手法と同様にして欠陥画素から出力されたデータに対して置換や補間等の処理がなされるように構成することも可能であり、適宜行われる。
【０３３６】
以上のように構成すれば、前述した放射線画像検出器１の各放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）取得の原理に従って、放射線画像撮影の前又は後にダーク読取を少なくとも１回行うだけで、放射線画像検出器１の各放射線検出素子（ｘ，ｙ）についてそれぞれ真値に近似する有効なオフセット補正値Ｏ（ｘ，ｙ）を取得することが可能となる。
【０３３７】
また、被写体の撮影に際し、ダーク読取の回数を少なくとも１回に抑えることができ、さらに、放射線画像検出器（ＦＰＤ）１から外部装置にデータを送信する場合に、画素（放射線検出素子（ｘ，ｙ））ごとの特性ばらつきが補正された最終的な画像データＦＯ（ｘ，ｙ）のみを１回送信するだけで済むため、ダーク読取値Ｄ（ｘ，ｙ）等を送信する必要がなく、データの送信に要する電力消費を低減させることが可能となる。
【０３３８】
特に、放射線画像検出器１がバッテリ内蔵型であり、画像データＦＯ（ｘ，ｙ）を無線通信するような場合でも、画像データＦＯ（ｘ，ｙ）の送信が１回で済むため、電力消費が低減されて内蔵バッテリ２１の消耗を防止することが可能となる。また、各放射線検出素子（ｘ，ｙ）について従来方式と同等の有効なオフセット補正値Ｏ（ｘ，ｙ）がそれぞれ取得され、良好な画像データＦＯ（ｘ，ｙ）を得ることが可能となる。
【０３３９】
また、欠陥画素判定ステップ（ステップＳ１）において上記の本実施形態に係る欠陥画素判定方法を用いることで、放射線検出素子が欠陥画素であるか否かが的確に判定され、欠陥画素であると判定されて登録されている放射線検出素子に対して置換や補間等の適切な処理が行われてオフセット補正値Ｏ（ｘ，ｙ）が算出されるため、有効なオフセット補正値Ｏ（ｘ，ｙ）を的確に取得することが可能となる。そのため、最終的に得られる画像データＦＯ（ｘ，ｙ）のＳＮ比を良好なものとすることが可能となる。
【０３４０】
なお、本実施形態の放射線画像検出器を用いた画像生成や、その前提となる上記の原理では、放射線画像撮影の前又は後のダーク読取で放射線検出素子（ｘ，ｙ）に対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の空間的平均値である温度補償変数Ｗ（ｘ，ｙ）として空間的統計値を算出する場合について説明した。また、過去のキャリブレーション時の複数回のダーク読取で放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）、及び当該放射線検出素子（ｘ，ｙ）に対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値ｄｍ（ｘ´，ｙ´）の空間的平均値ｗｍ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）として、それぞれ時間的統計値を算出する場合について説明した。
【０３４１】
しかし、空間的統計値や時間的統計値のとり方を、このように単純平均をとる算出法以外にも、種々のとり方を採用することが可能である。
【０３４２】
すなわち、空間的統計値のとり方としては、例えば、一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）のうち、当該一の放射線検出素子（ｘ，ｙ）に位置的に近い放射線検出素子（ｘ´，ｙ´）ほど出力された各ダーク読取値（ｘ´，ｙ´）に対する重みが大きくなるように重み付けして、その重み付け平均値を空間的統計値として算出するように構成することが可能である。
【０３４３】
また、例えば、一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値（ｘ´，ｙ´）の中央値（メディアン値）を空間的統計値として算出するように構成することも可能である。
【０３４４】
また、時間的統計値のとり方としては、例えば、過去のキャリブレーション時に複数回行われたダーク読取において、一の放射線検出素子（ｘ，ｙ）自身から出力された複数個のダーク読取値ｄｍ（ｘ，ｙ）や、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の複数個の空間的統計値である温度補償変数ｗｍ（ｘ´，ｙ´）のうち、時間的に現在に近いデータほどダーク読取値ｄｍ（ｘ，ｙ）や空間的統計値ｗｍ（ｘ´，ｙ´）に対する重みが大きくなるように重み付けして、その重み付け平均値を時間的統計値として算出するように構成することが可能である。
【０３４５】
また、例えば、一の放射線検出素子（ｘ，ｙ）自身から出力された複数個のダーク読取値ｄｍ（ｘ，ｙ）や、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の複数個の空間的統計値ｗｍ（ｘ´，ｙ´）の中央値を時間的統計値として算出するように構成することも可能である。
【０３４６】
なお、上記の欠陥画素判定方法において使用される空間的平均値ｗｍ（ｘ，ｙ）や空間的平均値Ｗ（ｘ，ｙ）の代わりに重み付け平均値や中央値を用いるように構成することが可能である。
【０３４７】
［異常画素判定システム、欠陥画素判定システム及び放射線画像生成システム］
前述したように、上述の［判定手法７］、［判定手法８］は、いずれもダーク読取値を使用するため、放射線の照射を伴わずに異常画素であるか否かを判定することができる。また、上述の［判定手法９］〜［判定手法１１］では、いずれも放射線を照射して得られた実写画像データを使用して異常画素であるか否かを判定することができる。さらに、上記のようにして異常画素であると判定された放射線検出素子（ｘ，ｙ）について上記の選択肢１）〜５）の方法のいずれかの方法で取り扱うことで、放射線画像検出器１において異常画素判定に基づく欠陥画素判定方法を実現することができる。
【０３４８】
そして、上記の［放射線画像検出器における異常画素判定や異常画素判定に基づく欠陥画素判定］では、この異常画素の判定や欠陥画素としての登録を放射線画像検出器１で行う場合について説明した。しかし、この異常画素の判定や欠陥画素としての登録を、放射線画像検出器１やコンソール等を含むシステムとして行うように構成することも可能である。以下に、本発明に係る異常画素判定システムや欠陥画素判定システムの実施形態の一例を示す。
【０３４９】
本実施形態では、異常画素判定システム３０や異常画素判定に基づいて欠陥画素判定を行う欠陥画素判定システムは、図２８に示すように、放射線画像検出器１（図１等参照）と、コンソール３１と、サーバ手段３９とを備えた放射線画像生成システムに構築されている。なお、以下では、コンソール３１が異常画素の判定を行う異常画素判定装置や異常画素判定に基づく欠陥画素判定を行う欠陥画素判定装置の機能を果たす場合について説明するが、サーバ手段３９等にそれらの機能を持たせるように構成することも可能である。
【０３５０】
放射線画像検出器１は、撮影室Ｒ１に設けられたブッキー装置３３の保持部３３ａに装填されて用いられるようになっている。ブッキー装置３３には、携帯情報端末様の小型の操作部３３ｂが設けられている。
【０３５１】
また、本実施形態では、放射線発生装置３４の１つの放射線源３４ａが、各ブッキー装置３３に共有されて用いられるようになっており、放射線画像検出器１をブッキー装置３３に装填して用いる場合には、対応する放射線源３４ａから放射線が照射される際に、放射線発生装置３４の放射線発生のタイミング制御が放射線画像検出器１の制御手段６と連動する仕組みになっており、放射線発生装置３４の放射線発生のタイミングに基づいて、放射線画像検出器１の各種制御が行われる。なお、各ブッキー装置３３に、それぞれ放射線発生装置３４の各放射線源３４ａを対応付けて設けるように構成することも可能である。
【０３５２】
放射線画像検出器１は、ブッキー装置３３に装填されずに、単独でフリーの状態で用いることもできるようになっている。その場合、例えばベッドタイプのブッキー装置３３の上に放射線画像検出器１を乗せた状態で使用したり、患者が放射線画像検出器１を手で抱えた状態で使用したりすることができる。
【０３５３】
放射線画像検出器１には、ダーク読取値Ｄ（ｘ，ｙ）等を無線方式で送信するための通信手段であるアンテナ装置３（図１参照）が設けられている。また、放射線画像検出器１には、放射線画像検出器１がブッキー装置３３の保持部３３ａに装填された際に保持部３３ａに設けられた図示しない電極と接続してダーク読取値Ｄ（ｘ，ｙ）等を有線方式で送信するための通信手段及び放射線画像検出器１への電源供給手段である端子１３（図２参照）が、放射線画像検出器１のアンテナ装置３や電源スイッチ１１の存在する筐体２の側面部分とは反対側（対面）の側面部分に設けられている。また、端子１３と同一側面部分には、バッテリ２１を充電するための電源供給手段である端子２２が設けられている。
【０３５４】
このように、放射線画像検出器１は、ブッキー装置３３に装填された際には端子１３と保持部３３ａの電極とが接続されてダーク読取値Ｄ（ｘ，ｙ）等がケーブル３２を介して有線方式で中継端末３５に送られ、中継端末３５を介してコンソール３１に送信されるようになっている。また、放射線画像検出器１は、ブッキー装置３３に装填された際には、ケーブル３２、端末１３を介して中継端末３５から放射線画像検出器１へ電力が供給されるように構成されている。
【０３５５】
また、ブッキー装置３３に装填されずに放射線画像検出器１が単独で用いられる場合には、アンテナ装置３を介してダーク読取値Ｄ（ｘ，ｙ）等を無線方式で送信するようになっている。撮影室Ｒ１には、放射線画像検出器１が無線方式でコンソール３１にダーク読取値Ｄ（ｘ，ｙ）等を通信する際に中継する無線アンテナ３６を備える中継端末３５が設けられている。そして、放射線画像検出器１のアンテナ装置３から無線方式で送信されたダーク読取値Ｄ（ｘ，ｙ）等の情報は、無線アンテナ３６で受信され、中継端末３５を介してコンソール３１に送信されるようになっている。
【０３５６】
また、放射線画像検出器１は、ブッキー装置３３に装填された際にはケーブル３２、端子１３を介して中継端末３５から電力の供給を受けるが、単独で用いられる場合には、内蔵のバッテリ２１の電力により動作するようになっている。なお、中継端末３５にはバッテリ充電用の図示しない電力供給手段である端子が設けてあり、放射線画像検出器１の端子２２を中継端末３５上の端子に接触させると、中継端末３５から図示しない端子、端子２２を介して放射線画像検出器１のバッテリ２１に電力が供給され、バッテリ２１が充電されるようになっている。しかし、この形態には限定されない。
【０３５７】
また、放射線画像検出器１内には、図示しないタグが内蔵されている。この場合、タグとして、いわゆるＲＦＩＤ（Radio Frequency IDentification）タグが用いられており、タグには、タグの各部を制御する制御回路や放射線画像検出器１のＩＤ等の固有情報を記憶する記憶部がコンパクトに内蔵されている。
【０３５８】
前室Ｒ２の入口の近傍には、放射線画像検出器１のＲＦＩＤタグを読み取るタグリーダ３７が設置されている。タグリーダ３７は、内蔵する図示しないアンテナを介して電波等に所定の指示情報を乗せて発信し、前室Ｒ２に入室し或いは退室する放射線画像検出器１を検出して、放射線画像検出器１のＩＤ等をコンソール３１に送信するようになっている。
【０３５９】
なお、放射線画像検出器１内にＲＦＩＤタグを内蔵させる場合の例を示したが、放射線画像検出器１内にＲＦＩＤタグを内蔵させる代わりに、放射線画像検出器１の筐体外側表面に放射線画像検出器１のＩＤ等を書き込んだバーコードを添付しておき、これをバーコードリーダーで読み取るようにしても良い。この場合、タグリーダ３７が、バーコードリーダーとなり、前室Ｒ２に入室し或いは退室する放射線画像検出器１のバーコードを読み取って、放射線画像検出器１のＩＤ等をコンソール３１に送信する。
【０３６０】
なお、タグリーダ３７やバーコードリーダーは前室Ｒ２の入口の近傍の代わりに、撮影室Ｒ１の入口近傍に設置し、撮影室Ｒ１への放射線画像検出器１の入退室を管理するようにしても良い。
【０３６１】
このように、ＲＦＩＤやバーコードを使って、前室Ｒ２や撮影室Ｒ１への放射線画像検出器１の入室或いは退室をコンソール３１に通知することで、放射線画像検出器１をどの放射線発生装置（本例の場合は放射線発生装置３４）と連動されれば良いかを、コンソール３１や放射線発生装置３４に自動的に知らせることができる。このようなＲＦＩＤやバーコードを使った放射線画像検出器１の撮影室等への入室或いは退室の管理は、撮影室や放射線発生装置が複数存在する施設で有効に機能する。なお、ＲＦＩＤやバーコードを使用しない場合は、放射線画像検出器１の撮影室等への入室或いは退室の情報を、使用者が直接コンソール３１へ入力するようにしても良い。
【０３６２】
前室Ｒ２には、異常画素判定システム（欠陥画素判定システム）３０や放射線画像生成システム全体の制御を行うコンソール３１が設けられており、コンソール３１には、前述した中継端末３５やタグリーダ３７、放射線発生装置３４の本体部３４ｂ等が接続されており、また、中継端末３５を介してブッキー装置３３等が接続されている。
【０３６３】
コンソール３１は、図示しないＣＰＵ（Central Processing Unit）やＲＯＭ、ＲＡＭ、ハードディスク等の記憶手段、入出力インターフェース等がバスに接続されたコンピュータで構成されている。
【０３６４】
また、コンソール３１には、ネットワークＮＷを介してコンピュータからなるサーバ手段３９が接続されている。また、サーバ手段３９は、ハードディスク等からなる記憶手段３８が接続されている。
【０３６５】
記憶手段３８及びコンソール３１自身の図示しない記憶手段には、放射線画像生成システムで使用可能な各放射線画像検出器１について、前述した空間的平均値を算出するために、各放射線検出素子（ｘ，ｙ）にどの放射線検出素子（ｘ´，ｙ´）を対応付けるかの情報が放射線画像検出器１のＩＤに対応付けられて予めそれぞれ記憶されている。
【０３６６】
また、コンソール３１は、放射線画像検出器１から、キャリブレーション時に行われたダーク読取で各放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）や、放射線画像撮影の直前又は直後に１回以上行われたダーク読取において得られたダーク読取値Ｄ（ｘ，ｙ）、Ｄｍ（ｘ，ｙ）が有線方式或いは無線方式が送信されてくると、それらを当該放射線画像検出器１のＩＤに対応付けて、自らの記憶手段やサーバ手段３９の記憶手段３８に記憶させるようになっている。
【０３６７】
さらに、サーバ手段３９は、コンソール３１から放射線画像検出器１のダーク読取値ｄｍ（ｘ，ｙ）等が送信されてくると、それらを記憶手段３８に保存するとともに、当該放射線検出器１の各放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）や、各放射線検出素子（ｘ，ｙ）に対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄｍ（ｘ´，ｙ´）の温度補償変数としての空間的平均値ｗｍ（ｘ´，ｙ´）の時間的統計値（時間的平均値）ω（ｘ，ｙ）を算出して記憶手段３８に保存するようになっている。
【０３６８】
なお、このダーク読取値ｄ（ｘ，ｙ）の保存及び時間的平均値δ（ｘ，ｙ）等の算出、保存を、コンソール３１やコンソール３１の記憶手段で行うように構成することも可能である。
【０３６９】
コンソール３１は、ＲＯＭやハードディスク等の記憶手段に格納されている各種処理の実行に必要なプログラムを読み出してＲＡＭの作業領域に展開して、プログラムに従って処理を実行するようになっている。
【０３７０】
また、ＲＯＭやハードディスク等の記憶手段には、前述した異常画素判定プログラムが格納されており、本実施形態では、異常画素判定装置としてのコンソール３１が、この異常画素判定プログラムに従って上記の本発明に係る異常画素判定方法を実行するようになっている。また、異常画素判定プログラムは、前述した本発明に係る異常画素判定に基づく欠陥画素判定方法をも実行させるようになっており、その場合、欠陥画素判定装置としてのコンソール３１が異常画素判定プログラムに従って本発明に係る異常画素判定に基づく欠陥画素判定方法を実行するようになっている。
【０３７１】
具体的には、コンソール３１は、例えば、異常画素判定プログラムに従って前述した異常画素判定方法のうち［判定手法７］で説明した異常画素判定を行う場合、まず、放射線画像撮影の前又は後に少なくとも１回ダーク読取を行った放射線画像検出器１から送信されてきた各放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）を入手する。或いは、前もって自らの記憶手段やサーバ手段３９の記憶手段３８（以下、単に記憶手段という。）に記憶されている当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）を読み出す。
【０３７２】
そして、コンソール３１は、当該放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力されたダーク読取値Ｄ（ｘ´，ｙ´）の空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）との差分Ｅ（ｘ，ｙ）を算出し（上記（３８）式参照）、その絶対値｜Ｅ（ｘ，ｙ）｜が予め設定された閾値Ｅthを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定する。
【０３７３】
また、コンソール３１は、例えば、異常画素判定プログラムに従って前述した異常画素判定方法のうち［判定手法８］で説明した異常画素判定を行う場合、［判定手法７］と同様に一の放射線検出素子（ｘ，ｙ）についてダーク読取値Ｄ（ｘ，ｙ）と空間的平均値Ｗ（ｘ，ｙ）の差分Ｅ（ｘ，ｙ）を（３８）式に従って第１差分として計算し、また、過去のキャリブレーション時に複数回行われたダーク読取において（１２）式に従って計算されたダーク読取値ｄｍ（ｘ，ｙ）の分布の時間的平均値δ（ｘ，ｙ）とダーク読取値ｄｍ（ｘ，ｙ）の空間的平均値ｗｍ（ｘ，ｙ）の分布の時間的平均値ω（ｘ，ｙ）との差分ε（ｘ，ｙ）を第２差分として計算し、これらの第１差分Ｅ（ｘ，ｙ）と第２差分ε（ｘ，ｙ）との差分Λ（ｘ，ｙ）（上記（４１）式参照）の絶対値｜Λ（ｘ，ｙ）｜が予め設定された閾値Λthを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定する。
【０３７４】
このように構成すれば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）が、それと対応付けられた複数の放射線検出素子から出力されたダーク読取値Ｄ（ｘ´，ｙ´）の空間的平均値Ｗ（ｘ，ｙ）等とかけ離れた異常な値である場合に、的確に当該放射線検出素子（ｘ，ｙ）を異常画素であると判定することが可能となる。
【０３７５】
そして、異常画素であると判定した当該一の放射線検出素子（ｘ，ｙ）の情報を記憶手段（コンソール３１の記憶手段やサーバ手段３９の記憶手段３８）に一時的に記憶させておけば、その画素が欠陥画素として欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素に起因するアーチファクト（画像欠陥）が発生することを有効に防止することが可能となる。
【０３７６】
また、コンソール３１は、例えば、異常画素判定プログラムに従って前述した異常画素判定方法のうち［判定手法９］や［判定手法１０］で説明した異常画素判定を行う場合、まず、被写体が存在する状態で放射線が照射された放射線画像検出器１から送信されてきた各放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を入手する。或いは、前もって放射線画像撮影が行われ記憶手段に記憶されている当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を読み出す。
【０３７７】
そして、コンソール３１は、当該放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）との差分Ｖ（ｘ，ｙ）（上記（４５）式参照）が、予め設定された閾値Ｖthを越えた場合、或いはその絶対値｜Ｖ（ｘ，ｙ）｜が予め設定された閾値Ｖ´thを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定する。
【０３７８】
また、コンソール３１は、例えば、異常画素判定プログラムに従って前述した異常画素判定方法のうち［判定手法１１］で説明した異常画素判定を行う場合、被写体が存在しない状態で放射線が照射された放射線画像検出器１から送信されてきた各放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を入手する。或いは、前もって被写体が存在しない状態で放射線画像撮影が行われ記憶手段に記憶されている当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を読み出す。
【０３７９】
そして、コンソール３１は、当該放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）と、この放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）との差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ，ｙ）｜が予め設定された閾値Ｖ´thを越えた場合に、当該一の放射線検出素子（ｘ，ｙ）が異常画素であると判定する。
【０３８０】
このように構成すれば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）が、それと対応付けられた複数の放射線検出素子から出力された実写画像データＦ（ｘ´，ｙ´）の空間的平均値Ａ（ｘ，ｙ）とかけ離れた突発的に異常な値である場合に、的確に当該放射線検出素子（ｘ，ｙ）を異常画素であると判定することが可能となる。
【０３８１】
そして、異常画素であると判定した当該一の放射線検出素子（ｘ，ｙ）の情報を記憶手段（コンソール３１の記憶手段やサーバ手段３９の記憶手段３８）に一時的に記憶させておけば、その画素が欠陥画素として欠陥画素マップに登録していなくとも、一時的に欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素が発生することを有効に防止することが可能となる。
【０３８２】
更に、撮影対象となる部位が、微小石灰化診断の為の乳房撮影等の画素単位の出力レベルが重要となる場合には、コンソール上で使用予定の放射線画像検出器IDを入力して、当該検出器の異常履歴を表示せしめ、過去に少なくとも一度、突発的な異常を呈した画素が放射線画像検出器のどの領域であるかを技師が確認できれば、当該放射線画像検出器により出力された最終画像データにおける微小石灰化病変部に対応する画像に対する信頼性を向上させることができる。なお、技師は、関心領域となる検出器中央部に、現在は正常であるが、過去に異常を呈したことのある画素が存在する検出器を使用せず、新たな検出器を選定することも可能であり、診断精度を維持し、誤診を予防することができる。
【０３８３】
また、関心領域内に当該異常画素が存在する場合には、被写体に対する放射線画像検出器の位置を変更して、照射野外（関心領域外）とすることで診断精度を維持することができる（乳房画像の読影時において、通常、照射野外領域は、眩惑防止の為、正常／異常に係らず一律低輝度化処理を行う）。
【０３８４】
さらに、コンソール３１で上記のようにして異常画素であると判定された放射線検出素子（ｘ，ｙ）について上記の選択肢１）〜５）の方法のいずれかの方法で取り扱うことで、異常画素判定システム（欠陥画素判定システム）で異常画素判定に基づく欠陥画素判定方法を実現することができる。
【０３８５】
具体的には、コンソール３１が、異常画素であると判定した放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）である異常な信号値の大きさに応じて、その信号値の大きさが予め設定された閾値を越える場合に欠陥画素であると判定したり、異常画素であると判定した放射線検出素子（ｘ，ｙ）から異常な信号値を出力する確率が予め設定された発生確率の閾値を越えた場合に欠陥画素であると判定するように構成することができる。
【０３８６】
そして、欠陥画素であると判定した放射線検出素子の情報を自らの記憶手段やサーバ手段３９の記憶手段３８に保存されている欠陥画素マップに登録するように構成することで、前述した定常的な欠陥画素と同程度に異常な信号値を出力する異常画素を欠陥画素として的確に判定して管理することが可能となるとともに、そのような異常画素（欠陥画素）について、周囲の画素値を用いて補間処理により異常画素自体の異常値の補正を行う補正処理等を施すことが可能となり、最終的な画像データＦＯ（ｘ，ｙ）の中に異常画素が発生することを有効に防止することが可能となる。
【０３８７】
なお、コンソール３１の記憶手段には、コンソール３１に前述した［判定手法１］〜［判定手法６］に従って欠陥画素判定を実行させるための欠陥画素判定プログラムが格納されており、コンソール３１は、上記の欠陥画素判定方法の［判定手法１］〜［判定手法６］のうちのいずれかの判定手法を実行し、或いはそれらの判定手法を組み合わせて実行するようになっている。また、それにより欠陥画素と判定された放射線画像検出器１の放射線検出素子（ｘ，ｙ）の情報は欠陥画素マップに登録される。
【０３８８】
一方、ＲＯＭやハードディスク等の記憶手段には、前述した実写画像データ取得ステップ（図２７のステップＳ２）や、ダーク読取ステップ（ステップＳ３）、温度補償変数算出ステップ（ステップＳ４）、時間的平均値等算出ステップ（ステップＳ５）、オフセット補正値算出ステップ（ステップＳ６）、画像補正ステップ（ステップＳ７）を実現するための画像生成プログラムが格納されており、コンソール３１は、この画像生成プログラムを読み出してＲＡＭの作業領域に展開して、この画像生成プログラムに従って必要な処理を実行し、また、装置各部を制御するようになっている。
【０３８９】
以下、システムが放射線画像生成システムとして機能する場合について説明する。コンソール３１は、放射線技師や医師等から放射線画像撮影の指示を受けると、撮影条件（管電流、管電圧等）や照射開始を制御する放射線発生装置３４を駆動する等して実写画像データ取得ステップ（図２７のステップＳ２）を実行する。また、放射線画像撮影の前又は後に、少なくとも１回のダーク読取を行うように放射線画像検出器１に指示信号を送信して、ダーク読取ステップ（ステップＳ３）を実行する。
【０３９０】
放射線画像検出器１は、自らのＩＤと、各放射線検出素子（ｘ，ｙ）について実写画像データ取得ステップ（ステップＳ２）で取得した放射線画像の実写画像データＦ（ｘ，ｙ）とダーク読取ステップ（ステップＳ３）で取得したダーク読取値Ｄ（ｘ，ｙ）とを有線方式又は無線方式でコンソール３１に送信する。コンソール３１は、放射線画像検出器１から実写画像データＦ（ｘ，ｙ）とダーク読取値Ｄ（ｘ，ｙ）とが送信されてくると、それらを自らの記憶手段に保存する。
【０３９１】
そして、取得した実写画像Ｆ（ｘ，ｙ）から、まず、１／１６等に圧縮した間引きデータを生成して表示し、再撮影の必要性有無を迅速に判断可能とする。なお放射線画像検出器が間引きデータを生成し送信する場合には、当該間引きデータを表示する。
この段階で、ポジショニング不良等が認められれば、再撮影を実行することとなり、オフセット補正値の生成処理や欠陥や異常画素有無の判定処理は行なわず、再撮影のデータ取得まで待機する。
【０３９２】
また、間引きデータ表示時に、ポジショニングと共に各画素出力が所定範囲内か否か、
言い換えると、飽和状態ではなく被写体を透過したＸ線量に比例した出力値を示しているか否かを知る為に、間引きデータに対しオフセット／ゲイン補正処理済の画像を表示することとしても良い。
【０３９３】
また、コンソール３１は、ＲＯＭ等の記憶手段に記憶されている異常画素判定プログラムに従って、前述した異常画素判定方法（［判定手法７］〜［判定手法１１］参照）を実行し、放射線画像検出器１から送信されてきた実写画像データＦ（ｘ，ｙ）やダーク読取値Ｄ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）が異常画素であるか否かを判定する。なお、［判定手法１１］を実行する場合には、被写体が存在しない状態で放射線画像検出器１に対して放射線発生装置３４の放射線源３４ａから放射線が照射される。
【０３９４】
そして、コンソール３１は、異常画素と判定できる放射線検出素子（ｘ，ｙ）があれば、その番号（ｘ，ｙ）を記憶手段に一時的に記憶させるとともに、上記（４８）式等に従って異常画素と判定された放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）を正常値に置き換えて記憶手段に記憶させる。
【０３９５】
コンソール３１は、続いて、温度補償変数算出ステップ（ステップＳ４）を実行する。コンソール３１は、放射線画像検出器１のＩＤを参照して、記憶手段から各放射線検出素子（ｘ，ｙ）に対応付ける複数の放射線検出素子（ｘ´，ｙ´）の情報を読み出し、それらから出力された各ダーク読取値Ｄ（ｘ´，ｙ´）を用いて、上記（３）式に従って各放射線検出素子（ｘ，ｙ）についてそれぞれ温度補償変数Ｗ（ｘ，ｙ）を算出する。
【０３９６】
その際、コンソール３１は、一の放射線検出素子（ｘ，ｙ）について温度補償変数Ｗ（ｘ，ｙ）を算出するにあたって、当該一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子の中に欠陥画素が含まれる場合には、前述した［算出手法１］〜［算出手法４］のいずれかの手法を用いて欠陥画素から出力されたダーク読取値Ｄ（ｘｓ，ｙｓ）をその近傍の放射線検出素子から出力されたダーク読取値で置換したり補間するなどして、当該一の放射線検出素子（ｘ，ｙ）の温度補償変数Ｗ（ｘ，ｙ）を算出する。
【０３９７】
また、コンソール３１は、予め過去のキャリブレーション時等に得られたデータに基づいて時間的平均値等算出ステップ（ステップＳ５）を実行するが、本実施形態では、当該放射線画像検出器１についての時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗｍ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）は、前述したように、予めサーバ手段３９で予め計算され、記憶手段３８に保存されている。コンソール３１がこの段階で自らそれらの値を算出するように構成してもよい。
【０３９８】
コンソール３１は、続いて、オフセット補正値算出ステップ（ステップＳ６）を実行する。その際、コンソール３１は、タグリーダ３７を介して送信されてきた放射線画像検出器１のＩＤをネットワークＮＷを介してサーバ手段３９に送信する。サーバ手段３９は、放射線画像検出器１のＩＤに基づいて、当該放射線画像検出器１の時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗｍ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）を記憶手段３８から読み出して、コンソール３１に送信する。
【０３９９】
コンソール３１は、当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）について、サーバ手段３９から入手したダーク読取値ｄｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗｍ（ｘ´，ｙ´）の時間的平均値ω（ｘ，ｙ）とから、それらの差分ε（ｘ，ｙ）（＝δ（ｘ，ｙ）−ω（ｘ，ｙ））を算出する。そして、算出したε（ｘ，ｙ）と、温度補償変数算出ステップ（ステップＳ４）で算出した当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）についての各温度補償変数Ｗ（ｘ，ｙ）とを加算して、各放射線検出素子（ｘ，ｙ）に対するオフセット補正値Ｏ（ｘ，ｙ）を算出する。
【０４００】
コンソール３１は、続いて、画像補正ステップ（ステップＳ７）を実行し、実写画像データ取得ステップ（ステップＳ２）で各放射線検出素子（ｘ，ｙ）から取得した各実写画像データＦ（ｘ，ｙ）から、オフセット補正値算出ステップ（ステップＳ６）で各放射線検出素子（ｘ，ｙ）について算出した各オフセット補正値Ｏ（ｘ，ｙ）をそれぞれ差し引き、その差分に各放射線検出素子（ｘ，ｙ）について予め算出されているゲイン補正値Ｇ（ｘ，ｙ）を乗算して、各放射線検出素子（ｘ，ｙ）について最終的な画像データＦＯ（ｘ，ｙ）を生成する。
【０４０１】
このように放射線画像生成システムを構成することで、上記の放射線画像検出器を用いた画像生成方法の効果と同様の効果が有効に発揮され、各放射線検出素子（ｘ，ｙ）についてそれぞれ有効なオフセット補正値Ｏ（ｘ，ｙ）が取得され、実写画像データ取得ステップ（ステップＳ２）で生成された実写画像データＦ（ｘ、ｙ）に対して画素（放射線検出素子１４）ごとの特性ばらつきが補正されて最終的な画像データＦＯ（ｘ，ｙ）が生成される。
【０４０２】
また、それとともに、放射線画像撮影の前又は後にダーク読取を少なくとも１回行うだけで済み、放射線画像の実写画像データＦ（ｘ，ｙ）と少なくとも１回分のダーク読取値Ｄ（ｘ，ｙ）を放射線画像検出器１から送信するだけで済むため、放射線画像撮影に要する電力消費を低減させることが可能となる。
【０４０３】
特に、放射線画像検出器１がバッテリ内蔵型である場合には、ダーク読取が少なくとも１回で済み、ダーク読取値Ｄ（ｘ，ｙ）の送信が少なくとも１回分だけで済むため、放射線画像撮影やダーク読取値Ｄ（ｘ，ｙ）の読取及び送信に要する電力消費が低減され、内蔵バッテリ２１の消耗を防止することが可能となる。また、放射線画像検出器１を有線方式で使用する場合であっても、ダーク読取回数やダーク読取値Ｄ（ｘ，ｙ）の送信回数が低減できるため、消費電力を抑えることができる。
【０４０４】
また、上記の本実施形態に係る欠陥画素判定方法や本発明に係る異常画素判定方法、或いは本発明に係る異常画素判定に基づく欠陥画素判定方法を用いて異常画素判定や欠陥画素判定を行うことで、放射線検出素子が異常画素や欠陥画素であるか否かを的確に判定し、欠陥画素であると判定して登録している放射線検出素子に対して置換や補間等の適切な処理を行い、異常画素と判定された放射線検出素子に対しても欠陥画素と同等に置換や補間等の適切な処理を行ってオフセット補正値Ｏ（ｘ，ｙ）を算出することが可能となり、有効なオフセット補正値Ｏ（ｘ，ｙ）を的確に取得することが可能となる。そのため、最終的に得られる画像データＦＯ（ｘ，ｙ）のＳＮ比を良好なものとすることが可能となる。
【０４０５】
なお、本実施例では、放射線画像撮影の前又は後にダーク読取を１回行う場合に特に有効な手法として説明を行ったが、ダーク読取を２回以上行う場合であっても同様な効果を得ることができる。すなわち、キャリブレーション時に実施したダーク読取り回数をＭ回、放射線画像撮影の前又は後に行うダーク読取り回数をＫ回とした場合に、Ｍ＞Ｋの関係が成立すれば、本発明で説明したのと同様の効果が得られることは言うまでもない。
【０４０６】
すなわち、Ｋ＞２の場合は、ダーク読取値をＤｋ（ｘ，ｙ）（ｋ=１〜Ｋ、Ｋ＞２）とおき、
【数１７】
2010074644.tif 000019
を本発明の実施例のＤ（ｘ，ｙ）にあてはめて考えれば良い。すなわち、
Ｄ（ｘ，ｙ）=Ｄｋave（ｘ，ｙ）    …（５６）
として考えれば良い。
【０４０７】
また、当該放射線画像検出器１の放射線検出素子（ｘ，ｙ）中に欠陥画素がある場合でも、上記の空間的統計値（空間的平均値）の算出手法を用いることで、欠陥画素から出力されるダーク読取値Ｄが適切に置換され或いは補間されて、有効な空間的統計値（空間的平均値）である温度補償変数Ｗ（ｘ，ｙ）やｗｍ（ｘ´，ｙ´）を算出することが可能となる。また、それにより、オフセット補正値Ｏ（ｘ，ｙ）を的確に算出することが可能となる。
【０４０８】
なお、本発明が上記の実施の形態に限定されず、適宜変更可能であることはいうまでもない。本発明の概念は、一の放射線検出素子と同じように温度変動する放射線検出素子（ｘ´，ｙ´）を予め選択して当該一の放射線検出素子（ｘ，ｙ）に予め対応付けるものである。そのため、放射線検出素子自体は、上記の実施形態のように放射線検出素子タイプである必要はなく、他の構造を有するものを２次元的に配置したセンサパネル部を備える放射線画像検出器にも適用可能であることはいうまでもない。
【０４０９】
また、本発明では、放射線検出素子の画素単位での温度変動を取り上げて説明を行ったが、温度変動以外の因子であっても、ある画素値が近傍の画素値と類似の変化をもたらすような変動であれば、同様の手法を適応することで同様の効果を得られることは言うまでもない。
【０４１０】
さらに、本実施形態では、データ処理をコンソール３１で行う場合について説明したが、実写画像データＦ（ｘ，ｙ）から各画素ごとに各オフセット補正値Ｏ（ｘ，ｙ）をそれぞれ差し引き、実写画像データＦ（ｘ、ｙ）を補正して最終的な画像データＦＯ（ｘ，ｙ）を生成するまでの全てのデータ処理を放射線画像検出器１で行うように構成することも可能である。
【０４１１】
さらに、例えば、上記の実施形態においては、放射線画像撮影時とは別に行われるキャリブレーション時にダーク読取を行うこととしたが、放射線画像撮影の前又は後に行うダーク読取値を保存しておき、それを用いるように構成してもよい。この場合には、例えば連続して行われた放射線画像撮影の複数回分のダーク読取値を正規化する等して保存しておき、当該複数回分のダーク読取値を用いて上記の実施形態と同様に時間的統計値（時間的平均値）δ（ｘ，ｙ）やω（ｘ，ｙ）を算出するようにしてもよい。
【０４１２】
また、本実施例では、オフセット補正を例に取り、説明を行ったが、ゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、放射線画像検出器に対して所定の条件で一様な放射線を照射し、この読取値（本実施例の中の実写画像データＦ（ｘ，ｙ）に相当）に対してもオフセット補正が適応される。このため、本実施例では、実写画像データＦ（ｘ，ｙ）に対してオフセット補正を行う場合を例に取り説明を行ったが、前記理由によりゲイン補正値Ｏ（ｘ，ｙ）を求める際にも同様な処置が実施可能であり、かつ同様の効果があることは言うまでもない。
【図面の簡単な説明】
【０４１３】
【図１】本実施形態に係る放射線画像検出器の外観構成を示す図である。
【図２】放射線画像検出器のセンサパネル部及び読取部の構成に示す等価回路図である。
【図３】放射線検出素子に割り当てられた番号を説明する図である。
【図４】複数回のダーク読取ごとに放射線検出素子から出力されるダーク読取値を説明する図である。
【図５】複数回のダーク読取で放射線検出素子から出力される複数のダーク読取値の時間的ゆらぎの分布を説明するグラフである。
【図６】放射線検出素子に対応付ける複数の放射線検出素子のとり方の一例を示す図である。
【図７】放射線検出素子が異なるとダーク読取値の時間的ゆらぎの分布は平均値や標準偏差が異なることを説明するグラフである。
【図８】空間的平均値Ｗ（ｘ，ｙ）の分布の標準偏差が小さい分布となることを説明するグラフである。
【図９】Ｄ（ｘ，ｙ）の時間的ゆらぎの分布とＤ（ｘ’，ｙ’）の時間的ゆらぎの分布の温度変化によるシフト度合いが異なるとブロードとなる温度補償変数の分布を示すグラフである。
【図１０】１回のダーク読取では放射線検出素子からオフセット補正値の真値を推定することが難しいことを説明するグラフである。
【図１１】複数回のダーク読取ごとに放射線検出素子から出力されるダーク読取値及び空間的平均値を説明する図である。
【図１２】複数回のダーク読取ごとに出力されるダーク読取値や空間的平均値が正規分布状に分布することを説明するグラフである。
【図１３】放射線検出素子の温度が時間的に変動することを説明するグラフである。
【図１４】放射線画像撮影の前又は後に行われる１回のダーク読取で放射線検出素子から出力されるダーク読取値や空間的平均値について推定される時間的ゆらぎの分布を説明するグラフである。
【図１５】複数回のダーク読取中に温度変動が起こるとダーク読取値等の分布はブロードになるが、温度補正済みダーク読取値の分布はブロードにならない様子を示すグラフである。
【図１６】（Ａ）〜（Ｃ）は一の放射線検出素子が正方領域の中心ではない位置に位置するように設定された領域の例を示す図である。
【図１７】所定個数の放射線検出素子ごとに１個ずつ接続される読み出しＩＣを説明する図である。
【図１８】領域に属する各放射線検出素子からの電気信号が隣接する２つの読み出しＩＣで別々に読み出される場合を説明する図である。
【図１９】（Ａ）、（Ｂ）は通常の仕方で複数の放射線検出素子が選択される場合、（Ｃ）、（Ｄ）は隣接する読み出しＩＣの境界の近傍をセンサパネル部の周縁部分と同様に扱う場合を説明する図である。
【図２０】放射線検出素子からのダーク読取値のゆらぎの分布の標準偏差が平均的な分布とそれよりも格段に大きい分布を表すグラフである。
【図２１】放射線検出素子からのダーク読取値のゆらぎの分布の平均値が閾値よりも大きい場合を表すグラフである。
【図２２】差分のゆらぎの分布の標準偏差が平均的な分布を表すグラフである。
【図２３】差分のゆらぎの分布の標準偏差が図２２に示した場合よりも大きい分布を表すグラフである。
【図２４】注目画素が放射線検出素子（４，４）であり、それに対応付けられた７×７個の放射線検出素子うち放射線検出素子（６，６）が欠陥画素である場合の例を説明する図である。
【図２５】図２４の例で領域外の位置に存在する放射線検出素子で欠陥画素を置換する例を示す図である。
【図２６】図２４の例で放射線検出素子（７，６）が欠陥画素の場合に領域外の位置に存在する放射線検出素子で欠陥画素を補間する例を説明する図である。
【図２７】放射線画像検出器を用いた画像生成の処理手順を示すフローチャートである。
【図２８】本実施形態に係る欠陥画像判定システムを含む放射線画像生成システムの構成を示す図である。
【符号の説明】
【０４１４】
１  放射線画像検出器
３  アンテナ装置（無線方式の通信手段）
４、５、６  センサパネル部、読取部及び制御手段（信号値取得手段）
６  制御手段（判定手段）
７  記憶手段
１３  端子（通信手段）
１４、（ｘ，ｙ）  放射線検出素子
２１  バッテリ
３０  異常画素判定システム、欠陥画素判定システム
３１  コンソール（異常画素判定装置、欠陥画素判定装置）
３８  記憶手段
３９  サーバ手段（異常画素判定装置、欠陥画素判定装置）
Ａ（ｘ，ｙ）  実写画像データの空間的平均値
Ｄ（ｘ，ｙ）  ダーク読取値
Ｅ（ｘ，ｙ）  差分、第１差分
Ｅth  閾値
Ｆ（ｘ，ｙ）  実写画像データ
Ｖ（ｘ，ｙ）  差分
Ｖth  閾値
Ｖ´th  閾値
Ｗ（ｘ，ｙ）  空間的平均値
（ｘｓ，ｙｓ）  欠陥画素
δ（ｘ，ｙ）  ダーク読取値の時間的平均値
ε（ｘ，ｙ）  第２差分
｜Λ（ｘ，ｙ）｜  第１差分と第２差分との差分の絶対値
Λth  閾値
ω（ｘ，ｙ）  ダーク読取値の空間的平均値の時間的平均値
【図１】
2010074644.tif 000020
【図２】
2010074644.tif 000021
【図３】
2010074644.tif 000022
【図４】
2010074644.tif 000023
【図５】
2010074644.tif 000024
【図６】
2010074644.tif 000025
【図７】
2010074644.tif 000026
【図８】
2010074644.tif 000027
【図９】
2010074644.tif 000028
【図１０】
2010074644.tif 000029
【図１１】
2010074644.tif 000030
【図１２】
2010074644.tif 000031
【図１３】
2010074644.tif 000032
【図１４】
2010074644.tif 000033
【図１５】
2010074644.tif 000034
【図１６】
2010074644.tif 000035
【図１７】
2010074644.tif 000036
【図１８】
2010074644.tif 000037
【図１９】
2010074644.tif 000038
【図２０】
2010074644.tif 000039
【図２１】
2010074644.tif 000040
【図２２】
2010074644.tif 000041
【図２３】
2010074644.tif 000042
【図２４】
2010074644.tif 000043
【図２５】
2010074644.tif 000044
【図２６】
2010074644.tif 000045
【図２７】
2010074644.tif 000046
【図２８】
2010074644.tif 000047