WO2010032498A1

WO2010032498A1 - 放射線画像生成システム

Info

Publication number: WO2010032498A1
Application number: PCT/JP2009/054033
Authority: WO
Inventors: 久米川
Original assignee: コニカミノルタエムジー株式会社
Priority date: 2008-09-19
Filing date: 2009-03-04
Publication date: 2010-03-25
Also published as: JP2011250810A

Abstract

　間引きデータを画面上に迅速に表示させることが可能な放射線画像生成システムを提供する。　放射線画像生成システム３０は、放射線画像検出器１と、コンソール３１と、オフセット補正値生成手段３１と、放射線発生装置３４とを備え、コンソール３１は、間引きデータｆ（ｘ，ｙ）又は実写画像データＦ（ｘ，ｙ）に基づいて作成した間引きデータに基づく画像、或いは、間引きデータ又は実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及びオフセット補正値に基づくオフセット補正処理を行って間引きデータ又は実写画像データに基づいて作成した間引きデータに対する画像処理後の画像、或いは、実写画像データに対してゲイン補正処理及びオフセット補正値に基づくオフセット補正処理を行った画像処理後の画像、のうちの少なくとも１つとを、前記表示部に表示させる。

Description

放射線画像生成システム

　本発明は、放射線画像生成システムに関するものである。

　通称フラットパネルディテクタ（Flat Panel Detector：ＦＰＤ）と呼ばれる固体撮像素子を２次元的に配置した放射線画像検出器には、検出素子として、ａ－Ｓｅ（アモルファスセレン）のような光導電物質を用いて放射線エネルギーを直接電荷に変換し、この電荷を２次元的に配置されたＴＦＴ（Thin Film Transistor：薄膜トランジスタ）等の信号読み出し用のスイッチ素子によって画素単位に電気信号として読み出す直接方式や、放射線エネルギーをシンチレータ等で光に変換し、この光を２次元的に配置されたフォトダイオード等の光電変換素子で電荷に変換してＴＦＴ等によって電気信号として読み出す間接方式等が良く知られている。

　そして、いずれの方式においても、被写体を透過してきた放射線を放射線画像検出器で検出して得られた実写画像データに対してゲイン補正やオフセット補正等を行い、実写画像データを補正する必要があることが知られている。

　一般的に、実写画像データの補正では、下記（１）式に示される通り、放射線画像検出器の各放射線検出素子（センサパネル部における座標は（ｘ，ｙ））から出力された実写画像データＦ（ｘ，ｙ）からオフセット補正値Ｏ（ｘ，ｙ）を差し引き、その差分にゲイン補正値Ｇ（ｘ，ｙ）を乗算することにより、最終的な画像データＦ_Ｏ（ｘ，ｙ）を得るようにして補正が行われる。なお、放射線画像検出器の１つの放射線検出素子から出力される実写画像データＦ（ｘ，ｙ）に基づいて生成される画像データＦ_Ｏ（ｘ，ｙ）は、撮影画像（すなわち放射線画像）における１画素分の画像データに相当するため、以下、放射線検出素子のことを画素と呼ぶ場合がある。
　　Ｆ_Ｏ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）－Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）　　…（１）

　このように、実写画像データの補正では、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）を得ることが必要となるため、放射線画像検出器に対してキャリブレーションを定期的に行い、経時的に特性が変動し得るオフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）を更新するのが一般的である。

　しかしながら、この方法では、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が一致していることが前提となるため、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が異なる場合は、特に温度依存性の高いオフセット補正値Ｏ（ｘ，ｙ）が適正な値からずれてしまい、最終的な画像データＦ_Ｏ（ｘ，ｙ）のＳＮ比を劣化させてしまうという問題がある。

　この問題を解決するために、放射線画像撮影ごとに、撮影の直前や直後に放射線を照射しない状態で各放射線検出素子からの出力値（以下、ダーク読取値Ｄ（ｘ，ｙ）という。また、ダーク読取値をダーク画像と呼ぶ場合がある。）を検出して、当該放射線画像撮影におけるオフセット補正値Ｏ（ｘ，ｙ）が算出される場合もある。これは、放射線画像撮影で実写画像データＦ（ｘ，ｙ）が得られた時点における放射線検出素子の温度特性とできるだけ同じ温度条件下でオフセット補正値Ｏ（ｘ，ｙ）を得るための処理である。

　しかし、実写画像データＦ（ｘ，ｙ）を取得する場合と同様に、ダーク読取値Ｄ（ｘ，ｙ）を取得する際にも、各種電気ノイズ、すなわち、フォトダイオードの暗電流ノイズ、ＴＦＴ過渡ノイズ、ＴＦＴサーマルノイズ、ＴＦＴリークノイズ、ＴＦＴから電荷を読み出すデータラインの寄生容量によって生じるサーマルノイズ、読み出し回路内部のアンプノイズ、Ａ／Ｄ変換によって生じる量子化ノイズなどが影響するため、仮に同一の温度条件下でダーク読取値Ｄ（ｘ，ｙ）を読み取ったとしても、ダーク読取値Ｄ（ｘ，ｙ）には、これら電気ノイズに起因する信号値のゆらぎ（ばらつき）が生じる。

　そのため、放射線画像撮影の直前や直後にダーク読取値Ｄ（ｘ，ｙ）を読み取るとしても、読み取られたダーク読取値Ｄ（ｘ，ｙ）の値がすなわちその温度条件等の撮影条件におけるオフセット補正値Ｏ（ｘ，ｙ）の真値であるとは必ずしも言えない。そこで、放射線画像撮影の直前や直後にダーク読取を複数回行って各ダーク読取値Ｄ（ｘ，ｙ）の平均値を算出し、その平均値をオフセット補正値Ｏ（ｘ，ｙ）として採用することがしばしば行われる（例えば特許文献１～３等参照）。

　これは、放射線画像撮影の直前や直後に、放射線画像撮影における放射線検出素子等の温度特性と同じ温度条件下で複数回読み出されたダーク読取値Ｄ（ｘ，ｙ）の平均値を算出すれば、ダーク読取値Ｄ（ｘ，ｙ）のゆらぎが緩和、若しくは相殺されるため、その平均値は、その撮影条件下におけるオフセット補正値Ｏ（ｘ，ｙ）の真値に等しく、或いは少なくともそれに近い値となるという考えに基づくものである。そして、その平均値であるオフセット補正値Ｏ（ｘ，ｙ）を用いて実写画像データＦ（ｘ，ｙ）を補正すれば、補正後の最終的な画像データＦ_Ｏ（ｘ，ｙ）のＳＮ比を良好なものとすることができる。

　一方、このようなオフセット補正をはじめとする補正処理等の画像処理は、従来から、放射線画像検出器等の撮影装置とは異なる画像処理プロセッサやコンソール等の処理装置で行われることも多い（例えば、特許文献４等参照）。また、近年、バッテリを内蔵し、ケーブルを介さずに無線方式により外部の処理装置等との間で実写画像データＦ（ｘ，ｙ）等の送受信を行うポータブルの放射線画像検出器が開発されている（例えば、特許文献５等参照）。

　また、放射線画像検出器では、製造過程において、正常な画素値を出力できない欠陥画素が発生することがある。例えば、欠陥画素の例としては、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）などがある。

　これら欠陥画素は、ＴＦＴ製造過程で発生する場合もあるし、ＴＦＴ上に光電変換素子を形成する際に発生する場合もある。放射線画像検出器では、照射した放射線量に比例して画素値が適正に変化することが求められるが、これら欠陥画素では、照射した放射線量に応じた正常な画素値の変化は期待できない。

　診断画像に対してこのような欠陥画素が発生すると、誤診の原因になったり、診断行為の妨げになるため、欠陥画素が存在する場合は、この欠陥画素を周囲の画素値を用いて補間処理により補正を行い、最終的な画像データＦ_Ｏ（ｘ，ｙ）中に欠陥画素が視認できないように処理することが求められる。

　このような補正処理を行うためには、予め欠陥画素の画素位置が分かっていなければならない。このような欠陥画素は、製造後の出荷検査や、定期的なキャリブレーション時に放射線照射を伴う検査を実施することで検出され、検出された欠陥画素の画素位置が欠陥画素マップに登録される（例えば、特許文献６参照）。この欠陥画素マップは放射線画像検出器毎に構築され、以後、個々の放射線画像検出器と対応付けて運用、管理される。

　前述したオフセット補正、ゲイン補正を行う際にも、この欠陥画素マップが参照され、欠陥画素の画素値をオフセット補正やゲイン補正に誤って使用しないように管理される。
米国特許第５４５２３３８号明細書米国特許第６２２２９０１号明細書米国特許第７０４１９５５号明細書特開平１１－１１３８８９号公報特開平７－１４０２５５号公報特開２００１－８１９８号公報

　欠陥画素の例として、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）を上げたが、これらの欠陥画素は非常に容易に検出することができる。

　しかしながら、上記のいずれにも属さないが、ある確率で異常な信号値を出力する画素（以下、異常画素と呼ぶ。）が存在する。これらの画素は欠陥画素マップに登録されていないため、ある確率で最終的な画像データＦ_Ｏ（ｘ，ｙ）中に異常な値を持つ画素を生成させ、診断の妨げとなる危険性があった。

　また、これまでの欠陥画素判別方法では、欠陥画素を判別する際、放射線の照射を行った画像情報を用いて欠陥画素を検出していたため、画像生成のために使用される時間の関係で、欠陥画素を判別するために使用する画像数が数枚からせいぜい十数枚に限定されてしまうという問題があった。従って、定常的に異常値を出力する欠陥画素は検出できてもある確率で発生する異常画素を捕まえられなかったり、標準偏差などの統計値を用いて欠陥画素の判別を行う際に、精度の良い統計値を求めることが困難であった。そのため、放射線の照射回数を増やそうとすると生産性が低下し、生産コストが増加するという問題点があった。

　また、このような欠陥画素や異常画素は、経時的に増加する傾向にある、そのため、工場出荷時に欠陥画素を登録しておけば安全ということではなく、ユーザー環境に設置後も、定期的にキャリブレーションを行い、新たな欠陥画素や異常画素が発見された場合には、その画素位置を欠陥画素マップに登録したり、異常画素として適切に処理することが好ましい。しかしながら、ユーザー先に設置後に行うキャリブレーション作業においても、定期的にユーザーに複数回の放射線照射を義務づけることになり、ユーザーの作業負荷が増大するという課題があった。

　また、キャリブレーション処理により欠陥画素や異常画素を判別する際に、放射線画像検出器内の各素子の温度変動を考慮した欠陥画素や異常画素を判別する方法が提案されていないため、温度変化の影響を受けずに精度良く欠陥画素や異常画素を判定することが困難であった。

　一方、バッテリを内蔵し、データ等を無線方式で外部装置に送信するポータブルの放射線画像検出器では、放射線画像撮影で撮影された実写画像データを外部装置に送信する際に比較的時間がかかりバッテリを消費してしまうため、放射線画像中で被写体の撮影位置等を確認して再撮影が必要か否かを判断するための画像として、放射線画像検出器の全画素分のデータである実写画像データの代わりに、実写画像データを間引いてデータ量を少なくした間引きデータを送信することができるようにシステムが構成されているものもある。

　このようなシステムでは、放射線画像検出器から短時間で間引きデータが送信され、コンソールの表示画面等に間引かれた画像データが瞬時に表示されるため、放射線画像検出器のバッテリの消耗の度合いが軽減されるとともに、放射線技師等が表示された間引きデータを確認して再撮影が必要か否かを即座に決定することができる。このように、間引きデータには画面上に迅速に表示されることが要請される。

　また、最終的には実写画像データの送信が必要となる為、間引きデータ送信後に、実写画像データを送信すると操作する技師にとって２度手間となるので、送信は実写画像データ１回のみとする運用も想定されるが、この運用においても、再撮影要否が迅速に確認できるように再撮影要否確認要画像を表示することが望まれる。

　更に、最終的には、異常・欠陥画素に対する補正を含み、各撮影部位に応じた諧調処理等を行った画像処理済の良好な診断画像を提供する必要がある。

　本発明は、以上のような課題を解決するためになされたものであり、再撮影要否確認用画像（一般的な運用では間引きデータ）を画面上に迅速に表示させることが可能な放射線画像生成システムを提供することを目的とする。また、最終的な診断画像データから異常画素や欠陥画素の影響を的確に排除することが可能な放射線画像生成システムを提供することを目的とする。

　前記の問題を解決するために、本発明の放射線画像生成システムは、
　２次元状に配置された複数の放射線検出素子と、
　放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
　前記実写画像データに基づいて間引きデータを作成可能な制御手段と、
　前記実写画像データ、前記ダーク読取値及び前記間引きデータを送信可能な通信手段と、
を有する放射線画像検出器と、
　前記放射線画像検出器から送信された前記実写画像データ及び前記間引きデータに対して画像処理を行う画像処理部と、
　前記実写画像データ、前記間引きデータに基づく画像、又はそれらに対して画像処理を行った後の画像を表示する表示部と、
を有するコンソールと、
　前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値に基づいて、前記放射線画像検出器における当該放射線画像撮影でのオフセット補正値を生成するオフセット補正値生成手段と、
　前記放射線画像検出器に対して放射線を照射する放射線発生装置と、
を備え、
　前記放射線画像検出器は、放射線画像撮影の前又は後に少なくとも１回行われたダーク読取における前記ダーク読取値と、放射線画像撮影に基づく前記実写画像データ及び前記実写画像データに基づく前記間引きデータのうちの少なくとも１つとを前記コンソールに送信し、
　前記コンソールは、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに基づく画像、或いは、前記間引きデータ、又は前記実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行って、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに対する画像処理後の画像、或いは、前記実写画像データに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行った画像処理後の画像、のうちの少なくとも１つを、前記表示部に表示させることを特徴とする。

　本発明のような方式の放射線画像生成システムによれば、施設の運用形態に応じた、再撮影要否確認の為の画像データを短時間で表示可能となる。また、再撮影要否確認後には、放射線検出素子の実写画像データを補正するためのオフセット補正値が、その放射線検出素子から出力されたダーク読取値と、その放射線検出素子に予め対応付けられた複数の放射線検出素子から出力されたダーク読取値とを用いて算出されるため、ダーク読取値の取得が放射線画像撮影の前又は後に１回程度行われたダーク読取であっても、十分に的確なオフセット補正値を得ることが可能となり、画素（放射線検出素子）ごとの特性ばらつきが補正されて、ＳＮ比が良好な最終的な診断画像データを生成することが可能となる。さらに、本発明特有の異常画素判定方法や欠陥画素判定方法を用いることで、最終的な画像データから異常画素や欠陥画素の影響を的確に排除して、良好な最終的な診断画像データを生成することが可能となる。

本実施形態に係る放射線画像検出器の外観構成を示す図である。放射線画像検出器のセンサパネル部及び読取部の構成に示す等価回路図である。放射線検出素子に割り当てられた番号を説明する図である。複数回のダーク読取ごとに放射線検出素子から出力されるダーク読取値を説明する図である。複数回のダーク読取で放射線検出素子から出力される複数のダーク読取値の時間的ゆらぎの分布を説明するグラフである。放射線検出素子に対応付ける複数の放射線検出素子のとり方の一例を示す図である。放射線検出素子が異なるとダーク読取値の時間的ゆらぎの分布は平均値や標準偏差が異なることを説明するグラフである。空間的平均値Ｗ（ｘ，ｙ）の分布の標準偏差が小さい分布となることを説明するグラフである。Ｄ（ｘ，ｙ）の時間的ゆらぎの分布とＤ（ｘ’，ｙ’）の時間的ゆらぎの分布の温度変化によるシフト度合いが異なるとブロードとなる温度補償変数の分布を示すグラフである。１回のダーク読取では放射線検出素子からオフセット補正値の真値を推定することが難しいことを説明するグラフである。複数回のダーク読取ごとに放射線検出素子から出力されるダーク読取値及び空間的平均値を説明する図である。複数回のダーク読取ごとに出力されるダーク読取値や空間的平均値が正規分布状に分布することを説明するグラフである。放射線検出素子の温度が時間的に変動することを説明するグラフである。放射線画像撮影の前又は後に行われる１回のダーク読取で放射線検出素子から出力されるダーク読取値や空間的平均値について推定される時間的ゆらぎの分布を説明するグラフである。複数回のダーク読取中に温度変動が起こるとダーク読取値等の分布はブロードになるが、温度補正済みダーク読取値の分布はブロードにならない様子を示すグラフである。一の放射線検出素子が正方領域の中心ではない位置に位置するように設定された領域の例を示す図である。一の放射線検出素子が正方領域の中心ではない位置に位置するように設定された領域の例を示す図である。一の放射線検出素子が正方領域の中心ではない位置に位置するように設定された領域の例を示す図である。所定個数の放射線検出素子ごとに１個ずつ接続される読み出しＩＣを説明する図である。領域に属する各放射線検出素子からの電気信号が隣接する２つの読み出しＩＣで別々に読み出される場合を説明する図である。通常の仕方で複数の放射線検出素子が選択される場合を説明する図である。通常の仕方で複数の放射線検出素子が選択される場合を説明する図である。隣接する読み出しＩＣの境界の近傍をセンサパネル部の周縁部分と同様に扱う場合を説明する図である。隣接する読み出しＩＣの境界の近傍をセンサパネル部の周縁部分と同様に扱う場合を説明する図である。放射線検出素子からのダーク読取値のゆらぎの分布の標準偏差が平均的な分布とそれよりも格段に大きい分布を表すグラフである。放射線検出素子からのダーク読取値のゆらぎの分布の平均値が閾値よりも大きい場合を表すグラフである。差分のゆらぎの分布の標準偏差が平均的な分布を表すグラフである。差分のゆらぎの分布の標準偏差が図２２に示した場合よりも大きい分布を表すグラフである。注目画素が放射線検出素子（４，４）であり、それに対応付けられた７×７個の放射線検出素子うち放射線検出素子（６，６）が欠陥画素である場合の例を説明する図である。図２４の例で領域外の位置に存在する放射線検出素子で欠陥画素を置換する例を示す図である。図２４の例で放射線検出素子（７，６）が欠陥画素の場合に領域外の位置に存在する放射線検出素子で欠陥画素を補間する例を説明する図である。実写画像データの例を示す写真である。図２７Ａの実施画像データに基づいて作成された間引きデータの例を示す写真である。放射線画像検出器を用いた画像生成の処理手順を示すフローチャートである。本実施形態に係る欠陥画像判定システムを含む放射線画像生成システムの構成を示す図である。

　以下、本発明に係る放射線画像生成システムの実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。

［放射線画像検出器の基本的な構成］
　まず、本発明の放射線画像生成システムに用いられる放射線画像検出器の基本的な構成について説明する。

　放射線画像検出器（フラットパネルディテクタ）１は、図１に示すように、内部を保護する筐体２を備えており、筐体２の放射線入射面Ｘの内側には、照射された放射線を光に変換する図示しないシンチレータ層が形成されている。シンチレータ層は、例えばＣｓＩ：ＴｌやＧｄ_２Ｏ_２Ｓ：Ｔｂ、ＺｎＳ：Ａｇ等の母体内に発光中心物質が付活された蛍光体を用いて形成されたものを用いることができる。

　シンチレータ層の放射線が入射する側の面とは反対側の面側には、図２の等価回路図に示すように、放射線検出素子として、シンチレータ層から出力された光を電気信号に変換する複数のフォトダイオード１４が２次元状に配置されたセンサパネル部４が設けられている。１つのフォトダイオード１４から出力される電荷（信号値）は１つの画素を形成する。また、後で詳しく説明するように、各フォトダイオード１４にはそれぞれ信号読み出し用のスイッチ素子であるＴＦＴ１５が接続されている。

　なお、以下では、上記のようにシンチレータ層で放射線を光に変換してフォトダイオード等の光電変換素子で検出する、いわゆる間接方式の放射線画像検出器１を用いる場合について説明するが、放射線画像検出器は、この他にも、前述したシンチレータ層を介さず検出素子で入射した放射線を直接電気信号に変換する、いわゆる直接方式の放射線画像検出器を用いることも可能であり、その場合にも本発明を適用することが可能である。

　また、以下、これらの各方式の放射線画像検出器に用いられる検出素子を、あわせて放射線検出素子という。すなわち、放射線検出素子は、例えば本実施形態のような間接方式の放射線画像検出器１では１個のフォトダイオード１４、それに接続されたＴＦＴ１５及びシンチレータ層の当該フォトダイオード１４に対応する部分で構成され、例えば直接方式の放射線画像検出器では検出素子とそれに接続されたＴＦＴ等のスイッチ素子とで構成される。

　放射線画像検出器１には、バッテリ２１（図２参照）が内蔵されている。また、図１に示すように、本実施形態では、放射線画像検出器１の筐体２の側面部分に設けられたバッテリ交換用の蓋部材１０には、無線通信手段であるアンテナ装置３が埋め込まれて設けられている。さらに、筐体２の側面部分には、放射線画像検出器１の電源スイッチ１１や各種の操作状況等を表示するインジケータ１２、ユーザーからの指示を入力するための入力ボタン２３等が設けられている。

　図２に示すように、センサパネル部４の近傍には、センサパネル部４の各放射線検出素子の出力値を読み取る読取部５が設けられている。読取部５は、マイクロコンピュータ等からなる制御手段６や、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）、フラッシュメモリ等からなる記憶手段７、走査駆動回路８、読み出し回路９等で構成されている。

　センサパネル部４に２次元状に配置された複数の放射線検出素子には、図３に示すように、それぞれセンサパネル部４における放射線検出素子の行方向の位置ｘと列方向の位置ｙとを各成分とする座標（ｘ，ｙ）が当該放射線検出素子の番号（ｘ，ｙ）として予め割り当てられている。以下、個々の放射線検出素子を特定する場合には、放射線検出素子（ｘ，ｙ）という。

　なお、図３では、８×１６個の放射線検出素子（ｘ，ｙ）が記載されているが、これは簡略化して表現したものであり、実際にはさらに多くの放射線検出素子（ｘ，ｙ）が２次元状に配置されていてそれぞれ番号が割り当てられている。また、座標（ｘ，ｙ）（すなわち放射線検出素子の番号（ｘ，ｙ））を画素の番号（ｘ，ｙ）もしくは画素位置（ｘ．ｙ）と標記する場合もある。

　センサパネル部４及び読取部５の構成についてさらに説明すると、図２の等価回路図に示すように、センサパネル部４の各放射線検出素子（ｘ，ｙ）の一方の電極にはそれぞれ信号読み出し用のスイッチ素子であるＴＦＴ１５のソース電極が接続されている。また、各放射線検出素子（ｘ，ｙ）の他方の電極にはバイアス線Ｌｂが接続されており、バイアス線Ｌｂはバイアス電源１６に接続されていて、バイアス電源１６から各放射線検出素子（ｘ，ｙ）にバイアス電圧が印加されるようになっている。

　各ＴＦＴ１５のゲート電極はそれぞれ走査駆動回路８から延びる走査線Ｌｌに接続されており、各ＴＦＴ１５のドレイン電極はそれぞれ信号線Ｌｒに接続されている。各信号線Ｌｒは、それぞれ読み出し回路９内の増幅回路１７に接続されており、各増幅回路１７の出力線はそれぞれサンプルホールド回路１８を経てアナログマルチプレクサ１９に接続されている。また、アナログマルチプレクサ１９にはＡ／Ｄ変換器２０が接続されており、読み出し回路９はＡ／Ｄ変換器２０を介して制御手段６に接続されている。制御手段６には、記憶手段７が接続されている。

　放射線画像検出器１では、図示しない被写体を撮影する放射線画像撮影において、被写体を透過した放射線がシンチレータ層に入射すると、シンチレータ層からセンサパネル部４に光が照射され、光の照射を受けた量に応じて、放射線検出素子（ｘ，ｙ）内に電荷が蓄積される。

　そして、放射線画像撮影を終了し、放射線画像検出器１から実写画像データを電気信号として読み出す際には、走査線ＬｌからＴＦＴ１５のゲート電極に読み出し電圧を印加して各ＴＦＴ１５のゲートを開き、放射線検出素子（ｘ，ｙ）からＴＦＴ１５を介して蓄積された電荷を電気信号として信号線Ｌｒに取り出す。そして、電気信号を増幅回路１７で増幅する等して、アナログマルチプレクサ１９から順次Ａ／Ｄ変換器２０を介して制御手段６に出力する。

　制御手段６は、放射線検出素子（ｘ，ｙ）から出力され増幅された電気信号を、前述した放射線検出素子（すなわち画素）の番号（ｘ，ｙ）と対応付け、実写画像データＦ（ｘ，ｙ）として記憶手段７に保存するようになっている。

　ＴＦＴ１５に読み出し電圧を印加する走査線Ｌｌを順次走査して上記の読み出し処理を走査線Ｌｌごとに行うことで、センサパネル部４の全放射線検出素子（ｘ，ｙ）から電気信号をそれぞれ読み出し、各電気信号にそれぞれ画素の番号（ｘ，ｙ）を対応付け、各実写画像データＦ（ｘ，ｙ）として記憶手段７に保存するようになっている。

　このように、放射線画像検出器１では、複数の放射線検出素子（ｘ，ｙ）が２次元状に配置されたセンサパネル部４や、制御手段６や走査駆動回路８、読み出し回路９等で構成された読取部５等で、実写画像データＦ（ｘ，ｙ）を取得する画像データ取得手段が形成されている。

　また、放射線画像検出器１の画像データ取得手段では、上記のような実写画像データＦ（ｘ，ｙ）の取得のみならず、ダーク読取も行われるようになっている。

　ダーク読取では、放射線画像検出器１の複数の放射線検出素子（ｘ，ｙ）をリセットして電荷を放出させた後、各ＴＦＴ１５のゲートを閉じて、放射線画像検出器１を放射線が照射されない状態に保つ。そして、所定時間経過後（通常は、放射線照射時の放射線検出素子（ｘ，ｙ）への電荷蓄積時間と同じ時間経過後）、走査線ＬｌからＴＦＴ１５のゲート電極に読み出し電圧を印加して各ＴＦＴ１５のゲートを開き、各放射線検出素子（ｘ，ｙ）に溜まった電荷（暗電荷）を信号線Ｌｒに取り出し、上記と同様に、出力値を増幅回路１７で増幅する等してアナログマルチプレクサ１９から順次Ａ／Ｄ変換器２０を介して制御手段６に出力する。

　このようにして、放射線が曝射されない各放射線検出素子（ｘ，ｙ）から出力される電気信号がダーク読取値である。制御手段６は、各放射線検出素子（ｘ，ｙ）から出力された各電気信号を各画素の番号（ｘ，ｙ）と対応付けてダーク読取値Ｄ（ｘ，ｙ）として記憶手段７に保存するようになっている。なお、ＴＦＴ１５に読み出し電圧を印加する走査線Ｌｌを順次走査して、全放射線検出素子（ｘ，ｙ）からダーク読取値Ｄ（ｘ，ｙ）が読み出される。

　なお、前述したように、放射線画像検出器１では、各放射線検出素子（ｘ，ｙ）からの出力値特性変動把握のために各放射線検出素子（ｘ，ｙ）のキャリブレーションが行われる。通常の場合、キャリブレーションでは、被写体が介在しない状態で放射線画像検出器１に放射線を照射し、各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）ごとにゲイン補正値が算出される。

　また、キャリブレーションでは複数回ダーク読取も行われ、ダーク読取では、放射線画像検出器１に放射線を照射せずに各放射線検出素子（ｘ，ｙ）からダーク読取値Ｄ（ｘ，ｙ）が出力され、出力されたダーク読取値Ｄ（ｘ，ｙ）に基づいて各放射線検出素子（ｘ，ｙ）ごとにオフセット補正値が算出される。

　さらに、放射線画像撮影で実写画像データＦ（ｘ，ｙ）が得られた時点における放射線検出素子（ｘ，ｙ）の温度特性等の撮影条件と同じ撮影条件下でオフセット補正値を得るために、放射線画像撮影ごとに、撮影の直前又は直後に放射線を照射しない状態で各放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）を複数回検出して、その平均値を当該放射線画像撮影におけるオフセット補正値として算出する場合もあることは前述した通りである。

［放射線画像検出器の各放射線検出素子のオフセット補正値取得の原理］
　さて、本発明に係る放射線画像生成システムの説明の前に、放射線画像検出器の各放射線検出素子（各画素）に対するオフセット補正値を取得する際の原理について説明する。

　各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）を上記（１）式に従って補正するために、オフセット補正値（以下、Ｏ（ｘ，ｙ）と表す。）の真値が必要となる。しかし、前述したように、ダーク読取を行って実際に得られるデータは、前述した電気ノイズ等の影響のためにゆらぐ（変動する、もしくは誤差を持つ）ダーク読取値Ｄ（ｘ，ｙ）であり、オフセット補正値Ｏ（ｘ，ｙ）の真値を直接的には得ることはできない。

　そこで、従来は、前述したように、放射線画像撮影の直前や直後に、時間的に連続して複数回のダーク読取を行って一の放射線検出素子（ｘ，ｙ）から出力される複数のダーク読取値Ｄ（ｘ，ｙ）の平均値を算出してゆらぎを緩和、若しくは、相殺し、その平均値をオフセット補正値Ｏ（ｘ，ｙ）とする手法を採用した。

　具体的には、ダーク読取をＫ回行うとし、図４に示すように、ｋ回目（ｋ＝１～Ｋ）のダーク読取で当該放射線検出素子（ｘ，ｙ）から出力されるダーク読取値をＤ_ｋ（ｘ，ｙ）とすると、例えば図５に示すように、各ダーク読取で当該放射線検出素子（ｘ，ｙ）から出力される各ダーク読取値Ｄ_ｋ（ｘ，ｙ）は平均値Ｄ_ｋave（ｘ，ｙ）、標準偏差σ_Ｄｋ（ｘ，ｙ）の正規分布を形成するようにゆらぐことが知られている。これは、当該一の放射線検出素子（ｘ，ｙ）に対して繰り返しダーク読取を行うと、取得されるダーク読取値の値が、同じ画素位置であるにも関わらず読取り毎にゆらぐ（変動する、もしくは誤差を持つ）ことから、このゆらぎを「時間的ゆらぎ」と定義する。ゆらぎの大きさを表す統計的指標には、通常、ゆらぎの分布の標準偏差（ダーク読取値Ｄ_ｋ（ｘ，ｙ）の場合は標準偏差σ_Ｄｋ（ｘ，ｙ））が使用される。

　また、平均値Ｄ_ｋave（ｘ，ｙ）は各ダーク読取値Ｄ_ｋ（ｘ，ｙ）の時間軸上の平均値であるから、この平均値を「時間的平均値」と定義する。この平均値Ｄ_ｋave（ｘ，ｙ）をこの放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）として使用すれば、各ダーク読取値Ｄ_ｋ（ｘ，ｙ）が持つ値の時間的ゆらぎを緩和、若しくは相殺することができるので、一般的には、各ダーク読取値Ｄ_ｋ（ｘ，ｙ）の時間的平均値をオフセット補正値Ｏ（ｘ，ｙ）として使用する場合が多い。オフセット補正値Ｏ（ｘ，ｙ）を下記（２）式に示す。

　ダーク読取値Ｄ_ｋ（ｘ，ｙ）は、標準偏差σ_Ｄｋ（ｘ，ｙ）でゆらぐが、これをＫ回で平均化すると、ゆらぎの大きさは、一般に１／√Ｋになる。すなわち、オフセット補正値Ｏ（ｘ，ｙ）は標準偏差{１／√Ｋ・σ_Ｄｋ（ｘ，ｙ）} のゆらぎを持つ分布となる。従って、上記（２）式は、ダーク読取値Ｄ_ｋ（ｘ，ｙ）そのものをオフセット補正値として使用するよりも、ゆらぎの大きさが１／√Ｋの値をオフセット補正値として使用することを示している。

　しかしながら、このようにしてオフセット補正値Ｏ（ｘ，ｙ）を算出すると、ダーク読取やダーク読取値の送信を複数回行わなければならないため、電力を消費してしまう等の問題があることは前述した通りである。

　また、この問題を回避するために、キャリブレーション時に複数のダーク読取を行っておき、上記（２）式を用いてあらかじめ算出されたオフセット補正値Ｏ（ｘ，ｙ）を使用すると、キャリブレーションを行った時の放射線画像検出器内の各素子の温度と、放射線画像撮影を実施した時の放射線画像検出器内の各素子の温度が異なる場合は、オフセット補正値Ｏ（ｘ，ｙ）が適正な値からずれてしまうという問題があることも前述した通りである。

　そこで、オフセット補正方法の実施例を通して、キャリブレーションの手法や温度補償の手法を交えて説明する。

　この説明の中で取り上げられるキャリブレーションの手法や温度補償の手法は、本発明における欠陥画素判定方法でも使用するため、その基本的考え方については、以下の説明の中で論ずるものとする。また、オフセット補正方法の説明には、放射線画像撮影の直前又は直後に行うダーク読取を少なくとも１回だけ行って、たとえダーク読み取りが１回だけであっても、精度の良いオフセット補正値Ｏ（ｘ，ｙ）を算出する手法を用いる。ただし、本発明における欠陥画素判定方法は、以下に説明するオフセット補正方法に限定されるものではない。

　本手法（放射線画像撮影の直前又は直後に行うダーク読取を少なくとも１回だけ行って、たとえダーク読み取りが１回だけであっても、精度の良いオフセット補正値Ｏ（ｘ，ｙ）を算出する手法）の基本的な考え方は以下の通りである。

　まず、キャリブレーション時に、各放射線検出素子（ｘ，ｙ）に対応する各画素位置（ｘ，ｙ）に対して、温度変化のない所定の温度条件下で複数のダーク画像ｄ_ｋ（ｘ，ｙ）（以下、上記のように放射線画像撮影の直前や直後に行われる１回のダーク読取で得られるダーク読取値Ｄ（ｘ，ｙ）と区別するためにダーク読取値ｄ_ｋ（ｘ，ｙ）と表す。）を取得し、取得された複数のダーク画像ｄ_ｋ（ｘ，ｙ）を用いて、オフセット補正値δ（ｘ，ｙ）（以下、最終的に得られるオフセット補正値Ｏ（ｘ，ｙ）と区別するためにオフセット補正値δ（ｘ，ｙ）と表す。）を求めておく。なお、ダーク画像ｄ_ｋ（ｘ，ｙ）を取得している間は、各画素位置（ｘ，ｙ）の温度特性は同一であると仮定する。

　この時、そのオフセット補正値δ（ｘ，ｙ）を求めた時の各画素位置（ｘ，ｙ）の信号値の温度変化に伴う変化量を代弁する温度補償変数を定義し、これを各画素位置（ｘ，ｙ）毎に計算して記憶しておく。

　ここで、温度補償変数の位置づけについて説明する。温度補償変数は、各画素位置（ｘ，ｙ）における放射線検出素子（ｘ，ｙ）の温度そのものを測定することを目的としていない。知りたいのは、各放射線検出素子（ｘ，ｙ）の温度そのものではなく、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したかであるため（例えば、各放射線検出素子（ｘ，ｙ）の温度そのものが変化しても、その放射線検出素子（ｘ，ｙ）に対応する画素の信号値が変化していなければ、その放射線検出素子（ｘ，ｙ）の温度変化はなかったものと見なしても問題ない。）、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したか、その変化量を定量的に示唆する変数であることが望ましい。

　次に、放射線画像撮影の直前や直後に、ダーク読取を１回だけ行って、ダーク読取値Ｄ（ｘ，ｙ）を取得する。この時、ダーク読取値Ｄ（ｘ，ｙ）に対しても、上述と同様の温度補償変数を各画素位置（ｘ，ｙ）毎に計算して記憶する。

　次に、オフセット補正値δ（ｘ，ｙ）を求めた時の各画素位置（ｘ，ｙ）における温度補償変数とダーク読取値Ｄ（ｘ，ｙ）の各画素位置（ｘ，ｙ）における温度補償変数を比較し、その結果に基づいて、あらかじめ求めてあったオフセット補正値δ（ｘ，ｙ）を実際に放射線撮影された時の温度環境下におけるオフセット補正値Ｏ（ｘ，ｙ）に変換して、最終的な画像データＦ_Ｏ（ｘ，ｙ）を求める演算に使用する。

　通常は、過去のキャリブレーション時に取得したオフセット補正値Ｏ（ｘ，ｙ）を使用すると、温度変動の影響を受けるため、放射線画像撮影の直前や直後に取得したダーク読取値Ｄ（ｘ，ｙ）そのものをオフセット補正値Ｏ（ｘ，ｙ）として使用することが多い（すなわち、Ｏ（ｘ，ｙ）＝Ｄ（ｘ，ｙ）とする）。一方、本手法では、放射線画像撮影毎に取得したダーク読取値Ｄ（ｘ，ｙ）はオフセット補正値そのものとして直接的に使用せず、放射線撮影された時の温度環境下における各画素位置の温度補償変数の算出に使用する点が従来手法と異なっている。そして、あらかじめ複数のダーク画像を用いて算出されたオフセット補正値δ（ｘ，ｙ）をそのまま使用するのではなく、上述の温度補償変数を用いて、実際に放射線画像撮影が実施された時点への温度補正を実施した後にオフセット補正値Ｏ（ｘ，ｙ）として使用するという点が、従来手法にはない点である。

　この方法の利点は以下の通りである。

　まず第１に、最終的な画像データＦ_Ｏ（ｘ，ｙ）を求めるのに使用されるオフセット補正値Ｏ（ｘ，ｙ）を算出するための基準となるオフセット補正値δ（ｘ，ｙ）が、温度変化がないと見なせる所定の温度条件下の複数のダーク画像ｄ_ｋ（ｘ，ｙ）から算出されているため、電気ノイズ等の誤差成分（信号値の時間的なゆらぎ）が相殺された、真のオフセット補正値に近い値になっていることである。従って、放射線画像撮影時に取得された１回のダーク読取値Ｄ（ｘ，ｙ）そのものをオフセット補正値として使用する場合に比べて、補正誤差が少なく、ＳＮ比が良好な最終的な画像データＦ_Ｏ（ｘ，ｙ）を求めることができる。

　第２に、各画素位置（ｘ，ｙ）毎に温度補償変数を計算することで、温度変化による画素値の変動を各画素位置（ｘ，ｙ）毎に個別に補正しているため、各画素位置（ｘ，ｙ）の信号値がどのように温度変化の影響を受けても、これを補正し、良好なＳＮ比の最終的な画像データＦ_Ｏ（ｘ，ｙ）が得られることである。そして、この画素単位の温度補正が、たった１回のダーク画像の取得で行えることである。

　ここで重要になるのが、各画素位置（ｘ，ｙ）毎の温度補償変数に何を使用するかである。各画素位置（ｘ，ｙ）にあらかじめ温度センサを内蔵できれば良いが、この方法は現実的ではない。また、現実的な個数の温度センサを放射線画像検出器１に内蔵しておくことも考えられるが、これでは、各画素位置（ｘ，ｙ）毎に正確な温度補正を実現することができない。また、求めたいのは、放射線検出素子（ｘ，ｙ）の温度そのものではなく、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって、各画素位置（ｘ，ｙ）における信号値がどのように変化したかである。

　上記より、温度補償変数の要件は、以下の４項目であると考えられる。
要件１）各画素単位（もしくは、画素に近い小領域単位毎に）に求められる変数であること。
要件２）温度変化に伴う信号値の変化を把握できる変数であること。
要件３）ダーク画像Ｄ（ｘ，ｙ）から求められる変数であること。
要件４）時間的なゆらぎが少ない変数であること。

　この４項目の要件を満足する温度補償変数として、本件発明者は、ダーク読取値Ｄ（ｘ，ｙ）の注目画素（ｘ，ｙ）に対する周辺画素の画素値の平均値（もしくはメディアン値等）に注目した。

　放射線検出素子（ｘ，ｙ）や電気回路などの各素子が温度変化の影響を受けた結果は、ダーク読取値Ｄ（ｘ，ｙ）に反映される。我々が知りたいのは、放射線検出素子（ｘ，ｙ）や電気回路などの各素子の温度変化によって各画素の信号値がどのように変化したかであるので、温度補償変数にダーク読取値Ｄ（ｘ，ｙ）を用いるのは理にかなっており、かつ、要件１～３を満足している。

　次に選択した温度補償変数が要件４を満足するか否かについて、図６を用いて説明する。

　図６に示すように、例えば複数の放射線検出素子（ｘ，ｙ）として当該一の放射線検出素子（ｘ，ｙ）すなわち、画素位置（ｘ，ｙ）を中心とする５×５の正方領域内に存在する複数の画素位置（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）の信号値を温度補償変数算出に使用する。なお、当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ，ｙ）の選択の仕方等については後で考察する。

　しかし、当該一の放射線検出素子（ｘ，ｙ）を含む２５個の各放射線検出素子（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）から出力されるダーク読取値Ｄ（ｘ－２，ｙ－２）～Ｄ（ｘ＋２，ｙ＋２）の信号値の時間的なゆらぎの分布は、全てが図５に示したような分布であるとは限らず、実際には、例えば図７に示すように、平均値μ_Ｄや標準偏差σ_Ｄが異なる分布となるのが一般的である。

　本手法では、下記（３）式に示すように、当該１回のダーク読取でそれらの放射線検出素子（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）から出力された全ダーク読取値Ｄ（ｘ－２，ｙ－２）～Ｄ（ｘ＋２，ｙ＋２）の平均値Ｗ（ｘ，ｙ）を算出して、これを温度補償変数として採用する。なお、この場合の平均値（温度補償変数）Ｗ（ｘ，ｙ）はダーク読取値Ｄ（ｘ，ｙ）の画素位置（ｘ，ｙ）に対する空間的平均値である。（３）式中のＮは、この場合２５（＝５×５）である。

　ここで、この温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）が要件４）を満足するか否かについて考察する。

　ダーク読取値Ｄ（ｘ－２，ｙ－２）～Ｄ（ｘ＋２，ｙ＋２）は、それぞれが独立変数であり、かつそのゆらぎは時間的にも空間的にもランダムであるから、ダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの時間的平均値をμ_Ｄ（ｘ，ｙ）、時間的ゆらぎの分布の標準偏差をσ_Ｄ（ｘ，ｙ）と表し、温度補償変数Ｗ（ｘ，ｙ）の時間的ゆらぎの時間的平均値をμ_Ｗ（ｘ，ｙ）、時間的ゆらぎの分布の標準偏差をσ_Ｗ（ｘ，ｙ）と表すと、

となる（図８参照）。

　ここで、σ_Ｄ（ｘ，ｙ）の中に極めて大きな値（異常値）を持つ欠陥画素や多くの異常画素が含まれていなければ、通常、
　　σ_Ｗ（ｘ，ｙ）＜σ_Ｄ（ｘ，ｙ）
の関係式が成立する。

　例えば、σ_Ｄ（ｘ－２，ｙ－２）～σ_Ｄ（ｘ＋２，ｙ＋２）のそれぞれの値は異なるが、平均を取るとσ_Ｄ（ｘ，ｙ）にほぼ等しいと仮定すると、上記（５）式より、

となり、σ_Ｗ（ｘ，ｙ）は、σ_Ｄ（ｘ，ｙ）の値の１／√Ｎとなる。Ｎ＝２５の場合は、１／√Ｎ＝１／５であるから、温度補償変数Ｗ（ｘ，ｙ）は、時間的ゆらぎの分布の広がり（標準偏差）が、ダーク読取値Ｄ（ｘ，ｙ）にくらべて１／５の分布となる。すなわち、図８に示すように、広がりが狭い（標準偏差が小さい）分布となり、要件４）を満足する。

　極端な例として、画素位置（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）の内、注目画素である（ｘ，ｙ）以外の（Ｎ?１）の画素位置におけるダーク読取値時間的ゆらぎの標準偏差値が、注目画素位置におけるダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの標準偏差値σ_Ｄ（ｘ，ｙ）に対して、Ｎ倍の大きさを持ったと仮定しても、Ｎ＞２であれば、

となり、
　　σ_Ｗ（ｘ，ｙ）＜σ_Ｄ（ｘ，ｙ）
が成立する。

　Ｎ＝２５の場合、
　　ＳＱＲＴ［｛１＋（Ｎ－１）・Ｎ｝／Ｎ^２］
　　　　　　　　＝ＳＱＲＴ｛（１＋２４×２５）／（２５×２５）｝
　　　　　　　　≒０．９８１
となるから、
　　σ_Ｗ（ｘ，ｙ）≒０．９８１・σ_Ｄ（ｘ，ｙ）＜σ_Ｄ（ｘ，ｙ）
通常、σ_Ｄ（ｘ?２，ｙ?２）～σ_Ｄ（ｘ＋２，ｙ＋２）のそれぞれの値は異なるが、平均を取るとσ_Ｄ（ｘ，ｙ）にほぼ等しい値を取ることから、上記方法は要件４）を満足すると言って良い。

　温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布は、当該一の放射線検出素子（ｘ，ｙ）を中心とする正方領域内の複数の放射線検出素子１４（以下、所定の領域Ｒに属する放射線検出素子１４という意味で放射線検出素子１４_Ｒと表記する。）の数が多くなればなるほど、すなわちＮの値が大きい程、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの標準偏差σ_Ｗ（ｘ，ｙ）が小さくなる。これは、すなわち、σ_Ｗ（ｘ，ｙ）が十分に小さい値を取るとき、
　　Ｗ（ｘ，ｙ）≒μ_Ｗ（ｘ，ｙ）　　…（６）
と近似できることを意味している。

　すなわち、放射線画像撮影の直前又は直後に行われる１回のダーク読取で放射線検出素子１４_Ｒから出力されるダーク読取値Ｄ（ｘ－２，ｙ－２）～Ｄ（ｘ＋２，ｙ＋２）（以下、所定の領域Ｒに属する放射線検出素子１４_Ｒから出力されるダーク読取値という意味でダーク読取値Ｄ_Ｒと表記する。）は、図７に例示したように各時間的ゆらぎの分布内でその時間的ゆらぎの平均値μ_Ｄ（ｘ－２，ｙ－２）～μ_Ｄ（ｘ＋２，ｙ＋２）を中心としてゆらぐが、上記（３）式に従って算出される空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）は、空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の分布の時間的平均値μ_Ｗ（ｘ，ｙ）にほぼ等しい値が算出される。

　従って、当該一の放射線検出素子（ｘ，ｙ）を含む各放射線検出素子１４_Ｒから出力されるダーク読取値Ｄ_Ｒの時間的ゆらぎの分布の平均値μ_Ｄや標準偏差σ_Ｄが異なるとしても（図７参照）、それらの空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出すれば、各放射線検出素子１４_Ｒから出力された各ダーク読取値Ｄ_Ｒの時間的ゆらぎが有効に相殺されて数値の変動が非常に小さくなり、時間的ゆらぎの分布の平均値μ_Ｗ（ｘ，ｙ）にほぼ等しい空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出することができる。

　言い換えれば、複数のダーク読取を実施しなければ求められない時間的ゆらぎの平均値を、１回のダーク読取で得られたダーク読取値Ｄ（ｘ，ｙ）の空間的平均値Ｗ（ｘ，ｙ）で、少ない誤差をもって代用できるということである。

　ここで、一の放射線検出素子（ｘ，ｙ）について温度補償変数Ｗ（ｘ，ｙ）を算出する際の当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）の選択の仕方について説明する。

　放射線検出素子（ｘ，ｙ）では、通常、放射線検出素子（ｘ，ｙ）や、この信号値を読み出すための読み出し回路の素子の温度が低いと当該放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）は小さな値となり、放射線検出素子（ｘ，ｙ）の温度Ｔが上昇するにつれて、出力されるダーク読取値Ｄ（ｘ，ｙ）の値も大きくなっていくことが知られている。そのため、例えば、図７に示したような各放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μ_Ｄ（ｘ，ｙ）も、温度が高くなると、同図のグラフ上で右方向にシフトしていく。

　また、例えば放射線画像検出器１の電源がＯＮされるとセンサパネル部４の温度が次第に上昇していくが、センサパネル部４上の各放射線検出素子（ｘ，ｙ）や読み出し回路の素子の温度変動は一様ではなく、例えば、周囲に多数の放射線検出素子（ｘ´，ｙ´）が存在するセンサパネル部４の中央部分の放射線検出素子（ｘ，ｙ）と、センサパネル部４の周縁部分の放射線検出素子（ｘ，ｙ）とでは、温度の変動の仕方が異なるし、読み出し回路を構成する読み出しＩＣが異なれば、温度の変動の仕方が異なる。また、同じ読み出しＩＣ内であっても、増幅回路１７が異なれば、温度の変動の仕方が異なる。

　このような状況において、一の放射線検出素子（ｘ，ｙ）について空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出する際の当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）として、例えばセンサパネル部４上で当該一の放射線検出素子（ｘ，ｙ）から遠く離れた放射線検出素子（ｘ´，ｙ´）を選ぶと、その温度変動が当該一の放射線検出素子（ｘ，ｙ）の温度変動と異なることが予測される。

　そのため、例えば、図７に示したような各放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の時間的ゆらぎの分布の温度による左右方向へのシフトの度合が、当該一の放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布のシフトの度合と異なるものとなり、グラフ上の相対的な位置がそれぞれの画素位置の温度に依存して変動する。

　このように、温度によって放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布と他の放射線検出素子（ｘ´，ｙ´）の時間的ゆらぎの分布とのグラフ上の位置が、それぞれの画素位置の温度に依存して相対的にずれてしまうと、図９に示すように、それらの空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の分布はブロードになり、図８に示したような標準偏差σ_Ｗ（ｘ，ｙ）が小さい分布が得られなくなる。そのため、上記（３）式に従って算出される空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）が、空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）の時間的平均値μ_Ｗ（ｘ，ｙ）にほぼ等しい値となるという有益な効果を活用することができなくなってしまう。

　そこで、本手法では、放射線画像撮影の直前又は直後に行われる１回のダーク読取の結果から一の放射線検出素子（ｘ，ｙ）について空間的平均値（温度補償変数）Ｗ（ｘ，ｙ）を算出するための当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）として、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動する放射線検出素子（ｘ´，ｙ´）が予め選択されて、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けるようになっている。

　その際、センサパネル部４上の全ての放射線検出素子（ｘ，ｙ）について同じ選択の仕方とされる必要はなく、放射線検出素子ごとに選択の仕方が異なってもよい。また、選択の仕方としては、例えば上記のように当該一の放射線検出素子（ｘ，ｙ）を中心とするｎ×ｎの正方領域内に存在する複数の放射線検出素子（ｘ，ｙ）_Ｒを用いることが可能であるが、後述する図１６Ａ～図１６Ｃや図１９Ｃ、図１９Ｄに示すように、正方領域は必ずしも当該一の放射線検出素子（ｘ，ｙ）が中心に存在するように設定されなくてもよい。

　また、正方領域のサイズは、本実施例では５×５や７×７の例を用いたが（７×７の例は後述する図２４～図２６の説明文を参照）、このサイズに限定するものではない。１×１よりも大きいサイズであれば、２×２や３×３であっても良いし、９×９もしくはそれ以上のサイズであっても良い。また、正方領域の代わりに長方形状の領域を設定したり、不定形状の領域とすることも可能である。

　さらに、当該一の放射線検出素子（ｘ，ｙ）以外の放射線検出素子（ｘ´，ｙ´）は当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動する放射線検出素子であればよく、上記のような当該一の放射線検出素子（ｘ，ｙ）を含む空間的に連続する領域内に存在する放射線検出素子（ｘ，ｙ）_Ｒでなく、センサパネル部４上に点在する放射線検出素子（ｘ，ｙ）であってもよい。なお、これらの同一の温度特性を示す検出素子群は、工場検査等において、予め把握することができる。また、施設設置後にも再設定することが可能である。

　次に、上記の温度補償変数として使用する空間的平均値Ｗ（ｘ，ｙ）と当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値との関係について説明する。

　放射線画像撮影の直前又は直後に１回だけ行われるダーク読取では、図１０に示したように、一の放射線検出素子（ｘ，ｙ）からは、時間的なゆらぎを含む１個のダーク読取値Ｄ（ｘ，ｙ）しか得られず、オフセット補正値Ｏ（ｘ，ｙ）の真値は分からない。また、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）から出力される１個ずつの各ダーク読取値Ｄ（ｘ－２，ｙ－２）～Ｄ（ｘ＋２，ｙ＋２）も時間的なゆらぎを含むものであり、それらに基づいて直接的に当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値を得ることはできない。

　そこで、本手法では、以下で詳しく説明するように、過去の（当回の被写体の放射線画像撮影以前に行われた）キャリブレーション時にダーク読取を複数回行って当該一の放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄ（ｘ，ｙ）を得ておき、まず、一の放射線検出素子（ｘ，ｙ）から出力された複数回分のダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）を算出する。

　また、温度補償変数として、過去のキャリブレーション時の各回のダーク読取で、当該一の放射線検出素子（ｘ，ｙ）を中心とする前記正方領域内の複数の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値ｄの空間的平均値（以下、上記の放射線画像撮影の直前や直後に行われる１回のダーク読取で得られる温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）と区別するためにキャリブレーション時に得られる温度補償変数としての空間的平均値をｗ（ｘ，ｙ）と表す。）を算出し、さらに各回の温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の複数回分の時間的平均値ω（ｘ，ｙ）を算出する。

　そして、上記の当該一の放射線検出素子（ｘ，ｙ）自身のダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）と、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）とを用いて、当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するようになっている。

　なお、本発明でキャリブレーションという場合、前述したような定期的に（例えば毎月）行われるゲイン補正値及びオフセット補正値の両方を更新するためのキャリブレーションだけでなく、ゲイン補正値に比べ、比較的温度変動が大きいオフセット補正値のみを対象とし、必要に応じて実行されるオフセットキャリブレーションも含まれる。

　すなわち、本発明では、キャリブレーションとは、放射線画像検出器１の工場出荷時、施設への導入時、或いは放射線が曝射されていない待機時に、必要に応じて、ダーク読取値Ｄ（ｘ，ｙ）を含む実写画像データＦ（ｘ，ｙ）を補正するためのデータの全部又は一部を予め取得して作成するための作業をいい、上記のように定期的に行われるキャリブレーションよりも広い概念である。

　従って、施設への導入後のキャリブレーションも、定期的なキャリブレーションよりも短い周期で（例えば毎日）行われる場合もある。また、１回のキャリブレーションやオフセットキャリブレーションでは、通常、複数回のダーク読取が行われる。

　なお、あまりに古いデータを用いて時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗ（ｘ，ｙ）、その時間的平均値ω（ｘ，ｙ）を算出するとそれらの値の信頼性が低下するため、それらの値の算出の元となるダーク読取値ｄ（ｘ，ｙ）は、現在から所定回数前までのダーク読取で得られた値や、最新の（すなわち前回の）キャリブレーションやオフセットキャリブレーションで行われた複数回のダーク読取で得られた値を用いることが好ましい。

　また、ダーク読取値ｄ（ｘ，ｙ）から時間的平均値δ（ｘ，ｙ）等を算出する処理は、放射線画像検出器１で行ってもよく、或いは放射線画像検出器１からコンソール３１（図２９参照）等の外部の装置にダーク読取値ｄ（ｘ，ｙ）を送信して外部装置で行うように構成することも可能である。

　さて、上記の過去の複数回のダーク読取で得られた当該一の放射線検出素子（ｘ，ｙ）自身の各ダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）は、下記（７）式に従って算出される。なお、下記（７）式で、Ｍは時間的平均値δ（ｘ，ｙ）等の算出に用いられるダーク読取値ｄ_ｍ（ｘ，ｙ）の数、すなわち前記所定回数を表す。

　図１１に示すように、ｍ回目（ｍ＝１～Ｍ）のダーク読取で当該一の放射線検出素子（ｘ，ｙ）からダーク読取値ｄ_ｍ（ｘ，ｙ）が出力され、それらの値をヒストグラムにまとめて表すと、図１２に示すように、時間的平均値δ（ｘ，ｙ）を中心とし標準偏差σ_ｄ（ｘ，ｙ）とする正規分布状に分布する。この分布は、ダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的ゆらぎを表している。

　また、過去のキャリブレーション時のｍ回目（ｍ＝１～Ｍ）のダーク読取で、当該一の放射線検出素子（ｘ，ｙ）を中心とする前記正方領域内の複数の放射線検出素子（ｘ，ｙ）_Ｒから出力される各ダーク読取値ｄ_ｍ（ｘ，ｙ）の空間的平均値ｗ_ｍ（ｘ，ｙ）は、上記（３）式と同様の下記（８）式に従って算出される。なお、この場合、（８）式中のＮも２５（＝５×５）である。

　そして、各回の空間的平均値ｗ_ｍ（ｘ，ｙ）のＭ回分の時間的平均値ω（ｘ，ｙ）は、下記（９）式に従って算出される。

　図１１に示すように、ｍ回目（ｍ＝１～Ｍ）のダーク読取でそれぞれ空間的平均値ｗ_ｍ（ｘ，ｙ）が算出され、それらの値をヒストグラムにまとめて表すと、図１２に示すように、時間的平均値ω（ｘ，ｙ）を中心とし、標準偏差σ_ｗ（ｘ，ｙ）の正規分布状に分布する。

　なお、上記（９）式に（８）式と（７）式を代入して変形すると、

となることから分かるように、過去のキャリブレーション時のＭ回のダーク読取での各放射線検出素子（ｘ，ｙ）についての時間的平均値δ（ｘ，ｙ）が算出されれば、各回ごとにそれぞれ空間的平均値ｗ_ｍ（ｘ，ｙ）を算出しなくても、各放射線検出素子（ｘ，ｙ）ごとの時間的平均値δ（ｘ，ｙ）を用いてダーク読取値ｄ_ｍ（ｘ，ｙ）の空間的平均値ｗ_ｍ（ｘ，ｙ）のＭ回分の時間的平均値ω（ｘ，ｙ）を算出することができる。

　ここで、過去のキャリブレーション時の複数回のダーク読取で当該一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）と、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値との関係について考察する。

　キャリブレーションは、一般的に、患者撮影を行わない始業準備中、又は、終業準備中に行われることが多く、キャリブレーション時のダーク読取は、通常、予め設定された同一の温度Ｔ_０の条件下で行われると見做すことができる。従って、過去のキャリブレーション時のダーク読取で得られる当該一の放射線検出素子（ｘ，ｙ）に関する時間的平均値δ（ｘ，ｙ）及び温度補償変数としての空間的平均値ｗ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）は、上記温度Ｔ_０における値である。

　一方、放射線画像撮影が行われる際の放射線検出素子（ｘ，ｙ）やセンサパネル部４上の当該放射線検出素子（ｘ，ｙ）の近傍の温度Ｔは、図１３に示すように、必ずしもキャリブレーション時の温度Ｔ_０ではない。

　特に、放射線画像検出器１では、電源ＯＮの状態でも、放射線画像撮影に使用されない場合には、電力の消費を抑えるために自動的又は手動の操作で電源消費状態がスリープモードとされるようになっているものも多く、図１３に示すように、撮影可能モードでは放射線検出素子（ｘ，ｙ）等の温度Ｔが上昇し、スリープモード（図中のＳの部分参照）ではセンサパネル部４への通電がされなくなるため、放射線検出素子（ｘ，ｙ）等の温度Ｔが低下する。

　また、前述したように、放射線検出素子（ｘ，ｙ）では、通常、放射線検出素子（ｘ，ｙ）等の温度Ｔが低いと出力されるダーク読取値Ｄ（ｘ，ｙ）は小さな値となり、放射線検出素子（ｘ，ｙ）等の温度Ｔが上昇するにつれて出力されるダーク読取値Ｄ（ｘ，ｙ）の値も大きくなる。

　そのため、図１４に示すように、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の中心値（平均値）や、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ，ｙ）_Ｒから出力されるダーク読取値Ｄ_Ｒの平均値である温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布も、温度Ｔによって大きくなったり小さくなったりする。

　すなわち、図１４中のダーク読取値Ｄ（ｘ，ｙ）の分布や、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布が、温度Ｔによって、図１４中の横軸上を左右に変動する。一般には、温度Ｔが上昇すると、分布は横軸上を右側にシフトする傾向にあり、温度Ｔが下降すると、分布は横軸上を左側にシフトする傾向にある。なお、今回の放射線画像撮影の直前や直後に行われた１回のダーク読取で一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）等のゆらぎの分布はあくまで推定されるものであるため、図１４では、破線で表示されている。

　しかし、前述したように、その場合でも、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子（ｘ´，ｙ´）_Ｒは同じように温度変動する。そのため、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）について推定されるゆらぎの分布の図１４のグラフ上の位置と、当該一の放射線検出素子（ｘ，ｙ）に対応付けられた放射線検出素子（ｘ，ｙ）_Ｒから出力されるダーク読取値Ｄ_Ｒの平均値である温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布のグラフ上の位置との相対的な位置関係（すなわち、それぞれの分布の平均値の距離）は温度が変動しても変わらないはずである。

　ところで、もし今回の放射線画像撮影が温度Ｔ_０の環境下で行われた場合、図１４に示したダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μ_Ｄ（ｘ，ｙ）や、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μ_Ｗ（ｘ，ｙ）は、図１２に示した過去のキャリブレーション時の複数回のダーク読取で得られた複数のダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的ゆらぎの分布の平均値δ（ｘ，ｙ）や、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子１４_Ｒから出力される各ダーク読取値ｄ_ｍ（ｘ，ｙ）の空間的平均値ｗ_ｍ（ｘ，ｙ）の時間的ゆらぎの分布の平均値ω（ｘ，ｙ）にほぼ等しくなる。

　また、上記のダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μ_Ｄ（ｘ，ｙ）と温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μ_Ｗ（ｘ，ｙ）との相対的な位置関係は上記のように、温度Ｔが変化しても変わらないため、各分布の平均値の差分、すなわちダーク読取値Ｄ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μ_Ｄ（ｘ，ｙ）と温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の時間的ゆらぎの分布の平均値μ_Ｗ（ｘ，ｙ）との差分も、温度Ｔが変化しても変わらない。

　そのため、今回の放射線画像撮影の直前又は直後に行われた１回のダーク読取で得られる当該一の放射線検出素子（ｘ，ｙ）のダーク読取値Ｄ（ｘ，ｙ）の分布の時間的平均値μ_Ｄ（ｘ，ｙ）と、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の分布の時間的平均値μ_Ｗ（ｘ，ｙ）との差分をε（ｘ，ｙ）と置くと（図１４参照）、
　　ε（ｘ，ｙ）＝μ_Ｄ（ｘ，ｙ）－μ_Ｗ（ｘ，ｙ）　　…（１１）
は、温度Ｔにかかわらず、図１２に示した過去のキャリブレーション時の複数回のダーク読取で得られた複数のダーク読取値ｄ_ｍ（ｘ，ｙ）の分布の時間的平均値δ（ｘ，ｙ）と、当該一の放射線検出素子（ｘ，ｙ）を含む放射線検出素子１４_Ｒから出力される各ダーク読取値ｄ_ｍ（ｘ，ｙ）の空間的平均値ｗ_ｍ（ｘ，ｙ）の分布の時間的平均値ω（ｘ，ｙ）との差分、
　　δ（ｘ，ｙ）－ω（ｘ，ｙ）
に等しくなる。

　すなわち、
　　ε（ｘ，ｙ）＝μ_Ｄ（ｘ，ｙ）－μ_Ｗ（ｘ，ｙ）
　　　　　　　　＝δ（ｘ，ｙ）－ω（ｘ，ｙ）　　…（１２）
が成立する。

　（１２）式は、
　　μ_Ｄ（ｘ，ｙ）＝ε（ｘ，ｙ）＋μ_Ｗ（ｘ，ｙ）　　…（１３）
これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）－ω（ｘ，ｙ）の値が分かっていれば、後は、μ_Ｗ（ｘ，ｙ）さえ分かれば、μ_Ｄ（ｘ，ｙ）を推定することができることを表している。

　ここで、前述したように、（６）式から
　　Ｗ（ｘ，ｙ）≒μ_Ｗ（ｘ，ｙ）　　…（６）
これを（１３）式に代入すると、
　　μ_Ｄ（ｘ，ｙ）≒ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）　　…（１４）
の関係が得られる。

　これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）－ω（ｘ，ｙ）の値が分かっていれば、後は、放射線画像撮影時に取得される１つのダーク画像から算出される温度補償変数Ｗ（ｘ，ｙ）さえ計算できれば、μ_Ｄ（ｘ，ｙ）を推定することができることを表している。

　そして、前述した従来の場合と同様に、ダーク読取値Ｄ（ｘ，ｙ）について推定される時間的ゆらぎの分布の平均値μ_Ｄ（ｘ，ｙ）は、まさに当該放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）の真値とみなすことができるから、結局、上記（１４）式は、
　　Ｏ（ｘ，ｙ）≒ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）　　…（１５）
と書き換えることができる。

　これは、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）－ω（ｘ，ｙ）の値が分かっていれば、後は、放射線画像撮影時に取得される１つのダーク画像から算出される温度補償変数Ｗ（ｘ，ｙ）さえ計算できれば、Ｏ（ｘ，ｙ）、すなわちオフセット補正値の真値を推定することができることを表している。言い換えれば、キャリブレーション時にあらかじめ算出できる変数、すなわち、ε（ｘ，ｙ）＝δ（ｘ，ｙ）－ω（ｘ，ｙ）の値が分かっていれば、放射線画像撮影時の放射線検出素子や電気回路などの各素子がその後どのように温度変化していても、たった１回のダーク画像における各画素位置毎の温度補償変数Ｗ（ｘ，ｙ）を計算すれば、放射線画像撮影時の放射線検出素子や電気回路などの各素子の温度状態における理想的なオフセット補正値Ｏ（ｘ，ｙ）が求められることを表している。

　そこで、本手法では、
　　Ｏ（ｘ，ｙ）＝ε（ｘ，ｙ）＋Ｗ（ｘ，ｙ）　　…（１６）
（ただし、ε（ｘ，ｙ）＝δ（ｘ，ｙ）－ω（ｘ，ｙ））
に従って放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するようになっている。

　これにより、前述した課題を解決するための骨格となる手法が提案できた。

　ところで、これまでの説明では、キャリブレーションを行った時の放射線画像検出器１内の各素子の温度は一定と見なして来たが、キヤリブレーション実施時、ダーク画像の取得数（上述のＭの値）が大きくなると、キャリブレーション中にも放射線画像検出器１内の各素子の温度が変動してしまう。

　上記（７）式～（１０）式からも明白なように、Ｍの値が大きいほど、ｄ_ｍ（ｘ，ｙ）、ｗ_ｍ（ｘ，ｙ）の時間的なゆらぎが相殺され、（１６）式の精度を左右する、ω（ｘ，ｙ）やδ（ｘ，ｙ）の精度が向上するが、Ｍの値が大きいほど、ω（ｘ，ｙ）やδ（ｘ，ｙ）の値は温度変動の影響を受けやすくなることを意味している。

　そこで、この問題に対しても、これまで説明してきた温度補正の考え方を適応する。

　（１２）式に（７）式、（９）式を代入すると、

　ここで、
　　λ_ｍ（ｘ，ｙ）＝ｄ_ｍ（ｘ，ｙ）－ｗ_ｍ（ｘ，ｙ）　　…（１８）
とおくと、

　ここで、λ_ｍ（ｘ，ｙ）は、キャリブレーション時のダーク読取値ｄ_ｍ（ｘ，ｙ）の温度変化を温度補償変数ｗ_ｍ（ｘ，ｙ）によって温度補正された、温度補正済みダーク読取値λ_ｍ（ｘ，ｙ）と考えることができる。

　キャリブレーション時に放射線画像検出器１内の各素子の温度が変動してしまうと図１５に示すように、ｄ_ｍ（ｘ，ｙ）の分布は非常にブロードな分布となり、従って、上記（７）式で求められるｄ_ｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）は、正しい値を示せなくなる。

　そこで、一般的にキャリブレーション時にダーク読取を複数回実施し、その平均値を求める場合（従来は、この方法でオフセット補正値を直接求めることは前述した通りである）、直接ｄ_ｍ（ｘ，ｙ）の平均値を計算するのではなく、上述のように、温度補正済みダーク読取値λ_ｍ（ｘ，ｙ）の平均値（時間的平均値）であるε（ｘ，ｙ）を求めるようにすれば、温度変動が補正された正確なダーク読取の時間的平均値を求めることができるはずである。すなわち、ε（ｘ，ｙ）は温度補正済みダーク読取の時間的平均値と呼ぶことが可能である。

　温度補正済みダーク読取値λ_ｍ（ｘ，ｙ）や、温度補正済みダーク読取の時間的平均値ε（ｘ，ｙ）は、キャリブレーション時、さらには、放射線画像撮影時においても、ダーク読取を複数回実施し、その平均（時間的平均）を取る作業全てに適応することができる。また、時間的平均回数であるＭの値が大きくとも、温度変動の影響を受けにくくすることができる。

　また、キャリブレーション時にゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、同様な温度補償を適応することが可能である。ゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、あらかじめ定めておいた所定の放射線量で、被写体を置かない状態での放射線画像データを取得し、これにオフセット補正を施してゲイン補正値Ｇ（ｘ，ｙ）を計算するため、この際に取得する放射線画像データを、本手法の実写画像データＦ（ｘ，ｙ）と見なせば、本手法と同様の温度補正のかかったゲイン補正データＧ（ｘ，ｙ）を得ることができる。また、ゲイン補正値Ｇ（ｘ，ｙ）を求める際に使用されるオフセット補正値に対しても、同様に前述の温度補正を行うことが好ましい。これにより、精度の高いゲイン補正値Ｇ（ｘ，ｙ）を求めることができ、結果として、ＳＮ比の高い良好な画像データを得ることが可能となる。

　本手法において、キャリブレーション時のＭ回のダーク読取で放射線検出素子（ｘ，ｙ）から出力された各ダーク読取値ｄ_ｍ（ｘ，ｙ）の分布（図１２参照）においては、時間的平均値δ（ｘ，ｙ）が算出されるが、さらに、分布の標準偏差σｄ（ｘ，ｙ）を算出し、
　　ｄ_ｍ（ｘ，ｙ）－δ（ｘ，ｙ）　　…（２０）
　　ｄ_ｍ（ｘ，ｙ）／σｄ（ｘ，ｙ）　　…（２１）
或いは、
　　（ｄ_ｍ（ｘ，ｙ）－δ（ｘ，ｙ））／σｄ（ｘ，ｙ）　　…（２２）
等の演算を行って各ダーク読取値ｄ_ｍ（ｘ，ｙ）を正規化し、正規化された各ダーク読取値に基づいて上記の処理を構成することも可能である。

［放射線検出素子に予め対応付ける他の放射線検出素子の選択方法の例］
　前述したように、一の放射線検出素子（ｘ，ｙ）について温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）を算出するために当該一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子（ｘ´，ｙ´）については、当該一の放射線検出素子（ｘ，ｙ）と同じように温度変動するものであればよく、それらの選択の仕方としては、上記のように種々の手法があり得る。ここでは、放射線画像検出器１の実際的な構成に即した選択の仕方の一例について説明する。

　例えば、図６では、放射線画像検出器１のセンサパネル部４上で一の放射線検出素子（ｘ，ｙ）の周囲の行方向や列方向に十分な数の他の放射線検出素子（ｘ´，ｙ´）がある場合に、当該一の放射線検出素子（ｘ，ｙ）を中心として例えば５×５の正方領域内に存在する他の放射線検出素子（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）を選択する場合を示した。しかし、当該一の放射線検出素子（ｘ，ｙ）がセンサパネル部４の周縁部分にある場合には、このように例えば５×５の正方領域の中心に当該一の放射線検出素子（ｘ，ｙ）を位置させることができない場合もある。

　従って、そのような場合には、例えば図１６Ａ～図１６Ｃに示すように、図中で濃く着色されて表される当該一の放射線検出素子（ｘ，ｙ）の位置などを選択すれば良い。

　なお、この場合、キャリブレーション時のダーク読取における空間的平均値ｗ_ｍ（ｘ，ｙ）の算出やそれらの時間的平均値ω（ｘ，ｙ）の算出においても、当該放射線検出素子（ｘ，ｙ）についての領域Ｒは同じ要領で設定される。また、図１６Ａ～図１６Ｃにおいて、格子状に表される各放射線検出素子の上方及び左方に記載されている数字は、それぞれ各放射線検出素子（ｘ，ｙ）の行方向及び列方向の位置を表す数字であり、左方の座標（ｘ，ｙ）は当該放射線検出素子（ｘ，ｙ）の座標を表すものである。

　一方、各放射線検出素子（ｘ，ｙ）から電気信号を読み出し増幅するための読み出し回路９（図２参照）は、図１７に示すように、通常、行方向に１２８画素や２５６画素ずつの放射線検出素子（ｘ，ｙ）ごとに１個ずつの読み出しＩＣ９１、９２、…が接続されており、それが必要な数だけ並設されて構成されている。そして、読み出し回路９側の読み出しＩＣ９１、９２、…ごとの温度特性やノイズ特性等が必ずしも同一ではないため、各放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）が読み出しＩＣ９１、９２、…ごとに異なる場合がある。

　そのため、図１８に示すように、放射線検出素子（ｘ，ｙ）のオフセット補正値Ｏ（ｘ，ｙ）を算出するために設定される領域Ｒに属する各放射線検出素子から電気信号を読み出す読み出しＩＣが、隣接する読み出しＩＣ９１、９２で別々に読み出される場合、その領域Ｒに属する各放射線検出素子（ｘ，ｙ）については温度補償変数Ｗ（ｘ，ｙ）やω（ｘ，ｙ）の値が誤差を持ってしまい、オフセット補正値Ｏ（ｘ，ｙ）を良好に算出できない場合がある。

　そこで、このような場合には、例えば図１９Ｃ、図１９Ｄに示すように、読み出しＩＣ９１に接続されている放射線検出素子（ｘ，ｙ）と読み出しＩＣ９２に接続されている放射線検出素子（ｘ，ｙ）とを分けて考える。そして、読み出しＩＣ９１、９２の境界Ｌの近傍に存在する放射線検出素子（ｘ，ｙ）を、上記のセンサパネル部４の周縁部分にある場合と同様に扱うことで、一の放射線検出素子（ｘ，ｙ）に予め対応付けられる複数の放射線検出素子を同一の読み出しＩＣ９１に接続される放射線検出素子の中から選択するように構成することが可能となる。

　なお、図１９Ａ、図１９Ｂには、図６等に示した通常の仕方で一の放射線検出素子（ｘ，ｙ）に予め対応付けられる複数の放射線検出素子が選択される場合が示されている。

［欠陥画素判定方法や異常画素判定方法等の特徴］
　以下、本発明における欠陥画素判定方法や異常画素判定方法等の特徴について説明する。

　まず、欠陥画素は、製造後の出荷検査や、定期的なキャリブレーション時に放射線照射を伴う検査を実施することで検出され、検出された欠陥画素の画素位置が欠陥画素マップに登録し、個々の放射線画像検出器と対応付けて運用、管理されることは前述した通りである。

　欠陥画素の例として、定常的に大きな画素値を出力する画素（例えば、信号値が常に飽和レベルまで達している画素）や、もしくは定常的に小さな画素値しか出力しない画素（例えば、常に信号値がゼロレベルの値しか出力しない画素）、また、常に一定の信号値しか出力しない画素（放射線を照射しても出力する画素値が変化しない画素）については、容易に判別することが可能であるため、このような欠陥画素については、既に判別されて欠陥画素マップに登録されているものとする。

　ここで解決しようとしているのは、ある確率で異常な信号値を出力する異常な画素についてであり、これらの異常な画素を放射線を照射せずとも精度良く検出し、欠陥画素であるか否かを判定する方法、もしくは、異常な画素を欠陥画素マップに登録する方法、もしくは異常な画素を欠陥画素マップに登録せずに異常画素として適切に処理する方法、もしくは、放射線を照射した画像データ中の異常画素について対処する方法についてである。

　そのため、本発明における主たる方法では、キャリブレーション時や放射線画像撮影の前後に取得するダーク読取値（放射線照射を必要としない読取信号）を使用する。また、欠陥画素や異常画素を判別する際に、放射線画像検出器内の各素子の温度変動を考慮した欠陥画素を判別方法についても説明を行う。また、放射線を照射して得られた画像データの中から異常画素を検出し適切に対処する方法についても説明を行う。

　まず、本発明で解決したい欠陥画素と異常画素についての説明を行う。（１）式に示したように、最終的な画像データＦ_Ｏ（ｘ，ｙ）を得るには、実写画像データＦ（ｘ，ｙ）からオフセット補正値Ｏ（ｘ，ｙ）を差し引き、その差分にゲイン補正値Ｇ（ｘ，ｙ）を乗算する。
　　Ｆ_Ｏ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）－Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）　　…（１）

　従って、実写画像データＦ（ｘ，ｙ）が正常な画素値を有していても、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）に異常があれば、その演算結果である最終的な画像データＦ_Ｏ（ｘ，ｙ）は異常値を示すため、これは異常な画素（本来補正すべき画素）ということになる。

　また、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）が正常であっても、実写画像データＦ（ｘ，ｙ）の値が異常であれば、その演算結果である最終的な画像データＦ_Ｏ（ｘ，ｙ）は異常値を示すため、これも異常な画素（本来補正すべき画素）ということになる。

　ある確率でオフセット補正値Ｏ（ｘ，ｙ）、ゲイン補正値Ｇ（ｘ，ｙ）、実写画像データＦ（ｘ，ｙ）で異常画素値を出力する画素位置（ｘ，ｙ）は、その異常値の強度（どれくらい異常かを示す指針）や異常値を出力する頻度に応じて分類され、許容できない、もしくは許容すべきでないと判断された異常な画素は、欠陥画素として欠陥画素マップに登録を行う。

　また、突発的に異常値を出力する画素であっても、その異常値の強度（どれくらい異常かを示す指針）や異常値を出力する頻度が許容されるレベルと判断される場合は、異常画素として欠陥画素と同様の処理（補正処理など）を実施するが、欠陥画素として欠陥画素マップに登録しないことも許容する。もちろん、異常画素を全て欠陥画素として欠陥画素マップに登録しても良い。

　ここで、ゲイン補正値Ｇ（ｘ，ｙ）については、キャリブレーション時に放射線を照射して算出するものであるが、オフセット補正値Ｏ（ｘ，ｙ）については、放射線を照射せずに算出できるパラメータである。

　本発明者らは、このオフセット補正値の算出に使用するダーク読取値Ｄ（ｘ，ｙ）もしくはダーク読取値ｄ_ｍ（ｘ，ｙ）に着目し、ある確率で発生する異常画素をダーク読取値Ｄ（ｘ，ｙ）もしくはダーク読取値ｄ_ｍ（ｘ，ｙ）から見出し、判別し、欠陥画素や異常画素として処理する（場合によっては異常画素を欠陥画素として欠陥画素マップに登録する）方法を見出した。また、放射線画像データ中の異常画素についても検出し、適切に処理する方法についても見出した。

［欠陥画素判定方法］
　以下、まず、本発明における欠陥画素判定方法の実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。

　前述した図５等に示したように、同じ温度Ｔの環境下にあったとしても、放射線画像検出器１の一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）にはゆらぎ（ばらつき）が生じ、ゆらぎの分布は平均値μ（ｘ，ｙ）を中心に標準偏差σ（ｘ，ｙ）を有する正規分布状の分布となる。

　しかし、放射線検出素子（ｘ，ｙ）によっては、分布の平均値μ（ｘ，ｙ）が他の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の分布の平均値μ（ｘ´，ｙ´）から大きくかけ離れていたり、或いはゆらぎが非常に大きい場合、すなわち分布の標準偏差σ（ｘ，ｙ）が他の放射線検出素子（ｘ´，ｙ´）から出力されるダーク読取値Ｄ（ｘ´，ｙ´）の分布の標準偏差σ（ｘ´，ｙ´）に比べて非常に大きいような場合がある。

　そのため、そのような放射線検出素子（ｘ，ｙ）では、最終的に得られる画像データＦ_Ｏ（ｘ，ｙ）の安定性も悪くなり、欠陥画素と見なすことができる程の異常な値となってしまう場合がある。以下、このように、欠陥画素と見なすことができる放射線検出素子（ｘ，ｙ）を欠陥画素（ｘｓ，ｙｓ）と表す。

　また、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）は正常な範囲に入るものの、時折異常な値を出力する画素が存在する。このような画素については、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）のみに注目していても検出することはできない。このように、ダーク読取値Ｄ（ｘ，ｙ）の標準偏差σ（ｘ，ｙ）は正常な範囲に入るものの、時折異常な値を出力する画素も、許容範囲を越えるものは異常画素と見なす。

　そこで、本発明における欠陥画素判定方法では、放射線検出素子（ｘ，ｙ）がある確率で異常画素となる場合、その発生確率や異常度合いから、欠陥画素（ｘｓ，ｙｓ）に登録するか否かを判定するようになっている。

　前述したように、過去のキャリブレーション時に複数回（Ｍ回）行われたダーク読取において放射線検出素子（ｘ，ｙ）自身から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）の複数回分（Ｍ回分）のデータが、各放射線検出素子（ｘ，ｙ）ごとに蓄積される。そして、それらのダーク読取値ｄ_ｍ（ｘ，ｙ）をヒストグラムにまとめて表すと、図１２に示したように、各放射線検出素子（ｘ，ｙ）ごとにダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的ゆらぎの分布の時間的平均値δ（ｘ，ｙ）と標準偏差σｄ（ｘ，ｙ）等の時間的統計値を算出することができる。

　そこで、このように、過去のキャリブレーション時に複数回行われたダーク読取において得られた、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）の値そのもの、もしくはダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的ゆらぎの分布における時間的平均値δ（ｘ，ｙ）や標準偏差σｄ（ｘ，ｙ）等の時間的統計値に基づいて、下記の各判定手法に従って、当該一の放射線検出素子（ｘ，ｙ）が欠陥画素か否かを判定するようになっている。

［判定手法１］
　まず、時間的統計値として、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的ゆらぎの分布における標準偏差σｄ（ｘ，ｙ）に着目し、標準偏差σｄ（ｘ，ｙ）が予め設定された閾値σｄ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　　…（４３）
が成り立つ場合には、当該一の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）の大きさは、他の放射線検出素子（ｘ^＊，ｙ^＊）から出力されるダーク読取値Ｄ（ｘ^＊，ｙ^＊）の大きさとの違いが大きく許容範囲を越えるものであり、異常画素であると判定される。

　なお、差分Λ（ｘ，ｙ）は、（４１）式に示したように、共に温度補償を行った変数である差分ε（ｘ，ｙ）と差分Ｅ（ｘ，ｙ）どうしの演算であるため、ダーク読取値Ｄ（ｘ，ｙ）の温度変化による誤差が相殺されている。

　なお、［判定手法７］、［判定手法８］で示した異常画素の取扱については、以下の選択肢１）～５）の方法のいずれかを採用することが可能である。

選択肢１）偶発的に発生した異常画素と見なし、欠陥画素としては登録しない。
選択肢２）欠陥画素（ｘｓ，ｙｓ）として、直ちに欠陥画素マップに登録する。
選択肢３）異常画素の値の大きさに応じて、欠陥画素として登録するか否かを決める。すなわち、異常画素の値が非常に大きい場合（別途定められた欠陥画素への登録を判定するための閾値を越えた場合）は、欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録する。
選択肢４）異常画素として、欠陥画素とは別に画素位置を記憶しておき、その後、同じ画素位置の信号値がある頻度で異常画素となれば、その時点で発生確率の判定を行い、予め定められた発生確率の閾値を越えた時点で欠陥画素（ｘｓ，ｙｓ）として欠陥画素マップに登録する。
選択肢５）上記、選択肢３）と選択肢４）を組み合わせて実施する。

　上記［判定手法４］～［判定手法８］については、温度補償変数である空間的平均値ｗ_ｍ（ｘ，ｙ）や空間的平均値Ｗ（ｘ，ｙ）、及びこれらの温度補償変数を使用して求められた新たな変数である差分ε（ｘ，ｙ）、差分ｅ_ｍ（ｘ，ｙ）、差分Ｅ（ｘ，ｙ）、絶対値｜Λ（ｘ，ｙ）｜などが使われている。ここで、空間的平均値ｗ_ｍ（ｘ，ｙ）や空間的平均値Ｗ（ｘ，ｙ）の中に、欠陥画素として登録されていない突発的な異常画素が存在したならば、温度補償係数の値の精度が低下してしまう懸念がある。

　しかしながら、温度補償係数の場合は、ダーク読取値（Ｄ（ｘ，ｙ）やＤ_ｍ（ｘ，ｙ）やｄ_ｍ（ｘ，ｙ））の注目画素位置（ｘ，ｙ）に対してその周辺のＮ個の画素値の空間的平均値であるため、仮にＮ個の周辺画素の中に欠陥画素として登録されていない突発的な異常画素が存在したとしても、その影響は１／Ｎに緩和されているので大きな問題とはならないが、別途手段を設けて突発的な異常画素を検出し、温度補償変数の計算に使用する画素から除外するようにしても良い。

　さて、ここまでは、欠陥画素や異常画素をダーク読取値を用いて判別する方法について述べてきた。しかしながら、前述したように、良好な最終的な画像データＦ_Ｏ（ｘ，ｙ）を得るには、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）のみならず、実写画像データＦ（ｘ，ｙ）についても良好であること求められる。

　すなわち、オフセット補正値Ｏ（ｘ，ｙ）やゲイン補正値Ｇ（ｘ，ｙ）に対して欠陥画素や異常画素の対処が適切になされていても、実写画像データＦ（ｘ，ｙ）の画素値の中に突発的な異常画素が含まれていれば、最終的な画像データＦ_Ｏ（ｘ，ｙ）にも異常画素が認識されることになるため好ましくない。従って、実写画像データＦ（ｘ，ｙ）についても、欠陥画素補正だけでなく、突発的に発生する異常画素に対して適切な処置を施すことが好ましい。

　しかしながら、実写画像データＦ（ｘ，ｙ）については、下記に示すような特徴があるため、これまで説明して来たダーク読取値ほど扱いが簡単ではない。

１）実写画像データＦ（ｘ，ｙ）は殆どの場合、１回しか放射線撮影されない実写画像データＦ（ｘ，ｙ）として扱う必要がある。同じ被写体に対して複数回撮影される場合もあるが、その場合、複数回の放射線画像撮影の間に被写体を動かしたり被写体の向きを変えたりするため、おのおのが独立した実写画像データＦ（ｘ，ｙ）であると考えなければならない。

２）実写画像データＦ（ｘ，ｙ）を取得するための放射線画像撮影とは、ゲイン補正データを取得するためのベタ画像撮影（ベタ画像撮影とは被写体を置かずに、一様な放射線を照射する撮影のこと）をのぞき、必ず被写体を置いての撮影となるため、実写画像データＦ（ｘ，ｙ）には被写体の画像が写り込んでいる。

　しかも、前述したように、実写画像データＦ（ｘ，ｙ）は１回しか放射線撮影されない独立したデータであるため、実写画像データＦ（ｘ，ｙ）の値が被写体情報を正確に反映しているのか否かを客観的に判断することが難しい。例えば、同じ被写体に対して同条件での放射線画像撮影が複数回行われれば（複数回の撮影の間、被写体は動いていないと仮定する）、所定の画素位置（ｘ，ｙ）について複数回の画像データのばらつきを比較することにより、統計的な見地から被写体情報をほぼ正確に反映しているか否かを判断できる。

　例えば、複数回の放射線画像撮影における画素位置（ｘ，ｙ）の信号値の平均値を計算し、その平均値からの各信号値のばらつきを測定した時、ばらつきが少なければ、それらの信号値は被写体情報をほぼ正確に反映していると判断できる。一方、複数回の放射線画像撮影における画素位置（ｘ，ｙ）の信号値の平均値に対して、ある信号値の値が平均値から大きく逸脱した値を有していれば、その信号値は被写体情報を正確に反映しておらず、ノイズ等の影響を受けた異常画素であると判断できる。

　このように、同じ被写体に対して同条件での放射線画像撮影が複数回行われれば、統計的な見地から被写体情報をほぼ正確に反映しているか否かを判断できるが、実写画像データＦ（ｘ，ｙ）は１回しか放射線撮影されない独立したデータであるため、このような統計的判断をすることができない。

　従って、各画素値に対して、本来どのような値であるべきか（各画素値のあるべき値、すなわち各画素値の真値）が分からない。例えば、ダーク読取値であれば、キャリブレーション時に取得した複数回のデータに対して今回のデータが正常であるか異常であるかは、上記のような統計的処理によって判断することができるが、１回しか放射線撮影されない実写画像データＦ（ｘ，ｙ）に対しては同様の判断を行う術がない。

　このような課題の中で、実写画像データＦ（ｘ，ｙ）中に発生した突発的な異常画素を判別し、適切な処置を施す手法について以下に説明を行う。

［判定手法９］
　まず、実写画像データＦ（ｘ，ｙ）の各画素値は、被写体の放射線吸収率によって様々に変化するが、それは画素サイズレベルのミクロの視点で見れば、非常になだらかに変化する信号値である（通常の画素サイズは１００～２００ミクロンの範疇で選択されることが多く、１００ミクロン以下の選択はあっても２００ミクロンを越える画素サイズの選択は行われない場合が多い）。特に特定画素だけが大きく変化するような画素は存在しないと言って良い。また、最も画素レベルに近いサイズを有する病変としては、乳房撮影における微小石灰化があるが、これは放射線を吸収する方向に作用するため、信号値が減衰する方向で信号値が変化する。

　ここで問題にしたい異常画素は、殆どが信号値が大きくなるように作用するため（信号値が小さくなるように作用する異常は既に欠陥画素として登録されている）、放射線を吸収する方向に作用する病変が異常画素と混同されることはない。

　一方、信号値が大きくなるように作用する病変で、かつ画素レベルのサイズを持つものは存在しないため、独立した画素で周囲の画素に対して大きな値を持つものは突発的な異常画素と見なして問題ない。本発明ではこの特性を利用して、実写画像データＦ（ｘ，ｙ）の突発的な異常画素の判別を行う。

　さらに突発的に発生する画素は常に孤立した１画素である。そもそも突発的に発生する異常画素は、発生頻度が低く（発生頻度が高いものは、既に欠陥画素として登録されている）、かつ、１枚のフラットパネルディテクタは通常５メガ画素前後の画素数を有するため、今回問題にしている突発的な異常画素が並んで、もしくはクラスター状の塊として発生する確率は天文学的に小さい。

　従って、上述したように、今回課題として取り上げている突発的な異常画素は常に孤立した１画素として発生すると見なして良い。もちろん、隣り合って発生する場合や、クラスター状の塊として発生する場合を考慮しても良いことは言うまでもないが、以下の説明では、孤立した１画素として発生すると見なして説明を行う。

　上記とまとめると、今回課題としている実写画像データＦ（ｘ，ｙ）の突発的な異常画素は以下の性質をもつものとして考えられる。

１）異常画素は必ず信号値が大きくなる（放射線が多く照射された画素のような挙動を示す）
２）異常画素は孤立した１点である。異常画素の周囲の画素は、欠陥画素として登録されている画素以外は必ず正常な画素である。
３）正常画素は、信号値が大きくなる方向にはなめらかに変化し、１画素単位で大きな画素変化を示さない。

　上記突発的異常画素を判別する方法の１つの実施例を、前述した図６を使って以下に示す。

　例えば複数の放射線検出素子（ｘ，ｙ）として当該一の放射線検出素子（ｘ，ｙ）すなわち、画素位置（ｘ，ｙ）を、実写画像データの注目画素位置とする。そして、画素位置（ｘ，ｙ）を中心とする５×５の正方領域内に存在する２５の画素位置（ｘ－２，ｙ－２）～（ｘ＋２，ｙ＋２）を比較領域として定義し、この比較領域の信号値から、注目画素である画素位置（ｘ，ｙ）の信号値を除く２４個の信号値の平均値をＡ（ｘ，ｙ）と表すと、

　ここでＬの値は、この場合２４（＝５×５?１）である。本例では、Ｌ＝２４とし、画素位置（ｘ，ｙ）中心とする５×５の正方領域を比較領域として説明を行うが、比較領域の選定方法はこれに限定するものでない。Ｌ＝８（＝３×３?１）として、画素位置（ｘ，ｙ）中心とする３×３の正方領域を比較領域としても良いし、例えば７×７もしくはそれ以上のサイズの正方領域を比較領域として採用しても本発明の効果を損なうものではない。また、比較領域として画素位置（ｘ，ｙ）囲む領域が選定されていれば、比較領域は必ずしも正方領域でなくとも良い。また、演算回数を減らすために、比較領域を画素位置（ｘ，ｙ）の上下、左右の４画素としたり、比較領域を画素位置（ｘ，ｙ）の上下もしくは左右の２画素としても良い。

　次に、（４５）式に従って、差分Ｖ（ｘ，ｙ）を計算する。
　　Ｖ（ｘ，ｙ）＝Ｆ（ｘ，ｙ）－Ａ（ｘ，ｙ）　　…（４５）

　この差分Ｖ（ｘ，ｙ）が、予め設定された閾値Ｖth以下であれば、すなわち、
　　Ｖ（ｘ，ｙ）≦Ｖth　　…（４６）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。

　一方、差分Ｖ（ｘ，ｙ）が、予め設定された閾値Ｖthよりも大きい場合、すなわち、
　　Ｖ（ｘ，ｙ）＞Ｖth　　…（４７）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）を異常画素と判定して、（４９）式に定義される平均値Ｂ（ｘ，ｙ）で実写画像データＦ（ｘ，ｙ）の値を置き換える（（４８）式参照）。
　　Ｆ（ｘ，ｙ）←Ｂ（ｘ，ｙ）　　…（４８）

　ここで、画素位置（ｘ，ｙ）を中心とする３×３の正方領域内に存在する９つの画素位置（ｘ－１，ｙ－１）～（ｘ＋１，ｙ＋１）を計算領域として定義し、この計算領域の信号値から、注目画素である画素位置（ｘ，ｙ）の信号値を除く８個の信号値の平均値を平均値Ｂ（ｘ，ｙ）とする。

　（４９）式において、Ｐの値は、この場合８（＝３×３?１）である。本例では、Ｐ＝８とし、画素位置（ｘ，ｙ）中心とする３×３の正方領域を計算領域として説明を行うが、計算較領域の選定方法はこれに限定するものでない。Ｐ＝２４（＝５×５?１）として、画素位置（ｘ，ｙ）中心とする５×５の正方領域を比較領域としても良いし、例えば７×７もしくはそれ以上のサイズの正方領域を比較領域として採用しても良い。また、計算領域として画素位置（ｘ，ｙ）囲む領域が選定されていれば、計算領域は必ずしも正方領域でなくとも良い。また、演算回数を減らすために、計算領域を画素位置（ｘ，ｙ）の上下、左右の４画素としたり、計算領域を画素位置（ｘ，ｙ）の上下もしくは左右の２画素としても良い。

　また、ここでは、比較領域Ａ（ｘ，ｙ）と計算領域Ｂ（ｘ，ｙ）を個別に定義したが（比較領域と計算領域は必ずしも一致している必要はない）、比較領域Ａ（ｘ，ｙ）と計算領域Ｂ（ｘ，ｙ）を同一のものとして扱っても良い。なお、比較領域や計算領域の中に欠陥画素が存在した場合は、欠陥画素に隣接する正常画素を使用して平均値Ａ（ｘ，ｙ）や平均値Ｂ（ｘ，ｙ）を計算すれば良い。

　このように、実写画像データＦ（ｘ，ｙ）の中に突発的に異常画素に対する対応を行えば、最終的な画像データＦ_Ｏ（ｘ，ｙ）の中に異常画素は発生することがなくなり、良好な診断画像を安定的に供給することが可能となる。なお、ここで検出された異常画素の取扱（欠陥画素として登録するか否か）については、上述の［判定手法８］の後に記述した選択肢１）～５）の方法のいずれかを採用することが可能である。

［判定手法１０］
　なお、上記では、実写画像データＦ（ｘ，ｙ）中の異常画素は必ず信号値が大きくなる（放射線が多く照射された画素のような挙動を示す）として説明を行ったが、異常画素が、信号値が小さくなる方向に変化する場合も考慮しても良い。すなわち、

１）異常画素は必ず周囲の画素値に対して大きくなる方向、もしくは小さくなる方向に異常値を持つ。
２）異常画素は孤立した１点である。異常画素の周囲の画素は、欠陥画素として登録されている画素以外は必ず正常な画素である。
３）正常画素は、はなめらかに変化し、１画素単位で大きな画素変化を示さない。
として判断を行っても良い。

　すなわち、（４５）式のＶ（ｘ、ｙ）に対する絶対値｜Ｖ（ｘ、ｙ）｜を定義し、この絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´th以下であれば、すなわち、
　　｜Ｖ（ｘ、ｙ）｜≦Ｖ´th　　…（５０）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。

　一方、差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´thよりも大きい場合、すなわち、
　　｜Ｖ（ｘ、ｙ）｜＞Ｖ´th　　…（５１）
であれば、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）を異常画素と判定して、（４９）式に定義される平均値Ｂ（ｘ，ｙ）で実写画像データＦ（ｘ，ｙ）の値を置き換えるようにしても良い（（４８）式参照）。以下、［判定手法９］での説明と同様である。

［判定手法１１］
　また、前述したように、乳房撮影における微小石灰化の病変については病変サイズが画素サイズレベルに近いため、乳房撮影を行って周囲の画素から出力された信号値と大きく異なる実写画像データＦ（ｘ，ｙ）を出力した放射線検出素子があったとしても、それが病変部に存在する現実の微小石灰化を撮影したものなのか、単に異常値が出力されたものなのかを判別することは難しく、乳房撮影を行って得られた実写画像データＦ（ｘ，ｙ）に基づいて放射線検出素子を異常画素と見なしてよいか否かの判定を行うことは困難である。

　しかし、その一方で、特に乳房撮影に用いられる放射線画像検出器１については、放射線が照射された際にどの画素（放射線検出素子）がどの程度異常な信号値を出力するかを予め把握しておかないと、微小石灰化が存在するのに異常値と誤判定したり、異常値を微小石灰化と誤診してしまうという問題が生じる。

　そこで、放射線画像検出器１に対して、被写体が存在しない状態で、放射線画像検出器１に一様に放射線を照射して、各放射線検出素子（ｘ，ｙ）から実写画像データＦ（ｘ，ｙ）を出力させ、その被写体が存在しない状態で取得された各実写画像データＦ（ｘ，ｙ）に基づいて、上記［判定手法１０］と同様にして異常画素の判定を行うように構成することができる。

　すなわち、被写体が存在しない状態で放射線画像検出器１の放射線入射面Ｘ（図１参照）に対してその全面にわたってほぼ一定の線量となるように放射線を照射して、各放射線検出素子（ｘ，ｙ）からほぼ同じ信号値の実写画像データＦ（ｘ，ｙ）が出力されるべき状況とする。

　そして、その状態で実際に各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）について、上記（４４）式及び（４５）式に従ってＶ（ｘ，ｙ）を算出してＶ（ｘ、ｙ）に対する絶対値｜Ｖ（ｘ、ｙ）｜を算出し、この絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´th以下であれば、すなわち、
　　｜Ｖ（ｘ、ｙ）｜≦Ｖ´th　　…（５０）
であれば、当該一の放射線検出素子は異常画素ではないと判定する。そして、実際の乳房撮影において、当該一の放射線検出素子から出力された実写画像データＦ（ｘ，ｙ）の値は変更せずにそのまま使用する。また、更に、ダーク読取値に基づく異常画素か否かの判定も行い、両方で異常画素では無いと判定された場合にのみ実写画像データをそのまま使用するとより好ましい。

　しかし、各放射線検出素子（ｘ，ｙ）からほぼ同じ信号値の実写画像データＦ（ｘ，ｙ）が出力されるべき状況であるにもかかわらず、差分Ｖ（ｘ，ｙ）の絶対値｜Ｖ（ｘ、ｙ）｜が、予め設定された閾値Ｖ´thよりも大きい場合、すなわち、
　　｜Ｖ（ｘ、ｙ）｜＞Ｖ´th　　…（５１）
であれば、当該一の放射線検出素子を異常画素と判定する。このような場合には、周辺画素は正常組織を検出した出力値となっているので、当該異常画素に対し、周辺画素の値を用いた補間処理した値を使用することはできず、（１）異常画素が存在することを告知し他の放射線画像検出器を使用するよう促す、（２）当該異常画素のアドレスを表示し、異常画素の存在する領域を照射野領域外（乳房領域外）として撮影を行うよう被写体に対する放射線画像検出器の位置を変更するよう促す、等を技師等の操作者に対してコンソールを介して告知する。

　さて、これまでは、異常画素の判定方法、補正方法をオフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に対して個別に実施する説明を行って来たが、異常画素の判定方法、補正方法をオフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に対して個別に実施せずに、それらの演算結果である最終的な画像データＦ_Ｏ（ｘ，ｙ）に対して実施する方法もある。

　最終的な画像データＦ_Ｏ（ｘ，ｙ）は（１）式
　　Ｆ_Ｏ（ｘ，ｙ）＝（Ｆ（ｘ，ｙ）－Ｏ（ｘ，ｙ））×Ｇ（ｘ，ｙ）　　…（１）
によって計算されるため、オフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）に異常画素が存在すれば、最終的な画像データＦ_Ｏ（ｘ，ｙ）も異常画素となる。

　異常画素は殆どの場合、オフセット補正値Ｏ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）の画素値が大きくなる方向に変化するため、上記（１）式の演算結果では、オフセット補正値Ｏ（ｘ，ｙ）に異常画素があれば最終的な画像データＦ_Ｏ（ｘ，ｙ）の信号値が小さくなる方向に変化し、実写画像データＦ（ｘ，ｙ）に異常画素があれば最終的な画像データＦ_Ｏ（ｘ，ｙ）の信号値が大きくなる方向に信号値が変化することになる。

　このため、判定方法としては、上記［判定手法１０］の実写画像データＦ（ｘ，ｙ）を最終的な画像データＦ_Ｏ（ｘ，ｙ）と置き換えて考えれば、［判定手法１０］で説明した方法と同様の方法で異常画素を判定することが可能である。

［欠陥画素がある場合の空間的平均値の算出手法］
　次に、一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子の中に、欠陥を有する放射線検出素子（以下、欠陥画素という。）が含まれる場合の当該一の放射線検出素子（ｘ，ｙ）の温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）の算出手法について説明する。

　以下の説明を分かり易くするために、ここでは、図２４に示すように、温度補償変数としての空間的平均値Ｗ（ｘ，ｙ）を算出する当該一の放射線検出素子（注目画素）が放射線検出素子（４，４）であり、放射線検出素子（４，４）にそれを中心とする７×７個の放射線検出素子（１，１）～（７，７）が対応付けられており、放射線検出素子（６，６）が欠陥画素である場合について説明する。

［算出手法１］
　最も単純な算出手法としては、欠陥画素（６，６）を含めて、当該一の放射線検出素子（４，４）に予め対応付けられた複数の放射線検出素子（１，１）～（７，７）から出力されるダーク読取値Ｄ（１，１）～Ｄ（７，７）を用いて、

により放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出することが考えられる。しかし、この算出手法では、Ｗ（４，４）が欠陥画素（６，６）のダーク読取値Ｄ（６，６）の影響を受けてしまう可能性が残る。

［算出手法２］
　当該一の放射線検出素子（４，４）に予め対応付けられた複数の放射線検出素子（１，１）～（７，７）のうち、欠陥画素（６，６）を除く放射線検出素子から出力されたダーク読取値Ｄ（１，１）～Ｄ（７，７）（ただしＤ（６，６）を除く。）のみを用いて、

により放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出することが可能である。

［算出手法３］
　欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子から出力されたダーク読取値に置き換えて空間的平均値Ｗ（４，４）の算出を行うように構成することも可能である。

　例えば、欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子（５，６）から出力されたダーク読取値Ｄ（５，６）を用いて、
　　Ｄ（６，６）←Ｄ（５，６）　　…（５４）
と置き換える。そして、上記（５２）式に従って放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出する。

　なお、この場合、欠陥画素（６，６）から出力されるダーク読取値Ｄ（６，６）を、領域Ｒの内部で欠陥画素（６，６）に近接する放射線検出素子（５，６）等から出力されるダーク読取値Ｄで置換する代わりに、例えば図２５に示す放射線検出素子（８，６）のように、領域Ｒ外の位置に存在する放射線検出素子であって欠陥画素（６，６）に近接する放射線検出素子から出力されるダーク読取値Ｄで欠陥画素（６，６）から出力されるダーク読取値Ｄ（６，６）を置換するように構成することも可能である。

［算出手法４］
　欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）を、欠陥画素（６，６）の近傍の複数の放射線検出素子から出力された各ダーク読取値Ｄで補間して空間的平均値Ｗ（４，４）の算出を行うように構成することも可能である。

　例えば、欠陥画素（６，６）から出力されたダーク読取値Ｄ（６，６）の値を、欠陥画素（６，６）の近傍の放射線検出素子（５，５）、（６，５）、（７，５）、（５，６）、（７，６）、（５，７）、（６，７）、（７，７）から出力された各ダーク読取値Ｄ（５，５）、Ｄ（６，５）、Ｄ（７，５）、Ｄ（５，６）、Ｄ（７，６）、Ｄ（５，７）、Ｄ（６，７）、Ｄ（７，７）を用いて、例えばそれらの値の平均値を算出してＤ（６，６）の値をその平均値で補間する。そして、上記（５２）式に従って放射線検出素子（４，４）の空間的平均値Ｗ（４，４）を算出する。

　なお、この場合、例えば欠陥画素が図２６に示すように、領域Ｒの端部に存在する放射線検出素子（７，６）であるような場合には、欠陥画素（７，６）から出力されたダーク読取値Ｄ（７，６）の値を、欠陥画素（７，６）の近傍の放射線検出素子（５，５）、（６，５）、（７，５）、（５，６）、（６，６）、（５，７）、（６，７）、（７，７）から出力された各ダーク読取値Ｄ（５，５）、Ｄ（６，５）、Ｄ（７，５）、Ｄ（５，６）、Ｄ（６，６）、Ｄ（５，７）、Ｄ（６，７）、Ｄ（７，７）を用いて補間するように構成することが可能である。

　また、図２６に示すように、領域Ｒ外の位置に存在する放射線検出素子であって欠陥画素（７，６）に近接する放射線検出素子（８，５）、（８，６）、（８，７）から出力されるダーク読取値Ｄ（８，５）、Ｄ（８，６）、Ｄ（８，７）をも用いて、欠陥画素（７，６）から出力されたダーク読取値Ｄ（７，６）の値を、放射線検出素子（６，５）、（７，５）、（８，５）、（６，６）、（８，６）、（６，７）、（７，７）、（８，７）から出力された各ダーク読取値Ｄ（６，５）、Ｄ（７，５）、Ｄ（８，５）、Ｄ（６，６）、Ｄ（８，６）、Ｄ（６，７）、Ｄ（７，７）、Ｄ（８，７）を用いて補間するように構成することも可能である。

　なお、上記の各算出手法では、領域Ｒ内に１つの欠陥画素、すなわち１つの点欠陥がある場合について説明したが、複数の点欠陥が存在する場合にも同様の算出手法を用いることができる。また、欠陥画素が線状に存在する線欠陥がある場合においても、同様の算出手法を用いることができる。

［間引きデータの作成］
　さて、本発明では、放射線画像検出器１の制御手段６や、後述する放射線画像生成システム３０のコンソール３１で、実写画像データＦ（ｘ，ｙ）に基づいて間引きデータを作成する。

　間引きデータの作成においては、放射線画像検出器１の制御手段６やコンソール３１は、放射線画像撮影装置１の各放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）の中から、所定の間引き率で実写画像データＦ（ｘ，ｙ）を間引くことで間引きデータを作成するようになっている。ここで、本実施形態において作成される間引き画像データの例について説明する。

　例えば、１７５μｍの画素サイズで読み取られた半切サイズ（１４×１７インチ）の実写画像データ（画素数は２０１０×２４００）を、縦横共に１／１５に圧縮する場合、送信される画像データの画素数は１３４×１６０となる。この場合、例えば、放射線画像検出器１の記憶手段７に一旦記憶された実写画像データＦ（ｘ，ｙ）や、放射線画像検出器１から送信されコンソール３１の記憶手段等に記憶された実写画像データＦ（ｘ，ｙ）が読み出される際に、縦横方向にそれぞれ１５画素おきに実写画像データＦ（ｘ，ｙ）を読み出すことで間引きデータが作成される。

　放射線画像検出器１の制御手段６で間引きデータを作成する場合、上記の例では、制御手段６は、例えば図２７Ａに例示するように、２０１０×２４００の画素による実写画像データＦ（ｘ，ｙ）を放射線検出素子（ｘ，ｙ）から読み取って生成するとともに、この実写画像データＦ（ｘ，ｙ）に基づいて、図２７Ｂに例示するような間引きデータｆ（ｘ，ｙ）を作成し、両者を対応付けて記憶手段７に保存しておく。

　なお、間引きデータｆ（ｘ，ｙ）についてはデータを作成すると同時にコンソール３１に送信するように構成することも可能である。また、上記のように縦横の間引き率を同一とすれば、画像のアスペクト比を原画像（元画像）と同一に保つことができる。なお、間引きデータｆ（ｘ，ｙ）の縮小化率（圧縮率）はここに例示したものに限定されない。

［放射線画像生成システム］
　以下に、本発明に係る放射線画像生成システムの実施形態の一例を示す。本実施形態では、放射線画像生成システム３０は、図２９に示すように、放射線画像検出器１（図１等参照）と、コンソール３１と、サーバ手段３９とを備えている。

　なお、以下では、コンソール３１が、放射線画像検出器１から送信されてきたダーク読取値に基づいて放射線画像検出器１における当該放射線画像撮影でのオフセット補正値を生成するオフセット補正値生成手段として機能する場合について説明するが、このオフセット補正値生成処理を放射線画像検出器１の制御手段６で行うように、すなわち放射線画像検出器１の制御手段６をオフセット補正値生成手段として構成することも可能である。

　放射線画像検出器１は、撮影室Ｒ１に設けられたブッキー装置３３の保持部３３ａに装填されて用いられるようになっている。ブッキー装置３３には、携帯情報端末様の小型の操作部３３ｂが設けられている。

　また、本実施形態では、放射線発生装置３４の１つの放射線源３４ａが、各ブッキー装置３３に共有されて用いられるようになっており、放射線画像検出器１をブッキー装置３３に装填して用いる場合には、対応する放射線源３４ａから放射線が照射される際に、放射線発生装置３４の放射線発生のタイミング制御が放射線画像検出器１の制御手段６と連動する仕組みになっており、放射線発生装置３４の放射線発生のタイミングに基づいて、放射線画像検出器１の各種制御が行われる。なお、各ブッキー装置３３に、それぞれ放射線発生装置３４の各放射線源３４ａを対応付けて設けるように構成することも可能である。

　放射線画像検出器１は、ブッキー装置３３に装填されずに、単独でフリーの状態で用いることもできるようになっている。その場合、例えばベッドタイプのブッキー装置３３の上に放射線画像検出器１を乗せた状態で使用したり、患者が放射線画像検出器１を手で抱えた状態で使用したりすることができる。

　放射線画像検出器１には、ダーク読取値Ｄ（ｘ，ｙ）等を無線方式で送信するための通信手段であるアンテナ装置３（図１参照）が設けられている。また、放射線画像検出器１には、放射線画像検出器１がブッキー装置３３の保持部３３ａに装填された際に保持部３３ａに設けられた図示しない電極と接続してダーク読取値Ｄ（ｘ，ｙ）等を有線方式で送信するための通信手段及び放射線画像検出器１への電源供給手段である端子１３（図２参照）が、放射線画像検出器１のアンテナ装置３や電源スイッチ１１の存在する筐体２の側面部分とは反対側（対面）の側面部分に設けられている。また、端子１３と同一側面部分には、バッテリ２１を充電するための電源供給手段である端子２２が設けられている。

　このように、放射線画像検出器１は、ブッキー装置３３に装填された際には端子１３と保持部３３ａの電極とが接続されてダーク読取値Ｄ（ｘ，ｙ）等がケーブル３２を介して有線方式で中継端末３５に送られ、中継端末３５を介してコンソール３１に送信されるようになっている。また、放射線画像検出器１は、ブッキー装置３３に装填された際には、ケーブル３２、端末１３を介して中継端末３５から放射線画像検出器１へ電力が供給されるように構成されている。

　また、ブッキー装置３３に装填されずに放射線画像検出器１が単独で用いられる場合には、アンテナ装置３を介してダーク読取値Ｄ（ｘ，ｙ）等を無線方式で送信するようになっている。撮影室Ｒ１には、放射線画像検出器１が無線方式でコンソール３１にダーク読取値Ｄ（ｘ，ｙ）等を通信する際に中継する無線アンテナ３６を備える中継端末３５が設けられている。そして、放射線画像検出器１のアンテナ装置３から無線方式で送信されたダーク読取値Ｄ（ｘ，ｙ）等の情報は、無線アンテナ３６で受信され、中継端末３５を介してコンソール３１に送信されるようになっている。

　また、放射線画像検出器１は、ブッキー装置３３に装填された際にはケーブル３２、端子１３を介して中継端末３５から電力の供給を受けるが、単独で用いられる場合には、内蔵のバッテリ２１の電力により動作するようになっている。なお、中継端末３５にはバッテリ充電用の図示しない電力供給手段である端子が設けてあり、放射線画像検出器１の端子２２を中継端末３５上の端子に接触させると、中継端末３５から図示しない端子、端子２２を介して放射線画像検出器１のバッテリ２１に電力が供給され、バッテリ２１が充電されるようになっている。しかし、この形態には限定されない。

　また、放射線画像検出器１内には、図示しないタグが内蔵されている。この場合、タグとして、いわゆるＲＦＩＤ（Radio Frequency IDentification）タグが用いられており、タグには、タグの各部を制御する制御回路や放射線画像検出器１のＩＤ等の固有情報を記憶する記憶部がコンパクトに内蔵されている。

　前室Ｒ２の入口の近傍には、放射線画像検出器１のＲＦＩＤタグを読み取るタグリーダ３７が設置されている。タグリーダ３７は、内蔵する図示しないアンテナを介して電波等に所定の指示情報を乗せて発信し、前室Ｒ２に入室し或いは退室する放射線画像検出器１を検出して、放射線画像検出器１のＩＤ等をコンソール３１に送信するようになっている。

　なお、放射線画像検出器１内にＲＦＩＤタグを内蔵させる場合の例を示したが、放射線画像検出器１内にＲＦＩＤタグを内蔵させる代わりに、放射線画像検出器１の筐体外側表面に放射線画像検出器１のＩＤ等を書き込んだバーコードを添付しておき、これをバーコードリーダーで読み取るようにしても良い。この場合、タグリーダ３７が、バーコードリーダーとなり、前室Ｒ２に入室し或いは退室する放射線画像検出器１のバーコードを読み取って、放射線画像検出器１のＩＤ等をコンソール３１に送信する。

　なお、タグリーダ３７やバーコードリーダーは前室Ｒ２の入口の近傍の代わりに、撮影室Ｒ１の入口近傍に設置し、撮影室Ｒ１への放射線画像検出器１の入退室を管理するようにしても良い。

　このように、ＲＦＩＤやバーコードを使って、前室Ｒ２や撮影室Ｒ１への放射線画像検出器１の入室或いは退室をコンソール３１に通知することで、放射線画像検出器１をどの放射線発生装置（本例の場合は放射線発生装置３４）と連動されれば良いかを、コンソール３１や放射線発生装置３４に自動的に知らせることができる。このようなＲＦＩＤやバーコードを使った放射線画像検出器１の撮影室等への入室或いは退室の管理は、撮影室や放射線発生装置が複数存在する施設で有効に機能する。なお、ＲＦＩＤやバーコードを使用しない場合は、放射線画像検出器１の撮影室等への入室或いは退室の情報を、使用者が直接コンソール３１へ入力するようにしても良い。

　前室Ｒ２には、欠陥画素判定システム３０や放射線画像生成システム全体の制御を行うコンソール３１が設けられており、コンソール３１には、前述した中継端末３５やタグリーダ３７、放射線発生装置３４の本体部３４ｂ等が接続されており、また、中継端末３５を介してブッキー装置３３等が接続されている。

　コンソール３１は、図示しないＣＰＵ（Central Processing Unit）やＲＯＭ、ＲＡＭ、ハードディスク等の記憶手段、入出力インターフェース等がバスに接続されたコンピュータで構成され、実写画像データＦ（ｘ、ｙ）や間引きデータｆ（ｘ，ｙ）に対して処理を行う図示しない画像処理部を備えている。また、コンソール３１には、ＣＲＴ（Cathode Ray Tube）やＬＣＤ（Liquid Crystal Display）等からなる表示部３１ａが設けられており、表示部３１ａの表示画面には、前述した間引きデータｆ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）、ゲイン及びオフセット補正処理や、欠陥や異常画素に対する補正処理や、撮影部位に応じた階調処理等の各種の画像処理が行われ最終的に生成された診断画像データＦ_Ｏ（ｘ，ｙ）等が表示されるようになっている。

　また、コンソール３１には、ネットワークＮＷを介してコンピュータからなるサーバ手段３９が接続されている。また、サーバ手段３９は、ハードディスク等からなる記憶手段３８が接続されている。

　記憶手段３８及びコンソール３１自身の図示しない記憶手段には、放射線画像生成システムで使用可能な各放射線画像検出器１について、前述した空間的平均値を算出するために、各放射線検出素子（ｘ，ｙ）にどの放射線検出素子（ｘ´，ｙ´）を対応付けるかの情報が放射線画像検出器１のＩＤに対応付けられて予めそれぞれ記憶されている。

　また、コンソール３１は、放射線画像検出器１から、キャリブレーション時に行われたダーク読取で各放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）や、放射線画像撮影の直前又は直後に１回以上行われたダーク読取において得られたダーク読取値Ｄ（ｘ，ｙ）、Ｄ_ｍ（ｘ，ｙ）が有線方式或いは無線方式が送信されてくると、それらを当該放射線画像検出器１のＩＤに対応付けて、自らの記憶手段やサーバ手段３９の記憶手段３８に記憶させるようになっている。

　さらに、コンソール３１は、当該放射線検出器１の各放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）や、各放射線検出素子（ｘ，ｙ）に対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力された各ダーク読取値ｄ_ｍ（ｘ´，ｙ´）の温度補償変数としての空間的平均値ｗ_ｍ（ｘ´，ｙ´）の時間的統計値（時間的平均値）ω（ｘ，ｙ）を算出してサーバ手段３９の記憶手段３８に保存するようになっている。

　なお、このダーク読取値ｄ（ｘ，ｙ）の保存及び時間的平均値δ（ｘ，ｙ）等の算出、保存を、コンソール３１やコンソール３１の記憶手段で行うように構成することも可能である。

　コンソール３１は、ＲＯＭやハードディスク等の記憶手段に格納されている各種処理の実行に必要なプログラムを読み出してＲＡＭの作業領域に展開して、プログラムに従って処理を実行するようになっている。また、本実施形態では、コンソール３１の記憶手段やサーバ手段３９の記憶手段３８（以下、単に記憶手段という。）には、放射線画像検出器１ごとに欠陥画素の画素位置（ｘ，ｙ）等の情報が登録された欠陥画素マップが記憶されている。

　また、ＲＯＭやハードディスク等の記憶手段には、前述した欠陥画素判定を実行するための欠陥画素判定プログラムや、異常画素判定を実行するための異常画素判定プログラム、間引きデータｆ（ｘ，ｙ）を表示させるための間引きデータ表示プログラム、放射線画像撮影で得られた実写画像データＦ（ｘ，ｙ）から最終的な診断画像データＦ_Ｏ（ｘ，ｙ）を生成するための放射線画像生成プログラム等が格納されている。

　コンソール３１は、キャリブレーション時などには、欠陥画素判定プログラムを読み出して、上記の欠陥画素判定方法を実行するようになっている。なお、欠陥画素判定プログラムにおいて、上記の欠陥画素判定方法の［判定手法１］～［判定手法６］のうち、いずれの判定手法が実行され、或いはそれらの判定手法が組み合わされて実行されるかは、適宜決められる。

　コンソール３１は、例えば、欠陥画素判定プログラムに従って前述した欠陥画素判定方法のうち［判定手法３］で説明した欠陥画素判定を行う場合、放射線画像検出器１から送信されてきたダーク読取値ｄ_ｍ（ｘ，ｙ）と予め定められた閾値ｄ_ｍthとを比較し、少なくとも１つのダーク読取値ｄ_ｍ（ｘ，ｙ）が閾値ｄ_ｍthよりも大きい場合（上記（２９）式参照）に、その放射線検出素子（ｘ，ｙ）を欠陥画素と判定して、記憶手段に保存されている欠陥画素マップに登録する。

　また、コンソール３１は、例えば、欠陥画素判定プログラムに従って前述した欠陥画素判定方法のうち［判定手法６］で説明した欠陥画素判定を行う場合、放射線画像検出器１から送信されてきたダーク読取値ｄ_ｍ（ｘ，ｙ）に基づいて、一の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値ｄ_ｍ（ｘ，ｙ）等と、当該一の放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子から出力された各ダーク読取値の空間的平均値ｗ_ｍ（ｘ，ｙ）との差分ｅ_ｍ（ｘ，ｙ）の絶対値｜ｅ_ｍ（ｘ，ｙ）｜を算出し、それと予め定められた閾値ｅ_ｍthとを比較し、差分の絶対値｜ｅ_ｍ（ｘ，ｙ）｜が閾値ｅ_ｍthよりも大きい場合（上記（３７）式参照）に、その放射線検出素子（ｘ，ｙ）を欠陥画素と判定して、記憶手段に保存されている欠陥画素マップに登録する。

　なお、上記の例では［判定手法３］や［判定手法６］を実施する場合について説明したが、欠陥画素判定方法の他の［判定手法１］や［判定手法２］、［判定手法４］、［判定手法５］を実施する場合も同様であり、また、それらを組み合わせた判定手法を実施する場合も同様である。

　また、前述したように、［判定手法７］～［判定手法１１］の異常画素判定方法においても、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）がひどく異常な値であったり、頻繁に異常値が出力される場合には、欠陥画素と見なされて記憶手段に保存されている欠陥画素マップに登録される。

　従って、欠陥画素マップには、放射線画像検出器１の製造後の出荷検査や、定期的なキャリブレーション時に検出される定常的な欠陥画素のみならず、各回ごとの放射線画像撮影の前又は後に１回以上行われるダーク読取で得られたダーク読取値Ｄ（ｘ，ｙ）や、実際の放射線画像撮影で取得された実写画像データＦ（ｘ，ｙ）、或いは被写体が存在しない状態で行われた放射線画像撮影で取得された実写画像データＦ（ｘ，ｙ）に基づいて欠陥画素と見なされた画素が含まれており、常時更新される。

　また、コンソール３１は、通常の放射線画像撮影時には、記憶手段から放射線画像生成プログラムのほか、間引きデータ表示プログラムや異常画素判定プログラム等を読み出して、同時並行で処理を行うようになっている。

　放射線画像生成プログラムは、最終的な画像データＦ_Ｏ（ｘ，ｙ）等を生成するためのプログラムであり、図２８に示すように、欠陥画素判定ステップ（ステップＳ１）や、実写画像データ取得ステップ（ステップＳ２）や、ダーク読取ステップ（ステップＳ３）、間引きデータ作成ステップ（ステップＳ４）、温度補償変数算出ステップ（ステップＳ５、Ｓ９）、時間的平均値等算出ステップ（ステップＳ６、Ｓ１０）、オフセット補正値算出ステップ（ステップＳ７、Ｓ１１）、間引き画像表示ステップ（ステップＳ８）、画像補正ステップ（ステップＳ１２）等からなる。コンソール３１は、この放射線画像生成プログラムを読み出してＲＡＭの作業領域に展開して、この放射線画像生成プログラムに従って必要な処理を実行し、また、装置各部を制御するようになっている。

　欠陥画素判定ステップ（ステップＳ１）では、前述したように欠陥画素判定プログラムに基づいて欠陥画素判定が行われ、欠陥画素と判定された放射線検出素子（ｘ，ｙ）が放射線画像検出器１に対応付けられて欠陥画素マップに記憶される。従って、欠陥画素判定ステップ（ステップＳ１）は、今回の放射線画像撮影の際に行われるというより、今回の放射線画像撮影の準備作業として行われる工程といえる。

　コンソール３１は、放射線技師や医師等から放射線画像撮影の指示を受けると、撮影条件（管電流、管電圧等）や照射開始を制御する放射線発生装置３４を駆動する等して実写画像データ取得ステップ（ステップＳ２）を実行して各放射線検出素子（ｘ，ｙ）から実写画像データＦ（ｘ，ｙ）を取得させる。また、放射線画像撮影の前又は後に、少なくとも１回のダーク読取を行うように放射線画像検出器１に指示信号を送信して、ダーク読取ステップ（ステップＳ３）を実行して各放射線検出素子（ｘ，ｙ）からダーク読取値Ｄ（ｘ，ｙ）を取得させる。また、間引きデータ作成ステップ（ステップＳ４）を実行して前述したように実写画像データＦ（ｘ，ｙ）から間引きデータｆ（ｘ，ｙ）を作成させる。

　なお、放射線画像撮影の前又は後にダーク読取を１回だけ行うように構成してもよく、また、２回以上行うように構成することも可能である。また、ダーク読取が放射線画像撮影の前に行われる場合には、実写画像データ取得ステップ（ステップＳ２）とダーク読取ステップ（ステップＳ３）の順番が逆になる。

　また、間引きデータ作成ステップ（ステップＳ４）は実写画像データ取得ステップ（ステップＳ２）の後、適宜のタイミングで行われる。さらに、放射線画像検出器１から実写画像データＦ（ｘ，ｙ）をコンソール３１に送信し、コンソール３１で間引きデータ作成ステップ（ステップＳ４）を実行して間引きデータｆ（ｘ，ｙ）を作成するように構成することも可能である。

　本実施形態では、放射線画像検出器１は、間引きデータ作成ステップ（ステップＳ４）を終了すると、自らのＩＤと、ダーク読取値Ｄ（ｘ，ｙ）と、間引きデータｆ（ｘ，ｙ）とを有線方式又は無線方式でコンソール３１に送信する。コンソール３１は、放射線画像検出器１からダーク読取値Ｄ（ｘ，ｙ）と間引きデータｆ（ｘ，ｙ）とが送信されてくると、それらを記憶手段に保存する。

　なお、この段階で、送信された間引きデータｆ（ｘ，ｙ）に基くサムネール画像（原画像サイズに対し１／１５のサイズの画像）を表示することとしても良い。この場合、オフセット／ゲイン補正や異常／欠陥画素補正は未だ実施されていないが、ポジショニングの確認には充分な画像であるので、最も早く再撮影要否を判断可能となる。

　コンソール３１は、続いて、温度補償変数算出ステップ（ステップＳ５）を実行する。コンソール３１は、放射線画像検出器１のＩＤを参照して、記憶手段から各放射線検出素子（ｘ，ｙ）に対応付ける複数の放射線検出素子（ｘ´，ｙ´）の情報を読み出し、それらから出力された各ダーク読取値Ｄ（ｘ´，ｙ´）を用いて、上記（３）式に従って各放射線検出素子（ｘ，ｙ）についてそれぞれ温度補償変数Ｗ（ｘ，ｙ）を算出する。

　その際、コンソール３１は、この段階では、間引きデータｆ（ｘ，ｙ）を出力した放射線検出素子（ｘ，ｙ）のみについて温度補償変数Ｗ（ｘ，ｙ）を算出する。また、一の放射線検出素子（ｘ，ｙ）について温度補償変数Ｗ（ｘ，ｙ）を算出するにあたって、当該一の放射線検出素子（ｘ，ｙ）に予め対応付ける他の放射線検出素子の中に欠陥画素が含まれる場合には、前述した［算出手法１］～［算出手法４］のいずれかの手法を用いて欠陥画素から出力されたダーク読取値Ｄ（ｘｓ，ｙｓ）をその近傍の放射線検出素子から出力されたダーク読取値で置換したり補間するなどして、当該一の放射線検出素子（ｘ，ｙ）の温度補償変数Ｗ（ｘ，ｙ）を算出する。

　また、コンソール３１は、予め過去のキャリブレーション時等に得られたデータに基づいて時間的平均値等算出ステップ（ステップＳ６）を実行するが、本実施形態では、当該放射線画像検出器１についての時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗ_ｍ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）は、前述したように、予めサーバ手段３９で予め計算され、記憶手段３８に保存されている。コンソール３１がこの段階で自らそれらの値を算出するように構成してもよい。

　コンソール３１は、続いて、オフセット補正値算出ステップ（ステップＳ７）を実行する。その際、コンソール３１は、放射線画像検出器１のＩＤと、間引きデータｆ（ｘ，ｙ）を出力した放射線検出素子（ｘ，ｙ）の画素番号（ｘ，ｙ）の情報をネットワークＮＷを介してサーバ手段３９に送信する。サーバ手段３９は、間引きデータｆ（ｘ，ｙ）を出力した放射線検出素子（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）や温度補償変数としての空間的平均値ｗ_ｍ（ｘ，ｙ）の時間的平均値ω（ｘ，ｙ）を記憶手段３８から読み出して、コンソール３１に送信する。

　コンソール３１は、間引きデータｆ（ｘ，ｙ）を出力した放射線検出素子（ｘ，ｙ）について、サーバ手段３９から入手したダーク読取値ｄ_ｍ（ｘ，ｙ）の時間的平均値δ（ｘ，ｙ）と、温度補償変数としての空間的平均値ｗ_ｍ（ｘ´，ｙ´）の時間的平均値ω（ｘ，ｙ）とから、それらの差分ε（ｘ，ｙ）（＝δ（ｘ，ｙ）－ω（ｘ，ｙ））を算出する。そして、算出したε（ｘ，ｙ）と、温度補償変数算出ステップ（ステップＳ５）で算出した当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）についての各温度補償変数Ｗ（ｘ，ｙ）とを加算して、間引きデータｆ（ｘ，ｙ）を出力した放射線検出素子（ｘ，ｙ）に対するオフセット補正値Ｏ（ｘ，ｙ）を算出する。

　また、コンソール３１は、当該放射線検出素子（ｘ，ｙ）のゲイン補正値Ｇ（ｘ，ｙ）を記憶手段から読み出し、上記（１）式に従ってゲイン補正処理及びオフセット補正処理を行って間引きデータｆ（ｘ，ｙ）に対する処理を行い、処理後の画像データｆ_Ｏ（ｘ，ｙ）を算出する。なお、この処理では、上記（１）式におけるＦ_Ｏ（ｘ，ｙ）及びＦ（ｘ，ｙ）をそれぞれｆ_Ｏ（ｘ，ｙ）、ｆ（ｘ，ｙ）と読み替えるものとする。

　続いて、コンソール３１は、間引き画像表示ステップ（ステップＳ８）を実行し、上記のように、間引きデータｆ（ｘ，ｙ）から算出した、補正処理後の画像データｆ_Ｏ（ｘ，ｙ）に基づく間引き画像を、表示部３１ａの表示画面に表示する。その際、表示部３１ａの表示画面には、サムネール画像（原画像に対し１／１５に縮小された間引き画像）が表示される。なお、必要に応じて、前記サムネール画像（間引き画像）に、撮影した部位に応じた階調処理や周波数強調処理等の画像処理を施して表示することとしても良い。

　この段階で表示される画像においては、異常・欠陥画素に対する補正処理は、未だ実行されていないが、オフセット／ゲイン補正（場合によっては諧調処理等を含む）されているので、各画素の出力値が適正出力範囲内にあるか否かを視認可能となる、この段階でも
画質的な（各画素の出力範囲的な）再撮影要否を判断できる。

　なお、この段階で表示する画像データは間引かれたデータであるが、各画素の出力範囲の視認性を向上させる為に、間引かれたデータを拡大補間して原画像データと同サイズで表示することとしても良いが、この場合、間引かれたデータに基く拡大補間なので、解像度は低下するのはやむをえない（図１３（Ｂ）参照）。

　この段階で、放射線技師や医師等は、コンソール３１の表示部３１ａの表示画面に表示された間引き画像を確認して、放射線画像中の被写体の撮影位置が適切な位置であるか否か等を判断し、再撮影が必要か否かを判断する。そして、再撮影が必要であると判断した場合には、実写画像データ取得ステップ（ステップＳ２）以降の処理を再度繰り返す。

　また、放射線技師や医師等は、コンソール３１の表示部３１ａの表示画面に表示された間引き画像を確認して、再撮影の必要はないと判断すると、コンソール３１に対して、最終的な診断画像データＦ_Ｏ（ｘ、ｙ）を得るべく、処理を続行するように指示を出す。

　コンソール３１は、放射線技師等から続行の指示を受けると、続いて、放射線画像検出器１から全放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を入手していなければ、当該放射線画像検出器１に送信要求信号を発信して、実写画像データＦ（ｘ，ｙ）を送信させる。

　そして、コンソール３１は、温度補償変数算出ステップ（ステップＳ９）を実行し、すでに記憶手段から読み出してある各放射線検出素子（ｘ，ｙ）に対応付ける複数の放射線検出素子（ｘ´，ｙ´）の情報に基づいて、それらから出力された各ダーク読取値Ｄ（ｘ´，ｙ´）を用いて、上記（３）式に従って各放射線検出素子（ｘ，ｙ）についてそれぞれ温度補償変数Ｗ（ｘ，ｙ）を算出する。この場合も、前述した［算出手法１］～［算出手法４］の手法を用いて置換や補間等の処理が適宜行われる。

　なお、この温度補償変数算出ステップ（ステップＳ９）からオフセット補正値算出ステップ（ステップＳ１１）までの処理は、間引きデータｆ（ｘ，ｙ）に対応する放射線検出素子（ｘ，ｙ）については既に行っているため、改めて行う必要はない。

　続いて、コンソール３１は、時間的平均値等算出ステップ（ステップＳ１０）を実行してサーバ手段３９の記憶手段３８から当該放射線画像検出器１の各放射線検出素子（ｘ，ｙ）についての時間的平均値δ（ｘ，ｙ）や温度補償変数ω（ｘ，ｙ）を入手し、オフセット補正値算出ステップ（ステップＳ１１）を実行して、上記と同様の手順で各放射線検出素子（ｘ，ｙ）についてオフセット補正値Ｏ（ｘ，ｙ）を算出する。

　そして、コンソール３１は、各放射線検出素子（ｘ，ｙ）のゲイン補正値Ｇ（ｘ，ｙ）を記憶手段から読み出し、上記（１）式に従ってゲイン補正処理及びオフセット補正処理を行って各実写画像データＦ（ｘ，ｙ）に対する処理を行い、処理後の画像データＦ_Ｏ（ｘ，ｙ）を算出する。

　本実施形態では、コンソール３１は、続いて、画像補正ステップ（ステップＳ１２）を実行するようになっている。

　画像補正ステップ（ステップＳ１２）では、前述した欠陥画素判定や異常画素判定が行われ、欠陥画素や異常画素と判定された放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）に基づく画像データＦ_Ｏ（ｘ，ｙ）に対して置換や補間等の画像補正がなされるようになっている。

　具体的には、コンソール３１は、例えば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）の中に、定常的な欠陥画素や上述の［判定手法１］～［判定手法６］に従って欠陥画素と判定された画素が存在する場合、欠陥画素マップに登録されたそれらの欠陥画素に対応する放射線検出素子（ｘｓ，ｙｓ）から出力された実写画像データＦ（ｘ，ｙ）に基づく画像データＦ_Ｏ（ｘ，ｙ）に対して置換や補間等の画像補正を行う。

　また、コンソール３１は、例えば、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）に基づいて上述の［判定手法７］や［判定手法８］に従って異常画素であると判定された画素が存在する場合や、放射線画像検出器１の放射線検出素子（ｘ，ｙ）から出力された実写画像データＦ（ｘ，ｙ）に基づいて上述の［判定手法９］や［判定手法１０］に従って異常画素であると判定された画素が存在する場合には、それらの放射線検出素子（ｘ，ｙ）の情報を記憶手段に一時的に記憶させる。そして、それと同時に、欠陥画素同等と見なして、欠陥画素と同様に周囲の画素値を用いて置換処理を行ったり補間処理を行って画像補正を行う。

　なお、その際、放射線検出素子（ｘ，ｙ）から出力されるダーク読取値Ｄ（ｘ，ｙ）や実写画像データＦ（ｘ，ｙ）がひどく異常な値であったり、頻繁に異常値が出力される場合には、当該放射線検出素子（ｘ，ｙ）は欠陥画素と見なされて、記憶手段に保存されている欠陥画素マップに登録される。

　以上のオフセット補正処理、ゲイン補正処理、欠陥画素補正処理、異常画素補正処理を行った後、撮影部位に応じた階調処理や周波数強調処理、粒状抑制処理等の画像処理を行い、最終的な診断用画像を生成する。

　以上のように、本発明に係る放射線画像生成システム３０によれば、コンソール３１は、間引きデータ、或いは、オフセット／ゲイン処理済の間引きデータに対し、欠陥や異常画素に対する補正処理を行わずに表示するので、迅速に再撮影要否確認用の間引き画像を表示でき、技師は迅速に再撮影要否を判断できる。

　そのため、放射線技師等が表示された間引きデータを確認して再撮影の要否を即座に判断することが可能となるとともに、患者も再撮影の要否の判断を長い間、待つ必要がなくなるため、患者の負担が軽減される。

　また、放射線検出素子（ｘ，ｙ）の実写画像データＦ（ｘ，ｙ）を補正するためのオフセット補正値Ｏ（ｘ，ｙ）が、その放射線検出素子（ｘ，ｙ）から出力されたダーク読取値Ｄ（ｘ，ｙ）と、その放射線検出素子（ｘ，ｙ）に予め対応付けられた複数の放射線検出素子（ｘ´，ｙ´）から出力されたダーク読取値Ｄ（ｘ´，ｙ´）とを用いて算出されるため、ダーク読取値Ｄ（ｘ，ｙ）の取得が放射線画像撮影の前又は後に１回程度行われたダーク読取であっても、十分に的確なオフセット補正値Ｏ（ｘ，ｙ）を得ることが可能となり、画素（放射線検出素子）ごとの特性ばらつきが補正されて、ＳＮ比が良好な最終的な診断画像データＦ_Ｏ（ｘ，ｙ）を生成することが可能となる。

　さらに、上記の異常画素判定方法や欠陥画素判定方法を用いることで、最終的な画像データＦ_Ｏ（ｘ，ｙ）から異常画素や欠陥画素の影響を的確に排除して、良好な最終的な画像データＦ_Ｏ（ｘ，ｙ）を生成することが可能となる。

　なお、オフセット補正処理、ゲイン補正処理が終了して最終的な画像データＦ_Ｏ（ｘ，ｙ）が算出されると、一般的には最終的な画像データＦ_Ｏ（ｘ，ｙ）に対数変換処理がかけられ、放射線強度の対数に対してリニアなデータに変換される。そして対数変換後の最終的な画像データをＬＯＧ（Ｆ_Ｏ（ｘ，ｙ））と表すと、ＬＯＧ（Ｆ_Ｏ（ｘ，ｙ））に対して各種画像処理が施される。

　従って、最終的な画像データＦ_Ｏ（ｘ，ｙ）の信号値の値が小さい部分ほど、対数変換によって信号値が大きく増幅され、最終的な画像データＦ_Ｏ（ｘ，ｙ）の信号値の値が大きい部分ほど、対数変換によって信号値が圧縮される。従って、［判別手法１］～［判定手法１１］で定められた各種閾値は、この対数変換処理によって信号値最終的にどのように増幅されたり圧縮されたりするのかを考慮の上決定することが好ましい。また、ダーク読取値や実写画像データ、最終的な画像データＦ_Ｏ（ｘ，ｙ）の値などに応じて、各種閾値を変化させたりすることも可能であり、適宜設定される。

　また、本実施例では、放射線画像撮影の前又は後にダーク読取を１回行う場合に特に有効な手法として説明を行ったが、ダーク読取を２回以上行う場合であっても同様な効果を得ることができる。すなわち、キャリブレーション時に実施したダーク読取り回数をＭ回、放射線画像撮影の前又は後に行うダーク読取り回数をＫ回とした場合に、Ｍ＞Ｋの関係が成立すれば、本発明で説明したのと同様の効果が得られることは言うまでもない。

　すなわち、Ｋ＞２の場合は、ダーク読取値をＤｋ（ｘ，ｙ）（ｋ=１～Ｋ、Ｋ＞２）とおき、

を本発明の実施例のＤ（ｘ，ｙ）にあてはめて考えれば良い。すなわち、
　　Ｄ（ｘ，ｙ）=Ｄｋ_ave（ｘ，ｙ）　　…（５６）
として考えれば良い。

　また、当該放射線画像検出器１の放射線検出素子（ｘ，ｙ）中に欠陥画素がある場合でも、上記の空間的統計値（空間的平均値）の算出手法を用いることで、欠陥画素から出力されるダーク読取値Ｄが適切に置換され或いは補間されて、有効な空間的統計値（空間的平均値）である温度補償変数Ｗ（ｘ，ｙ）やｗ_ｍ（ｘ´，ｙ´）を算出することが可能となる。また、それにより、オフセット補正値Ｏ（ｘ，ｙ）を的確に算出することが可能となる。

　なお、本発明が上記の実施の形態に限定されず、適宜変更可能であることはいうまでもない。本発明の概念は、一の放射線検出素子と同じように温度変動する放射線検出素子（ｘ´，ｙ´）を予め選択して当該一の放射線検出素子（ｘ，ｙ）に予め対応付けるものである。そのため、放射線検出素子自体は、上記の実施形態のように放射線検出素子タイプである必要はなく、他の構造を有するものを２次元的に配置したセンサパネル部を備える放射線画像検出器にも適用可能であることはいうまでもない。

　また、本発明では、放射線検出素子の画素単位での温度変動を取り上げて説明を行ったが、温度変動以外の因子であっても、ある画素値が近傍の画素値と類似の変化をもたらすような変動であれば、同様の手法を適応することで同様の効果を得られることは言うまでもない。

　さらに、本実施形態では、データ処理をコンソール３１で行う場合について説明したが、実写画像データＦ（ｘ，ｙ）から各画素ごとに各オフセット補正値Ｏ（ｘ，ｙ）をそれぞれ差し引き、実写画像データＦ（ｘ、ｙ）を補正して最終的な画像データＦ_Ｏ（ｘ，ｙ）を生成するまでの全てのデータ処理を放射線画像検出器１で行うように構成することも可能である。

　さらに、例えば、上記の実施形態においては、放射線画像撮影時とは別に行われるキャリブレーション時にダーク読取を行うこととしたが、放射線画像撮影の前又は後に行うダーク読取値を保存しておき、それを用いるように構成してもよい。この場合には、例えば連続して行われた放射線画像撮影の複数回分のダーク読取値を正規化する等して保存しておき、当該複数回分のダーク読取値を用いて上記の実施形態と同様に時間的統計値（時間的平均値）δ（ｘ，ｙ）やω（ｘ，ｙ）を算出するようにしてもよい。

　また、本実施例では、オフセット補正を例に取り、説明を行ったが、ゲイン補正値Ｇ（ｘ，ｙ）を求める際にも、放射線画像検出器に対して所定の条件で一様な放射線を照射し、この読取値（本実施例の中の実写画像データＦ（ｘ，ｙ）に相当）に対してもオフセット補正が適応される。このため、本実施例では、実写画像データＦ（ｘ，ｙ）に対してオフセット補正を行う場合を例に取り説明を行ったが、前記理由によりゲイン補正値Ｏ（ｘ，ｙ）を求める際にも同様な処置が実施可能であり、かつ同様の効果があることは言うまでもない。

　さらに、（ａ）間引きデータ又は実写画像データに基づいて作成した間引きデータに基づく画像と、（ｂ）間引きデータ又は実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及びオフセット補正値に基づく画像処理を行った画像処理後の間引きデータの画像と、（ｃ）実写画像データに対してゲイン補正処理及びオフセット補正値に基づく画像処理を行った画像処理後の画像との、いずれの画像をどのような順番でコンソール３１の表示部３１ａに表示するかを予め設定したり、或いはコンソール３１を操作することによって設定する等して、放射線画像生成システム３０の運用状況等に応じて操作者が適宜選択できるように構成することも可能である。

　（ａ）間引きデータ又は実写画像データに基づいて作成した間引きデータに対してオフセット／ゲイン補正や異常／欠陥画素補正を実施せずに表示すれば、表示されるまでの時間が最も短い時間で表示させることが可能となる。また、ポジショニングの確認を充分的確に行うことが可能となり、再撮影要否を即座に判断することが可能となる。

　その際、放射線画像は、従来、スクリーン／フィルムや輝尽性蛍光体シート等を用いて撮像されていた際には、各画素ごとに読み取られる画像データが、照射された放射線の線量の対数（log）に比例する形で読み取られることが多く、放射線技師等の操作者も対数表示に慣れている。そのため、間引きデータ又は実写画像データに基づいて作成した間引きデータを単にコンソール３１の表示部３１ａに表示する代わりに、各画素の画像データを対数変換処理して表示することが好ましい。

　また、（ｂ）間引きデータ又は実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及びオフセット補正値に基づく画像処理を行った後、画像処理後の間引きデータの画像をコンソール３１の表示部３１ａに表示する場合、間引きデータを単に、或いは対数変換処理した後に表示する場合に比べて、表示されるまでに時間を要する。しかし、例えば、オフセット／ゲイン補正し、さらに撮影部位に応じた階調処理等を行って表示することで、ポジショニングの確認のみならず、画像データが飽和しているか否かの確認を同時に行うことが可能となる。

　さらに、（ｃ）実写画像データに対してゲイン補正処理及びオフセット補正値に基づく画像処理や上記のような異常画素判定方法や欠陥画素判定方法を用いて適切に補正して画像処理後の画像をコンソール３１の表示部３１ａに表示する場合、上記の（ａ）、（ｂ）の場合よりも表示されるまでに時間を要する。

　しかし、上記の画像処理を施し、また、必要に応じて階調処理や周波数強調処理、粒状抑制処理等の撮影部位に応じた画像処理等を行うことで、前述したように、最終的な画像データから異常画素や欠陥画素の影響を的確に排除して、良好な最終的な画像データを生成することが可能となる。

　その際、上記の（ａ）～（ｃ）の各画像に対して画像処理に要する時間が異なり、コンソール３１の表示部３１ａに表示されるまでに要する時間が異なることを利用して、例えば、上記（ａ）～（ｃ）の各画像処理を同時に行い、表示部３１ａに表示する画像を、画像処理が終了した順に（ａ）、（ｂ）、（ｃ）の順番で順次切り替えて表示するように構成することも可能である。また、その際、操作者の操作により切り替えられて表示されなくなった画像を再表示できるように構成することも可能である。

　医療の分野において、診断用の放射線画像を撮影するための放射線画撮影システムに利用可能性がある。

符号の説明

１　放射線画像検出器
３　アンテナ装置（通信手段）
４、５、６　センサパネル部、読取部及び制御手段（画像データ取得手段）
６　制御手段（オフセット補正値生成手段）
１３　端子（通信手段）
１４、（ｘ，ｙ）　放射線検出素子
２１　バッテリ
３０　放射線画像生成システム
３１　コンソール（オフセット補正値生成手段）
３１ａ　表示部
３４　放射線発生装置
３８　記憶手段
Ａ（ｘ，ｙ）　実写画像データの空間的平均値
Ｄ（ｘ，ｙ）　ダーク読取値
Ｅ（ｘ，ｙ）　差分、第１差分
Ｅth　閾値
Ｆ（ｘ，ｙ）　実写画像データ
ｆ（ｘ，ｙ）　間引きデータ
Ｏ（ｘ，ｙ）　オフセット補正値
Ｖ（ｘ，ｙ）　差分
Ｖth　閾値
Ｖ´th　閾値
Ｗ（ｘ，ｙ）　空間的平均値
（ｘｓ，ｙｓ）　欠陥画素
δ（ｘ，ｙ）　ダーク読取値の時間的平均値
ε（ｘ，ｙ）　第２差分
｜Λ（ｘ，ｙ）｜　第１差分と第２差分との差分の絶対値
Λth　閾値
ω（ｘ，ｙ）　ダーク読取値の空間的平均値の時間的平均値

Claims

　２次元状に配置された複数の放射線検出素子と、
　放射線画像撮影では前記複数の放射線検出素子から実写画像データを取得し、放射線が照射されない状態で行われるダーク読取では前記複数の放射線検出素子からダーク読取値を取得する信号値取得手段と、
　前記実写画像データに基づいて間引きデータを作成可能な制御手段と、
　前記実写画像データ、前記ダーク読取値及び前記間引きデータを送信可能な通信手段と、
を有する放射線画像検出器と、
　前記放射線画像検出器から送信された前記実写画像データ及び前記間引きデータに対して画像処理を行う画像処理部と、
　前記実写画像データ、前記間引きデータに基づく画像、又はそれらに対して画像処理を行った後の画像を表示する表示部と、
を有するコンソールと、
　前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値に基づいて、前記放射線画像検出器における当該放射線画像撮影でのオフセット補正値を生成するオフセット補正値生成手段と、
　前記放射線画像検出器に対して放射線を照射する放射線発生装置と、
を備え、
　前記放射線画像検出器は、放射線画像撮影の前又は後に少なくとも１回行われたダーク読取における前記ダーク読取値と、放射線画像撮影に基づく前記実写画像データ及び前記実写画像データに基づく前記間引きデータのうちの少なくとも１つとを前記コンソールに送信し、
　前記コンソールは、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに基づく画像、或いは、前記間引きデータ、又は前記実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行って、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに対する画像処理後の画像、或いは、前記実写画像データに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行った画像処理後の画像、のうちの少なくとも１つを、前記表示部に表示させることを特徴とする放射線画像生成システム。
　前記オフセット補正値生成手段は、前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値とに基づいて、前記放射線画像検出器における当該放射線画像撮影でのオフセット補正値を生成することを特徴とする請求項１に記載の放射線画像生成システム。
　前記オフセット補正値生成手段は、前記コンソールに設けられていることを特徴とする請求項１又は請求項２に記載の放射線画像生成システム。
　前記コンソールは、
　前記放射線画像検出器における欠陥画素が登録された欠陥画素マップを記憶する記憶手段を備え、
　前記放射線画像検出器から前記実写画像データが送信されていなければ前記実写画像データを送信させるとともに、
　前記実写画像データの中に、前記欠陥画素マップに登録された前記欠陥画素に対応する放射線検出素子から出力されたデータが存在する場合には、前記実写画像データに対して前記ゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行った後の画像における当該放射線検出素子から出力されたデータの位置の前記画像処理後のデータを補正する処理を行うことを特徴とする請求項１から請求項３のいずれか一項に記載の放射線画像生成システム。
　前記コンソールは、
　前記放射線画像検出器から前記実写画像データが送信されていなければ前記実写画像データを送信させるとともに、
　前記実写画像データの中に、前記放射線画像検出器から送信された前記ダーク読取値に基づいて異常なダーク読取値を出力する異常画素であると判定される放射線検出素子から出力されたデータが存在する場合には、前記実写画像データに対して前記ゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行った後の画像における当該放射線検出素子から出力されたデータの位置の前記画像処理後のデータを補正する処理を行うことを特徴とする請求項１から請求項４のいずれか一項に記載の放射線画像生成システム。
　前記コンソールは、前記放射線画像検出器から送信された前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値とに基づいて、当該一の放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定することを特徴とする請求項５に記載の放射線画像生成システム。
　前記コンソールは、前記放射線画像検出器から送信された前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項６に記載の放射線画像生成システム。
　前記コンソールは、
　前記放射線画像検出器から送信された前記ダーク読取値に基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記ダーク読取値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値との第１差分を算出し、
　過去のキャリブレーション時に複数回行われた前記ダーク読取において当該一の放射線検出素子から出力された前記ダーク読取値の時間的平均値と、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記ダーク読取値の空間的平均値の時間的平均値との第２差分を算出し、
　前記第１差分と前記第２差分との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項６に記載の放射線画像生成システム。
　前記コンソールは、前記放射線画像検出器から前記実写画像データが送信された場合には、当該実写画像データに基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値とに基づいて、当該一の放射線検出素子が異常な信号値を出力する異常画素であるか否かを判定することを特徴とする請求項５から請求項８のいずれか一項に記載の放射線画像生成システム。
　前記コンソールは、前記実写画像データに基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項９に記載の放射線画像生成システム。
　前記コンソールは、前記実写画像データに基づいて、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分の絶対値が、予め設定された閾値を越えた場合に、当該一の放射線検出素子が異常な信号値を出力する異常画素であると判定することを特徴とする請求項９に記載の放射線画像生成システム。
　前記欠陥画素マップに登録される欠陥画素には、予め被写体が存在しない状態で行われた前記放射線画像撮影において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分の絶対値が予め設定された閾値を越えた場合に異常な実写画像データを出力する異常画素であると判定された前記放射線検出素子であって、出力される実写画像データの大きさに応じて欠陥画素であると判定された前記放射線検出素子が含まれることを特徴とする請求項４に記載の放射線画像生成システム。
　前記欠陥画素マップに登録される欠陥画素には、予め被写体が存在しない状態で行われた前記放射線画像撮影において、前記放射線画像検出器の一の前記放射線検出素子から出力された前記実写画像データと、前記一の放射線検出素子に予め対応付けられた複数の放射線検出素子から出力された前記実写画像データの空間的平均値との差分の絶対値が予め設定された閾値を越えた場合に異常な実写画像データを出力する異常画素であると判定された前記放射線検出素子であって、前記異常な実写画像データを出力する確率が予め設定された発生確率の閾値を越えて欠陥画素であると判定された前記放射線検出素子が含まれることを特徴とする請求項４に記載の放射線画像生成システム。
　前記コンソールは、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに基づく画像、或いは、前記間引きデータ、又は前記実写画像データに基づいて作成した間引きデータに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行って、前記間引きデータ又は前記実写画像データに基づいて作成した間引きデータに対する画像処理後の画像、或いは、前記間引きデータの画像処理後の画像に対し、欠陥或いは異常画素の補正処理が施された画像、或いは、前記実写画像データに対してゲイン補正処理及び前記オフセット補正値に基づくオフセット補正処理を行った画像処理後の画像、或いは、前記実写画像データに対する画像処理後の画像に対し、欠陥或いは異常画素の補正処理が施された画像、のうちの少なくとも１つを、前記表示部に表示させることを特徴とする請求項４から請求項１３のいずれか一項に記載の放射線画像生成システム。