JP4993144B2 - Ｐｅｔ検出器の光量出力を補正する方法 - Google Patents

Description

本発明は、医療用放射線撮像装置に取り付けて使用され、ＲＩ線源を含む薬剤が投与された被検体から発生する放射線を検出するＤＯＩ（Depth Of Interaction）検出器等のＰＥＴ検出器の光量出力を補正する方法に関する。

医療用放射線撮像装置の１つとして、ＰＥＴ（Positron Emission Tomography）装置が知られている。ＰＥＴ装置は、ＲＩ線源（陽電子放出核種）を含む薬剤が投与された被検体から発生する放射線（γ線）の対を同時計数することにより、ＲＩ線源の体内分布を画像化し、生体機能を見るための装置であり、生体や疾病患の臨床検査等に応用されている。

ＰＥＴ装置には放射線検出器（以下、「ＰＥＴ検出器」とする。）が備えられている。ＰＥＴ検出器は、通常、シンチレータ結晶の２次元的配列からなるシンチレータアレイと、このシンチレータアレイに光結合された少なくとも１本の光電子増倍管とから構成されている。そして、多数のＰＥＴ検出器が、それぞれのシンチレータアレイの検出面を被検体に向けた状態で、被検体の周囲を取り囲む円周上に配置されて１つの検出器リングが構成され、さらに、複数の検出器リングがその軸方向に重ねられて、ＰＥＴ装置のガントリに搭載され、ＰＥＴ装置の放射線検出部が構成される。

こうして、被検体にＲＩ線源を含む薬剤が投与されると、被検体から放射線（γ線）が反対方向に同時に放出されるイベントが発生し、この放射線（γ線）の対が放射線検出器のシンチレータ結晶に光電吸収され、それによって放射線（γ線）のエネルギーが光に変換され、さらに、この光は光電子増倍管によってパルス状の電気信号に変換され、この電気信号に基づいて放射線（γ線）の入射位置およびエネルギーが検出される。

この場合、放射線（γ線）の入射位置は、放射線（γ線）との相互作用により発光したシンチレータ結晶を識別することによって検出されるので、放射線（γ線）が放射線検出器の検出面に斜めに入射してシンチレータ結晶の深部で発光が生じた場合には、実際の放射線（γ線）の入射位置と放射線（γ線）の検出位置との間にずれが生じ、それによって放射線検出器の空間分解能が低下する。

そこで、微小なシンチレータ結晶の２次元的配列からなるシンチレータアレイを複数段積層してシンチレータブロックとし、それに光電子増倍管を光結合したＤＯＩ検出器が提案されている。ＤＯＩ検出器によれば、シンチレータ結晶の深さ方向の放射線検出情報が得られるので、放射線がＤＯＩ検出器の検出面に斜めに入射した場合でも、空間分解能の低下が抑制出来る（例えば、特許文献１、２参照）。

ところで、これらのＰＥＴ検出器においては、同じエネルギーを持つγ線に対してシンチレータ結晶毎に光量出力にばらつきがある。
そこで、各シンチレータ結晶の光量出力をエネルギーの次元に変換する。この補正係数はエネルギー補正係数と呼ばれ、放射線入射位置毎に算出されたエネルギー補正係数がエネルギー補正ルックアップテーブルに保持される。

この光量スペクトルは、光量出力毎のカウント数で表されるヒストグラムであり、従来のエネルギー出力補正法においては、エネルギー補正係数の算出精度は、ヒストグラムのビン数（チャンネル数）によって決定されていた。
したがって、従来法による補正精度はあまり高くなく、よって、従来法を用いてＤＯＩ検出器のエネルギー出力の補正を行うと、高いエネルギー分解能が得られないという問題を生じていた。また、ＤＯＩ検出器の場合は更にシンチレータ結晶の個数が多くなっているので、ＤＯＩ検出器以外の通常のＰＥＴ検出器と比べてルックアップテーブルを作る手間が増大していた。

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

したがって、本発明の課題は、ＰＥＴ検出器のエネルギー補正係数を高精度でまた高速に算出することができる方法を提供することにある。

上記課題を解決するため、第１発明は、シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段または複数段有するシンチレータブロックと、前記シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成されたＰＥＴ検出器の光量出力を補正する方法であって、
（１）基準線源から前記ＰＥＴ検出器の検出面に放射線を照射し、前記光電子増倍管からの光量に関する出力信号を収集するステップと、
（２）収集した前記出力信号に基づき放射線の入射位置を求めるための２次元位置マップを作成し、放射線と相互作用したシンチレータ結晶を識別するステップと、
（３）前記放射線と相互作用したシンチレータ結晶毎の光量に関する情報を蓄積した光量スペクトルを取得するステップと、
（４）前記光量スペクトルのピークチャンネルを小数点以下まで求め、その光量ピークが前記基準線源からの放射線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するためのエネルギー補正係数を算出し、前記エネルギー補正係数に基づいてエネルギー補正ルックアップテーブルを作成するステップと、を有することを特徴とする方法を構成したものである。

上記課題を解決するため、また、第２発明は、シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段または複数段有するシンチレータブロックと、前記シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成されたＰＥＴ検出器の光量出力を補正する方法であって、
（１）基準線源から前記ＰＥＴ検出器の検出面に放射線を照射し、前記光電子増倍管からの光量に関する出力信号を収集するステップと、
（２）収集した前記出力信号に基づき前記放射線の入射位置を求めるための２次元位置マップを作成し、放射線と相互作用したシンチレータ結晶を識別するステップと，
（３）前記放射線と相互作用したシンチレータ結晶毎の光量に関する情報を蓄積した光量スペクトルを取得するステップと、
（４）前記光量スペクトルのピークチャンネルを小数点以下まで求め、その光量ピークが前記基準線源からの放射線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するための第１のエネルギー補正係数を算出し、前記第１のエネルギー補正係数に基づいて第１のエネルギー補正ルックアップテーブルを作成するステップと、
（５）前記ステップ（１）〜（２）を順次実行した後、収集した前記ＰＥＴ検出器の前記出力信号に前記第１のエネルギー補正ルックアップテーブルを適用して、前記放射線と相互作用したシンチレータ結晶毎のエネルギースペクトルを取得するステップと、
（６）前記ステップ（５）で取得した、前記放射線と相互作用したシンチレータ結晶毎のエネルギースペクトルのエネルギーピークを小数点以下まで求め、前記エネルギーピークを前記基準線源からの放射線に固有のエネルギーピークに揃えるための第２のエネルギー補正係数を算出し、前記第２のエネルギー補正係数に基づいて第２のエネルギー補正ルックアップテーブルを作成するステップと、
（７）前記第２のエネルギー補正ルックアップテーブルから前記第１のエネルギー補正ルックアップテーブルに対する補正係数を求め、求めた補正係数を前記第１のエネルギー補正ルックアップテーブルに掛けることによって、前記第１のエネルギー補正ルックアップテーブルを補正するステップと、を有することを特徴とする方法を構成したものである。

第１発明によれば、エネルギー補正ルックアップテーブルを作成するにあたり、光量スペクトルのピークチャンネルを小数点以下まで求め、それに基づいてエネルギー補正係数を算出するようにしたので、エネルギー補正ルックアップテーブルを高精度で算出することができる。

また、第２発明によれば、先に取得された位置ルックアップテーブルとエネルギー補正ルックアップテーブルに基づいて収集したシンチレータ結晶毎の光量スペクトルから第２のエネルギー補正ルックアップテーブルを作成し、第２のエネルギー補正ルックアップテーブルを用いて元のエネルギー補正ルックアップテーブルを補正することによって、元のエネルギー補正ルックアップテーブルの精度を上げることができる。したがって、先に取得されたエネルギー補正ルックアップテーブルを適用して収集されたエネルギースペクトルを分析し、補正しきれなかった部分を補正するためのエネルギー補正ルックアップテーブルを算出し、それを用いて先のエネルギー補正ルックアップテーブルを再補正するようにすれば、より精度の高いエネルギー補正ルックアップテーブルを作成することができ、また、位置ルックアップテーブルを再作成せずに済むので、より高速にエネルギー補正ルックアップテーブルを作成することができる。

本発明の１実施例によるＤＯＩ検出器の光量出力を補正する方法のフロー図である。 本発明の方法が適用されるＤＯＩ検出器の構成を示す斜視図である。 シンチレータアレイの層毎の２次元位置マップの１例を示す平面図である。 本発明の別の実施例によるＤＯＩ検出器の光量出力を補正する方法のフロー図である。 ＧＳＯ層を２層備えた２個のＤＯＩ検出器の各層における１つのγ線入射位置のエネルギースペクトルを示すグラフである。 ＬＳＯ層およびＧＳＯ層の２層を備えた２個のＤＯＩ検出器の各層における１つのγ線入射位置のエネルギースペクトルを示すグラフである。

以下、添付図面を参照して本発明の好ましい実施例について説明する。図１は、本発明の１実施例によるＤＯＩ検出器の光量出力を補正する方法のフロー図である。図１を参照して、本発明によれば、まず、基準線源からＤＯＩ検出器の検出面に放射線が照射され、光電子増倍管からの光量に関する出力信号が収集される（ステップ１）。この実施例では、基準線源は、５１１ｋｅＶのγ線を放出する核種からなり、ＤＯＩ検出器は、図２に示すように、シンチレータ結晶の２次元的配列（検出素子マトリクスが９×１０）からなるシンチレータアレイ１を２段に積層したシンチレータブロック２と、シンチレータブロック２に光結合された４本の光電子増倍管３とから構成される。

次に、収集された各光電子増倍管３の出力信号を用いて重心演算が行われ、得られた重心位置が、シンチレータアレイの層毎に、図３に示すような９×１０領域に区分けした２次元位置マップ上に投影される。そして、２次元位置マップ上における各シンチレータ結晶に対応する位置ルックアップテーブルが作成され、それによって、γ線と相互作用したシンチレータ結晶が特定される（ステップ２）。

また、各シンチレータ結晶の光量出力を蓄積することで光量スペクトルが取得される（ステップ３）。そして、光量スペクトルのピーク位置が小数点以下まで求められ、そのピーク位置が、照射されたγ線のエネルギーに対応するとして、エネルギー補正係数が算出され、算出されたエネルギー補正係数を用いてエネルギー補正ルックアップテーブルが作成される（ステップ４）。

このステップ４について、具体的に説明する。
光量に関する出力が７ビットとすると、光量スペクトルはチャンネル毎に計数値がプロットされたものであるから、光量スペクトルのスケールは０〜１２７チャンネルとなる。ところが、従来法のように、光量スペクトルのピーク位置をチャンネルのスケールで算出すると、エネルギー補正係数は最大でチャンネルの数に対応する１２８通り（７ビット）の数値しかとることができない。その結果、エネルギー補正係数の算出精度、よってエネルギー出力の補正精度に限界が生じる。
したがって、本発明の方法では、光量スペクトルのピーク位置を７ビット（０〜１２７）のスケールではなく、小数点以下まで算出し、エネルギー補正係数を、例えば１０ビットの数値範囲で算出することによって、光量出力の補正精度を上げた。

光量スペクトルのピーク位置の算出は、次のようにして行う。まず、ピーク位置を算出する前に、光量スペクトルを滑らかな曲線によって近似する。この実施例では、光量スペクトルを、そのピークを中心とする３点を通る２次関数によって近似する。この場合、光量スペクトルをガウス分布によって近似してもよい。そして、この近似曲線に基づき、チャンネル数のスケール（７ビット）よりも細かいスケール（１０ビット）でピーク位置を算出する。それによって、従来は７ビットの整数値（チャンネル数）で求めていたピーク位置を、小数点以下まで求めることができる。
そして、エネルギー補正係数ＥＣを
ＥＣ＝６４［ｃｈ］／ピーク位置［ｃｈ］
×（５１２［ｋｅＶ］／γ線のエネルギー［ｋｅＶ］） （１）
によって求める。

こうして、シンチレータ結晶毎のエネルギー補正係数ＥＣを算出し、算出したエネルギー補正係数ＥＣに基づいてエネルギー補正ルックアップテーブルを作成する。

以上のように、本発明によれば、エネルギー補正ルックアップテーブルを作成するにあたり、光量スペクトルのピークチャンネルを小数点以下まで求め、光量出力が基準線源からのγ線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するためのエネルギー補正係数を算出し、エネルギー補正ルックアップテーブルを作成するので、補正精度が上がる。
また、上述の実施例では、本発明の方法をＤＯＩ検出器に適用したが、ＤＯＩ検出器以外のＰＥＴ検出器（シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段有するシンチレータブロックと、シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成される）にも適用することができ、それによって、ＰＥＴ検出器の光量出力の補正精度を上げることができる。

図４は、本発明の別の実施例による、ＰＥＴ検出器の光量出力をさらに高精度でかつ高速に補正する方法のフロー図である。図４を参照して、本発明によれば、まず、基準線源からＤＯＩ検出器の検出面に放射線が照射され、光電子増倍管からの光量に関する出力信号が収集される（ステップ１）。この実施例の場合も、図１の実施例と同様、基準線源は、５１１ｋｅＶのγ線を放出する核種からなり、ＤＯＩ検出器は、図２に示したものと同じ構成を有している。

次に、収集された出力信号を用いて重心演算が行われ、得られた重心位置が、シンチレータアレイの層毎に、２次元位置マップ上に投影される。そして、２次元位置マップ上における各シンチレータ結晶に対応する位置ルックアップテーブルが作成され、それによって、γ線と相互作用したシンチレータ結晶が識別される（ステップ２）。

また、シンチレータ結晶のそれぞれの光量出力を蓄積することで光量スペクトルが取得される（ステップ３）。そして、光量スペクトルのピークチャンネルが小数点以下まで求められ、その光量ピークが基準線源からのγ線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するための第１のエネルギー補正係数が算出され、それに基づいて第１のエネルギー補正ルックアップテーブルが作成される（ステップ４）。

さらに、ステップ（１）〜（２）が再び順次実行され、収集されたＤＯＩ検出器の出力信号に第１のエネルギー補正ルックアップテーブルが適用され、γ線と相互作用したシンチレータ結晶毎のエネルギースペクトルが取得される（ステップ５）。

そして、ステップ５で取得された、γ線と相互作用したシンチレータ結晶毎のエネルギースペクトルのエネルギーピークが小数点以下まで求められ、そのエネルギーピークを基準線源からのγ線に固有のエネルギーピークに揃えるための第２のエネルギー補正係数ＥＣ”が算出され、それに基づいて第２のエネルギー補正ルックアップテーブルＥＣＬＵＴ”が作成される（ステップ６）。

このステップ６について具体的に説明する。
エネルギースペクトルのピークは、第１のエネルギー補正係数の精度が低かったり、経時変化によってシンチレータの光量出力に変化があったりした場合は５１１ｋｅＶにならない。そこで、例えば、ピークを中心とする３点を通る２次関数を求め、この２次関数の極大値を求めることによってピーク位置が算出される。
そして、第２のエネルギー補正係数ＥＣ”が
ＥＣ”＝５１1［ｋｅＶ］／ピークエネルギー［ｋｅＶ］ （３）
で求められ、それに基づいて第２のエネルギー補正ルックアップテーブルＥＣＬＵＴ”が作成される。

第２のエネルギー補正ルックアップテーブルＥＣＬＵＴ”は、第１のエネルギー補正ルックアップテーブルＥＣＬＵＴ’に対する補正係数とすることができるので、
ＥＣＬＵＴ＝ＥＣＬＵＴ’×ＥＣＬＵＴ” （４）
とすることで、第１のエネルギー補正ルックアップテーブルＥＣＬＵＴ’が補正され、より精度の高いエネルギー補正ルックアップテーブルＥＣＬＵＴが作成される（ステップ７）。

本発明によれば、先に取得された位置ルックアップテーブルおよびエネルギー補正ルックアップテーブルに基づいて収集したシンチレータ結晶毎の光量スペクトルから第２のエネルギー補正ルックアップテーブルを作成し、第２のエネルギー補正ルックアップテーブルを用いて元のエネルギー補正ルックアップテーブルを補正することによって、元のエネルギー補正ルックアップテーブルの精度を上げることができる。
こうして、先に取得されたエネルギー補正ルックアップテーブルを適用して収集されたエネルギースペクトルを分析し、補正しきれなかった部分を補正するためのエネルギー補正ルックアップテーブルを算出し、それを用いて先のエネルギー補正ルックアップテーブルを再補正するようにすれば、より精度の高いエネルギー補正ルックアップテーブルを作成することができ、また、位置ルックアップテーブルを再作成せずに済むので、より高速にエネルギー補正ルックアップテーブルを作成することができる。

なお、上述の実施例では、本発明の方法をＤＯＩ検出器に適用したが、ＤＯＩ検出器以外のＰＥＴ検出器（シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段有するシンチレータブロックと、シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成される）にも適用することができる。

次に、ＥＣＬＵＴによれば、ＥＣＬＵＴ’に比べてどの程度補正の精度が上がるかを調べるべく実証実験を行った。この実験では、それぞれ発光減衰時間が異なるケイ酸ガドリニウム（ＧＳＯ）結晶から形成された２層のシンチレータアレイ（シンチレータアレイの層は、いずれも９×１０個のシンチレータ結晶から構成される）を備えたＤＯＩ検出器を２個（以下、「Ｇ−Ｇ検出器１」および「Ｇ−Ｇ検出器２」と呼ぶ。）と、ケイ酸ガドリニウム（ＧＳＯ）結晶から形成されたシンチレータアレイの層、およびケイ酸ルテチウム（ＬＳＯ）結晶から形成されたシンチレータアレイの層の２層（シンチレータアレイの層は、いずれも９×１０個のシンチレータ結晶から構成される）を備えたＤＯＩ検出器を２個（以下、「Ｌ−Ｇ検出器１」および「Ｌ−Ｇ検出器２」と呼ぶ。）作製した。

そして、以下の条件下において、４個のＤＯＩ検出器のそれぞれについて平均エネルギースペクトルを収集し（収集時間は３００秒）、解析した。
条件（１）：ＥＣＬＵＴ’を適用してエネルギースペクトルを収集した。
条件（２）：ＥＣＬＵＴを適用してエネルギースペクトルを収集した。
Ｇ−Ｇ検出器１およびＧ−Ｇ検出器２のそれぞれを構成する１つのシンチレータ結晶（γ線入射位置）について収集したエネルギースペクトルを図５のグラフに示し、Ｌ−Ｇ検出器１およびＬ−Ｇ検出器２のそれぞれの１つのシンチレータ結晶（γ線入射位置）について収集したエネルギースペクトルを図６のグラフに示す。図５および図６中、各検出器の欄において、上段のグラフは、当該検出器の２つのシンチレータアレイの層のうち、下側のシンチレータアレイの層（以下、「ＤＯＩ０」と呼ぶ。）のエネルギースペクトルを示し、下段のグラフは、上側のシンチレータアレイの層（以下、「ＤＯＩ１」と呼ぶ。）のエネルギースペクトルを示す。また、各エネルギースペクトルのチャンネル数は７ビットであり、そのピーク位置は１０ビットのスケールで求めてある。

さらに、収集したエネルギースペクトルを解析した。解析結果を下表に示す。

この表において、エネルギースペクトルのピーク位置が、６３．８７５チャンネルに近づくほど補正の精度は高いものと判定される。また、エネルギー分解能はエネルギースペクトルの半値幅で定義され、よって、エネルギー分解能が小さいほど、補正の精度は高いものと判定される。
また、表中、平均値は、それぞれ、シンチレータアレイの層毎に、それを構成する全シンチレータ結晶（９×１０個）について算出した平均値であり、ＥＷはピーク位置の補正後に設定したエネルギーウィンドウを表し、ΔＥＷを次式に従って算出した。
ΔＥＷ＝Δ｜ピーク位置−６３．８７５［ｃｈ］｜×８［ｋｅＶ／ｃｈ］
＋Δエネルギー分解能［％］×５１１［ｋｅＶ］ （５）
ここで、Δは、条件（１）の下で得られたデータと条件（２）の下で得られたデータとの差を表す。

この表からわかるように、ＥＣＬＵＴを適用することによって、Ｇ−Ｇ検出器では、ピーク位置で平均２．７２ｋｅＶ（３．８５−１．１３＝２．７２）改善され、エネルギー分解能で平均２．７５ｋｅＶ（９７．０９−９４．３４＝２．７５）改善され、合計で平均５．４７ｋｅＶの改善が実現し、また、Ｌ−Ｇ検出器では、ピーク位置が平均６．２０ｋｅＶ（７．５４−１．３４＝６．２０）改善され、エネルギー分解能が平均６．６６ｋｅＶ（９１．９８−８４．３２＝７．６６）改善され、合計で平均１３．８６ｋｅＶの改善が実現した。

本発明によれば、先に取得された位置ルックアップテーブルと第１のエネルギー補正ルックアップテーブルＥＣＬＵＴ’に基づいて収集したエネルギースペクトルから第２のエネルギー補正ルックアップテーブルＥＣＬＵＴ”を作成し、第２のエネルギー補正ルックアップテーブルＥＣＬＵＴ”を用いて元のエネルギー補正ルックアップテーブルＥＣＬＵＴ’を補正することによって、元のエネルギー補正ルックアップテーブルＥＣＬＵＴ’を補正し，更にその精度を上げることができる。つまり、先に取得された位置ルックアップテーブルとエネルギー補正ルックアップテーブルを適用して収集されたエネルギースペクトルを分析し、補正しきれなかった部分を補正するためのエネルギー補正ルックアップテーブルを算出し、それを用いて先のエネルギー補正ルックアップテーブルを再補正するようにすれば、より高速にかつ精度の高いエネルギー補正ルックアップテーブルを簡単に作成することができる。

１ シンチレータアレイ
２ シンチレータブロック
３ 光電子増倍管

Claims (2)

1. シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段または複数段有するシンチレータブロックと、前記シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成されたＰＥＴ検出器の光量出力を補正する方法であって、
   （１）基準線源から前記ＰＥＴ検出器の検出面に放射線を照射し、前記光電子増倍管からの光量に関する出力信号を収集するステップと、
   （２）収集した前記出力信号に基づき放射線の入射位置を求めるための２次元位置マップを作成し、放射線と相互作用したシンチレータ結晶を識別するステップと、
   （３）前記放射線と相互作用したシンチレータ結晶毎の光量に関する情報を蓄積した光量スペクトルを取得するステップと、
   （４）前記光量スペクトルのピークチャンネルを小数点以下まで求め、その光量ピークが前記基準線源からの放射線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するためのエネルギー補正係数を算出し、前記エネルギー補正係数に基づいてエネルギー補正ルックアップテーブルを作成するステップと、を有することを特徴とする方法。
2. シンチレータ結晶の２次元的配列からなるシンチレータアレイを１段または複数段有するシンチレータブロックと、前記シンチレータブロックに光結合された少なくとも１本の光電子増倍管とから構成されたＰＥＴ検出器の光量出力を補正する方法であって、
   （１）基準線源から前記ＰＥＴ検出器の検出面に放射線を照射し、前記光電子増倍管からの光量に関する出力信号を収集するステップと、
   （２）収集した前記出力信号に基づき前記放射線の入射位置を求めるための２次元位置マップを作成し、放射線と相互作用したシンチレータ結晶を識別するステップと，
   （３）前記放射線と相互作用したシンチレータ結晶毎の光量に関する情報を蓄積した光量スペクトルを取得するステップと、
   （４）前記光量スペクトルのピークチャンネルを小数点以下まで求め、その光量ピークが前記基準線源からの放射線に固有のエネルギーに対応するとして、光量出力をエネルギー出力に変換するための第１のエネルギー補正係数を算出し、前記第１のエネルギー補正係数に基づいて第１のエネルギー補正ルックアップテーブルを作成するステップと、
   （５）前記ステップ（１）〜（２）を順次実行した後、収集した前記ＰＥＴ検出器の前記出力信号に前記第１のエネルギー補正ルックアップテーブルを適用して、前記放射線と相互作用したシンチレータ結晶毎のエネルギースペクトルを取得するステップと、
   （６）前記ステップ（５）で取得した、前記放射線と相互作用したシンチレータ結晶毎のエネルギースペクトルのエネルギーピークを小数点以下まで求め、前記エネルギーピークを前記基準線源からの放射線に固有のエネルギーピークに揃えるための第２のエネルギー補正係数を算出し、前記第２のエネルギー補正係数に基づいて第２のエネルギー補正ルックアップテーブルを作成するステップと、
   （７）前記第２のエネルギー補正ルックアップテーブルから前記第１のエネルギー補正ルックアップテーブルに対する補正係数を求め、求めた補正係数を前記第１のエネルギー補正ルックアップテーブルに掛けることによって、前記第１のエネルギー補正ルックアップテーブルを補正するステップと、を有することを特徴とする方法。

Images