JP5984246B2 - 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法 - Google Patents

陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法 Download PDF

Info

Publication number
JP5984246B2
JP5984246B2 JP2012029889A JP2012029889A JP5984246B2 JP 5984246 B2 JP5984246 B2 JP 5984246B2 JP 2012029889 A JP2012029889 A JP 2012029889A JP 2012029889 A JP2012029889 A JP 2012029889A JP 5984246 B2 JP5984246 B2 JP 5984246B2
Authority
JP
Japan
Prior art keywords
signal value
value
detection position
unit
scintillator
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2012029889A
Other languages
English (en)
Other versions
JP2012189583A (ja
Inventor
ケント・シー・バー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Toshiba Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Medical Systems Corp filed Critical Toshiba Medical Systems Corp
Priority to CN201210061629.8A priority Critical patent/CN102670232B/zh
Publication of JP2012189583A publication Critical patent/JP2012189583A/ja
Application granted granted Critical
Publication of JP5984246B2 publication Critical patent/JP5984246B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/24—Measuring radiation intensity with semiconductor detectors
    • G01T1/248—Silicon photomultipliers [SiPM], e.g. an avalanche photodiode [APD] array on a common Si substrate

Landscapes

  • Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Radiation (AREA)
  • Nuclear Medicine (AREA)

Description

実施形態は、陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法に関する。
医用イメージングの分野においては、ガンマ線検出器、特に、陽電子放出コンピュータ断層撮影装置、すなわち、ＰＥＴ（Positron Emission Tomography）装置の利用が増えている。ＰＥＴ装置によるＰＥＴイメージングにおいて、放射性医薬品は、注入、吸入、又は食物摂取により、画像化される被検体に取り込まれる。放射性医薬品の投与後、薬剤の物理的及び生体分子的な特性により、薬剤は、被検体内の特定部位に集積する。薬剤の実際の空間分布、薬剤の蓄積領域の濃度、及び投与から最終排出までのプロセスの動態は全て、臨床的な重要性を持ち得る因子である。このプロセスにおいて、放射性薬剤に付着した陽電子放射体は、半減期、分岐比率などの同位元素の物理的性質に応じて、陽電子を放出する。
放射性核種は、陽電子を放出する。放出された陽電子が電子と衝突すると、対消滅イベントが発生し、陽電子及び電子は消滅する。多くの場合、対消滅イベントによって約１８０度反対方向に放出される２つのガンマ線（５１１ｋｅＶ）が発生する。
ＰＥＴ装置は、ほぼ同時に検出された２つのガンマ線を検出し、これらの検出位置間に、直線、すなわち、ＬＯＲ（Line−of−response）を引くことにより、高い確率で、対消滅したであろう位置を導き出すことができる。このＬＯＲを多数蓄積し、断層撮影再構成プロセスを実行することにより、本来の分布を推定することができる。２つの消滅イベントの位置情報に加え、正確な検出タイミング情報（数百ピコ秒以内）を利用することが可能な場合、ＴＯＦ（time−of−flight）の計算により、ＬＯＲに沿った消滅イベントの推定位置に関する情報を更に加えることができる。ＰＥＴ装置が有する時間分解能の限界によって、このＬＯＲに沿った位置決めの精度が決まることになる。また、本来の消滅イベントの位置を決定する際の限界によって、ＰＥＴ装置の最終的な空間分解能が決まることになる。一方で、同位元素の固有の特性（例えば、陽電子のエネルギー）もまた、（２つの消滅ガンマ線が発生するまでの陽電子の飛程や、２つの消滅ガンマ線間の角度に影響を与える因子なので）ＰＥＴ装置の空間分解能を決める要因となる。
W.W.Moses著，"Time of Flight in PET Revisited"，IEEE Transactions on Nuclear Science，Vol.50，No.5，pp.1325-1330
本発明が解決しようとする課題は、エネルギー分解能を向上することができる陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法を提供することである。
実施形態に係る陽電子放出コンピュータ断層撮影装置は、シンチレータアレイと、複数の光センサと、特定部と、格納部と、導出部とを備える。前記シンチレータアレイは、複数のシンチレータを有する。前記複数の光センサは、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する。前記特定部は、前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線のシンチレータとの相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する。前記格納部は、前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する。前記導出部は、前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値を、前記領域の単位毎に導出する。
図１は、ガンマ線検出器が様々な放射性同位元素からのガンマ線に曝された場合の、非線形なＳｉＰＭ（Silicon Photomultiplier）応答の例を示す。 図２Ａは、薄いライトガイドに連結された５×５のシンチレータアレイ（又は、シンチレータアレイ）から放出される光の分布の３つの例を示す。図２Ａに示す３つの例において、光の強度の積分値はいずれも同じになるが、光の広がりはそれぞれ大きく異なる。図２Ａは、識別された相互作用のシンチレータに基づく１つの非線形補正では不十分であることを示す。 図２Ｂ−１は、図２Ａに示す３つの例の内の１つの例において生成されるであろう概算の信号レベルの例を示す。 図２Ｂ−２は、図２Ａに示す３つの例の内の１つの例において生成されるであろう概算の信号レベルの例を示す。 図２Ｂ−３は、図２Ａに示す３つの例の内の１つの例において生成されるであろう概算の信号レベルの例を示す。図２Ｂ−１〜３に示すように、３つの例のいずれの場合も５１１ｋｅＶのエネルギー値が与えられているが、出力された信号レベル（出力信号の信号レベル）はそれぞれの場合で異なる。 図３は、実施形態に係るＰＥＴ装置の構成を説明するための図である。 図４は、実施形態に係るＰＥＴ装置内の中央処理装置の構成を説明するための図である。 図５は、光学系束ね処理及び４分割の光電子増倍管（ＰＭＴ（Photomultiplier Tube)）を用いたシンチレータアレイの読み取りから得られるフラッドヒストグラム（flood histogram）を示す。 図６は、アナログデータ束ね処理及び位置感応型アバランシェフォトダイオード（avalanche photodiode）を用いたシンチレータアレイの読み取りから得られるフラッドヒストグラムを示す。 図７Ａは、実施形態に係るＰＥＴ装置によって実行される方法のステップを示すフローチャートである。 図７Ｂは、別の実施形態に係るＰＥＴ装置によって実行される方法のステップを示すフローチャートである。 図８Ａは、ガンマ線がシンチレータアレイに当たることより発生した、光の入射に反応した光センサからの典型的なパルスを示す。 図８Ｂ−１は、図８Ａの時間軸を変更したものである。 図８Ｂ−２は、図８Ａに示すパルスに対応するエネルギーの情報を示すものである。 図９Ａは、実施形態に係るガンマ線検出器の構成を示す。図９Ａは、シンチレータアレイ、ライトガイド、及び光センサの位置関係を示す。 図９Ｂは、実施形態に係るガンマ線検出器の構成を示す。図９Ａは、シンチレータアレイ、ライトガイド、及び光センサの位置関係を示す。 図１０Ａは、フラッドヒストグラムを示す。 図１０Ｂは、フラッドヒストグラムに対応するシンチレータＩＤルックアップテーブルを示す。
以下、実施形態に係るＰＥＴ装置、ＰＥＴ装置に実行させるためのプログラム、及びＰＥＴ装置によって実行される方法を、図面を参照しながら説明する。
実施形態に係るＰＥＴ装置は、顕著な非線形性を示す光センサを有するガンマ線検出器を備える。その非線形性の程度は、光センサ表面上の光子束の空間分布に依存する。実施形態に係るＰＥＴ装置は、必要な電子チャネルの数を減らすために光学系の束ね処理又はアナログ信号の束ね処理が可能なガンマ線検出器を備え、非線形な補正を行う改良方法を提供する。
ＰＥＴ装置は、多数の対消滅ガンマ線の検出イベントを収集し、このイベントデータを基に画像再構成を行うことで、被検体の投与薬剤分布の画像を生成する。対消滅ガンマ線は、被検体に投与された薬剤から発生する陽電子が発生点の近傍の電子と対消滅することによって約１８０度の方向に発生する２本のガンマ線である。ＰＥＴ装置では、ガンマ線検出器の検出素子においてほぼ同時に検出されたイベント情報から、幾何学的な特性に応じてヒストグラム化が可能なＬＯＲが形成される。このＬＯＲを基に、再構成の対象データである投影データや、サイノグラムデータが規定される。また、ＬＯＲを画像データに１つずつダイレクトに書き加えることによっても、画像データを形成することができる。
従って、ＬＯＲはＰＥＴ装置において、投与薬剤の対消滅位置情報を与える基本的な要素であるが、この対消滅位置情報に関し、下記に示す更なる情報を得ることができる。まず、周知のこととして、ＰＥＴ装置がある点を描写する能力は、画像再構成領域全体にわたり空間的に変化するが、中心部でより高くなり、周辺に向かって徐々に低下する。ＰＳＦ（Point Spread Function）が、この特性を表現する典型的なものであり、近年、再構成プロセスに組み込まれている。次に、ＰＥＴ装置は、ＴＯＦ、すなわち、対消滅ガンマ線が検出器に到着する時間差情報を基に、ＬＯＲ上で対消滅イベントが発生した点の確率分布を計算することができる。
上記のプロセス（イメージングプロセス）は、多数の対消滅イベントを必要とする。十分なイメージング処理を可能にするために必要な対消滅イベントの数を決定するには、一つ一つのイメージング事例を解析しなければならないものの、典型的な長さ１００ｃｍのＦＤＧ（fluoro-deoxyglucose）での検査については、数億の回数を重ねる必要があるというのが、現在の研究での典型である。これだけのイベント数を得るのに要する時間は、薬剤の注入量及びＰＥＴ装置の感度と計数能力によって決まる。
ＰＥＴイメージングシステム、すなわちＰＥＴ装置は、被検体から放出されるガンマ線を検出するために、互いに対向する位置に配置されたガンマ線検出器を使用する。典型的には、ＰＥＴ装置は、各角度から飛来するガンマ線を検出するために、リング状に配置されたガンマ線検出器を使用する。従って、ＰＥＴ装置のスキャナは、典型的には、できるだけ多くの放射線を捕捉できるように略円筒形状であり、当然のことながら等方性でなければならない。また、スキャナは、円環の一部が欠けている形状のものでも良く、この場合では、欠けている角度を捕捉するために、ガンマ線検出器を回転させることも可能である。しかし、そのような手法は、ＰＥＴ装置の全体的な感度に重大な影響を与える。円筒形状の場合、１つの面に含まれる全てのガンマ線は、ガンマ線検出器に反応する可能性が高く、軸方向の寸法を大きくすると、放射線を捕捉する感度や性能に極めて有効である。従って、最良の構造は、全てのガンマ線が検出される可能性がある球体構造である。当然ながら、人体への適用では、球状構造は、極めて大きくなり、極めて高価にならざるを得ない。従って、現実的には、ガンマ線検出器の軸方向の長さが可変である円筒形状が、最新のＰＥＴ装置のスキャナの構造の基本である。
スキャナは多数の検出器モジュールから構成され、検出器モジュールは多数のシンチレータから構成される。ＰＥＴ装置の性能を良くするためには、できるだけ多くのシンチレータを配置し、できるだけ多くのガンマ線を止めて光に変換する必要がある。ＰＥＴ装置は、放射性同位体元素の時間的空間的分布を再構成するには、検出された各イベントのエネルギー値（すなわち、シンチレータ内で生成された光量）、位置、及びタイミングの情報を得なければならない。ほとんどの最新ＰＥＴスキャナは、数千個の個別のシンチレータで構成される。これらのシンチレータは、検出器モジュール内に配置され、対消滅イベントの位置を特定するために用いられている。典型的には、シンチレータ素子は、約４ｍｍ×４ｍｍの断面を有する。より小さな寸法又は大きな寸法、又は、正方形でない断面も可能である。シンチレータの長さ又は奥行きは、ガンマ線の阻止能を決定し、典型的には、１０〜３０ｍｍの範囲である。検出器モジュールは、スキャナの主要な構成要素である。
ＰＥＴ装置によるＰＥＴイメージングの性能は、高速かつ高輝度のシンチレータによる、ガンマ線から光への変換に依存している。ＰＥＴ装置は、検出器の中でシンチレーションが起こった位置（ガンマ線がシンチレータと相互作用を起こした位置）を特定し、個々のイベントを時間情報をもとにペアリングし（つまり、ある時間窓の中で検出された２つの検出イベントを結びつけ）、対消滅イベントの位置を計算することができる。これらの動作を行うには、非常に高速の検出器（検出器及び電子機器）が必要であり、優れた信号対ノイズ比も必要である。高品質の電子機器を用いれば、信号対ノイズ比は主に、検出プロセスに関係するポアソン分布によって決まる。つまり、より多くの光子を検出すれば信号対ノイズ比が向上し、従って空間分解能及び時間分解能が一層高まることになる。なお、検出プロセスでシンチレーション光を逃してしまった場合、検出器の設計及び電子機器の改善をもっても補償することはできない。シンチレータ内で発生した光に対して実際に検出できた光量の割合は、設計部品の効率を表すのにふさわしい指標である。そのため、光の検出量を最大化しようとする人は誰でも、光センサをできるだけシンチレータに接近させるとともに、反射及び他のエッジ効果を回避しようとするであろう。このようなことを行えば、ガンマ線検出器は、否応なく、シンチレータと光センサとの間の距離が短く、大きなシンチレータアレイを持つものになってしまうであろう。
上述したように、ＰＥＴ装置は、単なる計数器ではない。ＰＥＴ装置は、対消滅イベントの存在の検出の他に、対消滅イベントの位置の識別も行う必要がある。各相互作用の位置を識別できるように最も単純に設計すれば、おそらく、シンチレータ毎に、独立した光センサ及びＡ／Ｄ変換器を設けることになるであろう。光センサの物理的な大きさ、Ａ／Ｄ変換器に必要な電力、及びこれらのコストといった制約のために、光センサ数及び電子機器のチャネル数を削減する目的で、通常は何らかの手段によってデータを束ねている。
データ束ねの最も一般的な２つの形態は、光学系の束ね処理とアナログ信号の束ね処理である。光を４つの光電子増倍管に分散させるためにライトガイドを用いている方式が、光学系の束ね処理の例である。光がどのように複数の光センサに分配されるかを適切に記録することにより、いかなる光センサ応答の組み合わせに対しても、ガンマ線との相互作用イベントの位置を計算することができる。また、アナログ信号の束ね処理の例として、いわゆる位置感応型アバランシェフォトダイオード、位置感応型ＳｉＰＭがある。
光学系束ね処理又はアナログデータ束ね処理の設計では、各相互作用イベントの検出器ブロック内における相対位置が、アンガーロジック（重心計算）又は統計的手法を用いて計算される。ここで、「相対位置」とは、シンチレータとガンマ線の相互作用イベントが発生した実際の位置（例えば、被検体内の絶対的な位置）の相対的な位置であって、ガンマ線がシンチレータに入射した位置を計算することで特定される「検出位置」のことである。相対位置としてガンマ線が相互作用したシンチレータを特定するために、通常、ＰＥＴ装置は、結果として得られるフラッドヒストグラムを、シンチレータのルックアップテーブルに分割する。フラッドヒストグラム及びルックアップテーブルにおける位置に基づいて相互作用したシンチレータが識別された後、エネルギー値及びタイミングの補正が、シンチレータ毎に行われる。エネルギー値の補正として非線形補正も可能である。通常この補正は、多数の異なる同位体（例えば、５１１及び１２７５ｋｅＶの２２Ｎａ、６６２ｋｅＶの１３７Ｃｓ、３５６ｋｅＶの１３３Ｂａ、１２２ｋｅＶの５７Ｃｏ、６０ｋｅＶの２４１Ａｍ）を用いた計測に基づくものであり、シンチレーション毎に行われる。なお、後述するように、実施形態に係るＰＥＴ装置は、エネルギー値の補正を、シンチレータの単位ではなく、シンチレータの数より多い数に分割された領域の単位で行う。
図１は、ガンマ線検出器が様々な放射性同位元素からのガンマ線に曝された場合の、非線形なＳｉＰＭ応答の例を示す。ガンマ線検出器は、３ｍｍ×３ｍｍ×１０ｍｍのＬＹＳＯ（Lutetium Yttrium Oxyorthosilicate）結晶、及び、３６００個のセルを有する３ｍｍ×３ｍｍのＳｉＰＭ（浜松ホトニクス社製）を含む。
高エネルギー値のガンマ線を検出する場合（例えば、ＰＥＴ用の５１１ｋｅＶなど）、ガンマ線検出器内で複数のシンチレータでガンマ線がコンプトン散乱することがある。これによって、複数のシンチレータに対してガンマ線はエネルギーを付与する。一般に、５１１ｋｅＶのガンマ線の３０％超が、複数のシンチレータと相互作用を起こす。シンチレータに付与された総エネルギーが同じ場合であっても、光センサへの光の広がりは、相互作用の過程や、相互作用したシンチレータの数によって、かなり異なることがある。例えば、１つのシンチレータに５１１ｋｅＶのエネルギーが付与された場合と、２つの隣接するシンチレータに合計で５１１ｋｅＶのエネルギーが付与された場合とを比較すると、光センサでの光の広がりに大きな差があることが分かる。
この光の広がりの差は、光センサの非線形補正において問題となる。検出された光子の空間的分布によって、その非線形性の程度が異なるからである。これは、いわゆるＳｉＰＭを用いた場合である。ＳｉＰＭ、すなわち、シリコン光電子増倍管（固体光電子増倍管であるＳＳＰＭ（Solid−State Photomultiplier）、あるいは、ＧＡＰＤやＭＡＰＤとも呼ばれる）は、多数の独立したアバランシェフォトダイオードで構成される。また、ＳｉＰＭは、「マイクロセル」と呼ばれ、ガイガーモードで動作する。ガイガーモードにおいて、ＳｉＰＭは、マイクロセルが１個以上の光子を検出すると放電する。放電の間に放出される電荷は、マイクロセルの静電容量及び動作電圧によって決まる。放出される電荷は、放電を起こす光子の数に左右されない。例えば、シンチレータ内でのガンマ線の相互作用により発生するイベントによって、光のパルスがＳｉＰＭに当たると、多数のマイクロセルが放電し、電気パルスが生じる。この電気パルスの振幅は、放電するマイクロセルの数に比例する。光子の密度が非常に少ない場合、同じマイクロセルに複数の光子が当たる確率は非常に低いため、電気パルスの振幅は、光子数に対して直線的に変化する。光子の密度が増すほど、同じマイクロセルに複数の光子が当たる確率が増すため、信号の非線形性は強くなる。
図２Ａは、薄いライトガイドに連結された５×５のシンチレータアレイから放出される光の分布の３つの例を示す。図２Ａは、シンチレータアレイ内に５１１ｋｅＶのエネルギーを付与する、３つの異なるイベントで予測される光の強度の違いを示す。この例では、５×５のシンチレータアレイ内のシンチレータが光学的に分離され（すなわち、これらのシンチレータがそれぞれの間に反射材を有する）、薄いライトガイドに結合されて、更にアレイ状の光センサ又は位置感応型光センサに結合されている。図２Ａの格子線は、シンチレータの境界を示している。それぞれのケースで、積分した光の強度は同じである。第１（左端）のケースでは、エネルギーの全てが中央のシンチレータに付与され、ピークの光の強度（又は光子束）が最も高い。第２のケースでは、ガンマ線がエネルギーの６７％を中央のシンチレータに付与し（コンプトン散乱による）、残りのエネルギーの３３％を近傍のシンチレータに付与している。第３のケース（右端）では、ガンマ線のエネルギーは２つのシンチレータの間に均一に分布している。このケースでは、結果として得られるシンチレータあたりの光の最大強度が最も低い。光子の密度に依存して、光の強度と出力が変化するＳｉＰＭなどの光センサをこのガンマ線検出器の読み取りに使用すれば、積分した光の強度は同じだが、光の広がりが異なるため、光センサの出力は、ここに示される３つのケースで異なる。その結果、エネルギー分解能が劣る。例えば、図２Ａに示す例では、たとえ５１１ｋｅＶのエネルギーが各ケースでシンチレータに付与されるとしても、３つのケースのそれぞれで非線形補正された信号レベルは異なる（各ケースに同じ非線形補正を適用すると仮定する）。これがエネルギー分解能を悪化させるメカニズムの代表例である。
図２Ｂは、図２Ａに示された３つのケースで生成される、概算の信号レベルを示す。ここでは例として、２つのシンチレータを考える。光の広がりの違い及びＳｉＰＭの非線形性により、たとえ３つのケース全てが２つのシンチレータに合計で５１１ｋｅＶのエネルギーを付与しても、それぞれのケースで積分した信号レベルは異なる。
上記を詳細に説明する。例えば、１つのシンチレータの中央付近にガンマ線が入射した場合には、コンプトン散乱の効果は比較的小さくなり、光センサに向かう光の広がりは小さくなり易い。一方、例えば、１つのシンチレータの境界付近にガンマ線が入射した場合には、コンプトン散乱の効果は比較的大きくなり、光センサに向かう光の広がりは大きくなり易い。前者の場合、光の広がりは小さいので、光子は比較的少ない数のマイクロセル（光センサであるＳｉＰＭのマイクロセル）に当たり、後者の場合、光の広がりは大きいので、光子は比較的多い数のマイクロセルに当たると考えられる。ここで、ＳｉＰＭは、同じマイクロセルに連続して光子が当たった場合に、そのハードウェアの制限上、全ての光子を検出できない場合がある。一方で、マイクロセルが放電することによって生じる電気パルスの振幅は、放電するマイクロセルの数に比例する。すなわち、光子が比較的少ない数のマイクロセルに当たる前者の場合、同じマイクロセルに連続して光子が当たる確率も高くなるので、ＳｉＰＭが全ての光子を検出できない確率も高くなり、信号値は低い傾向となる。反対に、光子が比較的多い数のマイクロセルに当たる後者の場合、同じマイクロセルに連続して光子が当たる確率は低くなるので、ＳｉＰＭは全ての光子を検出できることになり、信号値は高い傾向となる。このように、１つのシンチレータ内のどの位置にガンマ線が入射したかが、信号値に影響を与える。これが、実施形態に係るＰＥＴ装置が、以下に説明するように、光センサによって出力された信号値に対して、シンチレータの単位よりも細かい単位で、非線形の補正を行う理由である。
本明細書で開示する実施形態は、非線形補正（non-linearity correction）を、サブピクセル（sub-pixel）単位又は連続的に適用する。ここで、サブピクセル単位とは、シンチレータの数より多い数に分割された領域の単位のことである。光学系束ね処理又はアナログデータ束ね処理を用いると、個々のシンチレータに対応するピークに関するフラッドヒストグラムの中のシンチレータとガンマ線の相互作用イベントの位置に基づいて、複数のシンチレータ間の相互作用についての追加情報を得ることができる。フラッドヒストグラム内のシンチレーションが起こった位置は、アンガーロジック（重心計算）に基づき計算されるので、２つのシンチレータ内にエネルギーを付与するイベントは、２つの個々のシンチレータに対応するピークの間に配置される。当業者には自明であるが、これらのケースで計算された相互作用位置は、必ずしも単一の物理的位置に対応するわけではなく、むしろ光センサにより生成される相対的信号レベルを示すものである。
更に、ノイズからの影響を無視すると、個々のシンチレータのピーク間の正確な位置は、２つのシンチレータ内に付与される相対的なエネルギーによって決まることになる。非線形性の変動も、２つのシンチレータ内に付与される相対的なエネルギーによって決まる。従って、フラッドヒストグラムをサブピクセル化し、異なる非線形補正を各サブピクセル領域に適用することにより、より良い非線形補正を行うことができる。あるいは、個々のシンチレータのピークの間で連続的に変化する非線形補正を適用することもできる。
図３は、実施形態に係るＰＥＴ装置の構成を説明するための図である。当業者には自明であるが、図３に示すガンマ線検出器システムは、ＰＥＴ装置やＴＯＦ型ＰＥＴ装置の一部を形成する。ＰＥＴ装置及びＴＯＦ型ＰＥＴ装置についての更なる解説は、簡潔にするため省略する。なお、ＴＯＦ型ＰＥＴ装置についての解説が非特許文献１にあり、その内容全体が参照により本明細書に組み込まれるものとする。
図３では、非線形の光センサ１３５及び１４０はライトガイド１３０の上に配置され、複数のシンチレータを有するシンチレータアレイ１０５がライトガイド１３０の下に配置されている。ＳｉＰＭやアレイ状のＳｉＰＭを含む任意の非線形光センサを使用しているガンマ線検出器に本実施形態が適用できることは、当業者には自明である。複数のシンチレータを有する第２のシンチレータアレイ３０５が、シンチレータアレイ１０５に対向して、ライトガイド３３０と、光センサ３３５及び３４０に重ねて配置される。
図３において、対消滅ガンマ線は、被検体（図示を省略）から放出されると、互いに約１８０度の反対方向に進む。対消滅ガンマ線は、シンチレータ１００とシンチレータ３００とで、ほぼ同時に検出される。そして、既定制限時間内に、シンチレータ１００とシンチレータ３００とでガンマ線が検出されると、相互作用イベント１１０が特定される。こうして、ガンマ線検出システムは、ガンマ線を、シンチレータ１００とシンチレータ３００とでほぼ同時に検出する。しかし簡略化のために、シンチレータ１００でのガンマ線検出のみをここで説明する。当業者には自明であるが、シンチレータ１００についての説明は、シンチレータ３００でのガンマ線検出にも同様に適用できる。
図３に戻り、各光センサ１３５、光センサ１４０、光センサ３３５、光センサ３４０は、データ収集部３５０又はデータ収集部３６０に接続されている。データ収集部３５０、データ収集部３６０は、シンチレーション光に応答する光センサ１４０、光センサ１３５、光センサ３４０、光センサ３３５で生成される対応波形を積分することにより、デジタル化した出力値を生成する。
データ収集部３５０、データ収集部３６０は、１ギガヘルツ〜５ギガヘルツのサンプリングレートで動作するシグマデルタ変換器等のアナログデジタル変換器を含む。あるいは、データ収集部３５０、データ収集部３６０は、一定のサンプリングレートではなく、電圧閾値をトリガに用いて光センサ波形をサンプリングする多重閾値サンプラを含んでもよい。当業者には自明であるが、本実施形態の範囲から逸脱することなく、他のサンプリング方法及びデータ収集部を含むことも可能である。例えば、エネルギー値とタイミングでそれぞれ別のチャネルを用いることも可能である。この場合、典型的にエネルギーチャネルには、整形フィルタ及びより低いサンプリングレートのアナログデジタル変換器を使用する。また、タイミングチャネルは、典型的には、複数の光センサからの信号を合計する。次に、合計されたタイミング信号が比較器に入力され、時間デジタル変換器で各イベントの到着毎にタイムスタンプが生成される。
出力値が得られると、出力値は演算部３７０に送られ、以下に説明する方法に基づきシンチレータび相互作用イベントのエネルギーレベルが決定される。次に、出力値及び到着時間が電子記憶装置３７５に記憶され、ディスプレイ３８５に表示が可能である。インタフェース３８０は、演算装置３７０の構成及び制御の両方又は一方を行うことと、演算部３７０に更なる命令を出すこととの両方又は一方のために、使用されてよい。
当業者には自明であるが、ディスプレイ３８５は、ＣＲＴ（Cathode Ray Tube）やＬＣＤ（Liquid Crystal Display）などであってよい。インタフェース３８０は、キーボード、マウス、トラックボール、マイクロホン、タッチスクリーンなど、中央処理装置を接続して機能させるための既知のデバイスであってよい。同様に当業者には自明であるが、電子記憶装置３７５は、ハードディスクドライブ、ＣＤ−ＲＯＭ、ＤＶＤディスク、フラッシュメモリ、あるいは別の中央処理装置であってよい。更に、電子記憶装置３７５は、演算部３７０から着脱可能又は分離可能であり、あるいはこれに取り付けられていてもよい。電子記憶装置３７５は、ネットワークを介して演算装置に接続されて、そのため別の部屋や建物等、演算部３７０関連の別の場所に設置されてもよい。
図４は、実施形態に係るＰＥＴ装置内の中央処理装置の構成を説明するための図である。演算部３７０は、メインメモリ４４０とＲＯＭ（Read Only Memory）４５０の両方又は一方の中に記憶されたデータと命令を処理する処理部４８０を含む。また処理部４８０は、ディスク４１０又はＣＤ−ＲＯＭ４２０に記憶した情報を処理してもよい。例示的な処理部４８０は、米国Ｉｎｔｅｌ社製のＸｅｎｏｎプロセッサ（登録商標）あるいは米国ＡＭＤ社製のＯｐｔｅｒｏｎプロセッサ（登録商標）であってよい。当業者には自明であるが、処理装置４８０はＰｅｎｔｉｕｍ（登録商標）プロセッサ（登録商標）やＣｏｒｅ２Ｄｕｏプロセッサ（登録商標）などであってもよい。こうして、ガンマ線検出器のための方法に対応する命令が、ディスク４１０、ＣＤ−ＲＯＭ４２０、メインメモリ４４０、あるいはＲＯＭ４５０のいずれか１つに記憶されてよい。
また演算部３７０は、インターネットやプライベートネットワークなどのネットワークとインタフェースで接続するために、米国Ｉｎｔｅｌ社製のＩｎｔｅｌ Ｅｔｈｅｒｎｅｔ（登録商標） ＰＲＯネットワークインタフェースカード（登録商標）などのネットワークインタフェース４７５を含んでもよい。ディスプレイ制御部４３０は、ディスプレイ３８５を接続して機能させるための米国ＮＶＩＤＩＡ社製のＮＶＩＤＩＡ Ｇ−Ｆｏｒｃｅ ＧＴＸグラフィックスアダプタ（登録商標）であってよい。また演算部３７０は、キーボード２９５、ポインティングデバイス２８５、又はマイクロホン、トラックボール、ジョイスティック、タッチスクリーンなどのその他の汎用インタフェースを接続して機能させるためのＩ／Ｏインタフェース４９０を含んでもよい。
ディスク制御部４６０が、ディスク４１０及びＣＤ−ＲＯＭ４２０又はＤＶＤドライブをバス４７０に相互接続する。ディスク４１０は、ハードディスクドライブ又はフラッシュメモリドライブであってよい。バス４７０は、業界標準アーキテクチャ（ＩＳＡ（Industry Standard Architecture））、拡張業界標準アーキテクチャ（ＥＩＳＡ（Extended Industry Standard Architecture））、ビデオ電子装置規格化協会（ＶＥＳＡ（Video Electronics Standards Association））、周辺部品相互接続（ＰＣＩ（Peripheral Component Interconnect））などの、演算部３７０の構成部品の全てを相互接続するものであってよい。演算部３７０の構成部品の一般的機能及び機能性については周知のものであるため、簡潔のため説明を省略する。もちろん、米国Ｆｒｅｅｓｃａｌｅ社製のＦｒｅｅｓｃａｌｅ ＣｏｌｄＦｉｒｅ、Ｉ．ＭＸ（登録商標）、及びＡＲＭプロセッサ（登録商標）など、当技術分野で知られている他の処理装置やハードウェアメーカー及び製品を使用してもよい。
また例示的な演算部３７０は、ＦＰＧＡ（Field Programmable Gate Array）、特ＡＳＩＣ（Application Specific Integrated Circuit）、マイクロコントローラ、ＰＬＤ（Programmable Logic Device）、又は光ディスクなどのコンピュータ可読媒体上で別々に実装されてもよい。更に、例示的な演算部３７０は、ＰＣ（Personal Computer）などの計算機器のハードウェアプラットフォームであり、処理部４８０は、Ｉｎｔｅｌ Ｐｅｎｔｉｕｍ（登録商標）プロセッサ（登録商標）など、当技術分野で知られている任意の処理装置である。メインメモリ４４０、ＲＯＭ４５０、ディスク４１０、又はＣＤ−ＲＯＭ４２０のいずれかに記憶されているコンピュータで読み取り可能な命令が、ユーティリティアプリケーション、バックグラウンドデーモン、又はオペレーティングシステムの構成要素、あるいはこれらの組み合わせとして提供され、処理部４８０、及びＭｉｃｒｏｓｏｆｔ Ｗｉｎｄｏｗｓ（登録商標） Ｖｉｓｔａ（登録商標）、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ（登録商標）、Ｌｉｎｕｘ（登録商標）、及びＡｐｐｌｅ Ｍａｃ ＯＳ（登録商標）等の当業者に知られているシステムと連動して実行される。
メインメモリ４４０及びＲＯＭ４５０の両方又は一方は、演算部３７０の記録及び同様の機能をサポートする。このためメインメモリ４４０は、ＲＡＭ（Random Access Memory）、フラッシュメモリ、ＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory）メモリなどであってよい。一方、ＲＯＭ４５０は、ＰＲＯＭなどの読み出し専用メモリである。更に、メインメモリ４４０及びＲＯＭ４５０の説明は、このようなメモリは周知のものであるため、簡潔のため省略する。
図５及び図６は、シンチレータ間のコンプトン散乱を明確に識別することができる、２つのフラッドヒストグラムを示す。図５は、光学系束ね処理として４分割の光電子増倍管を用いたシンチレータアレイの読み取りから得られるフラッドヒストグラムを示す。図５において、３６個のピークは、シンチレータアレイに含まれる３６個のシンチレータそれぞれを示す。また、シンチレータアレイ内の最近傍のものを結ぶ線として見える、ピーク間に位置する相互作用イベントは、シンチレータアレイ内でコンプトン散乱したイベントから生じている。
図６は、アナログデータ束ね処理として位置感応型アバランシェフォトダイオードを用いたシンチレータアレイの読み取りから得られるフラッドヒストグラムを示す。図６において、６４個のピークは、８×８のシンチレーションアレイに含まれる６４個のシンチレータそれぞれを示す。また、鋭いピーク間に位置する相互作用イベントは、シンチレータアレイ内でコンプトン散乱したイベントから生じている。図６において、最近傍のピークを結ぶ線のパターンがはっきりと見える。これは、最近傍のシンチレータのペアにエネルギーを付与するコンプトン相互作用によるものである。
具体例として、図５に示すケースを検討する。このケースでは、３６個の個別のシンチレータからなる６×６のシンチレータアレイが、４分割の光電子増倍管を用いて読み取られている。４分割の光電子増倍管を用いてシンチレータアレイを読み取る代わりに、４分割ＳｉＰＭアレイを用いて読み取るケースを想定する。従来の方法を用いると、図５のフラッドヒストグラムは３６個の領域に分割され、各領域が１つのシンチレータを表す。エネルギー補正係数は、非線形補正を含むかもしれないが、シンチレータ毎に適用されるであろう。一実施形態では、フラッドヒストグラムが多数のサブピクセル（例えば、９００個）に分割される。次に、異なる非線形補正がキャリブレーションデータ（calibration data）から求められ、各サブピクセル領域に適用される。
別の実施形態では、キャリブレーションデータから各サブ領域のための補正を求めた後に、非線形補正をフラッドヒストグラム内の位置の関数として連続的に変化させる数学関数が導かれる。この場合、数学関数は、キャリブレーションで使用する複数領域の中心間を補間する方法とみなされる。
以下に説明する実施形態に係るＰＥＴ装置は、シンチレータアレイと、複数の光センサと、特定部と、格納部と、導出部とを備える。シンチレータアレイは、複数のシンチレータを有する。複数の光センサは、所定エネルギー値のガンマ線がシンチレータと相互作用することで生成されたシンチレーション光を検出する。特定部は、シンチレーション光の拡散により一つ又は複数の光センサにて検出され、光センサから出力された信号値に基づいて、ガンマ線のシンチレータとの相互作用イベントの検出位置を、シンチレータの数より多い数に分割された領域の単位で特定する。格納部は、光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、検出位置に関連付けて記憶部に格納する。導出部は、検出位置に関連付けて記憶部に格納された合計信号値と、所定エネルギー値とに基づいて、相互作用イベントのエネルギー値を補正する補正値を、領域の単位毎に導出する。
また、実施形態に係るＰＥＴ装置において、特定部は、複数の相互作用イベントそれぞれについて、光センサから出力された信号値に基づいて検出位置を特定してもよい。格納部は、複数の相互作用イベントそれぞれについて合計信号値を導出し、導出した合計信号値を、検出位置に関連付けて記憶部に格納してもよい。導出部は、検出位置に関連付けて記憶部に格納された複数の合計信号値から、領域の単位毎の平均値を更に導出し、導出した平均値と所定エネルギー値とに基づいて、相互作用イベントのエネルギー値を補正する補正値を、領域の単位毎に導出してもよい。
また、実施形態に係るＰＥＴ装置において、特定部は、信号値を出力した光センサが複数である場合に、各光センサの位置を、各光センサにて出力される信号値に応じて加重平均することで、検出位置を特定してもよい。また、導出部は、所定エネルギー値と合計信号値との比を、補正値として、領域の単位毎に導出してもよい。また、光センサは、少なくとも一つのシリコン光電子増倍管を含んでもよい。また、ＰＥＴ装置は、ＴＯＦ型であってもよい。
なお、特定部、格納部、及び導出部などの各部は、例えば、図３を用いて説明した演算部３７０に備えられる。
また、別の実施形態に係るＰＥＴ装置は、シンチレータアレイと、複数の光センサと、特定部と、補正部とを備える。シンチレータアレイは、複数のシンチレータを有する。複数の光センサは、所定エネルギー値のガンマ線がシンチレータに入力することで生成されたシンチレーション光を検出する。特定部は、シンチレーション光の拡散により一つ又は複数の光センサにて検出され、光センサから出力された信号値に基づいて、ガンマ線の相互作用イベントの検出位置を、シンチレータの数より多い数に分割された領域の単位で特定する。補正部は、光センサから出力された信号値から合計信号値を導出し、導出した合計信号値と、特定部によって特定された検出位置とを用いてキャリブレーションデータを参照し、検出位置に対応する補正値に基づいて、合計信号値を補正する。キャリブレーションデータは、相互作用イベントのエネルギー値を補正する補正値が領域の単位毎に定義されたものである。
また、別の実施形態に係るＰＥＴ装置において、特定部は、信号値を出力した光センサが複数である場合に、各光センサの位置を、各光センサにて出力される信号値に応じて加重平均することで、検出位置を特定してもよい。また、補正部は、合計信号値に補正値を乗算することによって、合計信号値を補正してもよい。また、ＰＥＴ装置は、ＴＯＦ型であってもよい。
なお、特定部及び補正部などの各部は、例えば、図３を用いて説明した演算部３７０に備えられる。
また、別の実施形態に係るＰＥＴ装置は、シンチレータアレイと、複数の光センサと、特定部と、格納部と、導出部とを備える。シンチレータアレイは、複数のシンチレータを有する。複数の光センサは、所定エネルギー値のガンマ線がシンチレータに入力することで生成されたシンチレーション光を検出する。特定部は、シンチレーション光の拡散により一つ又は複数の光センサにて検出され、光センサから出力された信号値に基づいて、ガンマ線の相互作用イベントの検出位置を、シンチレータの数より多い数に分割された領域の単位で特定する。格納部は、光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、検出位置に関連付けて記憶部に格納する。導出部は、検出位置に関連付けて記憶部に格納された合計信号値と所定エネルギー値とに基づいて、相互作用イベントのエネルギー値を補正する補正値が、検出位置の関数として連続的に変化する数学的関数を導出する。
なお、特定部、格納部、及び導出部などの各部は、例えば、図３を用いて説明した演算部３７０に備えられる。
また、別の実施形態に係るＰＥＴ装置は、シンチレータアレイと、複数の光センサと、特定部と、格納部と、エネルギーウィンドウ導出部とを備える。シンチレータアレイは、複数のシンチレータを有する。複数の光センサは、所定エネルギー値のガンマ線がシンチレータに入力することで生成されたシンチレーション光を検出する。特定部は、シンチレーション光の拡散により一つ又は複数の光センサにて検出され、光センサから出力された信号値に基づいて、ガンマ線の相互作用イベントの検出位置を、シンチレータの数より多い数に分割された領域の単位で特定する。格納部は、光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、検出位置に関連付けて記憶部に格納する。エネルギーウィンドウ導出部は、検出位置に関連付けて記憶部に格納された合計信号値と所定エネルギー値とに基づいて、エネルギーウィンドウを、領域の単位毎に導出する。
図７Ａは、実施形態に係るＰＥＴ装置によって実行される方法のステップを示すフローチャートである。図７Ａは、補正値（「補正係数」とも称する）のセットを決定する方法を示す。この補正値は、ガンマ線検出器によって検出された相互作用イベントのエネルギー値を決定するために用いられる。また、このガンマ線検出器は、シンチレータ素子からなるシンチレータアレイ上に配置された、１つ以上の非線形光センサを有する。また、このガンマ線検出器は、光学系束ね処理又はアナログデータ束ね処理を用いる。図７Ａに示す方法によって、正確なエネルギー値を決定するための補正値のセットが決定され、補正値のセットが決定されることによって、ガンマ線検出器がキャリブレーションされる。言い換えると、図７Ａに示す方法の目的は、例えば図１に示すように、光センサで生成され、積算された信号と、相互作用イベントのエネルギー値との間の非線形関係を、相互作用イベントの相対位置に対応するサブピクセルのセル位置毎に決定することである。セル位置の数は、シンチレータアレイに含まれるシンチレータの数よりも多い。
ステップＳ７１０において、実施形態に係るＰＥＴ装置は、第１のエネルギー値（例えば、５１１ｋｅＶ）を有するガンマ線を発生する線源から、検出器がガンマ線を受ける。かかるガンマ線は、２２Ｎａ、６８Ｇｅ、又は１８Ｆの放射線源から発生し得るものであり、その線源の線量は、典型的には、数十マイクロキューリから数ミリキューリの範囲である。線源の線量は、十分に高い計数率が得られるように、すなわち、適切な時間内（数十分から数時間）に必要なデータが収集されるように、選択される。また、線源の線量は、計数率が高過ぎることで相互作用イベントが重なり合い、この結果として重大な誤差が生じることがないように、選択される。
ステップＳ７２０において、ガンマ線がシンチレータ素子からなるシンチレータアレイと相互作用し、シンチレーション光が生成され、生成されたシンチレーション光は、１つ以上の光センサによって検出される。実施形態に係るデータ収集部は、相互作用イベントを契機とするエネルギー値として光センサから出力された信号値（任意単位）を、光センサ毎に収集する。例えば、図８Ａは、５１１ｋｅＶのガンマ線がシンチレータアレイに当たることより発生した、光の入射に反応した光センサからの典型的なパルスを示す。また、図８Ｂ−１は、図８Ａの時間軸を変更したものであり、図８Ｂ−２は、図８Ａに示すパルスに対応するエネルギーの情報を示すものである。図８Ｂ−１に示すように、パルスの立ち上がりからタイミングの情報が得られ、また、パルスの積分からエネルギーの情報が得られる。
ステップＳ７３０において、実施形態に係るＰＥＴ装置は、各光センサから出力された信号値に基づいて、シンチレータイベントの相対位置及び合計信号値を決定する。図９Ａ及び図９Ｂは、実施形態に係るガンマ線検出器の構成を示す。また、図９Ａ及び図９Ｂは、相互作用イベントの相対位置を決定する１つの方法を示す。図９Ａ及び図９Ｂは、ガンマ線検出器の２つの図であり、シンチレータ、ライトガイド、及び光センサの位置関係を示す。
１つのシンチレータからのシンチレーション光は、ライトガイドによって複数の光センサの上に拡散する。実施形態に係るＰＥＴ装置は、相互作用イベントの相対位置を、異なる光センサによって受信される信号の割合に基づいて、以下のように計算することができる。
ｘ＝（Σｘi Signali ）／ΣSignali ・・・（１）
ｙ＝（Σｙi Signali ）／ΣSignali ・・・（２）
ここで、Signalｉは、ｉ番目の光センサから出力された信号値であり、ｘｉ及びｙｉは、ｉ番目の光センサの中心の位置である。なお、ＰＥＴ装置は、複数の光センサによって検出される相互作用イベントの相対位置を、他のアルゴリズムを用いて決定することもできる。また、実施形態に係るＰＥＴ装置は、相互作用イベントのシンチレーション光を受けた各光センサにおける信号値の合計（合計信号値）から、相互作用イベントのエネルギー値の総量を計算することができる。
一実施形態においては、相対位置は、所定数のセル位置の１つに対応するものとして、決定される。ここで、このセル（サブピクセル）位置の所定数は、シンチレータ素子からなるシンチレータアレイの中のシンチレータ素子の数より多い。
ステップＳ７４０において、実施形態に係るＰＥＴ装置は、相互作用イベントに関する合計信号値を、ステップＳ７３０において決定された相互作用イベントの相対位置と関連付けて記憶部に格納する。
ステップＳ７５０において、実施形態に係るＰＥＴ装置は、ステップＳ７１０〜Ｓ７４０を複数の相互作用イベントについて繰り返す。図１０Ａは、フラッドヒストグラムを示す。また、図１０Ｂは、フラッドヒストグラムに対応するシンチレータＩＤルックアップテーブルを示す。図１０Ａに示すように、従来、相互作用イベントの頻度を示すフラッドヒストグラムは、相互作用イベント毎に格納されたデータに基づいて生成されている。フラッドヒストグラムは、シンチレータ毎に個別のピークを示す。また、シンチレータや光センサなどの配置によって、独特な歪みパターンを生じさせる。図１０Ｂに示すように、ＰＥＴ装置は、フラッドヒストグラムを領域に分割することにより、相互作用イベントが起こったと考えられるシンチレータを、フラッドヒストグラムを用いて特定することができる。同じ領域に該当する全ての相互作用イベントが、同じシンチレータに割り当てられる。こうして、ルックアップテーブルが生成され、ＰＥＴ装置は、フラッドヒストグラム内の相互作用イベントの相対位置を、シンチレータの位置に変換することができる。
一実施形態において、図１０Ａに示すフラッドヒストグラム内の各領域は、複数のセル（サブピクセル）に再分割される。実施形態に係るＰＥＴ装置は、相互作用イベントのエネルギー値を補正する補正値をサブピクセル毎に決定する。
ステップＳ７６０において、実施形態に係るＰＥＴ装置は、記憶された相互作用イベントのデータに基づいて、複数の相互作用イベントに関する合計信号値の平均値を、相対位置（サブピクセルのセル位置）毎に決定する。また、ＰＥＴ装置は、合計信号値の平均値を、相互作用イベントを生成するために用いられたガンマ線の第１のエネルギー値と関連付けて記憶部に格納する。
ステップＳ７７０において、実施形態に係るＰＥＴ装置は、相対位置毎に決定された合計信号値の平均値と第１のエネルギー値とに基づいて、補正値を相対位置毎に決定する。例えば、ＰＥＴ装置は、所定のエネルギー値（例えば、５１１ｋｅＶ）と、サブピクセルのセル位置に対する合計信号値の平均値（任意単位）との比を、そのセル位置に関する補正値として設定する。
ステップＳ７８０において、実施形態に係るＰＥＴ装置は、所定のエネルギー値を有する更に別のガンマ線について、ステップＳ７１０〜Ｓ７７０を繰り返し、各所定のエネルギー値、各セル位置に対応する補正値を決定する。こうして、ステップＳ７７０の後、キャリブレーションプロセスの最後に、実施形態に係るＰＥＴ装置は、各セル（サブピクセル）位置毎に、例えば、データ値｛Ｓｉ、Ｅｉ｝のセットを持つようになる。ここで、Ｓｉは、合計信号値（任意単位）であり、Ｅｉは、エネルギー値（ｋｅＶ）である。ＰＥＴ装置は、このデータ値のセットを用いて、例えば図１に示すような、補正されたエネルギー値に信号値を関連付ける非線形曲線を、セル位置毎に生成する。そして、ＰＥＴ装置は、セル位置毎に生成された非線形曲線を用いて、相互作用イベントの信号値をエネルギー値に変換することができ、ＳｉＰＭ応答の非線形を補償することができる。
上述したように、別の実施形態に係るＰＥＴ装置は、セル位置毎にキャリブレーションデータから補正値を導いた後に、フラッドヒストグラム内の位置の関数であって連続的に変化する数学的関数を非線形補正として導く。この場合、数学的関数は、キャリブレーションで使用する複数の領域の中心同士の間を補間する方法とみなされる。
他の実施形態においては、あるエネルギーのガンマ線を持つ線源を使用し、その解析を行っている最中に、同時に別のエネルギーについても解析を行う。このとき、ＰＥＴ装置は、コンプトンエッジや後方散乱ピークといった、取得されたスペクトル内に既知のエネルギーを見出すことのできる特徴的な点を用いて、ガンマ線検出器の非線形応答をキャリブレーションすることもできる。
図７Ｂは、別の実施形態に係るＰＥＴ装置によって実行される方法のステップを示すフローチャートである。別の実施形態に係るＰＥＴ装置は、図７Ｂに示すように、図７Ａに示される方法により得られる補正値を用いて、相互作用イベントのエネルギー値を補正する。
ステップＳ７１５において、実施形態に係るＰＥＴ装置の検出器は、ある線源から発生するガンマ線を受ける。
ステップＳ７２５において、ガンマ線がシンチレータ素子からなるシンチレータアレイに当たると、シンチレーション光が生成され、生成されたシンチレーション光は、１つ以上の光センサによって検出される。実施形態に係るデータ収集部は、ガンマ線との相互作用イベントを契機とするエネルギー値として光センサから出力された信号値（任意単位）を、光センサ毎に収集する。
ステップＳ７３５において、実施形態に係るＰＥＴ装置は、各光センサから出力された信号値に基づいて、ガンマ線との相互作用イベントの相対位置及び合計信号値を、上述したように決定する。相対位置は、所定数のセル位置の１つに対応するものとして、決定される。相対位置は、キャリブレーションデータが取得された多数のセル位置の中の対応するセル位置に変換される。セル（サブピクセル）位置の所定数は、シンチレータ素子からなるシンチレータアレイの中のシンチレータ素子の数より多い。
ステップＳ７４５において、実施形態に係るＰＥＴ装置は、相対位置に対応するセル位置毎に記憶された合計信号値（ステップＳ７３５において決定されたもの）、及び、キャリブレーションデータに基づいて、相互作用イベントに対する補正エネルギー値を計算する。ここで、キャリブレーションデータは、セル位置毎の信号値とエネルギー値との非線形関係を定義するものである。例えば、キャリブレーションデータが、実質的に、信号値とエネルギー値との関係を示す非線形曲線を定義している。一実施形態では、実施形態に係るＰＥＴ装置は、ステップＳ７３５において決定された合計信号値に、非線形曲線を用いて補正値を掛けることで、補正されたエネルギー値を得る。当業者には自明であるが、ここに記載の実施形態を適用する前に、この実施形態の範囲を逸脱することなく、更なるゲイン又はオフセット補正を信号に適用することもできる。
従来のシステムより優れた点は、ここに記載の実施形態によれば、より良い非線形補正をすることができ、エネルギー分解能を改善することができることである。ＰＥＴにとって、エネルギー分解能の改善は、散乱フラクションの低減になると言い換えることができ、最終的に、画像品質を改善することができる。
なお、実施形態に係るＰＥＴ装置は、エネルギーウィンドウを領域の単位毎に導出するエネルギーウィンドウ導出部を更に備えてもよい。エネルギーウィンドウ導出部は、検出位置に関連付けて記憶部に格納された合計信号値と所定エネルギー値とに基づいて、エネルギーウィンドウを領域の単位毎に導出する。ここで、ＰＥＴ装置は、エネルギーウィンドウを設け、信号値がエネルギーウィンドウ内の信号のみを処理の対象としてもよい。このエネルギーウィンドウは、例えば固定の上限値及び下限値などが予め設定されることで運用される。１つのシンチレータ内のどの位置にガンマ線が入射したかが、信号値に影響を与えることは、上述した通りである。そこで、実施形態に係るＰＥＴ装置は、領域の単位毎に格納された合計信号値に基づいて、例えば、合計信号値を中心とする所定幅のエネルギーウィンドウを、領域の単位毎に決定してもよい。
なお、エネルギーウィンドウ導出部は、検出位置に関連付けて記憶部に格納された合計信号値に基づいて、エネルギー分布曲線を生成し、生成したエネルギー分布曲線下の中央の面積から全体面積の所定割合までを合計することによって、エネルギーウィンドウを導出してもよい。
上述したように、一実施形態では、イベントのエネルギーを決定するために用いられる補正係数を決定するための方法が提供される。イベントは、ガンマ線検出器によって検出される。ガンマ線検出器は、シンチレータ素子からなるシンチレータアレイの上に配置された、少なくとも１つの非線形光センサを有する。ガンマ線検出器は、光学系束ね処理又はアナログデータ束ね処理を用いる。上述の方法は、（ａ）第１の所定のエネルギー値を有するガンマ線を生成するステップを含む。また、上述の方法は、（ｂ）少なくとも１つの非線形光センサによって生成される対応信号値を取得するステップを含む。少なくとも１つの非線形光センサのそれぞれは、イベントの発生を示す、生成されたガンマ線の１つの到着に応答して、シンチレータ素子からなるアレイの中の少なくとも１つのシンチレータから放出されたシンチレーション光を受けたことに応答して、対応信号値を生成する。また、上述の方法は、（ｃ）少なくとも１つの非線形光センサのそれぞれから取得した信号値に基づき、イベントの相対位置と合計信号値とを決定するステップを含む。この相対位置は、所定数のセル位置の１つであり、セル位置の所定数は、シンチレータ素子からなるアレイ内のシンチレータ素子の数よりも多い。また、上述の方法は、（ｄ）合計信号値を、決定されたセル位置と関連付けて記憶するステップを含む。また、上述の方法は、（ｅ）記憶されたイベントデータを生成するために、複数のイベントに対して、取得ステップ、決定ステップ、及び記憶ステップを繰り返すステップを含む。また、上述の方法は、（ｆ）セル位置毎に、記憶されたイベントデータに基づき、第１の所定のエネルギー値に対する平均合計信号値を決定するステップを含む。また、上述の方法は、（ｇ）セル位置毎に、決定された平均合計信号値と第１の所定のエネルギー値とに基づき補正係数を決定するステップを含む。
上述した実施形態の一態様によれば、上述の方法は、更に、（１）第２の所定のエネルギー値を有する第２のガンマ線を生成するステップを含む。また、上述の方法は、（２）セル位置毎に、第２の所定のエネルギー値に対応する第２の補正係数を決定するために、ステップ（ａ）〜（ｇ）を繰り返すステップを含む。また、上述の方法は、（３）セル位置毎に、各決定された補正値を、対応する所定のエネルギー値と関連付けて記憶するステップを含む。
別の態様によれば、上述の方法は、更に、セル位置毎に、記憶された補正係数に基づき、信号値とエネルギー値との非線形関係を決定するステップを含む。
更に別の態様によれば、相対位置を決定する上述したステップは、相対位置に対応するために所定数のセル位置の１つを決定するステップを含む。
更に別の態様によれば、相対位置を決定する上述したステップは、（１）イベントに対応するシンチレーション光を受ける、少なくとも１つの非線形光センサのそれぞれのｘ−ｙ位置に基づき、加重平均を計算するステップを含む。平均は、イベントに対応するシンチレーション光を受けた、少なくとも１つの非線形光センサのそれぞれの、対応する取得した信号値によって加重値を与えられる。また、相対位置を決定する上述したステップは、（２）少なくとも１つの非線形光センサのそれぞれの、取得した信号値を合計することにより、合計信号値を決定するステップを含む。
更に別の態様によれば、補正値を決定する上述したステップは、セル位置毎の第１の所定のエネルギー値に対応する補正値を、第１の所定のエネルギー値と、セル位置に対して決定された合計信号値との比として求めるステップを含む。
更に別の態様によれば、少なくとも１つの非線形光センサは、少なくとも１つのシリコン光電子増倍管を含む。
別の実施形態では、コンピュータプログラムを格納するコンピュータ可読媒体が提供され、このプログラムはコンピュータによって実行される。コンピュータは、イベントのエネルギーを決定するために用いられる補正係数を決定する。イベントは、ガンマ線検出器によって検出される。ガンマ線検出器は、シンチレータ素子からなるシンチレータアレイの上に配置された、少なくとも１つの非線形光センサを有する。ガンマ線検出器は、光学系束ね処理又はアナログデータ束ね処理を用いる。上述したプログラムにより上述したコンピュータは、（１）少なくとも１つの非線形光センサによって生成される対応信号値を受信するステップを実行する。少なくとも１つの非線形光センサのそれぞれは、イベントの発生を示す、生成されたガンマ線の１つの到着に応答して、シンチレータ素子からなるアレイの中の少なくとも１つのシンチレータから放出されたシンチレーション光を受けたことに応答して、対応信号値を生成する。また、コンピュータは、（２）少なくとも１つの非線形光センサのそれぞれから取得した信号値に基づき、イベントの相対位置と合計信号値とを決定するステップを実行する。この相対位置は、所定数のセル位置の１つであり、セル位置の所定数は、シンチレータ素子からなるアレイ内のシンチレータ素子の数よりも多い。また、コンピュータは、（３）合計信号値を、決定されたセル位置と関連付けて記憶するステップを実行する。また、コンピュータは、（４）記憶されたイベントデータを生成するために、複数のイベントに対して、取得ステップ、決定ステップ、及び記憶ステップを繰り返すステップを実行する。また、コンピュータは、（５）セル位置毎に、記憶されたイベントデータに基づき、第１の所定のエネルギー値に対する平均合計信号値を決定するステップを実行する。また、コンピュータは、（６）セル位置毎に、決定された平均合計信号値と第１の所定のエネルギー値とに基づき補正係数を決定するステップを実行する。
別の実施形態では、イベントのエネルギーを決定するために用いられる補正係数を決定するための方法が提供される。イベントは、ガンマ線検出器によって検出される。ガンマ線検出器は、シンチレータ素子からなるシンチレータアレイの上に配置された、少なくとも１つの非線形光センサを有する。上述の方法は、（１）複数のガンマ線を生成するステップを含む。また、上述の方法は、（２）少なくとも１つの非線形光センサによって生成される対応信号値を取得するステップを含む。少なくとも１つの非線形光センサのそれぞれは、イベントの発生を示す、生成されたガンマ線の１つの到着に応答して、シンチレータ素子からなるアレイの中の少なくとも１つのシンチレータから放出されたシンチレーション光を受けたことに応答して、対応信号値を生成する。また、上述の方法は、（３）少なくとも１つの非線形光センサのそれぞれから取得した信号値に基づき、イベントの相対位置と合計信号値とを決定するステップを含む。この相対位置は、所定数のセル位置の１つであり、セル位置の所定数は、シンチレータ素子からなるアレイ内のシンチレータ素子の数よりも多い。また、上述の方法は、（４）決定された相対位置に対応するセル位置に対する合計信号値と、記憶されたキャリブレーションデータとに基づき、イベントに対する補正エネルギー値を計算するステップを含む。記憶されたキャリブレーションデータは、セル位置毎の信号値とエネルギー値との非線形関係を定義する。
別の実施形態では、コンピュータプログラムを格納するコンピュータ可読媒体が提供され、このプログラムはコンピュータによって実行される。コンピュータは、シンチレータ素子からなるシンチレータアレイの上に配置された、少なくとも１つの非線形光センサを有するガンマ線検出器によって検出されるイベントに対応する補正エネルギー値を決定する。上述したプログラムにより上述したコンピュータは、（１）少なくとも１つの非線形光センサによって生成される対応信号値を受信するステップを実行する。少なくとも１つの非線形光センサのそれぞれは、イベントの発生を示す、生成されたガンマ線の１つの到着に応答して、シンチレータ素子からなるアレイの中の少なくとも１つのシンチレータから放出されたシンチレーション光を受けたことに応答して、対応信号値を生成する。また、コンピュータは、（２）少なくとも１つの非線形光センサのそれぞれから取得した信号値に基づき、イベントの相対位置と合計信号値とを決定するステップを実行する。この相対位置は、所定数のセル位置の１つであり、セル位置の所定数は、シンチレータ素子からなるアレイ内のシンチレータ素子の数よりも多い。また、コンピュータは、（３）決定された相対位置に対応するセル位置に対する合計信号値と、記憶されたキャリブレーションデータとに基づき、イベントに対する補正エネルギー値を計算するステップを実行する。記憶されたキャリブレーションデータは、セル位置毎の信号値とエネルギー値との非線形関係を定義する。
別の実施形態では、イベントのエネルギーを決定するために用いられる補正係数を決定するための方法が提供される。イベントは、ガンマ線検出器によって検出される。ガンマ線検出器は、シンチレータ素子からなるシンチレータアレイの上に配置された、少なくとも１つの非線形光センサを有する。上述の方法は、（ａ）第１の所定のエネルギー値を有するガンマ線を生成するステップを含む。また、上述の方法は、（ｂ）少なくとも１つの非線形光センサによって生成される対応信号値を取得するステップを含む。少なくとも１つの非線形光センサのそれぞれは、イベントの発生を示す、生成されたガンマ線の１つの到着に応答して、シンチレータ素子からなるアレイの中の少なくとも１つのシンチレータから放出されたシンチレーション光を受けたことに応答して、対応信号値を生成する。また、上述の方法は、（ｃ）少なくとも１つの非線形光センサのそれぞれから取得した信号値に基づき、イベントの相対位置と合計信号値とを決定するステップを含む。この相対位置は、所定数のセル位置の１つであり、セル位置の所定数は、シンチレータ素子からなるアレイ内のシンチレータ素子の数よりも多い。また、上述の方法は、（ｄ）合計信号値を、決定されたセル位置と関連付けて記憶するステップを含む。また、上述の方法は、（ｅ）記憶されたイベントデータを生成するために、複数のイベントに対して、取得ステップ、決定ステップ、及び記憶ステップを繰り返すステップを含む。また、上述の方法は、（ｆ）記憶されたイベントデータと、第１の所定のエネルギー値とに基づき、連続的に変化する数学的補正関数のパラメータを、シンチレータ素子からなるアレイ内の位置の関数として決定するステップを含む。数学的補正関数は、取得したイベントの合計信号値における空間的変化を表す。
更に、この実施形態の上述した方法は、イベント毎に、イベントの決定された数学的補正関数と、決定された相対位置とを用いて、イベントのエネルギーを決定するステップを含む。
以上述べた少なくとも一つの実施形態の陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法によれば、エネルギー分解能を向上することが可能になる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
１００ シンチレータ
１０５ シンチレータアレイ
１３０ ライトガイド
１３５、１４０ 光センサ
３００ シンチレータ
３０５ シンチレータアレイ
３３０ ライトガイド
３３５、３４０ 光センサ

Claims (18)

  1. 複数のシンチレータを有するシンチレータアレイと、
    所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサと、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線のシンチレータとの相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定部と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納部と、
    前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値を、前記領域の単位毎に導出する導出部と
    を備えたことを特徴とする陽電子放出コンピュータ断層撮影装置。
  2. 前記特定部は、複数の前記相互作用イベントそれぞれについて、前記光センサから出力された信号値に基づいて前記検出位置を特定し、
    前記格納部は、複数の前記相互作用イベントそれぞれについて、前記合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納し、
    前記導出部は、前記検出位置に関連付けて記憶部に格納された複数の合計信号値から、前記領域の単位毎の平均値を更に導出し、導出した平均値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値を、前記領域の単位毎に導出することを特徴とする請求項１に記載の陽電子放出コンピュータ断層撮影装置。
  3. 前記特定部は、前記信号値を出力した光センサが複数である場合に、各光センサの位置を、各光センサにて出力される信号値に応じて加重平均することで、前記検出位置を特定することを特徴とする請求項１又は２に記載の陽電子放出コンピュータ断層撮影装置。
  4. 前記光センサは、少なくとも一つのシリコン光電子増倍管を含むことを特徴とする請求項１〜３のいずれか一つに記載の陽電子放出コンピュータ断層撮影装置。
  5. ＴＯＦ（Time Of Flight）型であることを特徴とする請求項１〜４のいずれか一つに記載の陽電子放出コンピュータ断層撮影装置。
  6. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置に実行させるためのプログラムであって、
    所定エネルギー値のガンマ線が前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値を収集する収集手順と、
    前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定手順と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納手順と、
    前記収集手順と前記特定手順と前記格納手順とを繰り返す繰り返し手順と、
    前記検出位置に関連付けて記憶部に格納された合計信号値から、前記領域の単位毎の平均値を導出する平均値導出手順と、
    前記平均値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値を、前記領域の単位毎に導出する補正値導出手順と
    を実行させるためのプログラム。
  7. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置によって実行される方法であって、
    所定エネルギー値のガンマ線が前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値を収集する収集工程と、
    前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定工程と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納工程と、
    前記収集工程と前記特定工程と前記格納工程とを繰り返す繰り返し工程と、
    前記検出位置に関連付けて記憶部に格納された合計信号値から、前記領域の単位毎の平均値を導出する平均値導出工程と、
    前記平均値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値を、前記領域の単位毎に導出する補正値導出工程と
    を含んだことを特徴とする方法。
  8. 前記所定エネルギー値である第１エネルギー値とは異なる第２エネルギー値について、前記収集工程と、前記特定工程と、前記格納工程と、前記繰り返し工程と、前記平均値導出工程と、前記補正値導出工程とを繰り返し、前記第２エネルギー値に対応する補正値を、前記領域の単位毎に導出する第２補正値導出工程を更に含んだことを特徴とする請求項７に記載の方法。
  9. 複数のシンチレータを有するシンチレータアレイと、
    所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサと、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定部と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値と、前記特定部によって特定された検出位置とを用いて、前記相互作用イベントのエネルギー値を補正する補正値が前記領域の単位毎に定義されたキャリブレーションデータを参照し、前記検出位置に対応する補正値に基づいて、前記合計信号値を補正する補正部と
    を備えたことを特徴とする陽電子放出コンピュータ断層撮影装置。
  10. 前記特定部は、前記信号値を出力した光センサが複数である場合に、各光センサの位置を、各光センサにて出力される信号値に応じて加重平均することで、前記検出位置を特定することを特徴とする請求項９に記載の陽電子放出コンピュータ断層撮影装置。
  11. 前記補正部は、前記合計信号値に前記補正値を乗算することによって、前記合計信号値を補正することを特徴とする請求項９又は１０に記載の陽電子放出コンピュータ断層撮影装置。
  12. ＴＯＦ型であることを特徴とする請求項９〜１１のいずれか一つに記載の陽電子放出コンピュータ断層撮影装置。
  13. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置に実行させるためのプログラムであって、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定手順と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値と、前記特定手順によって特定された検出位置とを用いて、前記相互作用イベントのエネルギー値を補正する補正値が前記領域の単位毎に定義されたキャリブレーションデータを参照し、前記検出位置に対応する補正値に基づいて、前記合計信号値を補正する補正手順と
    を実行させるためのプログラム。
  14. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置によって実行される方法であって、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定工程と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値と、前記特定工程によって特定された検出位置とを用いて、前記相互作用イベントのエネルギー値を補正する補正値が前記領域の単位毎に定義されたキャリブレーションデータを参照し、前記検出位置に対応する補正値に基づいて、前記合計信号値を補正する補正工程と
    を含んだことを特徴とする方法。
  15. 複数のシンチレータを有するシンチレータアレイと、
    所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサと、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定部と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納部と、
    前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値が、検出位置の関数として連続的に変化する数学的関数を、補間を行うことにより導出する導出部と
    を備えたことを特徴とする陽電子放出コンピュータ断層撮影装置。
  16. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置に実行させるためのプログラムであって、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定手順と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納手順と、
    前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値が、検出位置の関数として連続的に変化する数学的関数を、補間を行うことにより導出する導出手順と
    を実行させるためのプログラム。
  17. 複数のシンチレータを有するシンチレータアレイと、所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサとを含むガンマ線検出器を備えた陽電子放出コンピュータ断層撮影装置によって実行される方法であって、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定工程と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納工程と、
    前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値との比に基づいて、前記相互作用イベントのエネルギー値を補正する補正値が、検出位置の関数として連続的に変化する数学的関数を、補間を行うことにより導出する導出工程と
    を含んだことを特徴とする方法。
  18. 複数のシンチレータを有するシンチレータアレイと、
    所定エネルギー値のガンマ線が前記シンチレータに入力することで生成されたシンチレーション光を検出する複数の光センサと、
    前記シンチレーション光の拡散により一つ又は複数の前記光センサにて検出され、前記光センサから出力された信号値に基づいて、前記ガンマ線の相互作用イベントの検出位置を、前記シンチレータの数より多い数に分割された領域の単位で特定する特定部と、
    前記光センサから出力された信号値から合計信号値を導出し、導出した合計信号値を、前記検出位置に関連付けて記憶部に格納する格納部と、
    前記検出位置に関連付けて記憶部に格納された合計信号値と前記所定エネルギー値とに基づいて生成されたエネルギー分布曲線に基づいて、エネルギーウィンドウを、前記領域の単位毎に導出するエネルギーウィンドウ導出部と
    を備えたことを特徴とする陽電子放出コンピュータ断層撮影装置。
JP2012029889A 2011-03-11 2012-02-14 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法 Active JP5984246B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210061629.8A CN102670232B (zh) 2011-03-11 2012-03-09 正电子发射计算机断层摄影装置、以及通过它执行的方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/045,610 2011-03-11
US13/045,610 US8294110B2 (en) 2011-03-11 2011-03-11 Method for improved correction of SiPM non-linearity in multiplexed radiation detectors

Publications (2)

Publication Number Publication Date
JP2012189583A JP2012189583A (ja) 2012-10-04
JP5984246B2 true JP5984246B2 (ja) 2016-09-06

Family

ID=46794667

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012029889A Active JP5984246B2 (ja) 2011-03-11 2012-02-14 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法

Country Status (2)

Country Link
US (1) US8294110B2 (ja)
JP (1) JP5984246B2 (ja)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8884213B2 (en) * 2011-07-06 2014-11-11 Siemens Medical Solutions Usa, Inc. Energy correction for one-to-one coupled radiation detectors having non-linear sensors
US9063087B2 (en) * 2012-03-21 2015-06-23 Shimadzu Corporation Radiation detector signal processor and radiation detector provided therewith
JP6017916B2 (ja) * 2012-10-16 2016-11-02 株式会社豊田中央研究所 光検出器
US9164144B2 (en) 2012-12-27 2015-10-20 General Electric Company Characterization and calibration of large area solid state photomultiplier breakdown voltage and/or capacitance
JP6325650B2 (ja) * 2013-04-24 2018-05-16 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 光検出器画素アレイを有する放射線検出装置、ｐｅｔ撮像システム、光フォトン加算方法、及び光フォトン加算方法のステップを実行するためのコンピュータプログラム
JP2014241543A (ja) * 2013-06-12 2014-12-25 株式会社東芝 光検出装置およびｃｔ装置
CN105637386B (zh) * 2013-10-14 2019-01-15 皇家飞利浦有限公司 在正电子发射断层摄影（pet）能量直方图中的直方图平滑
KR101646651B1 (ko) * 2014-10-30 2016-08-08 서강대학교산학협력단 다중 문턱전압을 이용한 의료 영상기기의 신호처리 시스템 및 방법
US10527741B2 (en) 2015-04-07 2020-01-07 Siemens Medical Solutions Usa, Inc. Setup of SIPM based PET detector using LSO background radiation
EP3351970A4 (en) * 2015-09-17 2019-06-26 Shimadzu Corporation INFORMATION PROCESSING DEVICE, RADIATION DETECTOR, RADIATION IMAGING APPARATUS AND PROGRAM
JP6608241B2 (ja) * 2015-10-28 2019-11-20 浜松ホトニクス株式会社 放射線位置検出器、ｐｅｔ装置、プログラム及び記録媒体
US10371835B2 (en) 2016-01-11 2019-08-06 General Electric Company Microcell interconnection in silicon photomultipliers
US9817134B1 (en) 2016-09-12 2017-11-14 General Electric Company Universal readout for light-sharing detector
JP6737154B2 (ja) 2016-12-02 2020-08-05 株式会社島津製作所 放射線検出装置
US10466368B2 (en) 2017-02-02 2019-11-05 General Electric Company Systems and methods for detector output adjustment
WO2018202878A1 (en) * 2017-05-04 2018-11-08 Koninklijke Philips N.V. Positron emission tomography (pet) timing calibration using coincidences involving high-energy cascade gamma from prompt-gamma positron emitters
CN109212582B (zh) * 2018-08-14 2020-11-27 北京永新医疗设备有限公司 基于多层闪烁晶体探测器的***均匀性校正方法及装置
US10310098B1 (en) * 2018-10-05 2019-06-04 Canon Medical Systems Corporation Determine position of scattered events in pixelated gamma detector using inverse energy weighting
EP3864440A4 (en) * 2018-10-09 2022-05-25 Shenzhen Xpectvision Technology Co., Ltd. METHODS AND SYSTEMS FOR FORMING IMAGES BY RADIATION
US10962662B2 (en) 2018-12-06 2021-03-30 Canon Medical Systems Corporation Method and apparatus to calibrate timing of a detector using multi-channel events
KR20210119279A (ko) 2019-02-15 2021-10-05 더 리서치 파운데이션 포 더 스테이트 유니버시티 오브 뉴욕 유사 각기둥 광 가이드 어레이를 갖는 고분해능 깊이-인코딩 pet 검출기
US11255985B2 (en) 2019-05-31 2022-02-22 Canon Medical Systems Corporation Method and apparatus to use a broad-spectrum energy source to correct a nonlinear energy response of a gamma-ray detector
CN112014870A (zh) * 2019-05-31 2020-12-01 佳能医疗***株式会社 放射线检测装置、能量校正方法以及程序
US11789130B2 (en) * 2020-02-13 2023-10-17 Sensors Unlimited, Inc. Detection pixels and pixel systems

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4212061A (en) * 1977-12-21 1980-07-08 Medtronic, Inc. Radiation signal processing system
US4323977A (en) * 1979-11-20 1982-04-06 Siemens Gammasonics, Inc. Non-uniformity energy correction method and apparatus
JPH08160141A (ja) * 1994-12-09 1996-06-21 Toshiba Corp シンチレーションカメラ
US5585637A (en) * 1995-06-09 1996-12-17 Adac Laboratories Multi-head nuclear medicine camera for dual SPECT and PET imaging
JPH10197641A (ja) * 1996-12-30 1998-07-31 Shimadzu Corp 放射線カメラ
US7709805B2 (en) * 2003-08-01 2010-05-04 General Electric Company Method for generating optical anisotropy in scintillators using pulsed lasers
US20050023473A1 (en) * 2003-08-01 2005-02-03 Burr Kent Charles System and method for reducing optical crosstalk in multi-anode photomultiplier tube
US20060011853A1 (en) * 2004-07-06 2006-01-19 Konstantinos Spartiotis High energy, real time capable, direct radiation conversion X-ray imaging system for Cd-Te and Cd-Zn-Te based cameras
US7411197B2 (en) * 2005-09-30 2008-08-12 The Regents Of The University Of Michigan Three-dimensional, position-sensitive radiation detection
EP1949135A2 (en) * 2005-10-17 2008-07-30 Koninklijke Philips Electronics N.V. Pmt gain and energy calibrations using lutetium background radiation
US7193208B1 (en) * 2005-10-24 2007-03-20 General Electric Company Time-of-flight capable high resolution pet detector

Also Published As

Publication number Publication date
US8294110B2 (en) 2012-10-23
US20120228484A1 (en) 2012-09-13
JP2012189583A (ja) 2012-10-04

Similar Documents

Publication Publication Date Title
JP5984246B2 (ja) 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法
US8907290B2 (en) Methods and systems for gain calibration of gamma ray detectors
EP2867701B1 (en) Digital positron emission tomography (dpet) energy calibration method
JP6615432B2 (ja) イメージングシステムおよび校正方法
US9029786B2 (en) Nuclear medicine imaging apparatus, and nuclear medicine imaging method
CN102670232B (zh) 正电子发射计算机断层摄影装置、以及通过它执行的方法
US9474501B2 (en) Hybrid method based on simulation and experimental data to normalize pet data
JP2014158713A (ja) コンピュータ断層撮影装置、校正プログラム、及び光子数校正装置
US8558181B2 (en) Positron emission tomography system with hybrid detection geometries and sampling
JP6188418B2 (ja) 偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体
JP7317586B2 (ja) 医用画像処理装置、方法及びプログラム
US9091771B2 (en) System and method for improving detection of gamma interactions in a positron emission tomography system
JP6125309B2 (ja) 偶発同時計数推定方法及び偶発同時計数推定装置
JP6912304B2 (ja) 波高頻度分布取得装置、波高頻度分布取得方法、波高頻度分布取得プログラム及び放射線撮像装置
JP6054050B2 (ja) 核医学イメージング方法、核医学イメージング装置及び記憶媒体
JP4594855B2 (ja) 核医学診断装置及び放射線カメラ並びに核医学診断装置における放射線検出方法
US9207334B1 (en) Methods and systems for a light sensor in gamma ray detectors
US8288728B2 (en) Method and apparatus to facilitate crystal identification in a pet detector
JP5632221B2 (ja) 可変なライトガイド厚を有するガンマ線検出モジュールおよびガンマ線スキャナ・システム
JP2024515639A (ja) 陽電子放射断層撮影デバイスを修正するためのシステムおよび方法
JP6825706B2 (ja) 放射線検出装置およびそれを備えた核医学診断装置
JP5794765B2 (ja) 核医学イメージング装置および核医学イメージング装置の作動方法
Goertzen et al. PET instrumentation
WO2013168778A1 (ja) 偶発同時計数推定方法及び偶発同時計数推定装置
JP2015230194A (ja) 位置同定型検出器を用いた放射線撮像装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150205

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20151102

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160113

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160301

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160418

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20160513

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160705

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160729

R150 Certificate of patent or registration of utility model

Ref document number: 5984246

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350