(57)【要約】  画像再構成方法１００を実行するために、非一時的コンピュータ読取可能な媒体は、少なくとも１つの電子プロセッサ２０を含むワークステーション１８により読取可能且つ実行可能な命令を記憶する。該方法は、複数の放射線検出器１７により検出された撮像データのフレームにおいて該複数の放射線検出器のシングルス率を決定するステップ；各放射線検出器に対するエネルギ補正係数Ｎwgtを前記撮像データのフレームの取得の間に前記放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定するステップ；各放射線検出器に対するシングルス有効計数時間補正係数を前記シングルス率及び該放射線検出器に対して決定された前記エネルギ補正係数から決定するステップ；システム同時計数有効計数時間補正係数をシステムシングルス率から決定するステップ；放射線検出器の対を接続する複数の同時計数線（ＬＯＲ）の各ＬＯＲに対して、該ＬＯＲに対する有効計数時間補正係数を該ＬＯＲにより接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数及び前記決定されたシステム同時計数有効計数時間補正係数から決定するステップ；及び前記撮像データのフレームを前記決定されたＬＯＲに対する有効計数時間補正係数を用いて再構成するステップ；を含む。
2021511523.tif 000002
【特許請求の範囲】
【請求項１】
画像再構成方法を実行するために少なくとも１つの電子プロセッサを含むワークステーションにより読取可能且つ実行可能な命令を記憶した非一時的コンピュータ読取可能な媒体であって、前記画像再構成方法が、
複数の放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率を決定するステップと、
各放射線検出器に対するエネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定するステップと、
各放射線検出器に対するシングルス有効計数時間補正係数を、当該放射線検出器に対して決定された前記エネルギ補正係数及び前記シングルス率から決定するステップと、
放射線検出器の対を接続する複数の同時計数線の各同時計数線に対して、当該同時計数線に対する有効計数時間補正係数を、当該同時計数線により接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定するステップと、
前記撮像データのフレームを前記決定された同時計数線に対する有効計数時間補正係数を用いて再構成するステップと、
を有する、非一時的コンピュータ読取可能な媒体。
【請求項２】
前記エネルギ補正係数を決定するステップは、
エネルギＥのガンマ線により生成されるトリガの統計的計数である関数Ｎ(Ｅ)を決定するステップと、
各放射線検出器に対する前記エネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布にわたってＮ(Ｅ)を平均することにより決定するステップと、
を有する、請求項１に記載の非一時的コンピュータ読取可能な媒体。
【請求項３】
前記撮像データのフレームの収集の間において各放射線検出器に入射するガンマ線のエネルギスペクトル分布を、前記撮像データのフレームから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定するステップを更に含む、請求項２に記載の非一時的コンピュータ読取可能な媒体。
【請求項４】
前記複数の放射線検出器の各放射線検出器に対して前記エネルギ補正係数を決定するために同一の関数Ｎ(Ｅ)が使用される、請求項２又は請求項３に記載の非一時的コンピュータ読取可能な媒体。
【請求項５】
前記関数Ｎ(Ｅ)を決定するステップが、各エネルギＥにおけるＮ(Ｅ)の平均値を得るためにシミュレーションされた線源からの放射線検出器へのエネルギ蓄積のモンテカルロシミュレーションを実行するステップを有する、請求項２から４の何れか一項に記載の非一時的コンピュータ読取可能な媒体。
【請求項６】
低エネルギトリガ補正係数を、前記撮像データのフレームの収集の間において前記放射線検出器に入射するガンマ線のエネルギスペクトル分布の（ｉ）収集エネルギ窓に入る部分と、（ii）前記収集エネルギ窓及び下側トリガエネルギ窓を含む拡張エネルギ窓に入る部分との比に基づいて決定するステップ、
を更に有し、
前記各放射線検出器に対するシングルス有効計数時間補正係数が、更に、当該放射線検出器に対して決定された前記低エネルギトリガ補正係数から決定される、
請求項１から４の何れか一項に記載の非一時的コンピュータ読取可能な媒体。
【請求項７】
前記撮像データのフレームの収集の間において各放射線検出器に入射するガンマ線の前記エネルギスペクトル分布を、前記撮像データのフレームから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定するステップを更に含む、請求項６に記載の非一時的コンピュータ読取可能な媒体。
【請求項８】
前記各放射線検出器に対するシングルス有効計数時間補正係数が、更に、閾値より高いシングルス率に対する前記放射線検出器の処理限界による不感時間要因から決定される、請求項１から７の何れか一項に記載の非一時的コンピュータ読取可能な媒体。
【請求項９】
前記各同時計数線に対する有効計数時間補正係数が、当該同時計数線により接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定されると共に、更に前記撮像データのフレームに関するシステムシングルス率に依存する該同時計数線に対する同時計数係数から決定される、請求項１から８の何れか一項に記載の非一時的コンピュータ読取可能な媒体。
【請求項１０】
複数の放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率を決定するステップと、
各放射線検出器に対するエネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定するステップと、
低エネルギトリガ補正係数を、前記撮像データのフレームの収集の間において前記放射線検出器に入射するガンマ線のエネルギスペクトル分布の（ｉ）収集エネルギ窓に入る部分と、（ii）前記収集エネルギ窓及び下側トリガエネルギ窓を含む拡張エネルギ窓に入る部分との比に基づいて決定するステップと、
各放射線検出器に対するシングルス有効計数時間補正係数を、前記シングルス率、当該放射線検出器に対して決定された前記エネルギ補正係数及び前記低エネルギトリガ補正係数から決定するステップと、
放射線検出器の対を接続する複数の同時計数線の各同時計数線に対して、当該同時計数線に対する有効計数時間補正係数を、当該同時計数線により接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定するステップと、
前記撮像データのフレームを前記決定された同時計数線に対する有効計数時間補正係数を用いて再構成するステップと、
前記再構成された撮像データのフレームを表示するステップ及び該再構成された撮像データのフレームを記憶するステップの少なくとも一方と、
を有する、画像再構成方法。
【請求項１１】
前記エネルギ補正係数を決定するステップは、
エネルギＥのガンマ線により生成されるトリガの統計的計数である関数Ｎ(Ｅ)を決定するステップと、
各放射線検出器に対する前記エネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布にわたってＮ(Ｅ)を平均することにより決定するステップと、
を有する、請求項１０に記載の画像再構成方法。
【請求項１２】
前記撮像データのフレームの収集の間において各放射線検出器に入射するガンマ線のエネルギスペクトル分布を、前記撮像データのフレームから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定するステップを更に含む、請求項１１に記載の画像再構成方法。
【請求項１３】
前記複数の放射線検出器の各放射線検出器に対して前記エネルギ補正係数を決定するために同一の関数Ｎ(Ｅ)が使用される、請求項１１又は請求項１２に記載の画像再構成方法。
【請求項１４】
前記関数Ｎ(Ｅ)を決定するステップが、各エネルギＥにおけるＮ(Ｅ)の平均値を得るためにシミュレーションされた線源からの放射線検出器へのエネルギ蓄積のモンテカルロシミュレーションを実行するステップを有する、請求項１１から１３の何れか一項に記載の画像再構成方法。
【請求項１５】
前記撮像データのフレームの収集の間において各放射線検出器に入射するガンマ線の前記エネルギスペクトル分布を、前記撮像データのフレームから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定するステップを更に含む、請求項１０に記載の画像再構成方法。
【請求項１６】
前記各同時計数線に対する有効計数時間補正係数が、当該同時計数線により接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定されると共に、更に前記撮像データのフレームに関するシステムシングルス率に依存する該同時計数線に対する同時計数係数から決定される、請求項１０から１５の何れか一項に記載の画像再構成方法。
【請求項１７】
複数の放射線検出器を有する画像収集装置と、
少なくとも１つの電子プロセッサと、
を有する撮像システムであって、前記少なくとも１つの電子プロセッサが、
前記放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率を決定し、
低エネルギトリガ補正係数を、前記撮像データのフレームの収集の間において前記放射線検出器に入射するガンマ線のエネルギスペクトル分布の（ｉ）収集エネルギ窓に入る部分と、（ii）前記収集エネルギ窓及び下側トリガエネルギ窓を含む拡張エネルギ窓に入る部分との比に基づいて決定し、
各放射線検出器に対するシングルス有効計数時間補正係数を、前記シングルス率及び前記低エネルギトリガ補正係数から決定し、
放射線検出器の対を接続する複数の同時計数線の各同時計数線に対して、当該同時計数線に対する有効計数時間補正係数を、当該同時計数線により接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定し、
前記撮像データのフレームを前記決定された同時計数線に対する有効計数時間補正係数を用いて再構成する、
ようにプログラムされる、
撮像システム。
【請求項１８】
前記少なくとも１つの電子プロセッサが、
各放射線検出器に対するエネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定し、
各放射線検出器に対するシングルス有効計数時間補正係数を、更に、当該放射線検出器に対して決定された前記エネルギ補正係数から決定する、
ようにプログラムされる、
請求項１７に記載の撮像システム。
【請求項１９】
前記少なくとも１つの電子プロセッサが、
エネルギＥのガンマ線により生成されるトリガの統計的計数である関数Ｎ(Ｅ)を決定し、
各放射線検出器に対する前記エネルギ補正係数を、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布にわたってＮ(Ｅ)を平均することにより決定する、
ようプログラムされる、
請求項１８に記載の撮像システム。
【請求項２０】
前記少なくとも１つの電子プロセッサが、
前記撮像データのフレームの収集の間において各放射線検出器に入射するガンマ線のエネルギスペクトル分布を、前記撮像データのフレームから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定する、
ようにプログラムされる、請求項１９に記載の撮像システム。
【発明の詳細な説明】
【技術分野】
【０００１】
以下は、広くは、医療撮像技術、放射撮像技術、陽電子放射トモグラフィ（ＰＥＴ）撮像技術、単一光子放射コンピュータトモグラフィ（ＳＰＥＣＴ）撮像技術、医療画像判読技術、画像再構成技術及び関連する技術に関するものである。
【背景技術】
【０００２】
不感時間（dead-time）及びパイルアップ（pile-up；積み重なり）は、視野（ＦＯＶ）内の放射能に対して非線形応答を生じさせる陽電子放射トモグラフィ（ＰＥＴ）撮像システムの共通の問題である。不感時間は、２つの入力信号を別のイベントとして分解するための最小処理時間量として定義することができる。当該検出器が前のイベントを処理するために依然として忙しい間に新たな有効なイベントが到来した場合、該新たなイベントは逃され、計数し損ないにつながり得る。パイルアップの場合、２以上の到来光子が同一の検出器にエネルギを蓄積するが、単一のイベントとしてトリガされ、このことは不感時間損失と同様であり、読出エネルギを実際のエネルギよりも高くさせる。このことは、到来イベントのパイルアップシフトエネルギスペクトルを生じさせる。イベントのエネルギが、散乱光子のようにエネルギ窓の下側閾より元々低いものから該エネルギ窓内へシフトされると、このことは計数を得させる。イベントのエネルギをエネルギ窓内からエネルギ窓の高い方の閾より上側にシフトさせた場合、このことは計数損ないを生じさせる。不感時間及びパイルアップの全体としての効果は、“実効的”不感時間と見ることができる（Esther Vicente, Joaquin Lopez Herraiz他による文献“Improved dead-time correction for PET scanners: application to small-animal PET”, Physics in Medicine and Biology, March 2013参照）。一般的に、アクティビティが高いほど、不感時間による計数の損失（減損）の割合が高くなる。全体の不感時間損失割合が、被写体の大きさ、形状及び放射能分布と共に変化することも知られている（Yamamoto S, Amano M, Miura S他による文献“Deadtime correction method using random coincidence for PET”, J. Nucl. Med. 27 1925-8参照）。
【０００３】
定量的結果を実現するために、ＰＥＴデータは不感時間により発生する損失が補正されるべきである。
【０００４】
現在のＰＥＴシステムに配備される不感時間補正は、典型的に、画像再構成フェーズの間に実施される。簡単に言うと、このようなシステムは、不感時間による損失を全ての検出器の全システム平均シングルス率（個々の検出器のシングルス率（singles rate：シングルス計数率）から区別するためにシステムシングルス率と称する）の関数としてモデル化する。この方法は２つの仮定：即ち、当該システムシングルス率が寸法、形状又は放射能分布に無関係に同一である限り、（１）全ての応答線（同時計数線：ＬＯＲ）は不感時間による同じ量（又は割合）の計数損失を受け；（２）全てのＰＥＴ撮像調査は不感時間による同じ量（又は割合）の計数損失を受ける；という仮定に基づいている。更に詳細には、該方法は、１つのベッド位置における開始から終了まで平均化されたシステムシングルス率を採用し、不感時間補正係数vsシステムシングルス率のルックアップテーブルからスケーリング係数を検索し、次いで、この補間されたスケーリング係数を再構成の間における不感時間補正のために全てのイベント又は全てのＬＯＲに適用する。上記ルックアップテーブルは、均一に分布された放射能により充填され、当該ＰＥＴシステムのＦＯＶの中心に配置された円柱状ファントム（例えば、直径２０cmの）を用いる校正手順により事前に発生される。該ルックアップテーブルは、次いで、新たな校正により置換されるまで、被写体形状及び放射能分布の変化にも拘わらず全ての順次の調査のために使用される。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記方法は、幾つかの欠点を有している。個々の検出器（又は検出器ブロック）のシングルス計数率は被写体の活動及び減衰分布に依存して変化するので、不感時間による損失も各検出器及びＬＯＲにより変化する。上記方法は、この違いを無視し、全てのＬＯＲが不感時間による同じ割合の損失を共有すると仮定している。
【０００６】
更に、上記不感時間校正ルックアップテーブルはＦＯＶの中心に配置された特定の均一円柱状ファントムからのデータを用いて発生され、次いで、１つのパラメータ（システムシングルス計数率）のみに基づいて種々の患者／被写体の全ての順次の調査に適用される。実際には、例えば一層大きなファントム、線状ソース（線状線源）、点ソース（点線源）又は均一でない患者分布を用いる異なる設備の場合、実際の不感時間損失は、当該システムのシングルス計数率が同一であっても校正調査で観察されるものとは大幅に異なり得る。一般的に、ルックアップテーブルに記載される不感時間補正係数vsシステムシングルス計数率との間の関係は、この均一な円柱状ファントム設定のみに類似する被写体に対して有効である。他の設備に適用される場合、又は臨床的患者の場合、大きな偏り（バイアス）が生じ得る。事実、このことは、線状ソースが冷たい円柱状背景内に中心をずらして配置されるアメリカ電機工業会（ＮＥＭＡ）精度（相対計数率誤差）試験に対する主要な性能低下因子である。
【課題を解決するための手段】
【０００７】
以下は、これらの問題を克服する新規且つ改善されたシステム及び方法を開示する。
【０００８】
開示される一態様において、非一時的コンピュータ読取可能な媒体は、画像再構成方法を実行するために少なくとも１つの電子プロセッサを含むワークステーションにより読取可能且つ実行可能な命令を記憶する。上記方法は：複数の放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率（singles rates；シングルス計数率）を決定するステップと；各放射線検出器に対するエネルギ補正係数Ｎwgtを、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定するステップと；各放射線検出器に対するシングルス有効計数時間（singles live time）補正係数を、前記シングルス率及び当該放射線検出器に対して決定された前記エネルギ補正係数から決定するステップと；放射線検出器の対を接続する複数の同時計数線（応答線；ＬＯＲ）の各ＬＯＲに対して、当該ＬＯＲに対する有効計数時間補正係数を、当該ＬＯＲにより接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定するステップと；前記撮像データのフレームを前記決定されたＬＯＲに対する有効計数時間補正係数を用いて再構成するステップと；を有する。
【０００９】
開示される他の態様において、画像再構成方法は：複数の放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率を決定するステップと；各放射線検出器に対するエネルギ補正係数Ｎwgtを、前記撮像データのフレームの収集の間において当該放射線検出器に入射するガンマ線のエネルギスペクトル分布に基づいて決定するステップと；低エネルギトリガ補正係数Ｒを、前記撮像データのフレームの収集の間において前記放射線検出器に入射するガンマ線のエネルギスペクトル分布の（ｉ）収集エネルギ窓（In-Ewindow；Ｅ窓内）に入る部分と、（ii）前記収集エネルギ窓及び下側トリガエネルギ窓（Below-Ewindow；Ｅ窓未満）を含む拡張エネルギ窓（Full-Ewindow；全Ｅ窓）に入る部分との比に基づいて決定するステップと；各放射線検出器に対するシングルス有効計数時間補正係数を、前記シングルス率、当該放射線検出器に対して決定された前記エネルギ補正係数及び前記低エネルギトリガ補正係数から決定するステップと；放射線検出器の対を接続する複数の同時計数線（ＬＯＲ）の各ＬＯＲに対して、当該ＬＯＲに対する有効計数時間補正係数を、当該ＬＯＲにより接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定するステップと；前記撮像データのフレームを前記決定されたＬＯＲに対する有効計数時間補正係数を用いて再構成するステップと；前記再構成された撮像データのフレームを表示するステップ及び該再構成された撮像データのフレームを記憶するステップの少なくとも一方と；を含む。
【００１０】
開示される他の態様において、撮像システムは、複数の放射線検出器を有する画像収集装置を含む。少なくとも１つの電子プロセッサは：前記放射線検出器により検出される撮像データのフレームにおいて、該複数の放射線検出器のシングルス率を決定し；低エネルギトリガ補正係数Ｒを、前記撮像データのフレームの収集の間において前記放射線検出器に入射するガンマ線のエネルギスペクトル分布の（ｉ）収集エネルギ窓（Ｅ窓内）に入る部分と、（ii）前記収集エネルギ窓及び下側トリガエネルギ窓（Ｅ窓未満）を含む拡張エネルギ窓（全Ｅ窓）に入る部分との比に基づいて決定し；各放射線検出器に対するシングルス有効計数時間補正係数を、前記シングルス率及び前記低エネルギトリガ補正係数Ｒから決定し；放射線検出器の対を接続する複数の同時計数線（ＬＯＲ）の各ＬＯＲに対して、当該ＬＯＲに対する有効計数時間補正係数を、当該ＬＯＲにより接続される放射線検出器の対の前記決定されたシングルス有効計数時間補正係数から決定し；前記撮像データのフレームを前記決定されたＬＯＲに対する有効計数時間補正係数を用いて再構成する；ようにプログラムされる。
【００１１】
１つの利点は、改善された不感時間補正を伴う放射撮像システムを提供することにある。
【００１２】
他の利点は、個々のＬＯＲの間の及び撮像被写体毎の不感時間変動を補償することにある。
【００１３】
他の利点は、異なる撮像パラメータの間の不感時間損失の偏り（例えば、一層大きなファントム、線状線源（ソース）、点線源又は不均一な患者の分布等）を補償することにある。
【００１４】
所与の実施態様は、本開示を精読及び理解すれば当業者にとり明らかとなるように、上記利点の幾つか及び／又は他の利点を提供することができる。
【００１５】
本開示は、種々の構成要素及び構成要素の配置、並びに種々のステップ及びステップの配置の形態をとることができる。図面は、好ましい実施態様を図示するためだけのものであり、本開示を限定するものと見なされるべきではない。
【図面の簡単な説明】
【００１６】
【図１】図１は、一態様による画像再構成システムを概略的に示す。
【図２】図２は、図１のシステムの動作の例示的フローチャートを示す。
【図３】図３は、図１のシステムの動作の他の例示的フローチャートを示す。
【図４】図４は、図１のシステムに関するＮ(Ｅ)関係の一例を示す。
【図５】図５は、図１のシステムのシステムシングルス率λsystemの関数としてのＬcoincidenceの一例を示す。
【発明を実施するための形態】
【００１７】
不感時間は、高計数率において最も重要な補正である。この場合において、検出器がガンマ線を検出するために掛かる有限時間（即ち、“不感時間”）は、該検出器が後続のガンマ線を喪失し得る無視できない蓋然性を生じる。ＰＥＴにおいては応答線（同時計数線：ＬＯＲ）が検出されるので、関わる各検出器に対して不感時間を考慮しなければならない。対照的に、単一光子放射コンピュータトモグラフィ（ＳＰＥＣＴ）等の幾つかの他の撮像方式によれば、単一の光子が検出されるので、単一の検出器に対する不感時間のみが関係する。ＰＥＴにおける５１１keVガンマ線の検出の確率は各検出器に対する“有効計数時間（生き時間：live time）”の積を用いて計算可能であり、ここで、該有効計数時間は不感時間の逆である。個々の検出器の不感時間とは別に、ＬＯＲに沿う同時計数を発生することに関連する不感時間成分も存在し得る。この“同時検出”成分は、例えば、２つの時間的に一致する検出イベントを識別して、これらを同時計数として照合するための同時計数検出ユニット（ＣＤＵ）における処理時間を占め得る。ここでも、積の形成を可能にするために有効計数時間を用いると、これはＬＡｘＬＢｘＬcoinとなり、ここで、ＬＡ及びＬＢは別々に計算された当該２つの各検出器に対する有効計数時間であり、ＬcoinはＬＯＲ同時計数形成に関連する有効計数時間成分である。
【００１８】
不感時間を計算する目的のために利用可能なデータは、リストモードデータ、即ちタイムスタンプが付されたシングルスイベントのリストである。ここに開示される実施態様は、既存の不感時間計算を、リストモードデータで利用可能な（又はリストモードデータから推測可能な）低レベル情報を利用して追加の因数を考慮に入れることにより改善する。
【００１９】
第１に、不感時間はリストモードデータのみから推定されるものより大きくなり得る。何故なら、複数の検出器が単一のガンマ線によりトリガされることに起因し得る付加的不感時間が存在するからである。この場合において、“余分な”トリガはリストモードイベントとしては記録されない。何故なら、計算されるエネルギが“Ｅ窓内（In-Ewindow）”、即ちイベントが５１１keVガンマ線であるとされる大凡５１１keVを中心とするエネルギ窓、内にないからである。それでも、これらの余分なトリガは、該余分なトリガの間に到達する他のガンマ線が喪失され得る限り、不感時間に寄与するものである。これらの余分なトリガを考慮するために、形式Ｎ(Ｅ)の関係を採用する“1-to-N”補正が開発され、ここで、Ｅはガンマ線エネルギであり、Ｎ(Ｅ)はエネルギＥのガンマ線に対して統計的に発生するトリガの数である。例示的実施態様において、関数Ｎ(Ｅ)はモンテカルロシミュレーションにより推定される。これは一度限りのシミュレーションとして実行することができ、その後、1-to-N不感時間補正は該発生されたＮ(Ｅ)関数を使用する。
【００２０】
第２に、不感時間は低エネルギガンマ線でのトリガにより影響を受け得る。これらのガンマ線は、リストモードイベントとして記録されるには低過ぎるが、全エネルギが低過ぎると判定すべく信号を積分するために時間が掛けられるほど高い“Ｅ窓未満（Below-Ewindow）”範囲内に入る。例えば、Ewindowの下側閾値が４８０keVである場合におけるエネルギ４５０keVのイベントを考察する。この場合、該４５０keVのガンマ線は検出シーケンスをトリガし、該シーケンスはエネルギが４５０keVであると決定する。このエネルギはリストモードイベントのための下側閾値（４８０keV）より低いので、破棄される。しかしながら、この処理は、該４５０keVのガンマ線の検出の処理の間に到達する他のガンマ線が喪失され得る限りにおいて、不感時間に寄与する。このような効果を考慮するべく、各検出器ユニットの記録された（In-Ewindow）計数に対する全窓（即ち、In-Ewindow + Below-Ewindow）計数の比を推定するためにモンテカルロシミュレーションが採用される。
【００２１】
第３の効果は、ハードウェアの限界による喪失計数（数え落とし）である。例えば、幾つかのＰＥＴシステムは、メモリのオーバーフローにより、計数率が極端に高いと限界に達する。ここに開示される幾つかの実施態様においては、このハードウェア限界を考慮するために追加の不感時間が設けられる。
【００２２】
オプションとして考慮される第４の効果は、光電ピークをずらすことにより不感時間に影響を与え得る温度である。
【００２３】
ここに開示される幾つかの実施態様において、ＬＯＲ補正Ｌcoinはファントム測定から決定される。これは、この場合も、１つのタイプの検出システムに対して１回限りの校正であり得る。
【００２４】
ここではＰＥＴ撮像システムに関して説明されるが、開示される方法は、コンピュータトモグラフィ（ＣＴ）撮像システム（光子計数モードで動作する）、単一光子放射コンピュータトモグラフィ（ＳＰＥＣＴ）撮像システム、又はＰＥＴ／ＣＴ若しくはＳＰＥＣＴ／ＣＴ撮像システムのようなハイブリッドシステム等の他の撮像方式と関連して使用することもできる。撮像関連イベントがシングルス（ＰＥＴにおけるようにＬＯＲに沿う同時計数ではなく）であるＳＰＥＣＴ又はＣＴ等の方式において、同時計数態様に関係する当該開示される不感時間補正係数は組み込まれない。
【００２５】
本明細書で使用される場合、“Ｅ窓内（In-Ewindow）”なる用語（及び、その変形）は、蓄積されて検出されたエネルギが収集エネルギ窓内に入る光子を参照する。言い換えると、エネルギがIn-Ewindow内に入るイベントはリストモードデータに含まれるが、エネルギがIn-Ewindowの外側となるイベントはリストモードデータに含まれない。
【００２６】
本明細書で使用される場合、“Ｅ窓未満（Below-Ewindow）”なる用語（及び、その変形）は、エネルギが収集エネルギ窓の下側閾値より低いがＰＥＴ検出器をトリガする程十分に高い光子を参照する。
【００２７】
本明細書で使用される場合、“全Ｅ窓（Full-Ewindow）”なる用語（及び、その変形）は、In-Ewindow及びBelow-Ewindowの組み合わせを指す。
【００２８】
本明細書で使用される場合、“全入射シングルス（All incidence singles）” なる用語（及び、その変形）は、エネルギが検出処理をトリガする程十分に高い入射光子を参照する。全入射シングルスは、エネルギ窓評価基準を満たす光子及びエネルギがＥ窓未満である故に拒絶される光子の両方を含む。
【００２９】
本明細書で使用される場合、“検出器（Detector）”なる用語（及び、その変形）は、不感時間に関して同一の処理ユニットに属する個々の放射線検出器を指す。該検出器は、ＰＥＴシステムに依存して単一の検出器（クリスタル）又は一群の検出器であり得る。ＳＰＥＣＴの場合、単一の検出器は、適切には、ハニカムコリメータ又は他のタイプの放射線コリメータの開口に位置合わせされた検出器である。
【００３０】
本明細書で使用される場合、“有効計数時間（生き時間：Live-time）”なる用語（及び、その変形）は、入力イベントのうちの不感時間による計数の損失を伴う出力イベントの割合（又は等価的に確率）を指す。例えば、有効計数時間は０〜１の論理範囲を有する。１なる有効計数時間は不感時間無しに対応するであろう。０.８なる有効計数時間は、統計的に言って、収集時間の８０％は利用可能である一方、検出トリガを処理することにより占有されることにより当該時間の２０％は利用可能でないことを示すであろう。これらは、解説的な例に過ぎない。
【００３１】
本明細書で使用される場合、“測定されたシングルス率”なる用語（及び、その変形）は、検出システムから出力された、エネルギ窓の基準を満たす（即ち、“Ｅ窓内”の）シングルス率を示す。該シングルス率は、収集（利用可能なら）から直接的に得ることができ、又は遅延同時計数率等の他の情報から間接的に導出することができる。両者は測定されたシングルス率として扱われる。
【００３２】
本明細書で使用される場合、“システムシングルス率（λsystem）”なる用語（及び、その変形）は、全検出器にわたる測定されたシングルス率の和を示す。このシングルス率は、個々の検出器の測定されたシングルス率と区別するためのものである。
【００３３】
本明細書で使用される場合、“フレーム”なる用語（例えば、“撮像データのフレーム”）は、対応する再構成画像（即ち、“フレーム画像”）を生成するために再構成される、時間間隔にわたり収集された撮像データセットを示す。フレームは、軸方向に拡張された最終的再構成画像を生成するために、斯様な各フレームが対応するフレーム画像を一緒に連結されたフレーム画像と共に発生するように再構成される複ベッドＰＥＴ撮像の前後関係において固定されたベッド位置で収集することができる。他の例として、“単ベッド”ＰＥＴ撮像の前後関係において、フレームは、（単一の）対応するＰＥＴ画像を生成するために再構成される、単一のベッド位置における患者で収集される全データセットとすることもできる。フレームは時間的に延長することができ、特定の時間間隔の後に収集され、動的な又はゲートされたＰＥＴとすることもできる。
【００３４】
図１を参照すると、例示的医療撮像システム１０が図示されている。図１に示されるように、システム１０は画像収集装置１２を含んでいる。一例において、画像収集装置１２は、コンピュータトモグラフィ（ＣＴ）ガントリ１３を更に含む、ＰＥＴ／ＣＴ撮像システムのＰＥＴガントリを有することができる。他の例において、画像収集装置１２は、ＣＴ部品を備えない単独型ＰＥＴスキャナとすることができる。患者台（又はベッド）１４は患者をＰＥＴガントリ１２又はＣＴガントリ１３の検査領域１６内に装填するように配置される。ＰＥＴガントリ１２は放射線検出器１７のアレイ（図１には概略的に示されており、典型的には、ＰＥＴガントリ１２の放射線検出器は軸方向のＦＯＶに跨がるように配置された一連のＰＥＴ検出器リングとして構成される）を含む。
【００３５】
システム１０は、少なくとも１つの電子プロセッサ２０、少なくとも１つのユーザ入力装置（例えば、マウス、キーボード及び／又はトラックボール等）２２及び表示装置２４等の典型的な構成部品を備えるコンピュータ、ワークステーション又は他の電子データ処理装置１８も含む。幾つか態様において、表示装置２４はコンピュータ１８とは別個の部品であり得る。ワークステーション１８は１以上の非一時的記憶媒体２６（磁気ディスク、ＲＡＩＤ若しくは他の磁気記憶媒体；固体ドライブ、フラッシュドライブ、電子的に消去可能なリードオンリメモリ（EEROM）若しくは他の電子メモリ；光ディスク若しくは他の光記憶部；又はこれらの種々の組み合わせ等）も含むことができる。表示装置２４は、ユーザ入力装置２２からユーザ入力を受けるための１以上のフィールドを含むグラフィックユーザインターフェース（ＧＵＩ）２８を表示するように構成される。
【００３６】
当該少なくとも１つの電子プロセッサ２０は１以上の非一時的記憶媒体２６に動作的に接続され、該記憶媒体は、画像再構成方法又は処理１００を実行することを含む開示された処理を実行するために該少なくとも１つの電子プロセッサ２０により読み出し可能且つ実行可能な命令を記憶する。幾つかの例において、当該画像再構成方法又は処理１００は、少なくとも部分的にクラウド処理により実行することができる。
【００３７】
図２を参照すると、画像再構成方法１００の例示的実施態様がフローチャートとして概略的に示されている。当該処理を開始するために、前記画像収集装置１２（例えば、ＰＥＴ撮像装置）は、放射線検出器１７により検出されるＰＥＴ撮像データを収集するように構成され、又は前記少なくとも１つの電子プロセッサ２０により制御される。１０２において、当該少なくとも１つの電子プロセッサ２０は、放射線検出器により検出された撮像データのフレームにおいて、複数の放射線検出器１７のシングルス計数率を決定するようにプログラムされる。これは、個々のＰＥＴ検出器のシングルス率をリアルタイムに記録するＰＥＴシステムに関してはハードウェアを用いて実行することができ、又は収集されたリストモードＰＥＴ撮像データから計算することができる。例えば、一方法において、各ＬＯＲに関するシングルス率は遅延同時計数をカウントすることにより決定され（即ち、ＰＥＴ撮像スキャナの利用可能な同時計数検出能力を“遅延”同時計数時間窓を用いることにより活用する）、所与のＰＥＴ検出器に関するシングルス率は該ＰＥＴ検出器を含むＬＯＲに関する遅延同時計数データに基づいて決定される。
【００３８】
１０４において、当該少なくとも１つの電子プロセッサ２０は、各放射線検出器１７に対するエネルギ補正係数（Ｎwgt）を、当該フレームの撮像データの収集の間に該放射線検出器に入射したガンマ線のエネルギスペクトル分布に基づいて決定するようにプログラムされる。このようにするために、該少なくとも１つの電子プロセッサ２０は関数Ｎ(Ｅ)の値を決定するようにプログラムされ、ここで、Ｎ(Ｅ)はエネルギＥのガンマ線により生成されたトリガの統計的計数である。幾つか態様において、Ｎ(Ｅ)は各エネルギＥにおけるＮの平均値を求めるためにシミュレーションされたソースからの放射線検出器におけるエネルギ蓄積のモンテカルロシミュレーションを用いて事前に決定され、結果的Ｎ(Ｅ)関数は経験関数又はルックアップテーブル等として記憶される。次いで、当該少なくとも１つの電子プロセッサ２０は、各放射線検出器に対するエネルギ補正係数（Ｎwgt）を、Ｎ(Ｅ)を当該フレームの撮像データの収集の間に放射線検出器に入射したガンマ線のエネルギスペクトル分布にわたって平均することにより決定するようプログラムされる。当該複数の放射線検出器の各放射線検出器に対してエネルギ補正係数（Ｎwgt）を決定するために同一の関数Ｎ(Ｅ)が用いられる。
【００３９】
幾つか実施態様において、当該少なくとも１つの電子プロセッサ２０は、当該フレームの撮像データの収集の間に各放射線検出器に入射するガンマ線のエネルギスペクトル分布を、該フレームの撮像データから再構成される初期画像に対して実行されるモンテカルロシミュレーション処理を用いて決定するようにプログラムされる。
【００４０】
１０６において、当該少なくとも１つの電子プロセッサ２０は、低エネルギトリガ補正係数（Ｒ）を、当該フレームの撮像データの収集の間に放射線検出器に入射するガンマ線のうちの（ｉ）収集エネルギ窓（Ｅ窓内）内と、（ii）該収集エネルギ窓及び低トリガエネルギ窓（Ｅ窓未満）を含む拡張エネルギ窓（全Ｅ窓）内に入るエネルギスペクトル分布の部分との比に基づいて決定するようプログラムされる。幾つか例において、当該少なくとも１つの電子プロセッサ２０は、当該フレームの撮像データの収集の間に各放射線検出器に入射するガンマ線のエネルギスペクトル分布を、該フレームの撮像データから再構成される初期画像に対して実行されるモンテカルロシミュレーションを用いて決定するようにプログラムされる。単一のモンテカルロシミュレーションが、ステップ１０４及び処理１０６による比Ｒに関して各放射線検出器に入射するガンマ線のエネルギスペクトル分布を決定するための情報を生成することができると注記することができ、更に、処理１０４、１０６の順序を逆にすることができ、及び／又は共通のモンテカルロシミュレーション処理要素を用いて統合することができることが理解される。
【００４１】
１０８において、当該少なくとも１つの電子プロセッサ２０は、各放射線検出器に対するシングルス有効計数時間補正係数を、前記シングルス率（１０２からの）、当該放射線検出器に関して決定されたエネルギ補正係数（１０４からの）及び低エネルギトリガ補正係数（１０６からの）から決定するようプログラムされる。幾つかの例において、各放射線検出器に対する該シングルス有効計数時間補正係数は、更に、閾値より高いシングルス率に対する当該放射線検出器の処理限界による不感時間係数から決定される。
【００４２】
１１０において、当該少なくとも１つの電子プロセッサは、放射線検出器の対を接続する複数のＬＯＲの各ＬＯＲに対して、当該ＬＯＲに対する有効計数時間補正係数を、該ＬＯＲにより接続される放射線検出器の対に関する前記決定されたシングルス有効計数時間補正係数から決定するようプログラムされる。他の例において、各放射線検出器に対するシングルス有効計数時間補正係数は、更に、温度依存性エネルギピーク光電シフトを考慮に入れることにより温度センサの温度値から決定される。
【００４３】
１１２において、当該少なくとも１つの電子プロセッサ２０は、当該フレームの撮像データを前記決定されたＬＯＲに対する有効計数時間補正係数を用いて再構成するようプログラムされる。該少なくとも１つの電子プロセッサ２０は、次いで、表示装置２４を該再構成されたフレーム画像を表示するように制御し及び／又は非一時的記憶媒体２６を適切なデータベース（例えば、画像保存通信システム（Picture Archiving and Communications System））に当該臨床画像を記憶するように制御するようプログラムされる。多局撮像の場合、図２の処理は各フレームに関して適切に繰り返され、結果としてのフレーム画像は処理１１２において一緒に結合され、軸方向に拡張された再構成ＰＥＴ画像として表示される。
【００４４】
図３を参照すると、画像再構成方法２００の更に詳細な例示的実施態様がフローチャートとして図式的に示されている。２０２において、再構成されるべき撮像データが収集される（例えば、画像収集装置１２により）。２０４において、画像収集装置１２のシングルス率ＳＴが決定される。２０６において、シングルス率ＳＴは同時計数有効計数時間係数を決定するために同時計数有効計数時間モデルに入力される。２０８において、当該撮像データの遅延同時計数が決定される。２１０において、該遅延同時計数は各検出器１７の測定されたシングルス率を決定するためにシングルス率モデルに入力される。代わりに、各検出器１７のシングルス率は、画像収集装置１２のハードウェアから得られ得る。２１２において、当該撮像データの初期再構成が、中間再構成画像及び減衰マップを生成するために実行される。２１３において、該中間再構成画像及び減衰マップは、全対内比（full-to-in-ratio）Ｒ２１４及びエネルギスペクトルを決定するためにモンテカルロシミュレータに入力される。２１６において、上記エネルギスペクトルはＮ値を決定するために１→Ｎモデルに入力され、ここで、Ｎはガンマ線のエネルギにより生成されるトリガの統計的計数である。２１８において、測定された検出器シングルス率（２１０からの）、全対内比Ｒ（２１４からの）及びＮ値（２１６からの）は、シングルス有効計数時間補正係数を決定するためにシングルス有効計数時間ソルバに入力される。２２０において、同時計数有効計数時間補正係数（２０６からの）及びシングルス有効計数時間補正係数（２１８からの）は、最終的画像再構成処理に使用されるべきＬＯＲに対する有効計数時間係数を決定するために３要素コンバイナに入力される。
【００４５】
（実施例１）
以下の実施例においては、図３の処理の種々の演算の例示的構成が更に詳細に説明される。
【００４６】
ここに記載されるように、各ＬＯＲ（対応するＬＯＲを定める２つの検出器Ｌ（Ａ，Ｂ）を用いて示される）の有効計数時間係数は、式（１）に示されるように、３つの成分：即ち、検出器Ａにおけるシングルス検出の有効計数時間、Ｂにおけるシングルス検出の有効計数時間、及び同時検出形成部分の有効計数時間の積としてモデル化される。
【数１】
2021511523.tif 000003
この式において、λm,A及びλm,Bは、各々、検出器Ａ及びＢにおいて測定されたシングルス率であり；ＥＡ及びＥＢは、各々、検出器Ａ及びＢにおける入射シングルスのエネルギスペクトルであり；ＲＡ及びＲＢは、各々、Ａ及びＢにおける入射全Ｅ窓シングルス率の入射Ｅ窓内シングルス率に対する比であり；Ｌsingleはλｍ、Ｅ及びＲの関数としての有効計数時間であり、この時点における検出器センサの温度等の他の要因は無視する。Ｌcoincidenceは、測定されたシステムシングルス率λsystemの関数として近似される同時計数形成部分の有効計数時間である。
【００４７】
各ＬＯＲに関して導出されたシングルス有効計数時間は、当該イベントを検出する対応する放射線検出器Ａ，Ｂにおける全ピクセルに適用される。
【００４８】
Ｌsingleを導出するために、各検出器におけるシングルス検出の有効計数時間は、典型的に、実際のシステムに依存して非無力化可能（non-paralyzable；式２）、無力化可能（式３）又は混合された非無力化可能及び無力化可能（式４）を用いてモデル化される。これらの式において、τは各トリガされたイベントに対する平均での実効処理時間であり、λeventsは検出器毎の全トリガイベント率であり、θは前記混合モデルにおける非無力化可能成分の部分である。システムが完全に無力化可能なら、θはゼロである。
【数２】
2021511523.tif 000004
【数３】
2021511523.tif 000005
【数４】
2021511523.tif 000006
【００４９】
殆どの調査において、λeventsτは相対的に小さいので、式（２）〜（４）はテイラー展開を用いて（５）〜（７）に近似することができる。
【数５】
2021511523.tif 000007
【数６】
2021511523.tif 000008
【数７】
2021511523.tif 000009
【００５０】
１）Ｅ窓未満光子の寄与及び２）１→Ｎトリガ効果により、λeventsは、
【数８】
2021511523.tif 000010
と扱うことができ、ここで、λinci,in-Ewinは不感時間損失前の同時計数Ｅ窓内シングルス率であり、Ｒは入射全Ｅ窓シングルス率の入射Ｅ窓内シングルス率に対する比であり、Ｎは１→Ｎトリガ効果による係数であり、ｂは背景率（background rate）である。定義により、
【数９】
2021511523.tif 000011
とし、ここで、λMeasは、背景率ｂが既に減算された測定されたシングルス率である。式（８）及び（９）を式（７）に代入すれば、式（１０）を得る。
【数１０】
2021511523.tif 000012
【００５１】
この式は、
【数１１】
2021511523.tif 000013
に再整理することができる。
【００５２】
この式において、λMeas、Ｒ、Ｎ及びｂは、再構成の間における各検出器に対する定数として扱うことができる。この場合、式（１１）は未知の変数λinci,in-Ewinの三次方程式であり、更に、
【数１２】
2021511523.tif 000014
に再整理することができる。
【００５３】
各検出器に関して式（１２）を解くと、同時計数Ｅ窓内シングルス率λinci,in-Ewinが得られる。三次方程式は３つの根、即ち、１つの実根及び２つの共役複素根、又は３つの実根を有する。当該不感時間問題のシナリオにおいては、典型的に１つの実根及び２つの共役複素根のケースとなる。１つの実根だけが、λinci,in-Ewinに対する解となる。一旦λinci,in-Ewinが見付かったなら、シングルス検出の有効計数時間Ｌsingleは、上記式（９）を用いて容易に得ることができる。低計数率調査において、シングルス検出の有効計数時間は１.０に近付く。低計数率調査においてノイズに起因するバイアスを回避するために、有効計数時間は実際には式（１３）を用いて計算することができる。
【数１３】
2021511523.tif 000015
【００５４】
このアイデアに従えば、シングルス検出の有効計数時間を決定するための要点は、各検出に対して三次方程式（１２）を構築することである。該方程式において、背景率ｂ及び実効的イベント処理時間τは、所与のシステムに関して典型的に安定している。これらは、校正ダイナミック調査から適合させることができ、全ての検出器に対してシステムレベルの定数として扱うことができる。他の３つの変数、即ち、λMeas、Ｒ及びＮは、調査に依存して検出器により変化し、以下のように解決しなければならない。
【００５５】
前述したように、λMeasは、もし利用可能なら収集の間に検出器ハードウェアから直接的に発生することができ、又は遅延同時計数から導出することができる。λMeasは既に平均背景率ｂが減算されていることに注意されたい。
【００５６】
Ｎ及びＲを決定するために、所与の検出器システムモデルに関してＮの関係が入射エネルギの関数Ｎ(Ｅ)として構築される。検出器に衝突するガンマ光子は、自身のエネルギを複数の位置又は隣接する検出器ユニットに蓄積し得るので、該光子は１以上（Ｎ）のイベントをトリガし得る。Ｎの平均は、検出器の特性及び入射光子エネルギに依存する。従って、所与のＰＥＴ検出器モデルに関して、空気中の種々の放射エネルギ（５１１、５００、４７５、４５０、…、１００、５０、２５keV等）の点線源からの検出器へのエネルギ蓄積の一連のモンテカルロシミュレーションを実行し、各放射エネルギにおける一次光子当たりの平均のＮを得、次いで、これから入射エネルギの関数としてのＮの関係Ｎ(Ｅ)を構築する。これは、所与の検出器システムモデルに対する１回限りのシミュレーション作業である。Ｎ(Ｅ)関係の一例が図４に示されている。
【００５７】
如何なる再構成においても、放射能分布の中間画像を生成するために初期再構成が実行される。
【００５８】
各検出器に対して２つの事柄、即ち、入射γ光子のエネルギスペクトル分布、Ｈ(Ｅｋ)；及び入射Ｅ窓内シングルス計数に対する入射全Ｅ窓シングルス計数の比；を発生するために、オンザフライでの簡易モンテカルロシミュレーションが中間画像及び減衰マップを用いて実行される。上記比は、式（１２）に必要とされるＲ値である。
【００５９】
エネルギスペクトル分布Ｈ(Ｅｋ)から、正規化されたエネルギスペクトル密度関数Ｐ(Ｅｋ)が計算され、
【数１４】
2021511523.tif 000016
ここで、Ｍはエネルギスペクトルビンの数である。
【００６０】
ステップ１からのＮ(Ｅ)関係及びステップ３からのエネルギスペクトル密度関数Ｐ(Ｅｋ)を用いて、各検出器に対する重み付けＮ値を計算する。
【数１５】
2021511523.tif 000017
Ｎwgtは、式（１２）に必要とされるＮ値である。
【００６１】
各検出器に関してλMeas、Ｒ、Ｎ、τ及びｂ値が式（１２）に代入され、当該三次方程式をλinci,in-Ewinの実根のみに関して解く。
【００６２】
極端に高い計数率においてハードウェア限界に影響を与え得る（メモリバッファのオーバーフロー、データ転送帯域幅によるパッケージ欠落等）ＰＥＴシステムの場合、斯かるデータ損失は典型的な非無力化可能、無力化可能又は混合モデルを用いてモデル化することはできない。このような状況においては、測定された前記Ｅ窓内率λMeas及びトリガされたイベントの導出率（λinci,in-Ewin・Ｒ・Ｎ）の両方が特定の高閾値より高い場合、経験的モデルに基づいてλinci,in-Ewinに追加の調整を適用する。
【００６３】
シングルス検出の有効計数時間Ｌsingleは、非常に低い計数率調査に関しては式（９）又は式（１３）を用いて計算される。
【００６４】
Ｌcoincidenceを導出するために、同時計数形成成分の有効計数時間Ｌcoincidenceは、システムシングルス率λsystemの関数としてモデル化される。所与のＰＥＴカメラモデルに対して、Ｌcoincidenceとλsystemとの間の応答関係は、均一な円柱状ファントムの動的調査からシングルス検出不感時間のみの補正を伴う結果を既知の真実と比較することにより導出することができる。ここでも、これは所与のＰＥＴカメラモデルに対して１回限りの作業であり、全ての調査に適用可能である。Ｌcoincidence(λsystem)の一例が図５に示されている。
【００６５】
如何なる他の調査の再構成の間においても、Ｌcoincidenceは、実際のシステムシングルス率による前記ステップ１からの応答関数に基づいて決定され、全てのイベント及びＬＯＲに当てはまる。
【００６６】
Ｌsingle及びＬcoincidenceの上記導出は、主要な寄与要因のみを考慮に入れている。幾つかの他の要因も有効計数時間に影響を与え得る。例えば、センサの温度の変化も、感度の変化及び光電ピークのシフトにつながり得る。精度の更なる向上のために、これらの要因を当該モデルに含めることができ、斯かる情報が利用可能なら調整として適用することができる。
【００６７】
以上、本開示を好ましい実施態様に関連して説明した。当業者であれば、上記詳細な説明を精読し理解すれば変形例及び代替例を思い付くことができる。本発明は、添付請求項及びその均等物の範囲内に入る限りにおいて、全ての斯様な変形例及び代替例を含むと見なされることを意図するものである。
【図１】
2021511523.tif 000018
【図２】
2021511523.tif 000019
【図３】
2021511523.tif 000020
【図４】
2021511523.tif 000021
【図５】
2021511523.tif 000022
【国際調査報告】
2021511523.tif 000023
2021511523.tif 000024