JP2017009504A - Ｘ線データ処理装置、その方法およびプログラム - Google Patents

Abstract

【課題】比較器等を用いることなく、ピクセルアレイ型のＸ線検出器で検出され、チャージシェアの影響を受けたＸ線強度を補正することができるＸ線データ処理装置を提供する。【解決手段】フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理装置であって、検出部分ごとの検出値を受信し管理する管理部２１０と、検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を検出部分の有効面積率として算出する有効面積率算出部２３０と、管理されている計数値を算出された有効面積率で補正して真値を推定する補正部２５０と、を備える。【選択図】図４

Description

本発明は、特定線源のＸ線強度について計数値から真値を推定するＸ線データ処理装置、その方法およびプログラムに関する。

フォトンカウンティング方式のピクセルアレイ検出器においては、一般的にモノリシックなセンサを使用するため、センサの内部では、ピクセル間に明確な境界は存在しない。読み出しチップと接続する側に形成される読み出しパッド部分だけが、ピクセル状に形成されているのが一般的である（特許文献１参照）。

そのため、隣接するピクセルの境界部分でキャリアが拡散する場合、それらのピクセル間で１フォトン分の電荷が共有されるチャージシェアという現象が生じる。このチャージシェアが原因となり、入射しているＸ線の強度が一定であっても、各ピクセルに設定されている閾値に依存して計数が本来よりも多く計数される場合、少なく計数される場合が存在する。

従来、チャージシェアの影響を無視して測定を行うか、閾値を適当な値に設定することで影響を最小限にして測定が行われてきたが、このような方法ではチャージシェアの効果を取り除くことはできない。一方、チャージシェアによる影響そのものをキャンセルする回路を読み出しチップのピクセル間に実装する研究もなされている（非特許文献１参照）。

特開２０１２−２４２１１１号公報 特開２０１４−１５９９７３号公報

H.-E. Nilsson, B. Norlin, Frojdh, L. Tlustos, "Charge sharing suppression using pixel-to-pixel communication in photon counting X-ray imaging systems", Nuclear Instruments & Methods in Physics Research, Available online 6 February, A576, 2007, 243-247

しかしながら、チャージシェアによる影響をキャンセルする回路を実装する方法は、ピクセル間で信号の波高をそれぞれ比較する必要がある。そのため、比較器の実装により相応の面積が占められ、比較器間で比較するための閾値にばらつきが生じる。また、ピクセルサイズが小さいときに複数ピクセルに跨るチャージシェアに対応することが困難である。

本発明は、このような事情に鑑みてなされたものであり、比較器等を用いることなく、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器で検出され、チャージシェアの影響を受けたＸ線の計数値を補正することができるＸ線データ処理装置を提供することを目的とする。

（１）上記の目的を達成するため、本発明のＸ線データ処理装置は、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理装置であって、検出部分ごとの計数値を受信し管理する管理部と、前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出する有効面積率算出部と、前記管理されている計数値を前記算出された有効面積率で補正して真値を推定する補正部と、を備えることを特徴としている。これにより、比較器等を用いることなく、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器で検出され、チャージシェアの影響を受けたＸ線の計数値を補正することができる。

（２）また、本発明のＸ線データ処理装置は、前記有効面積率算出部が、各線源および各検出エネルギーの閾値に対して前記有効面積率を算出し、前記補正部は、前記算出された有効面積率を連立方程式の係数として表し、前記連立方程式を用いて、各検出エネルギーの閾値に対する計数値を各線源に対する補正値に１次変換することを特徴としている。これにより、多波長のＸ線が検出部分に入射する場合でも、波長依存性を有するチャージシェアの影響を小さくすることができる。

（３）また、本発明のＸ線データ処理装置は、前記検出部分に関するデータが、前記検出部分に応じて予め記憶された、前記検出部分のサイズおよび前記検出部分内の電荷の拡がりの分布を表すデータであることを特徴としている。これにより、検出部分の特性や形状に応じて、チャージシェアの影響を見積もり、有効面積率を算出することができる。

（４）また、本発明の方法は、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理の方法であって、検出部分ごとの計数値を受信し管理するステップと、前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出するステップと、前記管理されている計数値を前記算出された有効面積率で補正して真値を推定するステップと、を含むことを特徴としている。これにより、チャージシェアの影響を補正することができる。

（５）また、本発明のプログラムは、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理のプログラムであって、検出部分ごとの計数値を受信し管理する処理と、前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出する処理と、前記管理されている計数値を前記算出された有効面積率で補正して真値を推定する処理と、を含む一連の処理をコンピュータに実行させることを特徴としている。これにより、チャージシェアの影響を補正することができる。

本発明によれば、比較器等を用いることなく、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器で検出され、チャージシェアの影響を受けたＸ線の計数値を補正することができる。

本発明のＸ線回折システムの構成を示す概略図である。 主にＸ線検出器およびＸ線データ処理装置の構成を示すブロック図である。 （ａ）、（ｂ）それぞれチャージシェアが生じない場合と生じる場合のＸ線の検出を示す模式図である。 主にＸ線データ処理装置の構成を示すブロック図である。 ＲＯＩＣのゲインのばらつきを示すグラフである。 ピクセルサイズおよび閾値に応じたチャージシェアの影響を示す表である。 実測のＸ線プロファイルおよび各線源の有効面積率曲線を示すグラフである。 （ａ）、（ｂ）補正をしない場合および有効面積率を用いて補正した場合それぞれの画像を表す図である。 有効面積率を用いた補正をしない場合およびした場合の強度分布を表すグラフである。 蛍光Ｘ線低減モードで測定されたＸ線強度および有効面積率を用いた補正で蛍光Ｘ線を除去したＸ線強度を表すグラフである。 閾値を５ｋｅＶおよび７ｋｅＶのそれぞれとして測定された画像（上段）と有効面積率を用いて補正されたＣｕ光源画像およびＦｅ光源画像（下段）とを示す図である。

次に、本発明の実施の形態について、図面を参照しながら説明する。説明の理解を容易にするため、各図面において同一の構成要素に対しては同一の参照番号を付し、重複する説明は省略する。

（全体構成）
図１は、Ｘ線測定システム１０の構成の一例を示す概略図である。図１に示すようにＸ線測定システム１０は、Ｘ線源２０、試料Ｓ、Ｘ線検出器１００およびＸ線データ処理装置２００で構成されている。

Ｘ線源２０は、例えば、陰極であるフィラメントから放射された電子束を対陰極であるロータターゲットに衝突させてＸ線を発生させる。Ｘ線源２０から放射されるＸ線は、断面形状が円形または矩形のドット状の、いわゆるポイントフォーカスのＸ線ビームである。

ロータターゲットの外周面に、原子番号が互いに異なる複数種類の金属（例えば、ＭｏおよびＣｕ）を設けることで、Ｘ線源２０を多波長の線源にすることもできる。Ｃｕターゲットに電子が衝突したとき、特性線であるＣｕＫα線（波長１．５４２Å）を含むＸ線が放射され、Ｍｏターゲットに電子が衝突したとき、特性線であるＭｏＫα線（波長０．７１１Å）を含むＸ線が放射される。ロータターゲットから出射されるＸ線には、互いに違ったターゲット素材の特性Ｘ線であるＣｕＫα線およびＭｏＫα線が混在している。

試料Ｓは、試料支持装置により支持されている。試料支持装置は、試料Ｓの特性および測定の種類に応じて決められる。Ｘ線検出器１００は、試料Ｓで回折された回折Ｘ線および蛍光Ｘ線を検出する。Ｘ線データ処理装置２００は、測定された回折Ｘ線および蛍光Ｘ線のＸ線データを処理し、測定結果を表示させる。Ｘ線検出器１００およびＸ線データ処理装置２００の詳細については、後述する。

（Ｘ線検出器およびＸ線データ処理装置の構成）
図２は、主にＸ線検出器１００およびＸ線データ処理装置２００の構成を示すブロック図である。Ｘ線検出器１００は、複数のＸ線受光用のピクセル１１０（検出部分）を有しており、例えば２次元半導体検出器である。複数のピクセル１１０は、２次元的にアレイ化されており、規則的に配列されている。なお、検出器は、２次元半導体検出器に限られず、１次元半導体検出器であってもよい。

分別回路１２０は、複数のピクセル１１０の個々に接続されており、さらにはカウンタ部１３０が、分別回路１２０の個々に接続されている。カウンタ読み出し回路１５０は、各カウンタ部１３０に接続されている。

分別回路１２０は、ピクセル１１０のパルス信号をＸ線波長ごとに分別して出力する。カウンタ部１３０は、分別回路１２０によって波長毎に分別された信号のそれぞれの個数を計数する。カウンタ部１３０は、例えば、分別回路１２０によって分別された数のパルス信号をそれぞれカウントできるように分別数と同じ数のカウンタ回路を内蔵する。カウンタ読み出し回路１５０の出力信号は、エネルギーの閾値により分離されたＸ線データとしてＸ線データ処理装置２００に通信線を通して伝送される。

Ｘ線データ処理装置２００は、例えば、パーソナルコンピュータである。パーソナルコンピュータは、例えば、演算制御するためのＣＰＵ、データを記憶するためのメモリ、メモリ内の所定領域に記憶されたシステムソフト、およびメモリ内の他の所定領域に記憶されたアプリケーションプログラムソフト、等によって構成されている。

Ｘ線データ処理装置２００には、ユーザの入力を受け付ける入力部３００としてキーボード等が接続されている。ユーザは、入力部３００を介して測定結果の表示や補正の指示等を行うことができる。また、入力部３００は、線源および検出エネルギーの閾値に関するデータの入力を受け付ける。なお、線源には、試料への照射Ｘ線源だけでなく、蛍光Ｘ線源も含まれる。また、Ｘ線データ処理装置２００には、ディスプレイやプリンタ等の出力部４００が接続されている。出力部４００は、Ｘ線データ処理装置２００からの指示に従って測定結果を出力する。

図３（ａ）、（ｂ）は、それぞれチャージシェアが生じない場合と生じる場合のＸ線の検出を示す模式図である。図３（ａ）に示すように、Ｘ線が単一のピクセル１１０にのみ入射する場合には、チャージシェアが生じず、正確な測定が可能である。しかし、図３（ｂ）に示すようにＸ線検出器１００の表面近傍でできた電荷の雲が電極に到達するまでの間に拡がってチャージシェアが起こる。２つのピクセル１１０にわたって生じたチャージシェアにより、一方のピクセル１１０の検出したピークが低くなり、他方のピクセル１１０にも低いピークが検出されることになる。

（Ｘ線データ処理装置の構成）
図４は、主にＸ線データ処理装置２００の構成を示すブロック図である。Ｘ線データ処理装置２００は、管理部２１０、データ記憶部２２０、有効面積率算出部２３０、補正部２５０を備えており、Ｘ線の計数値から真値を推定する。なお、Ｘ線データ処理装置２００は、複数種類の光源の回折Ｘ線データを分離する場合に用いてもよいし、回折Ｘ線データから蛍光Ｘ線を分離する場合に用いてもよい。

管理部２１０は、Ｘ線検出器１００でピクセルごとに検出された計数値を受け取り、管理する。例えば、管理部２１０は、ピクセル１１０の番地（ｉ，ｊ）に関連付けて波長毎の回折Ｘ線強度を決定し、その結果のデータを記憶する。管理部２１０は、ユーザの指示に応じて、例えば記憶されたＣｕ線源の回折Ｘ線像とＭｏ線源の回折Ｘ線像の両方の回折像のデータを、出力部４００に表示させることができる。いずれか一方の回折像を表示させることもできるし、両方の像を同時に表示させることもできる。

データ記憶部２２０は、Ｘ線検出器１００のセンサの材質、構造およびピクセルに関するデータならびにＸ線源２０および検出エネルギーの閾値に関するデータを記憶している。ピクセルに関するデータは、ピクセルに応じて予め記憶された、ピクセルのサイズ、形状およびセンサ内の電荷の拡がりの分布を表すデータである。また、Ｘ線源２０および検出エネルギーの閾値に関するデータは、Ｘ線測定システム１０が用いられる際の条件を表すデータであり、ユーザにより入力される。これにより、ピクセルの特性や形状に応じて、チャージシェアの影響を見積もり、有効面積率を算出することができる。

有効面積率算出部２３０は、ピクセルに関するデータおよび入力されたデータを用いて、ピクセルの本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合をピクセルの有効面積率として算出する。特に、複数の線源および複数の閾値に対して有効面積率算出部２３０は、各線源および各検出エネルギーの閾値に対して有効面積率を算出する。

補正部２５０は、管理されている計数値を算出された有効面積率で補正して真値を推定する。このようにして、フォトンカウンティング方式のピクセルアレイ型のＸ線検出器で検出され、チャージシェアの影響を受けたＸ線の計数値を補正することができる。具体的には、複数の線源および複数の閾値に対して、補正部２５０は、算出された有効面積率を連立方程式の係数として表し、連立方程式を用いて、各検出エネルギーの閾値に対する計数値を各線源に対する補正値に１次変換する。これにより、多波長のＸ線がピクセルに入射する場合でも、波長依存性を有するチャージシェアの影響を小さくすることができる。

（補正および補正に用いられる有効面積率の算出）
チャージシェアによる有効面積率は閾値波長と、入射Ｘ線の波長に対する依存性が存在する。この波長依存性と、複数の閾値によって測定された計数値を用いて、波長毎に、ピクセルに到達したＸ線の計数の推定を行う。

異なる２つの波長をＡ、Ｂとし、異なる２つの適当な閾値をＬ、Ｈとする。また、波長Ａ、Ｂの光源から各ピクセルに到達するＸ線の真の光子数をそれぞれＩ_Ａ、Ｉ_Ｂとし、閾値Ｌ、Ｈで観測される計数をそれぞれＩ_Ｌ、Ｉ_Ｈとする。閾値Ｌと波長Ａによって決まるピクセル有効面積率をｐ_ＬＡのように表せば、以下の数式（１）ように観測される計数と真の光子数の関係が得られる。

ここで有効面積率からなる行列をＰとして、その逆行列を求めれば、以下の数式（２）となることから、波長Ａ、Ｂによる信号を区別せずに観測した計数を元に、線源Ａ、Ｂによる計数を独立に得ることが可能となる。

半導体センサ内での電荷の拡がりが正規分布に従い、その標準偏差がσであると仮定し、ピクセルの一辺の長さをｄとすると、有効面積率ｐは次の数式（３）のように近似できる（∵ｄ≫σ）。
ここで、λは対象としているＸ線の波長、λ_ＴＨは検出回路の閾値波長である。ｄ≫σが成立しない場合にはピクセルの４つの角による効果を考慮する必要がある。

上記のＸ線測定システム１０は、多色光源の分離に用いることができる。例えば、高エネルギー側の閾値と低エネルギー側の閾値によりそれぞれ分離されたＸ線の計数値Ｉ_ＨおよびＩ_Ｌから、以下の数式（４）によりＣｕ光源およびＭｏ光源の回折Ｘ線強度Ｉ_ＣｕおよびＩ_Ｍｏを求めることができる。

また、上記のＸ線測定システム１０は、温度補正後の一様性補正テーブルを仮想一様照射無しで再作成することにも応用できる。

（実施例１）
上記のＸ線測定システム１０を用いて測定を行った。Ｘ線検出器２０に用いた標準読み出し集積回路（ＲＯＩＣ）は、ゲインのばらつきを有する。例えば、６ｋｅＶに閾値を設定しても実際に検出するエネルギーは、ピクセルごとにばらつく。図５は、ＲＯＩＣのゲインのばらつきを示すグラフである。図５に示す例では、ゲインは、平均値μ＝６．０ｋｅＶ、標準偏差σ＝０．２ｋｅＶの分布を有している。

Ｘ線の検出には、一辺１００μｍの矩形で、センサ厚が３２０μｍ、キャリアの広がりがσ＝５．０μｍのガウス分布のピクセルを表面に有するＸ線検出器を用いた。ピクセルサイズの大小に応じてチャージシェアの影響の大小が異なり、有効面積率も異なる。図６は、ピクセルサイズおよび閾値に応じたチャージシェアの影響を示す表である。ＲＯＩＣの性能由来で設定値に対して１ｋｅＶ程度のばらつきが生じると、Ｘ線の計数値が８％程度ばらつく。

上記のようなＸ線検出器を用いて閾値に対するＣｕＫのＸ線強度を測定するとともに、閾値に対するＭｏ、Ｃｕ、Ｆｅ、Ｃｒ光源のそれぞれのＸ線検出の有効面積率を算出した。図７は、実測のＸ線プロファイルおよび各線源の有効面積率曲線を示すグラフである。算出において、有効面積率曲線は、ノイズピークの中心を０ｋｅＶとして、閾値を計算した。また、８．０４ｋｅＶの単色光について計算した有効面積−閾値曲線にＦＷＨＭ２４％（一例としてσ＝８．０４ｋｅＶで８３０ｅＶ）のガウス分布を畳み込んだ。その結果、実測のＣｕＫのＸ線プロファイルとＣｕ光源のＸ線検出の有効面積率の曲線とがほぼ一致した。

（実施例２）
上記のＸ線測定システムを用いて得られた実験データを、有効面積率を用いて補正する例を説明する。Ｃｕ光源を用いて、以下の（５）の通りに実験結果が得られたものとする。

このような実験結果に対して、単純な有効面積率による補正のみを行うと以下の（６）のようにＸ線の計数値が補正される。

上記の実験結果が、鉄製の試料にＣｕ光源のＸ線を照射して得られたものであれば、蛍光Ｘ線による低減分を反映して、以下の（７）のようにＣｕ光源由来のＸ線強度Ｉ_Ｃｕおよび蛍光Ｘ線のＸ線強度Ｉ_Ｆｅを計算できる。

また、入射するＸ線が単色であることが明らかであれば、バックグラウンドと有効面積１００％での強度が得られる。

（実施例３）
一辺１００μｍの矩形のピクセルのＸ線検出器を備えたＸ線測定システムにより、実際に検出されたＸ線強度の有効面積率による補正を行った。センサに並べて接続されたＲＯＩＣ間の計数値の補正には単純平均した値を採用した。図８は、（ａ）補正をしない場合および（ｂ）有効面積率を用いて補正した場合、それぞれの画像を表す図である。補正しない場合は、全体的に計数値が小さく、一方、補正した場合は、全体的に計数値が大きくなっている。

図９は、有効面積率を用いた補正をしない場合およびした場合の強度分布を表すグラフである。補正しない場合の曲線は、全体的に計数値が小さく、ピクセル間の計数値のばらつきが大きくなっている。一方、補正した場合の曲線は、全体的に計数値が大きく、ピクセル間の計数値のばらつきが小さくなっている。

（実施例４）
上記の実施例と同様の条件で、Ｃｕ光源を用い鉄製の粉末試料にＸ線を照射したときの回折パターンを撮影した。その際に、まず、Ｆｅの蛍光Ｘ線を低減するために閾値を７ｋｅＶに設定して撮影を行った（蛍光低減モード）。次に、閾値を５ｋｅＶおよび７ｋｅＶのそれぞれとしてＸ線画像を撮影した。そして、撮影された画像に対して各閾値とＣｕ光源（８．０４ｋｅＶ）およびＦｅ光源（６．４ｋｅＶ）のそれぞれに対して決まる有効面積率を係数として用いた連立方程式を解くことで、画像を補正し、各光源由来の画像に分離した（蛍光分離モード）。

図１０は、蛍光Ｘ線低減モードで測定されたＸ線強度および有効面積率を用いた補正で蛍光Ｘ線を分離したＸ線強度を表すグラフである。図１０は、蛍光低減モードの測定データ画像と蛍光除去画像をｘ＝３６０〜４１１の範囲で、平均計数プロファイルとしてプロットしたものである。蛍光分離モードのグラフは、蛍光低減モードのグラフに対し、ＳＢＲ（Signal to Background Ratio）で約３倍の向上が見られた。

このようにして、分離されたＣｕ光源(８．０４ｋｅＶ)由来の画像は、閾値を４．１ｋｅＶとして撮影した場合に得られるＣｕ光源由来の画像に相当する画像として得られた。また、Ｆｅ光源(６．４ｋｅＶ)由来の画像は、閾値３．３ｋｅＶとして撮影した場合に得られるＦｅ光源由来の画像に相当する画像として得られた。

図１１は、閾値を５ｋｅＶおよび７ｋｅＶのそれぞれとして測定された画像（上段）と有効面積率を用いて補正されたＣｕ光源画像およびＦｅ光源画像（下段）とを示す図である。閾値７ｋｅＶの画像（上段左）に比べ、有効面積率を用いて補正されたＣｕ光源画像（下段左）では、回折線が明確になり、蛍光Ｘ線をバックグラウンドとして分離、除去できているのが分かる。

１０ Ｘ線測定システム
２０ Ｘ線源
１００ Ｘ線検出器
１１０ ピクセル（検出部分）
１２０ 分別回路
１３０ カウンタ部
１５０ カウンタ読み出し回路
２００ Ｘ線データ処理装置
２１０ 管理部
２２０ データ記憶部
２３０ 有効面積率算出部
２５０ 補正部
３００ 入力部
４００ 出力部
Ｓ 試料

Claims (5)

1. フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理装置であって、
   検出部分ごとの計数値を受信し管理する管理部と、
   前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出する有効面積率算出部と、
   前記管理されている計数値を前記算出された有効面積率で補正して真値を推定する補正部と、を備えることを特徴とするＸ線データ処理装置。
2. 前記有効面積率算出部は、各線源および各検出エネルギーの閾値に対して前記有効面積率を算出し、
   前記補正部は、前記算出された有効面積率を連立方程式の係数として表し、前記連立方程式を用いて、各検出エネルギーの閾値に対する計数値を各線源に対する補正値に１次変換することを特徴とする請求項１記載のＸ線データ処理装置。
3. 前記検出部分に関するデータは、前記検出部分に応じて予め記憶された、前記検出部分のサイズおよび前記検出部分内の電荷の拡がりの分布を表すデータであることを特徴とする請求項１または請求項２記載のＸ線データ処理装置。
4. フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理の方法であって、
   検出部分ごとの計数値を受信し管理するステップと、
   前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出するステップと、
   前記管理されている計数値を前記算出された有効面積率で補正して真値を推定するステップと、を含むことを特徴とする方法。
5. フォトンカウンティング方式のピクセルアレイ型のＸ線検出器により検出されたＸ線の計数値から真値を推定するＸ線データ処理のプログラムであって、
   検出部分ごとの計数値を受信し管理する処理と、
   前記検出部分に関するデータおよび線源および検出エネルギーの閾値に関するデータを用いて、前記検出部分の本来の検出能力に対するチャージシェアの影響を受けた検出能力の割合を前記検出部分の有効面積率として算出する処理と、
   前記管理されている計数値を前記算出された有効面積率で補正して真値を推定する処理と、を含む一連の処理をコンピュータに実行させることを特徴とするプログラム。