(57)【要約】【課題】マルチチャネル検出におけるクロストークに起因するエネルギー測定値の劣化を高精度の補正すること。【解決手段】医用画像処理装置は、イベント取得部、補正係数取得部、選択部及び補正部を有する。イベント取得部は、複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すデータを複数の検出イベント毎に取得する。補正係数取得部は、複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得する。選択部は、データに基づいて複数の検出イベントの中からマルチチャネルイベントを選択する。補正部は、補正係数を前記マルチチャネルイベントに属する検出イベントのエネルギーに適用し、補正エネルギーを生成する。【選択図】      図１
2020064045.tif 000002
【特許請求の範囲】
【請求項１】
複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すデータを複数の検出イベント毎に取得するイベント取得部と、
前記複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得する補正係数取得部と、
前記データに基づいて前記複数の検出イベントの中からマルチチャネルイベントを選択する選択部と、
前記補正係数を前記マルチチャネルイベントに属する検出イベントのエネルギーに適用し、補正エネルギーを生成する補正部と、
を具備する医用画像処理装置。
【請求項２】
前記補正部は、前記マルチチャネルイベント毎に、前記マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーと前記補正係数との加算に基づいて前記補正エネルギーを生成する、請求項１記載の医用画像処理装置。
【請求項３】
前記補正部は、前記マルチチャネルイベント毎に、前記マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーと前記補正係数との乗算に基づいて前記補正エネルギーを生成する、請求項１記載の医用画像処理装置。
【請求項４】
前記補正部は、前記補正エネルギーを生成するため、
エネルギーウィンドウのうちの最大エネルギー及び最小エネルギーから前記補正係数を減算することによって前記エネルギーウィンドウをシフトし、
前記マルチチャネルイベント毎に、前記マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーを生成し、
前記シフトされたエネルギーウィンドウ内に合計エネルギーを有するマルチチャネルイベントを選択する、
請求項１記載の医用画像処理装置。
【請求項５】
前記補正部は、前記補正エネルギーを生成するため、
エネルギーウィンドウのうちの最大エネルギー及び最小エネルギーを前記補正係数で除算することによって前記エネルギーウィンドウをシフトし、
前記マルチチャネルイベント毎に、前記マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーを生成し、
前記シフトされたエネルギーウィンドウ内に合計エネルギーを有するマルチチャネルイベントを選択する、
請求項１記載の医用画像処理装置。
【請求項６】
前記データは、更に前記ガンマ線の検出時刻を表し、
前記選択部は、
前記複数の検出イベントのうちの最小エネルギー閾値を上回るエネルギーを有する２以上の特定の検出イベントを選択し、
前記選択された２以上の特定の検出イベントの任意の組合せ各々の中から、入射位置と検出時刻とエネルギーとの少なくとも一要素に基づいて、同一のマルチチャネルイベントに属する２以上の検出イベントを選択する、
請求項１記載の医用画像処理装置。
【請求項７】
前記補正係数は、前記複数の検出器要素内でのガンマ線のヒット数に応じて複数の補正係数を有し、
前記選択部は、前記マルチチャネルイベントを前記ヒット数に応じて分類し、分類後の各ヒット数のマルチチャネルイベントのエネルギーに、前記複数の補正係数のうちの対応するヒット数の補正係数を適用する、
請求項１記載の医用画像処理装置。
【請求項８】
前記補正係数は、前記複数の検出器要素各々の空間座標に応じた複数の補正係数を有し、
前記選択部は、前記マルチチャネルイベントのエネルギーに、前記複数の補正係数のうちの、前記マルチチャネルイベントに属する複数の検出イベントにそれぞれ対応する複数の入射位置の重み付き平均に対応する補正係数を適用する、
請求項１記載の医用画像処理装置。
【請求項９】
前記重み付き平均は、前記複数の検出イベントにそれぞれ対応する複数の入射位置の正規化加重和に基づいて計算され、
前記重みは、前記検出イベントのエネルギーを、非ゼロの実数で累乗した値である、
請求項８記載の医用画像処理装置。
【請求項１０】
前記空間座標は、前記複数の補正係数各々に対応するブロックにセグメント化され、
前記ブロックのサイズは、予め定義された空間分解能基準及び前記複数の補正係数の精度又は信頼性に関する予め定義されたカウント基準のうちの１つ又は複数を満たすように選択される、
請求項８記載の医用画像処理装置。
【請求項１１】
前記空間座標は、前記複数の補正係数各々に対応するブロックにセグメント化され、
前記ブロックのサイズは、マルチチャネルイベントのヒット数に応じて異なる、
請求項８記載の医用画像処理装置。
【請求項１２】
前記補正部は、前記マルチチャネルイベントに属する検出イベントのエネルギーに、前記補正係数を適用する前において、初期の非線形エネルギー補正を適用する、請求項１記載の医用画像処理装置。
【請求項１３】
前記補正部は、前記初期の非線形エネルギー補正として、検出されたガンマ線のエネルギー測定値と実際のエネルギーとの間の非線形関係を補正する、請求項１２記載の医用画像処理装置。
【請求項１４】
前記イベント取得部は、較正段階において、前記ガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表す較正データを複数の較正検出イベント毎に取得し、
前記選択部は、前記較正データに基づいて前記複数の較正検出イベントの中から較正マルチチャネルイベントを選択し、
前記補正係数取得部は、前記較正マルチチャネルイベント毎に、前記マルチチャネルイベントに属する２以上の検出イベントのエネルギーの代表エネルギーと、予め定義されたエネルギーとのシフト量に基づいて前記補正係数を生成する、
請求項１記載の医用画像処理装置。
【請求項１５】
前記補正係数は、乗法補正係数又は加法補正係数を含む、請求項１４記載の医用画像処理装置。
【請求項１６】
前記ガンマ線検出器は、各々が複数の検出器要素を含む複数の検出器要素ブロックを有し、
前記代表エネルギーは、前記複数の検出器要素ブロック各々によるマルチチャネルイベントのカウントのヒストグラムにおける中央値、モード及び平均値のうちの１つに対応するエネルギーである、
請求項１４記載の医用画像処理装置。
【請求項１７】
前記補正エネルギーに基づいて画像を再構成する再構成部を更に備える、請求項１記載の医用画像処理装置。
【請求項１８】
複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すデータを複数の検出イベント毎に取得し、
前記複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得し、
前記データに基づいて前記複数の検出イベントの中からマルチチャネルイベントを選択し、
前記補正係数を前記マルチチャネルイベントに属する検出イベントのエネルギーに適用して補正エネルギーを生成する、
ことを具備する医用画像処理方法。
【請求項１９】
コンピュータに、
複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すデータを複数の検出イベント毎に取得し、
前記複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得し、
前記データに基づいて前記複数の検出イベントの中からマルチチャネルイベントを選択し、
前記補正係数を前記マルチチャネルイベントに属する検出イベントのエネルギーに適用して補正エネルギーを生成する、
ことを実現させる医用画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医用画像処理装置、方法及びプログラムに関する。
【背景技術】
【０００２】
陽電子放射トモグラフィ（ＰＥＴ）において、放射性トレーサ剤は、患者に導入され、放射性トレーサ剤の物理的性質及び生体分子の性質により患者の体内の特定の位置に濃縮する。放射性トレーサ剤は、陽電子を放出する。陽電子が電子と衝突すると、消滅イベントが発生し、略１８０度離れて進む５１１ｋｅＶの２つのガンマ線を生成する。
【０００３】
ＰＥＴシステムは、ガンマ線の同時計数ペアを検出するために患者の周りに配置された検出器を使用する。検出器のリングは、各角度から来るガンマ線を検出するために使用される。従ってＰＥＴシステムの架台は、等方性放射の捕捉を最大化するために実質的に円筒である。検出器は、数千個の個々の結晶（例えば、オルトケイ酸ルテチウム（ＬＹＳＯ）又は他のシンチレータ結晶）から構成される。これらシンチレータ結晶は、シンチレーションイベントからの光パルスを測定するために光検出器を有するモジュール中にパッケージされた２次元シンチレータアレイを構成する。例えば、シンチレータ結晶からの光は、複数の光電子増倍管（ＰＭＴ）の間で共有されるか、又はシンチレータ結晶との１対１の対応を有するシリコン光電子増倍管（ＳｉＰＭ）によって検出される。
【０００４】
断層撮影の再構成の原理を介して放射性トレーサ剤の空間的及び時間的な分布を再構成するために、各検出されたイベントは、それのエネルギー（すなわち、生成された光の量）、それの位置、及びそれのタイミングについて特徴づけられる。２つのガンマ線を検出し、それらの検出位置間に応答線（ＬＯＲ：line-of-response）を引くことによって、対消滅の可能性のある位置を決定することができる。タイミング情報は、また、２つのガンマ線の飛行時間（ＴＯＦ：time-of-flight）情報に基づいて、ＬＯＲに沿う対消滅の統計的分布を決定するために使用される。多数のＬＯＲを蓄積することによって、患者内の放射能（例えば、放射性トレーサ剤の密度）の空間分布の立体画像を再構成することができる。
【０００５】
単光子放出計算機トモグラフィ（ＳＰＥＣＴ）は、各検出器要素（例えば、シンチレータ結晶）上でのガンマ線の入射の立体角を制限するためにコリメータが使用されること以外はＰＥＴと同様であり、ＬＯＲを決定するために一致を必要とするのとは反対に、単一のガンマ線検出イベントを使用して再構成を可能にする。
【０００６】
位置情報（例えば、ＬＯＲ）及びタイミング情報（例えば、ＴＯＦ）に加えて、ＰＥＴシステム及びＳＰＥＣＴシステムの検出器は、画像再構成プロセスにおいてエネルギー情報を取得及び使用する。エネルギー測定値に関しては、測定プロセスにおける非線形性及び／又は実施上の考慮すべき問題点が存在する。当該問題点としては、例えば、ガンマ線エネルギーがコンプトン散乱に起因して複数の検出器又はチャネルに吸収されるマルチチャネルのガンマ線検出におけるチャネル間の光又は電荷共有（クロストーク）がある。これらに起因してエネルギー測定値は理想的な線形応答から逸れる。従ってピクセル型ガンマ線検出器においてエネルギー測定値を補正する改善された技法が望まれている。
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明が解決しようとする課題は、マルチチャネル検出におけるクロストークに起因するエネルギー測定値の劣化を高精度の補正することである。
【課題を解決するための手段】
【０００８】
実施形態に係る医用画像処理装置は、複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すデータを複数の検出イベント毎に取得するイベント取得部と、前記複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得する補正係数取得部と、前記データに基づいて前記複数の検出イベントの中からマルチチャネルイベントを選択する選択部と、前記補正係数を前記マルチチャネルイベントに属する検出イベントのエネルギーに適用し、補正エネルギーを生成する補正部と、を具備する。
【図面の簡単な説明】
【０００９】
【図１】一実装形態による、エネルギー較正を生成する方法及びエネルギー測定値を補正するためにエネルギー較正を使用する方法の流れ図の例を示す図。
【図２Ａ】一実装形態による、陽電子放射トモグラフィ（ＰＥＴ）スキャナの透視図。
【図２Ｂ】一実装形態による、ＰＥＴシステムの概略図。
【図３Ａ】一実装形態による、単結晶及び光電子増倍管（ＰＭＴ）を有するガンマ線検出器モジュールの図。
【図３Ｂ】一実装形態による、アレイとして構成され、光検出器としてＰＭＴを使用する複数のシンチレータ結晶要素を有するガンマ線検出器モジュールの図。
【図３Ｃ】一実装形態による、アレイとして構成され、光検出器としてシリコン光電子増倍管（ＳｉＰＭ）を使用する複数のシンチレータ結晶を有するガンマ線検出器モジュールの図。
【図４Ａ】一実装形態による、ＰＭＴを使用するガンマ線検出器（ＧＲＤ）モジュールの上面図。
【図４Ｂ】一実装形態による、ＰＭＴ−ＧＲＤモジュールの側面図。
【図４Ｃ】一実装形態による、一連のシンチレータ結晶アレイの斜視図。
【図４Ｄ】一実装形態による、アンガーロジックを使用して計算された位置に応じたカウントのフラッドマップ較正の上面図。
【図５】本開示の幾つかの態様による、単一の１次ガンマ線から２つのエネルギー信号を生じるコンプトン散乱の一例を示す図。
【図６Ａ】一実装形態による、閾値超過時間（ＴＯＴ：time-over-threshold）測定のプロットを示す図。
【図６Ｂ】一実装形態による、所与のパルス形状について、閾値に対する比率ピーク高さに応じたＴＯＴのプロットを示す図。
【図７Ａ】シングルチャネル検出イベントについて、ピクセル間で信号共有がない場合のピクセル（すなわち、シンチレータ結晶）間のエネルギー分布の一例を示す図。
【図７Ｂ】シングルチャネル検出イベントについて、ピクセル間で５％の光共有がある場合のピクセル間のエネルギー分布の一例を示す図。
【図８Ａ】マルチチャネル検出イベントについて、ピクセル間で信号共有がない場合のピクセル（すなわち、シンチレータ結晶）間のエネルギー分布の一例を示す図。
【図８Ｂ】マルチャネル検出イベントについて、ピクセル間で５％の光共有がある場合のピクセル間のエネルギー分布の一例を示す図。
【図９】一実装形態による、互いに隣接する第１のシンチレータ結晶（７８）において測定されたＴＯＴ値と第２のシンチレータ結晶（６６）において測定されたＴＯＴ値とに応じたマルチチャネルイベントについてのカウント密度のプロットを示す図。
【図１０Ａ】シングルチャネル検出イベント（１ヒット）、２チャネル検出イベント（２ヒット）、３チャネル検出イベント（３ヒット）、及び４チャネル検出イベント（４ヒット）についての、エネルギー補正なしの合計エネルギーのヒストグラムのプロットを示す図。
【図１０Ｂ】シングルチャネル検出イベント（１ヒット）、２チャネル検出イベント（２ヒット）、３チャネル検出イベント（３ヒット）、及び４チャネル検出イベント（４ヒット）についての、エネルギー補正ありの合計エネルギーのヒストグラムのプロットを示す図。
【発明を実施するための形態】
【００１０】
陽電子放射トモグラフィ（ＰＥＴ）及び単光子放出計算機トモグラフィ（ＳＰＥＣＴ）のための画像再構成は、位置、時間及び／又はエネルギー情報を使用して実行される。エネルギー測定値は、様々な実施上の考慮すべき点により理想的な線形応答から逸脱する。例えば、非線形性は、エネルギーを表すために代わりの量として閾値超過時間を使用することの結果として測定プロセス中に起こり得る。更に、ガンマ線が、単一のガンマ線について２つ（以上）の信号を生じるコンプトン散乱を受けるときなどに発生するガンマ線のマルチチャネル検出は、単一のガンマの２つ（以上）の信号のエネルギー値を合計することによって幾分かは補正される。とはいえ、光共有及び／又は電荷共有により、隣接チャネルがマルチチャネルイベントに関与するとき、合計エネルギーは、ガンマ線がシングルチャネルイベントとして（例えば、コンプトン散乱なしで）検出された場合に測定されるエネルギーから逸脱する。そこで、本実施形態に係る方法及び装置は、実質的に線形な補正エネルギー値を生成する。
【００１１】
大部分の商用のＰＥＴシステムは、陽電子と電子との対消滅によって生成される５１１ｋｅＶのエネルギーを有するガンマ線を検出するためにシンチレータ結晶のアレイを使用する。ＰＥＴシステムは、エネルギー及び位置情報を与えることができる。ガンマ線のエネルギーは、シンチレータ結晶間での散乱（例えば、コンプトン散乱）により２つのシンチレータ結晶によって共有される場合がある。単一のガンマ線からのエネルギーが２つ以上のシンチレータ結晶又はチャネルにおいて検出されるとき、検出イベントは「マルチチャネル検出」と呼ばれる。散乱イベントの合計エネルギーは、全ての関係するシンチレータ結晶からの個々のエネルギーを加算することによって復元される（各シンチレータ結晶とそれらの対応する光検出器、増幅器及び電子回路は「チャネル」とも呼ばれる）。加算は、読出し及び電子回路の選定に応じてアナログ又はデジタル方法のいずれかを通して実現される。
【００１２】
シンチレータ結晶中に残留されるエネルギーは電気信号に変換され、当該電気信号はデジタル化される。このデジタル化プロセスは、様々な方法によって実行され得る。ガンマ線のエネルギー測定値をデジタル化するための方法として、閾値超過時間（ＴＯＴ：time-over-threshold）法がある。ＴＯＴ法は、コスト面で優れており、高いチャネル密度を必要とする実施例に容易に適用される。ＴＯＴ値は所与のチャネルに吸収される実際のエネルギーに単調に関係するが、ＴＯＴと実際のエネルギーとの関係性は完全な線形から逸脱することがある。本実施形態においては、このような線形性からのずれが補正される。
【００１３】
本明細書で使用する「エネルギー」という用語は、実際の又は真のエネルギーに直接的に関係する較正エネルギーを意味するように限定されない。概して、本明細書で使用する「エネルギー」という用語は、実際の又は真のエネルギーを表し、それに単調に関係するエネルギー座標を指定する。従って、「エネルギー」という用語は、コンテキストが別段に明確に示さない限り必ずしも実際の又は真のエネルギーを指すとは限らない。例えば、エネルギーの合計について本明細書で説明するとき、この合計は、実際のエネルギーに直接的に関係する較正値ではなく「エネルギー座標」に対して実行される。測定された／未加工のエネルギーＥraw（すなわち、「エネルギー座標」）の関係が、非線形関数Ｅraw＝ｆ（Ｅtrue）によって真のエネルギーＥtrueに関係するので、２チャネル検出からの２つの測定されたエネルギーｆ（Ｅ1）及びｆ（Ｅ2）の和（すなわち、Ｅ1＋Ｅ2＝ＥTotal、ＥTotalは入射の真のエネルギー、例えば、５１１ｋｅＶ）は、等価なシングルチャネル検出のための測定された／未加工のエネルギーと等しくならない、すなわち、ｆ（Ｅ1）＋ｆ（Ｅ2）≠ｆ（ＥTotal）＝ｆ（Ｅ1＋Ｅ2）である。従って、信号チャネル検出とマルチチャネル検出のエネルギーを正確に比較するために、エネルギー較正及び補正が、マルチチャネル検出の合計されたエネルギーに適用される。
【００１４】
読出しチャネルの数が極めて高いとき、超高回数レートでより良い性能が達成され、不感時間及びパイルアップ効果を低減し得る。一例は、シンチレータ結晶の１対１の読出しであり、各シンチレータ結晶が単一の光検出器に直接結合される。この場合、シンチレータ結晶は、シンチレーション光子に変換されたガンマ線エネルギーがガンマ線エネルギーを吸収したシンチレータ結晶内に維持されるように光学的に隔離される。光隔離は不完全であり、シンチレータ結晶間の最良の遮光材料の場合でも、隣接する結晶間に少量の光クロストークが生じてしまう。従って、本実施形態に係る装置及び方法は、チャネル間のクロストークを示すマルチチャネル検出中に測定されたエネルギーの差を考慮するエネルギー較正を与える。
【００１５】
同一の参照番号は、幾つかの図全体に亘り同じ又は対応する部分を指定するものとする。図１は、ＰＥＴデータ１０５中のエネルギー測定値を補正するための方法１００の流れ図の非限定的な例を示し、ＰＥＴデータ１０５中のエネルギー測定値を補正するために方法１００において使用されるべきエネルギー較正を決定するための方法１６０の流れ図の非限定的な例を示す。図２Ａ及び図２Ｂに、方法１００及び１６０を実行するために使用され得るＰＥＴスキャナの非限定的な例を示す。
【００１６】
図１に、ＰＥＴデータ１０５を補正するためにマルチチャネルイベントを考慮するためにＰＥＴデータのエネルギー較正を実行し、較正されたＰＥＴデータからＰＥＴ画像１５５を再構成する方法１００及び１６０のフローチャートを示す。方法１００及び１６０は、図２Ａ及び図２Ｂに示すＰＥＴシステム２００により実行される。
【００１７】
ステップ１１０においてプロセッサ２７０は、ガンマ線検出器からＰＥＴデータ１０５を取得し、前処理する。具体的には、プロセッサ２７０は、複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すＰＥＴデータ１０５を複数の検出イベント毎に取得する。検出イベントは、より詳細には、対消滅により発生した一対のガンマ線の同時計数イベントである。プロセッサ２７０は、ガンマ線検出器からデータ収集システム２７６を介してＰＥＴデータ１０５を取得する。本実施形態において検出器要素は、チャネルと同義である。検出イベントは、複数の検出器要素及びモジュールにおいて検出される。検出器モジュールの各々は、２つ以上のチャネル（すなわち、検出器要素）を含む。
【００１８】
ステップ１２０においてプロセッサ２７０は、ＰＥＴデータ１０５に基づいて、複数の検出イベント各々を、ヒット数に応じたグループに分類する。ヒット数は、シンチレータ結晶内でのガンマ線の相互作用の回数を意味する。相互作用は、シンチレータ結晶によるガンマ線の吸収を意味し、典型的には、光電吸収及びコンプトン散乱であり、更に電子対生成が加味されてもよい。グループは、シングルチャネルイベントに対応するグループとマルチチャネルイベントに対応するグループとに大別される。シングルチャネルイベントは、１ヒットに対応し、１次ガンマ線の検出イベントである。ヒット数は、本実施形態において多重度に近似することとなる。多重度は、一の入射ガンマ線を検出した光検出器又はチャネルの個数である。
【００１９】
まずプロセッサ２７０は、ＰＥＴデータ１０５に基づいて、複数の検出イベントの中からマルチチャネルイベントを選択する。次にプロセッサ２７０は、ＰＥＴデータ１０５から、当該マルチチャネル検出イベントを構成する複数の検出イベントのＰＥＴデータをグループ化する。すなわち、マルチチャネルイベントは、ヒット数に応じた、更に細かなグループに分類される。マルチチャネルイベントは、２ヒット以上に対応し、例えば、１回以上のコンプトン効果及び／又は電子対生成と１回の光電効果とを意味する。マルチチャネルイベントのグループは、１次散乱に対応するグループ、２次散乱に対応するグループ、更に高次の散乱に対応するグループ等がある。１次散乱に対応するグループは、２ヒット、すなわち、１次ガンマ線の一部を吸収し散乱ガンマ線として残りのエネルギーを放出する第１のシンチレータ結晶のヒットと散乱ガンマ線を吸収した第２のシンチレータ結晶のヒットとを含む。同様に、２次散乱に対応するグループは、３ヒット、すなわち、１次ガンマ線についての１ヒットと２つの散乱ガンマ線についての２ヒット）を含む。より高次散乱に対応するグループ、例えば、４ヒットを含む３次散乱に対応するグループも同様に定義可能である。マルチチャネルイベントは、例えば、検出信号の時間的及び空間的な近接度に基づいて、信号のエネルギーの和、及び合計されたエネルギーと検出信号の時間的及び空間的な近接度との任意の組合せに基づいて判別される。例えば、ガンマ線源が周知のエネルギー（例えば、陽電子消滅からのガンマ線について５１１ｋｅＶ）を有する場合、信号の合計が既知のエネルギーに近くなればなるほど、信号が同じマルチチャネルイベントに対応する可能性が高くなる。更に、時間的に互いにより近接して発生した信号は、同じマルチチャネルイベントに対応する可能性が高く、空間的に互いにより近接して発生した信号は、同じマルチチャネルイベントに対応する可能性が高くなる。更に、上記の３つの状態（すなわち、エネルギー、時間及び空間）の全てが充足されるとき、信号は、同じマルチチャネルイベントに対応する可能性が一層高くなる。プロセッサ２７０は、ＰＥＴデータに対する多変量統計解析を使用して、ＰＥＴデータからマルチチャネルイベントを選択する。
【００２０】
プロセッサ２７０は、例えば、２次以上のコンプトン散乱が発生したとき、特定の検出器モジュールの２つ以上のチャネルが３つ以上のチャネルを含むのか否かを決定する。幾つかの実装形態では、１次散乱のみがタイミング較正のために使用され、高次散乱についてのマルチチャネルイベントは廃棄される。他の実装形態では、１次散乱と高次散乱との両方のためのマルチチャネルイベントがタイミング較正のために使用される。従って、プロセッサ２７０は、対応する２チャネルイベントに２つ以上のチャネルのデータを分解するか、あるいはマルチチャネルイベントを完全に廃棄してもよい。
【００２１】
ステップ１３０においてプロセッサ２７０は、マルチチャネルイベントについて位置データを決定される。本明細書では、「位置データ」又はより単純に「位置」という用語は、３次元空間における絶対位置、所与の検出器モジュール内の２次元座標を指すか、又は検出器アレイモジュール内の個々のシンチレータ結晶を一意に識別するインデックスなどの結晶識別情報（ＩＤ）などを指す。例えば、結晶ＩＤは、３次元のユークリッド空間における絶対位置ではないが、検出イベントが発生した物理的位置を示す。従って、「位置データ」及び「位置」は、直線で囲まれた空間に対応するように較正された絶対位置を意味することに限定されず、むしろ、「位置データ」及び「位置」という用語は、別個の結晶要素を空間的に指定／識別する座標である。
【００２２】
１次散乱の場合において、検出器要素における２ヒットのグループはマルチチャネルイベントであるものとして選択される。ヒットの各々は、位置（ｘ）及びエネルギー（Ｅ）及び幾つかの実装形態では時間（ｔ）の量を含む検出器チャネルの検出信号である。従って、２つのヒットがある場合、第１及び第２の時間（すなわち、ｔ1及びｔ2）、第１及び第２の位置（すなわち、ｘ1及びｘ2）、第１及び第２のエネルギー（すなわち、Ｅ1及びＥ2）があることになる。マルチチャネルイベントの位置を決定するために任意の方法が使用される。例えば、エネルギー重み付けされた２次元座標が、例えば、次式を使用して各マルチチャネルイベントについて計算される。
【００２３】
【数１】
2020064045.tif 000003
【００２４】
ここで、ｘｉ→はチャネルｉの中央座標であり、Ｅiはチャネルｉのエネルギーである。例えば、Ｅiは、非線形較正されたエネルギーでもよいし、未加工のエネルギー値でもよい。ｗは、エネルギー重みＥｉＷを得るためにエネルギーＥiに適用される冪指数（exponent又はpower）であり、ｗは、０を除く任意の値である。後述するように、ｉ番目のチャネルのための中央座標ｘｉ→は、結晶識別情報（ＩＤ）に基づいて取得される。
【００２５】
ステップ１４０においてプロセッサ２７０は、補正されたＰＥＴデータを生成するために、エネルギー較正１８５をＰＥＴデータに適用する。エネルギー較正１８５は、エネルギーに関して閾値超過時間（ＴＯＴ）の値を線形化する補正、シングルチャネルイベントとマルチチャネルイベントとの間のエネルギー差の補正、又はこれらの補正の組合せを含むことができる。
【００２６】
幾つかの実装形態では、ＰＥＴデータ１０５のエネルギー値は、画像再構成のために、予め定義されたウィンドウ（例えば、５１１ｋｅＶ周辺のエネルギーウィンドウ）内にエネルギーを有する同時計数カウントのみを選択するために使用される。この場合、マルチチャネルエネルギーのスケーリングを適用し、シフト／スケーリングされたエネルギーに、５１１ｋｅＶ周辺のウィンドウを適用するのではなく、逆スケーリングがウィンドウに適用され、シフトされたウィンドウが、シフトされていないマルチチャネルエネルギーに適用される。このようにして、ウィンドウは、スケーリングのみを必要とし、所与のブロックについてシフトされたウィンドウ及び散乱／多重度の順序は、そのブロック及び多重度中のマルチチャネルイベントの各々に使用される。異なる較正は、異なるブロックに、及び所与のブロック内の異なるチャネル多様性に適用することができる。多重度は、マルチチャネル検出に関与するチャンネルの数を指す。例えば、１次コンプトン散乱は、２の多重度に対応する２チャネルで検出される。２次コンプトン散乱は、３の多重度に対応する３チャネルで検出され、以下同様である。
【００２７】
異なるエネルギー較正又はシフトが、各ブロックに、及び所与のブロック内の各多重度に適用される。また、これらの各エネルギー較正又はシフトは、ルックアップテーブル中に記憶される。多くのブロック及び多重度がありえても、所与のブロック中の多重度ごとに多くのマルチチャネルイベントがある。従って、計算の数は、所与の多重度及びブロックに対応する多くのマルチチャネルイベントの各々にエネルギーシフトを適用するのではなく、所与の多重度及びブロックについてのウィンドウにエネルギーシフトを適用することによって低減される。
【００２８】
ステップ１５０においてプロセッサ２７０は、任意の既知の再構成方法を使用して補正後のＰＥＴデータからＰＥＴ画像１５５を再構成する。例えば、ＰＥＴデータは、ボクセル位置に応じた放射能レベル（例えば、放射線トレーサ剤密度）の画像を再構成するために使用される。画像再構成は、逆投影法、フィルタ補正逆投影法、フーリエ変換ベースの画像再構成法、逐次画像再構成法、行列反転画像再構成法、統計画像再構成法、リストモード法、又は他の再構成法、あるいはそれらの組合せを使用して実行され得る。例えば、初期ＰＥＴ画像は、ＦＢＰ再構成されたＰＥＴ画像で初期化された順序付きサブセット期待値最大化（ＯＳ−ＥＭ：ordered subset expectation maximization）アルゴリズムを使用して再構成される。
【００２９】
較正段階においてプロセッサ２７０は、方法１６０の実行により、較正データ１６５からエネルギー較正１８５を生成する。較正データ１６５は、ＰＥＴシステム２００のガンマ線検出器の検出器要素により生成されたＰＥＴデータである。較正データ１６５は、陽電子放射線源からの同時計数カウントとして生成されるが、必ずしもそれから生成される必要はない。代替的に、較正データ１６５は、異なるガンマ線エネルギーを生成する不対のガンマ線（「不対放射線源」と呼ばれる）の１つ又は複数の陽電子放射線源を使用して生成されてもよい。プロセッサ２７０は、較正段階において、ガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表す較正データを複数の較正検出イベント毎に取得する。
【００３０】
例えば、較正データ１６５は、約３０年の半減期を有し、６６２ｋｅＶのエネルギーをもつガンマ線を生成するセシウム同位体１３７（Ｃｓ−１３７）の１つ又は複数の部分を有する不対放射線源を使用して生成される。Ｃｓ−１３７の当該部分は、ビーム内でほぼ均一なガンマ線の束密度を実現するように構成される。幾つかの実装形態では、ビームは、コーンビームでもよいし、又は全球である４πステラジアンに同位体的に放射されてもよい。Ｃｓ−１３７によって生成された６６２ｋｅＶのエネルギーを有するガンマ線は、陽電子放射の放射線源によって生成された５１１ｋｅＶのガンマ線よりも高い。複数の部分の代わりに、実質的に均一に分布された放射線源が使用されてもよい（一般に「フラッドソース」として知られる）。使用される他のガンマ線源は、例えば、以下を含む。
【００３１】
（ｉ）コバルト同位体６０（５．３年の半減期をもち、１．１７ＭｅＶ及び１．３３ＭｅＶのガンマ線エネルギーをもつＣｏ−６０）
（ｉｉ）ゲルマニウム同位体６８（０．７４年の半減期をもち、５１１ｋｅＶのガンマ線エネルギーをもつＧｅ−６８）
（ｉｉｉ）ナトリウム同位体２２（２．６年の半減期をもち、５１１ｋｅＶ及び１．２７５ＭｅＶのガンマ線エネルギーをもつＮａ−２２）
所与の放射線源によって生成されたガンマ線のエネルギーが、陽電子放射によって生成された５１１ｋｅＶのエネルギーよりも大きいとき、これらのより高いガンマ線エネルギーは、５１１ｋｅＶの大小双方を含む大きいエネルギー範囲をカバーするコンプトン散乱において様々な検出器要素に残留エネルギーを生じる。これは、ＴＯＴ値からエネルギーへのマッピングの較正をより高精度に行うことが可能になる。
【００３２】
ステップ１７０においてプロセッサ２７０は、較正データに基づいて、複数の較正検出イベントの中から較正マルチチャネルイベントを選択する。すなわち、プロセッサ２７０は、まず、較正データ１６５を、マルチチャネルイベントとシングルチャネルイベントとに分類する。較正データ１６５の総セットからのマルチチャネルイベントのこのソート及びフィルタ処理は、方法１００のステップ１２０において説明した方法と同様の方法を使用して実行される。
【００３３】
ステップ１８０においてプロセッサ２７０は、較正データ１６５に基づいてエネルギー較正１８５を生成する。上記で説明したように、エネルギー較正１８５は、ＴＯＴ非線形性補正、マルチチャネル補正又はその両方を含む。例えば、エネルギー較正１８５は、非線形補正を表す式のパラメータを取得するために各検出器要素の位置又は識別情報（ＩＤ）によってインデックス付けされたルックアップテーブルである。従って、エネルギー較正のパラメータ化は、検出器要素ベースによって検出器要素に対して実行される。例えば、以下で説明するように、電荷共有によるマルチチャネル検出に関して発生したエネルギーシフトΔＥは、次式によってパラメータ化される。
【００３４】
【数２】
2020064045.tif 000004
【００３５】
ここで、ｘは、シンチレータ結晶１中に残されたエネルギーの割合であり、δ1及びδ2は、隣接するシンチレータ結晶との２つの結晶の光クロストークである。このパラメータ化では、エネルギーシフトは、光クロストークのレベルに比例する。幾つかのシンチレータ結晶アレイでは、δ1及びδ2はまた、共有の方向に依存し得る（例えば、光共有は非対称であり得る）。従って、クロストークのレベルがシンチレータ結晶ごとに変動するとき、エネルギーシフトΔＥは、機能エネルギー共有として変化する。幾つかの実装形態では、光クロストークのレベル及びエネルギーシフトの量ΔＥは、以下で説明するように、検出器モジュール内のシンチレータ結晶の位置に依存し得る。
【００３６】
次に、環形状に構成された検出器モジュール（すなわち、ガンマ線検出器（ＧＲＤ））で構成されたＰＥＴシステム２００の非限定的な例の説明を与える。検出器モジュールの各々は、検出器要素の幾つかのアレイを含むことができる。ＧＲＤは、（例えば、光、赤外、及び紫外波長で）ガンマ線をシンチレーション光子に変換するためのシンチレータ結晶アレイを含み、シンチレーション光子は、光検出器によって検出される。図２Ａ及び図２Ｂに示す非限定的な例では、光検出器は、各シンチレータ結晶よりも大きい光電子増倍管（ＰＭＴ）である。好ましい一実施形態では、光検出器は、個々のシンチレータ結晶の断面積を近似する検出断面積を有することができるシリコン光電子増倍管（ＳｉＰＭｓ）であり、シンチレータ結晶と光検出器との間で１対１の対応を形成する。この実施形態を、図３Ｃ及び図５に示す非限定的な例によって示す。光検出器がシンチレータ結晶よりも大きく、従って、複数のシンチレータ結晶からの光信号を検出するために単一の光検出器が使用される場合、位置を決定するためにアンガーロジックが使用される。但し、シンチレータ結晶と光検出器との間に１対１の対応があるとき、アンガーロジックが必ずしも必要とされない。
【００３７】
図２Ａ及び図２Ｂは、方法１００及び１６０を実装することができるＰＥＴシステム２００の非限定的な例を示す。ＰＥＴシステム２００は、矩形の検出器モジュールを構成する複数のガンマ線検出器（ＧＲＤ）（例えば、ＧＲＤ１、ＧＲＤ２〜ＧＲＤＮ）を含む。例えば、検出器リングは、４０個のガンマ線検出器ＧＲＤを含む。なお、検出器リングが有するガンマ線検出器ＧＲＤの個数はこれに限定されず、４０個より多くても少なくてもよい。例えば、検出器リングは、４８個のガンマ線検出器ＧＲＤを有してもよいし、ＰＥＴシステム２００のよりボア径を大きくするため、より多くの数のガンマ線検出器ＧＲＤを有してもよい。
【００３８】
各ガンマ線検出器は、ガンマ放射線を吸収し、シンチレーション光子を放出する個々の検出器結晶の２次元アレイを含む。シンチレーション光子は、同じくガンマ線検出器に構成される光電子増倍管（ＰＭＴ）の２次元アレイによって検出される。検出器結晶のアレイとＰＭＴとの間に光ガイドが配設される。
【００３９】
代替的に、シンチレーション光子は、シリコン光電子増倍管（ＳｉＰＭ）アレイによって検出されてもよい。この場合、検出器結晶各々はＳｉＰＭを有することとなる。
【００４０】
各光検出器（例えば、ＰＭＴ又はＳｉＰＭ）は、シンチレーションイベントがいつ発生したのかと検出イベントを生成するガンマ線のエネルギーとを示すアナログ信号を生成する。更に、１つの検出器結晶から放出された光子が２つ以上の光検出器によって検出され、各光検出器において生成されたアナログ信号に基づいて、検出イベントに対応する検出器結晶は、例えば、アンガーロジック及び結晶復号を使用して決定される。
【００４１】
図２Ｂは、対象物ＯＢＪから放出されたガンマ線を検出するように構成されたガンマ線光子計数検出器（ＧＲＤ）を有するＰＥＴシステムの概略図を示す。ＧＲＤは、各ガンマ線検出に対応するタイミング、位置及びエネルギーを測定することができる。一実装形態では、ガンマ線検出器は、図２Ａ及び図２Ｂに示すように、リングに構成される。検出器結晶は、２次元アレイに構成された個々のシンチレータ材を有するシンチレータ結晶である。シンチレータ材は、任意の材料でよい。ＰＭＴは、シンチレーションイベントのアンガーロジック及び結晶復号を可能にするために各シンチレータ結晶からの光が複数のＰＭＴによって検出されるように構成される。
【００４２】
図２Ｂに、ＰＥＴシステム２００の構成の一例を示し、ここで、撮像されるべき対象物ＯＢＪが、天板２１６に載置される。ガンマ線検出器モジュールＧＲＤ１〜ＧＲＤＮは、対象物ＯＢＪ及び天板２１６を囲むように円周に沿って配置される。ガンマ線検出器は、架台２４０に固定的に接続されている円形構成要素２２０に固定的に接続され得る。架台２４０は、対象物ＯＢＪと天板２１６とが通ることができる開口を含む。消滅イベントにより対象物ＯＢＪから互いに反対方向に放出された一対のガンマ線は、ガンマ線検出器によって検出される。タイミング及びエネルギー情報は、ガンマ線ペアの同時計数を決定するために使用される。
【００４３】
図２Ｂは、ＰＥＴシステム２００に含まれる、ガンマ線検出データを取得、記憶、処理及び分配するための回路及びハードウェアを示す。回路及びハードウェアは、プロセッサ２７０、ネットワークコントローラ２７４、メモリ２７８及びデータ収集システム（ＤＡＳ）２７６を含む。ＰＥＴシステム２００はまた、ＤＡＳ２７６、プロセッサ２７０、メモリ２７８及びネットワークコントローラ２７４にガンマ線検出器により生成された検出測定値をルーティングするデータチャネルを含む。データ収集システム２７６は、ガンマ線検出器からの検出データの収集、デジタル化及びルーティングを制御することができる。一実装形態においてＤＡＳ２７６は、寝台による天板２１６の移動を制御する。プロセッサ２７０は、本明細書で説明するように、検出データから画像を再構成すること、検出データの前再構成処理及び画像データの後再構成処理を含む機能を実行する。
【００４４】
プロセッサ２７０は、本明細書で説明する方法１００及び１６０の様々なステップ及びそれらの変形形態を実行するように構成される。プロセッサ２７０は、個別論理ゲート、特定用途向け集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）又は他の複雑なプログラマブル論理デバイス（ＣＰＬＤ）として実装されるＣＰＵを含むことができる。ＦＰＧＡ又はＣＰＬＤ実装形態は、ＶＨＤＬ、Ｖｅｒｉｌｏｇ、又は任意の他のハードウェア記述言語でコーディングされ、コードは、直接ＦＰＧＡ又はＣＰＬＤ内で又は別個の電子メモリとして電子メモリ中に記憶される。更に、メモリは、ＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ（登録商標）又はフラッシュメモリなどの不揮発性でもよい。メモリはまた、スタティック又はダイナミックＲＡＭなどの揮発性であり、マイクロコントローラ又はマイクロプロセッサなどのプロセッサは、電子メモリならびにＦＰＧＡ又はＣＰＬＤとメモリとの間の相互作用を管理するために与えられてもよい。
【００４５】
代替的に、プロセッサ２７０中のＣＰＵは、方法１００及び１６０の様々なステップを実行するコンピュータ可読命令のセットを含むコンピュータプログラムを実行することができ、プログラムは、上記で説明した非一時的電子メモリ及び／あるいはハードディスクドライブ、ＣＤ、ＤＶＤ、フラッシュドライブ又は任意の他の既知の記憶媒体のいずれかに記憶される。更に、コンピュータ可読命令は、米国のＩｎｔｅｌからのＸｅｏｎプロセッサ又は米国のＡＭＤからのＯｐｔｅｒｏｎプロセッサなどのプロセッサ、ならびにマイクロソフト  ＶＩＳＴＡ、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ、ＬＩＮＵＸ（登録商標）、Ａｐｐｌｅ、ＭＡＣ−ＯＳ及び当業者に既知の他のオペレーティングシステムなどのオペレーティングシステムと併せて実行するユーティリティアプリケーション、バックグラウンドデーモン、又はオペレーティングシステムの構成要素、あるいはそれらの組合せとして与えられる。更に、ＣＰＵは、命令を実行するために並行して協働的に動作する複数のプロセッサとして実装される。
【００４６】
メモリ２７８は、ハードディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、フラッシュドライブ、ＲＡＭ、ＲＯＭ又は当技術分野において既知の任意の他の電子ストレージでよい。
【００４７】
米国のＩｎｔｅｌ  ＣｏｒｐｏｒａｔｉｏｎからのＩｎｔｅｌイーサネット（登録商標）ＰＲＯネットワークインターフェースカードなどのネットワークコントローラ２７４は、ＰＥＴイメージャの様々な部分の間をインターフェースすることができる。更に、ネットワークコントローラ２７４はまた、外部ネットワークとインターフェースすることができる。諒解され得るように、外部ネットワークは、インターネットなどの公衆ネットワークあるいはＬＡＮ又はＷＡＮネットワークなどのプライベートネットワーク、又はそれらの任意の組合せでよく、ＰＳＴＮ又はＩＳＤＮサブネットワークをも含んでもよい。外部ネットワークはまた、イーサネット（登録商標）ネットワークなどのワイヤードであるか、又はＥＤＧＥ、３Ｇ及び４Ｇワイヤレスセルラーシステムを含むセルラーネットワークなどのワイヤレスでもよい。ワイヤレスネットワークはまた、ＷｉＦｉ、Ｂｌｕｅｔｏｏｔｈ（登録商標）、又は既知の通信の任意の他のワイヤレス形態でもよい。
【００４８】
図３Ａ、図３Ｂ及び図３Ｃは、ガンマ線検出器の幾つかの非限定的な実装形態を示す。なお、ガンマ線検出器はガンマ線カメラとも呼ばれる。図３Ａでは、シンチレータ結晶は、単一のモノリシックブロックである。ここでは光電子増倍管（ＰＭＴ）として示される光検出器のアレイによって検出される２次光子にガンマ線光子を変換するシンチレーションイベントの位置が図示されている。シンチレーションイベントの位置は、アンガーロジックを使用して決定される。
【００４９】
図３Ｂでは、周期配列された個別のシンチレータ結晶に切り出され、反射材により光学的に分断されたシンチレータ結晶アレイが示される。ブロック中のシンチレータ結晶間の光隔離は、不完全であり、隣接するシンチレータ結晶との間で光共有が生じうる。光検出器がＰＭＴである場合、隣接するシンチレータ結晶間の光共有は、シンチレータ結晶を出た後に発生する光共有と比較して小さい。この場合、ＰＥＴデータにアンガーロジックを使用して、シンチレーションイベントの位置（ガンマ線の入射位置）が近似的に決定される。次いで、計算された入射位置に、フラッドマップに基づき較正されたルックアップテーブルを適用し、シンチレータ結晶アレイを構成する各シンチレータ結晶のインデックスを決定する。
【００５０】
図３Ｃでは、各シンチレータ結晶からの光が、各シリコン光電子増倍管（ＳｉＰＭ）によって検出される。各シンチレータ結晶が光検出器に一対一接続している場合、光検出器間での光共有が低減される。更に、この場合、各シンチレータ結晶は、単一の検出モジュール内の異なるシンチレータ結晶で発生する同時シンチレーションイベント間の識別を可能にすること（例えば、隣接するシンチレータ結晶間でコンプトン散乱を識別すること）によって分解能を拡張し得る。
【００５１】
長年にわたり、ＰＥＴについて最も一般に使用される光検出器はＰＭＴである。ＰＭＴは、入射光子を光電子に変換するための仕事関数を有する光電陰極材料を包含する真空管である。光電子は、印加された電場により加速され、ダイノードのカスケードと相互作用することによって増幅される。得られた電流は、初期シンチレーション光子の数、すなわち、ガンマ線によってシンチレーション結晶内に残留するエネルギーに比例する。
【００５２】
例えば、多くの小型のＳｉＰＭ等の位置高感度型のＰＭＴの性質を活用して、例えば、アンガーロジックを使用してシンチレータブロックをセグメント化することによって、光子検出の位置が決定される。図３Ｂ及び図３Ｃでは、サイズが数ミリメートルである小さい個々のシンチレータ結晶がブロック又はモジュールに稠密に配置されている。各シンチレータ結晶は光検出器に結合される。アンガーロジックでは、対消滅光子の相互作用位置は、光検出器からの相対信号を比較することによって決定される。光検出器からの相対信号は、分散したシンチレーション光子信号に起因する相対的な光共有によって決定される。次いで、決定された相互作用位置に基づいて、フラッドマップに基づき較正されたルックアップテーブルを利用して、ガンマ線が入射したシンチレータ結晶が決定される。
【００５３】
図４Ａ及び図４Ｂは、光検出器としてＰＭＴを使用するガンマ線検出器モジュールＧＲＤ１の上面図及び側面図をそれぞれ示す。図４Ａ及び図４Ｂに示すガンマ線検出器は、２つの異なるサイズのＰＭＴ（２５ｍｍのＰＭＴと３８ｍｍのＰＭＴ）を含む。２つの異なるサイズのＰＭＴを含むことにより、全て同じサイズのＰＭＴを含む場合に比して、シンチレータ結晶アレイをカバーする面積を広げることができる。アレイ１及びアレイ５は、１１画素×１６画素であり、アレイ２、アレイ３及びアレイ４は、８画素×１６画素である。図４Ｃは、５つの結晶アレイの斜視図を示し、図４Ｄは、アンガーロジックを使用して計算された概略位置に応じたカウントのヒストグラムのフラッドマップ較正を示す。アンガーロジックにより概略位置が算出されるので、ヒストグラムの極大値は、各シンチレータ結晶の中心に対応するものとして識別される。次いでヒストグラムの極大値は、位置情報（例えば、チャネルｉの中央座標のｘｉ→）を決定するときに使用されるルックアップテーブルにおける位置較正として記憶される。
【００５４】
図５は、第１のシンチレータ結晶中でコンプトン散乱を受け、その結果、残りのエネルギーが第２のシンチレータ結晶に残留される入射ガンマ線の非限定的な例を示す。第１及び第２のシンチレータ結晶は、次いで、それぞれ光電吸収を介して２次光子を生成し、その結果、対応するＳｉＰＭが、第１及び第２のエネルギー信号を生成する。散乱角度は小さく、第１のエネルギー信号中のエネルギーは、第２のエネルギー信号中のエネルギーよりも少ない。第１及び第２のパルスの形状が同じで大きさが異なる場合でも、図５に示すように、パルスの立ち上りエッジが所定値を超えた事に基づく到着時間に生じるタイムウォークオフセットに起因して、２つの信号が異なる時間に発生したと登録される。
【００５５】
図６Ａ及び図６Ｂは、ステップ１８０において較正され、ステップ１４０において適用されるＴＯＴ非線形補正についてのプロットを示す。図６Ａは、時間経過に応じて変化する、ガンマ線検出からのパルスの電圧のプロットを示す。図６Ａには、約１．２ミリボルトにおける予め定義された閾値が示される。パルスがこの閾値を超える持続時間がＴＯＴ値である。図６Ｂに示すように、ＴＯＴ値は、検出されたガンマ線のエネルギーに単調に関係する。この関係は、パルス形状がパルスの大きさの変化に関わらず歪まないとき、パルスの曲線下面積（Area Under the Curve）又はパルス振幅により表される。閾値を下回る信号では、ヒットは登録されない。
【００５６】
図６Ｂのプロットをモデル化する関数は、ステップ１８０において任意の方法に基づいて較正データから生成される。例えば、シングルチャネルイベントのみを使用して、ガンマ線について既知のエネルギー値を有する異なるエネルギー源が使用され、当該異なるエネルギー値に対応するＴＯＴ値が決定される。関数（例えば、シフトされた平方根）がこれらのポイントにカーブフィッティングされる。また、コンプトン散乱による散乱ガンマ線の測定データを使用して、より低いガンマ線エネルギーを補われる。それらは、入射ガンマ線よりも少ないエネルギーを有する。コンプトン散乱検出のエネルギーは、曲線下面積又はパルス振幅に基づいて決定される。この較正が、一度実行され、メモリに記憶され、次いで、ＰＥＴスキャンが実行又は処理されているときにメモリから呼び戻される。上記の較正の方法は、非限定的な例であり、ＴＯＴ値からエネルギーへのマッピングを較正する他の方法が、本明細書で説明する方法及び装置の趣旨から逸脱することなく使用されてもよい。
【００５７】
閾値超過時間（ＴＯＴ）法は、所与の閾値を上回る信号の総時間を測定することによって信号の振幅を推定する。測定された時間は、次いで、非線形変換（例えば、図６Ｂに示すマッピング関数）を通して振幅に変換して戻される。閾値の存在により、閾値を下回る全ての振幅は、０、すなわち、「ヒットなし」として判別される。ＴＯＴ法が、無視できない光クロストークを有する検出器アレイに適用された場合、入射ガンマ線のエネルギーが散乱により２つのシンチレータ結晶により共有されるとき、測定された合計エネルギーはシフトすることになる。
【００５８】
図７Ａは、シンチレータ結晶アレイ中の画素（すなわち、シンチレータ結晶）間にクロストークが発生しないシングルチャネル検出イベントの理想的な場合を示す。すなわち、検出されたエネルギーの１００％がガンマ線を吸収した画素中に残留する。図７Ｂは、シングルチャネルイベントにおける５％のクロストークの場合を示す。すなわち、光共有が５％であるとき、ガンマ線が入射した画素にエネルギーの８０％が残留し、各隣接画素はエネルギーの５％を受ける。
【００５９】
図８Ａ及び図８Ｂは、それぞれ、２ヒットマルチチャネル検出イベント（例えば、１つの１次ガンマ線及び１つの散乱ガンマ線）についてのクロストークなし及び５％のクロストークの理想的な場合を示す。この非限定的な例は、隣接シンチレータ結晶間の５％の光クロストークと１０％のＴＯＴ閾値とを仮定する。図８Ａでは、合計エネルギーは、シングルチャネルの場合と同じである６０＋４０＝１００となる。ただし、図８Ｂでは、合計エネルギーは、シングルチャネルの場合よりも５％多くなる（４８＋２）＋（３２＋３）＝８５となる。マルチチャネル検出では、閾値を上回るチャネルへのクロストークを介して合計エネルギーの３％及び２％が結合されるが、シングルチャネル検出では、クロストークを介してエネルギーが結合されるチャネルの全てが閾値を下回るので、図８Ａと図８Ｂとの比較により現れる差が生じる。すなわち、シングルチャネル検出の場合では、何れのクロストークエネルギーも測定エネルギー（すなわち、閾値を上回るエネルギー）に含まれないが、２チャネル検出の場合では、クロストークエネルギーの５％が測定エネルギーに取り戻される。従って検出器要素が完全な線形応答を有するときでも、クロストークは、シングルチャネル検出の測定エネルギーから逸れているマルチチャネル検出から測定された合計エネルギーに生じる。
【００６０】
クロストークの場合、シングルチャネル及び２チャネル検出について測定された合計エネルギーは、次式によって与えられる。
【００６１】
【数３】
2020064045.tif 000005
【００６２】
ここで、ｆは、第１のシンチレータ結晶（すなわち、１次ガンマ線からのエネルギーを吸収するシンチレータ結晶）に残留されたエネルギーの割合である。δ1は第２のシンチレータ結晶に対する第１のシンチレータ結晶の光クロストークであり、δ2は第１のシンチレータ結晶に対する第２のシンチレータ結晶の光クロストークである。幾つかの実装形態では、異なる光クロストークが隣接するシンチレータ結晶毎に使用され、個々の検出器は、エネルギーに対して完全な線形応答を有するとみなされる。マルチチャネルの場合、次式が得られる。
【００６３】
【数４】
2020064045.tif 000006
【００６４】
ここで、Ｅ0は、クロストークがない場合に測定されるエネルギーである。右側の第１の項Ｅ0ｆ（１−４δ1）は、第１のシンチレータ結晶中のエネルギーからクロストークによる損失を減じたものである。右側の第２の項Ｅ0ｆδ1は、第１のシンチレータ結晶からのクロストークにより第２のシンチレータ結晶により得られたエネルギーである。第３の項Ｅ0（１−ｆ）（１−４δ2）は、第２のシンチレータ結晶中のエネルギーからクロストークによる損失を減じたものである。右側の第４の項Ｅ0（１−ｆ）δ2は、第２のシンチレータ結晶からのクロストークにより第１のシンチレータ結晶により得られたエネルギーである。例えば、図８Ｂでは、ｆ＝０．６であり、δ1＝δ2＝０．０５であり、その結果、第１の項は、Ｅ0ｆ（１−４δ1）＝０．４８Ｅ0となり、第２の項は、Ｅ0ｆδ１＝０．０３Ｅ0となり、第３の項は、Ｅ0（１−ｆ）（１−４δ2）＝０．３２Ｅ0となり、第４の項は、Ｅ0（１−ｆ）δ2＝０．０２Ｅ0となる。上記の式は、２ヒットマルチチャネルの場合、次式に簡略化することができる。
【００６５】
【数５】
2020064045.tif 000007
【００６６】
シングルチャネルの場合、ｆ＝１であり、エネルギーは次式によって与えられる。
【００６７】
【数６】
2020064045.tif 000008
【００６８】
エネルギー差は、次式によって与えられる。
【００６９】
【数７】
2020064045.tif 000009
【００７０】
３ヒットや４ヒット等に対応するマルチチャネルイベントについて同様の式を導出することができる。
【００７１】
マルチチャネルのうちの２つのチャネルが隣接しないとき（例えば、図５に示すように、２つのチャネル間にシンチレータ結晶があるとき）、クロストークエネルギーは、閾値を超えるチャネルに結合するのではない。合計エネルギーはシングルチャネルイベント中に測定されるのと同じである。
【００７２】
上記の説明に鑑みて、所与のエネルギー測定値について、ガンマ線のエネルギーが複数のチャネルの間で共有される場合、較正されていない合計エネルギーはクロストークのレベルに依存することが観察される。例えば、ΔＥ（ｆ，δ1，δ2）についての上記の式は、光クロストークδ1及びδ2のレベルに比例するエネルギーシフトを示す。更に、クロストークがシンチレータ結晶ごとに変動する場合、エネルギーシフトは、エネルギー共有の関数として変化する。従って、エネルギーシフトの量は、空間依存のクロストークパラメータδ1（Ｘ→）及びδ２（Ｘ→）によってパラメータ化される。空間依存のクロストークパラメータδ1（Ｘ→）及びδ２（Ｘ→）は、位置Ｘ→によってインデックス付けされたルックアップテーブルを使用して取得される。すなわち、エネルギー較正１８５は、位置ｘ→における補正ルックアップテーブルを含み、これは、シンチレータ結晶間で共有するエネルギーの大きさを表す。マルチチャネルイベントごとに、エネルギー較正１８５は、ステップ１４０において、マルチチャネルイベントについて位置Ｘ→を計算し、位置Ｘ→を使用して空間依存のクロストークパラメータδ1（Ｘ→）及びδ２（Ｘ→）をルックアップし、次いで、２ヒットマルチチャネルイベントについて上記の補正係数（δ1＋３ｆ（δ2−δ1））等の補正係数を適用してマルチチャネルについて合計エネルギーを補正することによって適用される。
【００７３】
クロストークによる非線形性は、図９に示される。図９は、シンチレータ結晶間（すなわち、結晶ＩＤ６６及び結晶ＩＤ７８）のマルチチャネルイベントの測定カウントをグレイスケールで表す。クロストークがない場合に観察される線形応答に対応する線が図９に示されている。図９は、非線形エネルギー測定値について、ガンマ線のエネルギーが複数のチャネルによって共有される場合、較正されていない合計エネルギーは、エネルギーの分布態様に依存することを示す。
【００７４】
個々のチャネルでのエネルギー測定値の非線形性を近似的に較正するために初期の粗いエネルギー較正が使用されてもよい。例えば、初期エネルギー較正は、ＴＯＴ値とエネルギーとの間の非線形関係を補正することができる。初期の粗いエネルギー較正は、シングルチャネルイベントのみを使用して実行されてもよい。
【００７５】
上記で説明したように、エネルギー補正は、マルチチャネルイベント座標Ｘ→を使用してインデックス付けされるルックアップテーブルに空間的に依存し、それに基づき得る。次に、座標Ｘ→を決定する非限定的な方法が与えられる。マルチチャネルイベントは、上記で説明した初期の粗い較正を使用するか、あるいは既知のエネルギーを有するガンマ線源を使用した測定値から補正された較正（ＰＥＴ）データから収集される。次いで、合計エネルギーＥｓｕｍ＝ΣｉｎＥｉＷ及びエネルギー重み付けされた２次元座標Ｘ→＝Σｉｎｘi→ＥｉＷ／Ｅｓｕｍがマルチチャネルイベントごとに計算される。上記で説明したように、ｘi→はチャネルｉの中央座標であり、Ｅiは、エネルギー値（例えば、非線形較正されたエネルギー値又は未加工のエネルギー値）であり、ｗは、０以外の任意の値であり得るエネルギーの重みである。
【００７６】
エネルギーウィンドウは、信頼できるとみなすには低すぎるかあるいは高すぎるエネルギーを有するヒットをフィルタで除去するために使用されてもよい。
【００７７】
ステップ１８０におけるエネルギー較正の生成は、完全な２次元空間座標Ｘ→を複数のブロックにセグメント化することを含んでもよい。各ブロック中で、合計エネルギーの関数としてのカウント数のヒストグラムが、図１０Ａに示すように生成される。例えば、図１０Ａでは、「１ヒット」のヒストグラムは、シングルチャネルイベントに対応する。「２ヒット」のヒストグラムは、２つのヒットに関するマルチチャネルイベント（例えば、１次コンプトン散乱）に対応する。「３ヒット」のヒストグラムは、３つのヒットに関するマルチチャネルイベント（例えば、１つの１次ガンマ線と２つの散乱ガンマ線とを有する２次コンプトン散乱）に対応し、以下同様である。図１０Ａでは、マルチチャネルイベントに対応するヒストグラムのピークは、陽電子放射についての５１１ｋｅＶの既知のエネルギーで発生しない。従って、ヒストグラムのピークのエネルギーが決定され、決定されたピークを既知のエネルギー（例えば、陽電子放射源が使用されるときの５１１ｋｅＶ）にシフトするための補正パラメータが決定される。
【００７８】
補正パラメータは、乗法定数又は加法定数である。例えば、補正パラメータＣ（Ｘ→）が乗法定数であるとき、所与のブロックについての補正パラメータは、ヒストグラムＥｈｉｓｔ＝ｍｏｄｅ（Ｅｓｕｍ（Ｘ→））において、中央エネルギー（例えば、ピークエネルギー）に対する入射ガンマエネルギーＥ0（例えば、５１１ｋｅＶ）の比として計算される。
【００７９】
【数８】
2020064045.tif 000010
【００８０】
ここで、演算子（・）は、ヒストグラム中の最も高頻度の値を返す。既知のエネルギーにシフトされるべきである中央エネルギーとしてモードを使用する代わりに、ヒストグラムの中央又はピークエネルギーは、中央値又は平均ヒストグラムによって与えられるか、又はガウスなどの関数形態にヒストグラムを適合させることによって決定されてもよい。本明細書で説明する方法の趣旨から逸脱することなく、ヒストグラムＥhistの中央エネルギー（例えば、幾何又は算術平均など）を決定するために任意の方法が使用されればよい。
【００８１】
補正パラメータが加法定数であるとき、所与のブロックについて補正パラメータＣ（Ｘ→）は、入射ガンマエネルギーＥ0（例えば、５１１ｋｅＶ）と中央又はピークヒストグラムエネルギーＥhistとの差として次式のように計算される。
【００８２】
【数９】
2020064045.tif 000011
【００８３】
所与のブロックが信頼できる統計解析を実行するのに十分なカウントがないと決定され、ピークエネルギーＥhistを正確に決定することができない場合、補正パラメータＣ（Ｘ→）にデフォルト値が割り当てられる。デフォルト値は、例えば、補正パラメータＣ（Ｘ→）が乗法であるとき１の値、補正パラメータＣ（Ｘ→）が加法であるとき０の値に設定される。
【００８４】
エネルギー分解能をより良好にするため、多重度（すなわち、ヒット数）に応じた補正ルックアップテーブルと補正パラメータＣ（Ｘ→）とが生成及び使用されてもよい。例えば、多重度の様々なレベルごとに異なる補正パラメータＣ（Ｘ→）、例えば、２結晶イベント（すなわち、２つのヒット）のためにＣ２（Ｘ→）、３結晶イベント（すなわち、２つのヒット）のためにＣ３（Ｘ→）などが生成される。
【００８５】
更に、ブロックのサイズは、異なる補正パラメータＣ２（Ｘ→）又はＣ３（Ｘ→）等毎に異なる値に設定されてもよい。例えば、異なる多重度レベルに対応する座標空間Ｘ→のそれぞれのセグメント化（例えば、ブロックサイズ及び粗視化の程度）は、分解能について予め定義された目標及び予め定義された統計要件（例えば、所望の信号対雑音比ＳＮＲ）等の予め定義された基準に依存することができる。例えば、所与の多重度におけるブロックサイズ（すなわち、マルチチャネルイベントごとのヒット数）は、ブロックごとのカウント数が予め定義された閾値を超えることを保証することに基づき、その結果、多重度が高くなるとブロックサイズがより大きくなる。すなわち、多重度が増加するとカウントレートが減少し得る。更に、ブロックサイズが予め定義された最大分解能の限界を超える多重度は、破棄され、再構成から省略され得る。
【００８６】
例えば、図１０Ａに示すように、多重度「４ヒット」について、ヒストグラムのピーク又は中央エネルギーを決定することを困難にする低カウントレートが観察される。従って、これらのカウント数は、破棄され、断層撮影画像再構成中に使用されないことがある。この場合、最大多重度レベルは「３ヒット」になる。従って、最大多重度レベルと収集されたデータの必要とされる統計とは、目標とする分解能に応じて予め決定される。
【００８７】
ステップ１４０において、エネルギー較正１８５を適用し、ＰＥＴデータ１０５を補正するために、ステップ１２０において判別されたＰＥＴデータ１０５のマルチチャネルイベントの各々について、ＰＥＴデータ１０５から合計エネルギーＥsumが計算される。更に、ステップ１３０において計算された座標Ｘ→は、エネルギー較正１８５において、マルチチャネルイベントの各々について、補正係数Ｃｎ（Ｘ→）を探索するために使用される。ここでｎは多重度を指す。次に、補正パラメータ／係数Ｃｎ（Ｘ→）は、合計エネルギーＥsumに適用される。例えば、補正係数Ｃ（Ｘ→）が乗法であるとき、補正エネルギーは次式によって与えられる。
【００８８】
【数１０】
2020064045.tif 000012
【００８９】
補正係数Ｃ（Ｘ→）が加法であるとき、補正エネルギーは次式によって与えられる。
【００９０】
【数１１】
2020064045.tif 000013
【００９１】
幾つかの適用例では、再構成は、エネルギーがエネルギーウィンドウの範囲内にあるか否かのみに依存する。例えば、エネルギーウィンドウは、画像再構成において使用するのに十分信頼できる同時計数カウントを選択するために使用される。この場合、実際のエネルギー値は有意ではなく、より単純且つ高速な手法は、合計エネルギー値をシフトするのではなく、エネルギーウィンドウをシフトすることである。例えば、エネルギーウィンドウが［ＥMin，ＥMax］で、補正係数Ｃ（Ｘ→）が乗法であるとき、シフトされたエネルギーウィンドウは、［ＥMin／Ｃｎ（Ｘ→）、ＥMax／Ｃｎ（Ｘ→）］によって与えられる。同様に、補正係数Ｃ（Ｘ→）が加法であるとき、シフトされたエネルギーウィンドウは［ＥMin−Ｃｎ（Ｘ→）、ＥMax−Ｃｎ（Ｘ→）］によって与えられる。従って、多重度及びブロックごとに、補正係数Ｃ（Ｘ→）は、所与のブロック内の合計エネルギーＥsumの各々に繰り返し適用されるのではなく、エネルギーウィンドウの各端部ＥMin及びＥMaxに１回適用される。
【００９２】
図１０Ｂは、補正ＰＥＴデータが使用されることを除いて、図１０Ａと同様のヒストグラムを示す。図１０Ｂに示すように、２ヒットマルチチャネルイベントに関するエネルギーピークは、既知のエネルギー（すなわち、５１１ｋｅＶ）により厳密に一致するようにシフトされている。同様に、３ヒットマルチチャネルイベントに関するエネルギーピークは、既知のエネルギーに向かってシフトされている。
【００９３】
上記に鑑みて、本明細書で説明するエネルギー較正方法は、幾つかの利点を有する。第１に、これらの方法を使用して正確なエネルギー測定値が復元され、マルチチャネル／結晶イベントについて測定されたエネルギーがシングルチャネル／結晶イベントと同等になる。第２に、ＴＯＴ測定値の初期の粗いエネルギー較正の後に残る残余の非線形性が緩和され得る。第３に、本実施形態に係る方法は、計算量的に効率的であり、ルックアップテーブルを使用して容易に実装され得る。第４に、本実施形態に係る方法は、非物質的な結果に対してロバストであり得る。例えば、補正されたエネルギーは負のエネルギー値を生成しない。
【００９４】
本実施形態に係る方法の様々な実装形態は、任意の非線形エネルギー測定値及び／又はピクセル型ガンマ検出器のチャネル間のクロストークによりマルチチャネルイベントで起こる合計エネルギーのシフトを補正することを含むことができる。これらの方法は、（ｉ）ガンマ線について既知のエネルギーを有するガンマ線源を使用して較正データを取得することと、（ｉｉ）マルチチャネルイベントのエネルギー重み付けされた位置を計算することと、（ｉｉｉ）ガンマ線の既知のエネルギーを有する合計エネルギーを比較することによって補正ルックアップテーブルを生成することと、（ｉｖ）ＰＥＴスキャンからのＰＥＴデータ中のマルチチャネルイベントにルックアップテーブル中の補正係数を適用することを含むことができる。
【００９５】
幾つかの実装形態では、ガンマ検出器は、別個の結晶要素のアレイにセグメント化されるシンチレータ結晶を含むことができ、又は、代替的に、幾つかの他の実装形態では、ガンマ検出器は、別個の結晶要素のアレイにセグメント化されないシンチレータ結晶を含むことができる。更に、シンチレータ結晶は、光検出器に一対一で結合されてもよく、また、光検出器に一対一で結合されなくてもよい。
【００９６】
個々のエネルギー読取値に初期の粗い非線形補正を適用するオプションが設けられており、それによって、エネルギー較正の性能を改善する。幾つかの実装形態では、この初期の粗い非線形補正は、ルックアップテーブル又はパラメータ化された式を使用して実行される。パラメータ化された式は、第一原理に基づき、経験的又は実験的に決定される。
【００９７】
マルチチャネルイベントは、予め定義された時間ウィンドウ内で認識され、全て非ゼロエネルギーであると判別する１つのマルチチャネルを選択又はフィルタ処理され、予め定義された時間ウィンドウ内で非ゼロエネルギーであると判別するチャネル又はシンチレータ結晶は、予め定義された距離内にある。
【００９８】
幾つかの実装形態では、マルチチャネルイベントごとに、非線形補正又は未加工のいずれかのエネルギー読取値が、エネルギー重み付けされた２次元座標を計算するために使用される。エネルギーは、予め定義された非ゼロの冪指数ｗによって累乗される。
【００９９】
幾つかの実装形態では、較正データからのマルチチャネルイベントは、それらの２次元座標及びそれらの多重度（すなわち、マルチチャネルイベント中のヒット／チャネルの数）に従ってブロックにグループ化／クラスタリングされる。各グループ中で、ヒストグラムの中央又はピークエネルギーが決定され、補正係数（例えば、乗法補正係数あるいは加法補正係数）がガンマ線源の既知のエネルギーに中央又はピークエネルギーをシフトするために計算される。ヒストグラムピークの中央又はピークエネルギーを決定するのに十分なカウントがない場合、補正係数はデフォルト値に設定される。デフォルト値は、例えば、補正係数が乗法補正係数の場合、１に設定され、又は補正係数が加法補正係数の場合、０に設定される。
【０１００】
ピクセル型ガンマ検出器の２次元座標中のブロックへのセグメント化は、補正されるべきエネルギーシフトの信頼できる推定値を達成するのに必要である分解能の指定された目的／基準ならびに／あるいはブロックごとのカウントに基づいて各多重度でブロックサイズを生成するように、実行される。例えば、ブロックサイズは、エネルギー分解能の改善のバランスをとることと較正のために必要とされるリソースとの間のトレードオフに影響を及ぼすよう選択される。
【０１０１】
多重度各々について、補正係数とエネルギー重み付けされた２次元座標の位置とが関連付けられた補正ルックアップテーブル（ＬＵＴ）が形成されてもよい。
【０１０２】
多重度でのセグメント化に依存して、１つのＬＵＴあるいは複数のＬＵＴが生成されてもよい。複数のＬＵＴが生成された場合、異なるＬＵＴ中での２次元座標のセグメント化において使用されるブロックサイズは同じである必要がない。むしろ、異なるＬＵＴを生成するために異なる多重度のセグメント化において使用されるブロックサイズは、カウント密度、エネルギー補正の精度及び／又はエネルギー補正の大きさによって決定される。ＬＵＴは、オフライン補正のために、又はデータ収集中の（例えば、リアルタイムでの）オンライン補正のために使用される。
【０１０３】
全てのＬＵＴが較正されると、イベントのエネルギー重み付けされた座標及びヒットの回数を使用してＬＵＴ中の補正係数を探索することによってマルチチャネルイベントの合計エネルギーが得られる。
【０１０４】
あるエネルギーウィンドウ内のイベントを選択する場合、合計エネルギー補正は、イベント選択のために使用されるエネルギーウィンドウに直接的に適用される。また、補正されたエネルギーウィンドウは、次いで、イベントを保持し、画像の再構成のためにイベントを使用すべきか否かを決定するために較正されていないデータに適用されることになる。この場合、補正されたＰＥＴデータは、シフトされたエネルギーウィンドウの外部にイベント／カウントのない、シフトがイベントの記録されたエネルギー値に適用されることのないＰＥＴデータである。
【０１０５】
上記の説明において、本実施形態に係る撮像システムは、ＰＥＴシステムに適用されるとしたが、例えば、ＳＰＥＣＴシステム等の他の撮像システムに適用されてもよい。
【０１０６】
上記の説明の通り、本実施形態に係る医用画像処理装置は、プロセッサ２７０を有する。プロセッサ２７０は、イベント取得部、補正係数取得部、選択部及び補正部を有する。イベント取得部は、複数の検出器要素を有するガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すＰＥＴデータを複数の検出イベント毎に取得する。補正係数取得部は、複数の検出器要素のうちの一次ガンマ線のエネルギーの一部を吸収する検出器要素とは異なる検出器要素において散乱ガンマ線のエネルギーが吸収されるマルチチャネル検出中に生じる、クロストークに起因するエネルギーシフトを補正するための補正係数を取得する。補正係数は、エネルギー較正の一部である。選択部は、ＰＥＴデータに基づいて複数の検出イベントの中からマルチチャネルイベントを選択する。補正部は、補正係数をマルチチャネルイベントに属する検出イベントのエネルギーに適用し、補正エネルギーを生成する。
【０１０７】
補正係数の種類としては、加法補正係数と乗法補正係数とがある。加法補正係数の場合、補正部は、マルチチャネルイベント毎に、マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーと補正係数との加算に基づいて補正エネルギーを生成する。乗法補正係数の場合、補正部は、マルチチャネルイベント毎に、マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーと補正係数との乗算に基づいて補正エネルギーを生成する。
【０１０８】
補正エネルギーの生成方法としては、合計エネルギーに補正係数を適用する方法のみに限定されない。補正エネルギーを生成するため、補正部は、エネルギーウィンドウのうちの最大エネルギー及び最小エネルギーから補正係数を減算することによってエネルギーウィンドウをシフトし、マルチチャネルイベント毎に、マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーを生成し、シフトされたエネルギーウィンドウ内に合計エネルギーを有するマルチチャネルイベントを選択してもよい。あるいは、補正部は、補正エネルギーを生成するため、エネルギーウィンドウのうちの最大エネルギー及び最小エネルギーを補正係数で除算することによってエネルギーウィンドウをシフトし、マルチチャネルイベント毎に、マルチチャネルイベントに属する２以上の検出イベントのエネルギーの合計エネルギーを生成し、シフトされたエネルギーウィンドウ内に合計エネルギーを有するマルチチャネルイベントを選択してもよい。
【０１０９】
マルチチャネルイベントの選択の具体例は以下の通りである。ＰＥＴデータは、更にガンマ線の検出時刻を表している。選択部は、複数の検出イベントのうちの最小エネルギー閾値を上回るエネルギーを有する２以上の特定の検出イベントを選択し、選択された２以上の特定の検出イベントの任意の組合せ各々の中から、入射位置と検出時刻とエネルギーとの少なくとも一要素に基づいて、同一のマルチチャネルイベントに属する２以上の検出イベントを選択する。
【０１１０】
上記の通り、補正係数は、複数の検出器要素内でのガンマ線のヒット数に依らず共通の第１の補正係数でもよいし、ヒット数に応じて複数の第２の補正係数でもよい。第２の補正係数を用いる場合、選択部は、マルチチャネルイベントをヒット数に応じて分類し、分類後の各ヒット数のマルチチャネルイベントのエネルギーに、複数の第２の補正係数のうちの対応するヒット数の第２の補正係数を適用する。
【０１１１】
補正係数は、ヒット数ではなく、複数の検出器要素各々の空間座標に応じて複数の第３の補正係数を有してもよい。この場合、選択部は、前記マルチチャネルイベントのエネルギーに、複数の第３の補正係数のうちの、マルチチャネルイベントに属する複数の検出イベントにそれぞれ対応する複数の入射位置の重み付き平均に対応する第３の補正係数を適用する。重み付き平均は、複数の検出イベントにそれぞれ対応する複数の入射位置の正規化加重和に基づいて計算される。当該重みは、検出イベントのエネルギーを、非ゼロの実数で累乗した値である。空間座標は、複数の第３の補正係数各々に対応するブロックにセグメント化されてもよい。ブロックのサイズは、予め定義された空間分解能基準及び前記複数の補正係数の精度又は信頼性に関する予め定義されたカウント基準のうちの１つ又は複数を満たすように選択される。あるいは、ブロックのサイズは、マルチチャネルイベントのヒット数に応じて異なってもよい。
【０１１２】
上記の通り、補正部は、マルチチャネルイベントに属する検出イベントのエネルギーに、補正係数を適用する前において、初期の非線形エネルギー補正を適用してもよい。補正部は、初期の非線形エネルギー補正として、検出されたガンマ線のエネルギー測定値と実際のエネルギーとの間の非線形関係を補正する。
【０１１３】
上記の通り、ＰＥＴ撮像の前段階においてイベント取得部による較正データの取得（較正段階）が行われる。較正段階においてイベント取得部は、ガンマ線検出器により検出されたガンマ線の入射位置とエネルギーとを表すＰＥＴデータ（較正データ）を複数の検出イベント（較正検出イベント）毎に取得する。選択部は、較正データに基づいて、複数の較正検出イベントの中から較正マルチチャネルイベントを選択する。補正係数取得部は、較正マルチチャネルイベント毎に、マルチチャネルイベントに属する２以上の検出イベントのエネルギーの代表エネルギーと、予め定義されたエネルギーとのシフト量に基づいて補正係数を生成する。補正係数は、乗法補正係数又は加法補正係数を含む。ガンマ線検出器は、各々が複数の検出器要素を含む複数の検出器要素ブロックを有する。代表エネルギーは、複数の検出器要素ブロック各々によるマルチチャネルイベントのカウントのヒストグラムにおける中央値、モード及び平均値のうちの１つに対応するエネルギーである。
【０１１４】
本実施形態に係る医用画像処理装置は、ＰＥＴシステム又はＳＰＥＣＴシステムに組み込まれても良い。この場合、医用画像処理装置は、補正エネルギーに基づいて画像を再構成する再構成部を更に備える。
【０１１５】
以上説明した少なくとも１つの実施形態によれば、マルチチャネル検出におけるクロストークに起因するエネルギー測定値の劣化を高精度の補正することができる。
【０１１６】
本発明の幾つかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１１７】
２００  ＰＥＴシステム
２１６  寝台
２１６  天板
２２０  円形構成要素
２４０  架台
２７０  プロセッサ
２７４  ネットワークコントローラ
２７６  データ収集システム（ＤＡＳ）
２７６  データ収集システム
２７８  メモリ
【図１】
2020064045.tif 000014
【図２Ａ】
2020064045.tif 000015
【図２Ｂ】
2020064045.tif 000016
【図３Ａ】
2020064045.tif 000017
【図３Ｂ】
2020064045.tif 000018
【図３Ｃ】
2020064045.tif 000019
【図４Ａ】
2020064045.tif 000020
【図４Ｂ】
2020064045.tif 000021
【図４Ｃ】
2020064045.tif 000022
【図４Ｄ】
2020064045.tif 000023
【図５】
2020064045.tif 000024
【図６Ａ】
2020064045.tif 000025
【図６Ｂ】
2020064045.tif 000026
【図７Ａ】
2020064045.tif 000027
【図７Ｂ】
2020064045.tif 000028
【図８Ａ】
2020064045.tif 000029
【図８Ｂ】
2020064045.tif 000030
【図９】
2020064045.tif 000031
【図１０Ａ】
2020064045.tif 000032
【図１０Ｂ】
2020064045.tif 000033