JP2017521676A - 電離放射線のピクセル化検出器からの信号の解析 - Google Patents
電離放射線のピクセル化検出器からの信号の解析 Download PDFInfo
- Publication number
- JP2017521676A JP2017521676A JP2017513366A JP2017513366A JP2017521676A JP 2017521676 A JP2017521676 A JP 2017521676A JP 2017513366 A JP2017513366 A JP 2017513366A JP 2017513366 A JP2017513366 A JP 2017513366A JP 2017521676 A JP2017521676 A JP 2017521676A
- Authority
- JP
- Japan
- Prior art keywords
- pixel
- crosstalk
- signal
- pixels
- compensation
- 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.)
- Pending
Links
- 230000005865 ionizing radiation Effects 0.000 title claims abstract description 19
- 238000004458 analytical method Methods 0.000 title claims description 29
- 239000002245 particle Substances 0.000 claims abstract description 40
- 239000004065 semiconductor Substances 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims description 37
- 238000001228 spectrum Methods 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 7
- 230000001360 synchronised effect Effects 0.000 claims description 5
- 238000004846 x-ray emission Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 14
- 238000012545 processing Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 7
- 230000004907 flux Effects 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 229910052732 germanium Inorganic materials 0.000 description 4
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 description 4
- 238000004876 x-ray fluorescence Methods 0.000 description 4
- 229910003460 diamond Inorganic materials 0.000 description 3
- 239000010432 diamond Substances 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 229910052750 molybdenum Inorganic materials 0.000 description 2
- 239000011733 molybdenum Substances 0.000 description 2
- 238000000441 X-ray spectroscopy Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002800 charge carrier Substances 0.000 description 1
- 229910052729 chemical element Inorganic materials 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000192 extended X-ray absorption fine structure spectroscopy Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000001506 fluorescence spectroscopy Methods 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 238000009681 x-ray fluorescence measurement Methods 0.000 description 1
Images
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/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
- G01T1/366—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry with semi-conductor detectors
-
- 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/247—Detector read-out circuitry
-
- 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/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2921—Static instruments for imaging the distribution of radioactivity in one or two dimensions; Radio-isotope cameras
- G01T1/2928—Static instruments for imaging the distribution of radioactivity in one or two dimensions; Radio-isotope cameras using solid state detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
- G01N23/223—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material by irradiating the sample with X-rays or gamma-rays and by measuring X-ray fluorescence
-
- 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/243—Modular detectors, e.g. arrays formed from self contained units
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measurement Of Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本発明は、半導体検出器のピクセルに入射する光子または電離放射線の粒子のエネルギーを表す信号に関する。異なるピクセルからの信号間のクロストークは、クロストークを生じさせるイベントのタイミングに従って、クロストーク信号にアライメントされて適用される、タイム・ドメイン級数または関数の形式のクロストーク補償シグネチャを用いて補償される。【選択図】図4
Description
本発明は、X線等の電離放射線を検出する半導体検出器から受信した、特に、処理信号に複数のピクセル(画素)を有するような半導体検出器から受信した信号の処理に関する。本発明はまた、上記のような処理を実行するための装置、及び上記のような検出器とそのような装置とを組み合わせた構成に関する。単なる例として、本発明は、シンクロトロン光源を使用するX線分光法の分野で使用されてもよく、検出器は、モノリシック多元素高純度ゲルマニウム検出器であってもよい。
モノリシックでセグメント化された高純度ゲルマニウム(HPGe)検出器は、高分解能分光法のためのシンクロトロン光源で使用される。この種の光源の増加する輝度を活用するために、ビームラインは、X線のフラックス(線束)が増加するにつれて、飽和状態にしない、またはエネルギー分解能を失う検出器を必要とする。フラックスの増加に追随するために、検出器システムの性能に使用される2つの主なアプローチがある。それらは、与えられた立体角内のピクセル数を増加させ、性能低下なしに、与えられた数のピクセルがピクセルカウントレート当たり、より高く働くことを可能にするものである。
例えば、英国に拠点を置くダイヤモンドライトソース(Diamond Light Source)は、各ピクセルからの信号を解析し、入射するX線のスペクトルを導出可能なデジタル信号処理エレクトロニクスに結合される社内設計の高速低ノイズプリアンプを装備した、多数のピクセルを備えたモノリシックセグメント化HPGe検出器を使用して、これら両方のアプローチを主要な分光ビームラインに組み合わせる。
X線蛍光分光法の重要な仕様の2つ−エネルギー分解能と、ピーク対バックグラウンド比−は、1ピクセル当たりのカウント速度が増加するにつれて著しく低下し、エネルギーが増加するにつれて顕著になることが観察されている。
本発明は、関連する先行技術のこれら及び他の問題に対処することを求めている。
本発明者らは、ピクセル化半導体検出器における異なるピクセルからの電気信号の精密な検査が、隣接するピクセル上で生じるX線相互作用から他のピクセルへの考慮すべきクロストークを示し、性能低下につながることを観察した。したがって、本発明は、モノリシックセグメント化HPGe検出器及び他の電離放射線の検出器のピクセル間で発生するピクセル間クロストークを補正する信号処理アルゴリズムを提供する。本発明者らは、ピクセル間クロストークが、例えば高輝度分光ビームラインなどの蛍光分光検出器でのレート分解能の低下を招くことを発見した。これにより、本発明者らは、デジタルパルスプロセッサ(DPP)によってデジタル化されたHPGeピクセル検出器からX線蛍光信号を取得した。これらの信号は、同じDPPのエミュレーションに適用された。エネルギー正規化されたイベントクロストークシグネチャまたは各隣接ピクセルペアに固有のプロファイルを生成することを可能にする、ピクセル間クロストークを識別し、特徴付けるために、低い放射線発生率の実験データが使用された。
特に、これらのクロストークシグネチャまたはプロファイルは、メインピクセルにおけるイベントのタイミングに一致するようにアライメント(調整)または同期された後に、加算または減算などの重ね合わせを用いて隣接ピクセルからの信号に適用されるべく、タイム・ドメイン関数(時間的に連続的であってもよく、時間的に離散的であってもよい)として生成することができる。クロストークシグネチャはまた、イベントの高さ及び/または信号の他のプロパティ、特にメインピクセルにおけるイベントのかかるプロパティに応じてスケーリングされてもよい(そして、メインイベントの時間的に近位にある他のそのようなイベントをメインピクセルに任意に含む)。
言い換えれば、各クロストークシグネチャは、別のピクセルにおけるイベントによるクロストークの影響を受ける1つ以上のピクセルの信号において予想されるクロストーク信号の形状またはプロファイルを表す。この形状またはプロファイルは、イベントのプロパティ(及び任意選択的に他のイベントのプロパティ、例えば同一または他のピクセルの時間的に近似した他の検出されたイベント)に従ってスケーリングまたは他の方法で修正され、クロストークに対して補償するためにクロストークに影響されたピクセルに適用され得る。
広範囲の計数率をカバーするさらなる検出器信号は、その後、DPPエミュレーションを通過し、最初に、粗いイベント高さ推定値及びピクセルごとにイベント当たりのタイムスタンプを取得する。イベントのタイムスタンプを用いて、未処理の隣接ピクセルデータの遅延バージョンからのこれら固有の再スケーリングされた信号の減算、加算または他の組み合わせまたは重ね合わせを同期させる前に、イベントの高さは、ピクセルの正規化されたクロストークシグネチャをその近傍にスケーリングするために使用された。これにより、1回のパスで単一のピクセル上に多くのクロストークアーチファクトが効果的に排除された。その後、信号解析をさらに最適化できる、クロストークからファントムイベントトリガを低減するクロストーク補正データの第2パス処理が用いられ得る。
クロストーク補償(補正)されたデータは、追加の処理アルゴリズムによってさらに処理された。クロストーク補正なしで処理された同一の信号と比較すると、このエミュレーションの出力は、典型的にエネルギー分解能が35%向上することを示し、ピーク対バックグラウンド比が倍増した。
したがって、本発明は、電離放射線の半導体検出器から信号を受信するように構成されたアナライザを提供し、半導体検出器は、複数のピクセルを備え、各ピクセルは前記ピクセルに入射する放射線の電離粒子(イオン化粒子)または光子に対するピクセルの応答を表す一連のイベントフィーチャを含む信号を出力するように構成され、前記アナライザによって受信された信号は、相互間にクロストークを呈し、それによって前記イベントフィーチャまたは少なくとも一部のイベントフィーチャが、少なくとも1つの他のピクセル、例えば隣接ピクセルからの信号に対応するクロストークフィーチャを生じさせる。
典型的には、イベントフィーチャは、特徴的な形状またはプロファイルを有することになり、例えば、それらはピクセル信号におけるステップフィーチャであってもよく、各ステップのサイズまたは高さは、光子または電離放射線の粒子のエネルギーが前記ステップフィーチャを生じさせることを示す。このようにして、前記ピクセル信号は入射粒子に関するスペクトル情報を含む。前記検出器及び他の効果に依存して、この種のステップフィーチャは、典型的には、図4に見られるような信号の後続ステップ移動のような特徴の周囲の特性を含む、持続時間が数十から数百ナノ秒のオーダーである。次いで、他のピクセルからの電気信号におけるクロストーク信号は、同様の時間範囲である可能性が高く、単一のピークまたは2つ以上のピークの振動のような、様々な形状、プロファイルまたは形状をとることができる。
アナライザは、1つ以上のクロストーク補償またはキャンセルシグネチャ(クロストーク補償プロファイルとも呼ばれることがある)を格納するか、含むか、または備えられており、各ピクセルについて、そのピクセルからの信号におけるイベントフィーチャを検出するように構成され、1つ以上のクロストーク補償またはキャンセルシグネチャまたはプロファイルのうちの選択されたものを使用して、1つ以上の他のピクセルからのコンカレント(並行)信号の解析におけるクロストークを補償またはキャンセルするように構成される。これは、クロストークのための他のピクセル信号の初期補正、それに続く信号の必要な特性の検出のための解析、必要な特性を検出するための解析の一部としてのシグネチャの適用、これらの組み合わせ、または 他の方法によって、達成することができる。
クロストーク補償/キャンセルシグネチャは、例えば、データポイントの時系列として、解析関数として、または以下に説明する他の方法で提供され得る。各シグネチャは、別のピクセルにおけるイベントによるクロストークの予想される形状またはプロファイルに関連する形状またはプロファイルを有する。言い換えれば、クロストーク補償シグネチャは、典型的には、時間ベースの関数またはタイム・ドメイン関数であり、クロストークのキャンセルまたは低減を求めるために、隣接または他のピクセルにおける信号との(例えば重ね合わせによる)組み合わせを目的として、クロストークを引き起こすイベントフィーチャとアライメントまたは同期させることができる。このアライメントは、イベントのタイミングの検出と、そのイベントへのシグネチャの単純なアライメントによって、あるいは例えばクロストーク信号への補償シグネチャの最良適合を求めることによって、他のまたは追加の方法で実行され得る。
典型的には、クロストーク補償シグネチャまたはプロファイルは、イベントフィーチャの1つ以上のプロパティに関して正規化された使用形態で提供され、アナライザは、前記検出された各イベントフィーチャの前記1つ以上のプロパティを測定するように構成されるとともに、前記並行信号の補償解析で用いる前記測定されたプロパティに従って前記一つ以上のクロストーク補償シグネチャの前記選択されたものをスケーリングするように構成される。典型的には、スケーリングは、イベントフィーチャの大きさまたは高さを用いて実行され得、イベントフィーチャを生じる入射光子または粒子のエネルギーに通常は対応するかまたは直接関連してもよい。しかし、補償すべきイベントのプロパティ、及び/またはそのイベントに時間的に近い他のイベントをかかるスケーリングに使用することができる。例えば、補償すべきイベントに近接する複数のイベントの数及び/または近接度の尺度を、かかるスケーリングに使用することができる。
補償シグネチャのスケーリングは、シグネチャと補償すべき信号とのコンビネーション(重ね合わせなど)のプロセスの一部として実行されてもよく、またはこの組み合わせの前に行われてもよい。
クロストーク補償/キャンセルシグネチャは、1つ、より多くの、またはすべての隣接するピクセル、または他のピクセルに関して使用するために、各ピクセルに対して少なくとも1つのかかるシグネチャを含むことができ、より具体的には、それぞれの異なる隣接ピクセルまたは他のピクセルのコンカレント(並行)電気信号におけるクロストークの補償/キャンセルに使用するための異なるシグネチャを含む、各ピクセルのシグネチャのセットを含むことができる。
補償/キャンセルされるクロストークは、半導体検出器それ自体の内部の相互作用から生じることがあり、この場合、クロストークは、半導体検出器において隣接するピクセルにおいてイベントが検出されたピクセルに主に見られると予想されるが、一部のクロストークはまた、非隣接のピクセルにおけるかかる影響からも見られる。本発明は、主として隣接ピクセル、または隣接ピクセルに関してのみ、及び/またはすべての隣接ピクセルに関して適用されるクロストークシグネチャおよび補償を提供することができ、または他の非隣接ピクセルにも同様に適用することができる。しかしながら、補償されるクロストークは、半導体センサとアナライザとの間の電気経路、例えば、プリアンプ、導体、及び/またはこれらの経路の他の部分の間における電気経路にも生じる可能性がある。電気経路のかかる態様の幾何学的配置は、検出器におけるピクセルのレイアウトに対応しない場合があるので、これらの経路間のクロストークは、半導体検出器内で隣接するピクセルに関して同様に支配的ではないか、または単独で見られることはない。したがって、本発明はまた、隣接していてもいなくてもよいピクセル対に一般に適用されるクロストークシグネチャおよび補償を提供する。
各クロストーク補償シグネチャは、特定のピクセル対に適用されてもよく、特に、シグネチャが、ペアの第1要素の信号における検出されたイベントのために、ペアの第2要素の信号を補償するために使用されるという意味で、有向対に適用されてもよい。これらの目的のために、第1ピクセルと呼ぶことができる検出器の各ピクセルに対して、かかるペアのセットが提供されてもよく、セットの各ペアは、第1ピクセルの検出されたイベントに関して異なる第2ピクセルからの信号のクロストーク補償のために使用され、前記セットは、前記第1ピクセルに隣接する第2ピクセルのみを含むことができ、前記第1ピクセルに隣接する前記第2ピクセルのすべて、または他の第2ピクセルの組み合わせを含むことができる。
典型的には、クロストークの大きさは密接に関連しており、例えば、検出されたイベントの大きさまたは高さにほぼ線形に関連している。ここで、検出されたイベントは通常、検出された粒子または光子のエネルギーを表す。したがって、時間関数クロストーク補償シグネチャは、イベント高さまたは対応する粒子エネルギーなど、イベントフィーチャの1つ以上のプロパティに関して正規化された使用のための形式で提供されてもよい。次いで、アナライザは、クロストークを引き起こす並行(コンカレント)検出されたイベントフィーチャの1つ以上のプロパティに従ってクロストークシグネチャをスケーリングし、スケーリングされたシグネチャを他のピクセルの信号(これは、上述したように、半導体検出器において隣接するピクセルであってもよい)に適用するように構成することができ、例えば、加算、減算、または他の形式の組み合わせ、重ねあわせ、または他の数学的手段によって、同期ピクセル信号の対応する時系列とスケーリングされた補償シグネチャとを結合し、それによって、クロストークを少なくとも部分的にキャンセルする。典型的には、各検出されたイベントフィーチャは、対応するピクセルによって出力された信号のステップを含み、ステップの高さは、検出されたイベントフィーチャを生じさせる電離粒子または光子のエネルギーを表す。
本発明はまた、前記信号を前記請求項のいずれかのアナライザに出力するように構成された前記半導体検出器を含む装置を提供する。半導体検出器は、モノリシックであってもよく、セグメント化及び/またはピクセル化されていてもよく、ゲルマニウム検出器であってもよく、HPGe検出器であってもよく、これらの特性と他の特性の異なる組み合わせを含むことができる。
本発明は、蛍光、散乱、回折などの機構を介して検出器に到達するX線、ガンマ線、高エネルギー陽子、中性子及び他の粒子のような様々な異なる物理的機構から生じる様々な異なるタイプの電離放射線を検出するのに適した半導体検出器と共に使用することができる。
本発明はまた、上記の装置を含み、試料から蛍光的に散乱または放出されたX線を検出するように構成されたシンクロトロンビームラインを含む。
本発明はまた、上記に対応する方法、例えば電離放射線を検出するための半導体検出器から受信した信号を解析する方法を含み、それぞれ前記ピクセルに入射する電離粒子または光子に対する前記ピクセルの応答を表す一連のイベントフィーチャを含む前記信号を出力するように構成された複数のピクセルを備えて、該方法は、1つ以上のクロストーク補償またはキャンセルシグネチャを提供するステップと、各ピクセルについて、そのピクセルからの信号におけるイベントフィーチャを検出するステップと、各ピクセルについて、クロストークのための他の、例えば隣接するピクセルの信号を補償するか、または1つ以上のクロストークシグネチャの選択されたものを使用して検出されたイベントフィーチャからこのクロストークをキャンセルするステップと、を含む。
入射放射線の各光子または粒子に対するピクセルの応答、したがって結果としてのイベントフィーチャは、典型的には、粒子または光子のエネルギーに依存するので、半導体検出器は、スペクトル検出器とされることができ、イベントフィーチャを含む上述の信号は、スペクトル情報を含むイベントフィーチャを有するスペクトル信号である。
上述したように、補償シグネチャは、典型的には、クロストーク自体との同期にも対応するクロストークを引き起こすイベントとのアライメントまたは同期を使用する例えば重ね合わせにより、他のピクセルの信号と合成するためのタイム・ドメイン関数として提供される。補償シグネチャと補償されるクロストークフィーチャとの間の同期化は、クロストークを引き起こすイベントのタイミングを直接検出することによって達成することができるし、あるいは代わりに、時間整合(タイミングアライメント)を提供するために、補償シグネチャのマッチングまたはアライメントを他のピクセルのクロストークフィーチャに使用することもできる。
本発明はまた、例えば、各隣接ピクセルにおける信号への、または検出器におけるピクセルのすべてまたは一部のそれぞれにおいて、各ピクセルにおける検出されたイベントフィーチャに起因するクロストークの多数のインスタンス(すなわち、時系列またはタイム・ドメイン関数)を測定することによって、そして、クロストークの前記多数のインスタンスを組み合わせて、前記各ピクセルから各異なる隣接ピクセルまたは前記他のピクセルへのクロストークのための異なる前記補償/キャンセルシグネチャを形成することによって、1つ以上のクロストークシグネチャを導出または生成する方法を提供する。
本発明はまた、本書に記載のように、クロストーク補償またはキャンセルを用いて検出器信号から電離粒子または放射線の光子のエネルギースペクトルを導出する方法、およびこの方法でエネルギースペクトルを導き出すことを含むX線蛍光分光法を用いて試料を特徴付ける方法、を含む。
本発明はまた、適切なコンピュータ装置上で実行されるときに本書に記載のデータ処理方法を実施するように構成された、コンピュータプログラムコードおよびそのようなコードを持つコンピュータ可読媒体を提供する。
本発明はまた、リアルタイムで半導体検出器からの信号を解析(分析)するように構成されたアナライザ、およびアナライザがクロストークの補償を行うように、対応してリアルタイムで動作する方法であり、好ましくは、例えば、半導体検出器からのイベントフィーチャの出力と同じ速度で(かつ実質的に同時に)、各イベントのためのエネルギーを導出するためのさらなる解析も行うようにする方法を提供する。このために、アナライザは、例えば、フィールドプログラマブルゲートアレイファームウェアまたは他のデジタル信号プロセッサハードウェア要素を使用して実装することができる。
本発明の実施形態は、単なる例として、次の添付の図面を参照して、これより説明される。
図1は、ピクセル化された半導体検出器を用いて電離放射線を検出し、得られた信号を解析してエネルギースペクトルなどの電離放射線の特性(プロパティ)を決定するための装置を概略的に示す。
図2は、検出器の異なるピクセルからの信号間のクロストークを補償するための設備を含む、図1のアナライザ30のさらなる詳細を示す。
図3は、クロストーク補償の前の検出器のピクセルからの信号をグラフ化し、各ステップの高さがピクセルに入射する電離粒子または光子のエネルギーに近似する、ステップの形をとるイベントフィーチャを示す。
図4は、イベントフィーチャが検出されたピクセルに隣接するピクセルにおけるクロストークを補償するためのスケーラブルなクロストーク補償シグネチャを示す。
図5は、説明されたクロストーク補償技術を用いて、及び用いずに得られたモリブデン試料からのX線エネルギースペクトルをグラフ化している。
図1は、電離放射線8の粒子または光子を検出するための装置10の概略図を示す。この装置は、X線またはガンマ線のような電離放射線の粒子または光子が半導体検出器18に入射し突き当てることを可能にする、窓14を有する真空チャンバ12を備える。半導体検出器は、複数のピクセル20を含み、各ピクセルは、一連のイベントフィーチャを含む電気信号を出力するように構成され、各イベントフィーチャは、そのピクセルに当たる電離放射線の粒子または光子に対するピクセルの応答を表す。この種の各イベントフィーチャは、ピクセルからの時間変化信号のインターバルとして、例えば連続的な電気信号として、またはそのような電気信号からサンプリングされた時系列の電圧または電流値として存在する。ピクセルは、2次元の直線アレイまたは他の方法で配列され得る。
装置10はまた、電離放射線を検出器18に到達するよりも前にコリメートするためのコリメータ16と、典型的には各ピクセル毎に1つずつ、各ピクセルによって出力された電気信号を増幅する複数のプリアンプ22と、液体窒素の貯蔵部24を使用して、検出器18及び装置の他の部分を冷却するためのアレンジメントを含む。
ピクセルからの信号は、さらなる処理のためにアナライザ30に送られ、アナライザ30からの出力は、検出器18に当たる電離放射線8を特徴付けるために、例えばコンピュータ5で調査するか、またはさらに解析することができる。典型的には、検出器18及び/またはプリアンプ22からの電気信号は、アナライザで、またはアナライザで受信される前に、デジタル化され、信号の解析の少なくとも一部がデジタル化されたデータで実行される。
検出器18は、好ましくは、各入射粒子または光子のエネルギーを検出するスペクトル検出器である。特に、検出器は通常、半導体検出器である。典型的には、スペクトル半導体検出器は、入射粒子または光子がイオン化を引き起こす実質的に空乏化した領域を有するp−i−n構造を有することによって動作する。得られた電荷キャリアは電極に掃引され、収集される電荷の量は、通常、入射粒子または光子8のエネルギーに比例する。次いで、収集された電荷は、半導体検出器上の各ピクセルに形成された回路を利用して、出力電流または他の表現に変換され得る。各入射粒子または光子は、ピクセルに保持された電荷の増加と出力表示における結果的なステップをもたらし、各ステップの高さは、対応する入射粒子光子のエネルギーを表す。ピクセルは、飽和を避けるために定期的にリセットされなければならず、これは、グローバルリセット方式を使用して、または各ピクセルに対して独立して、同時に全ピクセル20の固定または制御可能な繰り返し速度で行われ得る。
典型的には、このような構成を有効にするために、各ピクセルに電荷感応型プリアンプが装備されている。ピクセルによって収集された電荷は、電荷を電圧信号に変換する電荷感応型プリアンプに入る。電荷はプリアンプの帰還コンデンサに蓄積され、このコンデンサは飽和を避けるために定期的にリセットする必要がある。飽和を回避する別の方法は、帰還コンデンサに並列に抵抗を用い、コンデンサが実際にリークしてからそれが連続的に放電されるようにする方法である。信号に過度のノイズを加えないと同時に、特定の用途に許容可能なコンデンサ放電レートを達成するため、抵抗器の値は適切に選択される必要がある。
高い粒子率または光子入射率では、電気信号のイベントフィーチャは特定のp−i−n構造から出力され、したがってピクセル20はより頻繁に相互に重なり、それによって出力データの品質が低下する。この困難に対処するために、同様の立体検出角をカバーしながら、複数のより小さい検出器のいずれかを一緒にパッケージングするか、または検出器をそれぞれ独立した出力電気信号を有する複数のセグメントまたはピクセル20を有するように構成することができる。この目的のために、検出器は、複数のピクセル20を形成するようにセグメント化されたゲルマニウムなどの半導体の単一片または結晶からなるモノリシックアレイ検出器であってもよい。そのような検出器は、当該分野で周知の様々なベンダ(例えばキャンベラインダストリーズ(Canberra Industries)など)から入手可能である(http://www.canberra.com参照)。
複数のピクセル20を含む電離放射線用の半導体検出器の難点は、アナライザによって受信された信号がクロストークを呈する可能性があることであり、それによって、電気信号における1つのイベントフィーチャは、1つ以上の他のピクセルからの信号における対応するクロストークフィーチャを生じさせるか、またはそれに対応する。この効果は、2つのピクセルが共通の長さの境界を共有する場合、すなわち隣接する場合に、特に顕著になる可能性がある。規則的に、正方形ピクセルの長方形アレイでは、2つのピクセルが隣接する各境界は同じ長さとなるが、正方形でないピクセルおよび他のアレイの場合は、2つの隣接ピクセル間の共通境界の長さを形成し、ピクセル間の正確な近接度は、検出器内で変化し得る。クロストークは、2つのピクセルが互いに接近しているか、より長い共通の境界を共有している場合に大きくなる可能性があるが、2つのピクセルが隣接していない場合には、1つ以上の他の介在ピクセルによって分離される。
この問題に対処するために、図1のアナライザ30は、このクロストークを補償するように構成されている。アナライザ30には補償シグネチャデータ32が提供され、任意の特定のピクセルについて、アナライザは、検出器18から受信した信号におけるイベントフィーチャを検出し、補償シグネチャデータ32を用いて、他の、例えば隣接するピクセルのうちの1つ以上からコンカレント(並行)信号を補償する。補償シグネチャデータは、典型的には、補償されるべきイベントのタイミングに基づく同期またはアライメントを用いて、補償されるべき信号との重ね合わせを使用して結合されるタイム・ドメイン信号を表す。
アナライザ30がこのクロストーク補償を実行するように構成され得る例示的な方法が図2に示されている。検出器18のピクセル1…nから受信された信号S1…Snは、それぞれが検出するように構成されたイベント検出器エレメント34−1…34−nを含むイベント検出器34に送られる。上記信号の1つでは、関連ピクセル20での粒子または光子8の入射によってイベントフィーチャが引き起こされる。そのようなイベントフィーチャは、例えば、信号の最低の高さ及び/または最大の長さのような所定の基準を満たしている信号のステップフィーチャとすることができる。図2において、信号S1は、高さE(入射粒子または光子のエネルギーに対応する)及びタイミングtのようなステップフィーチャ35を有するものとして示されている。フィーチャ35のタイミングは、典型的には、アナライザ30の内部クロックを参照することによって可能で、第1の信号のイベントフィーチャを、他の信号の対応するクロストークフィーチャと正確に位置合わせできるようにすべきである。このようなクロストークフィーチャ37は、ピクセル2…nから受信された信号の範囲内で図2に示されている。
アナライザ30はまた、信号S1…Snのうちの1つにおいて、信号S1…Snのうちの別のイベントフィーチャによって引き起こされるか、またはそれに関連するクロストークを補償するようにそれぞれ準備されたクロストーク補償器エレメント38−1〜38−nを備えるクロストーク補償器38を含む。例えば図2に示す信号S1のような信号の1つにおいてイベントフィーチャが検出されると、関連するイベント検出器エレメント(図2のエレメント34−1)は、他の検出器ピクセル、例えばイベントが検出された信号に隣接するピクセルに受信された信号を処理する少なくともそれらのクロストーク補償器エレメントにトリガ信号36を送る。このトリガ信号は、典型的には、検出されたイベントフィーチャの高さE及びタイミングtの一方または両方を含むことができる。トリガ信号36は、明示的または暗示的に、イベントフィーチャが生じたピクセルの識別子Pも運ぶ。
関連するトリガ信号、すなわち、別のピクセルにおける検出されたイベントフィーチャに関連するトリガ信号を受信すると、クロストーク補償器エレメント38−iは、クロストーク補償データ32を使用することによって、トリガ信号を発生させるクロストークの影響について信号Siを補償する。これは様々な方法で達成できる。図2の例では、タイム・ドメインクロストークシグネチャの別々のセット40が各ピクセルP1…Pnに対して格納され、各セット40は、そのセットが関連付けられているピクセルに隣接する各ピクセルのシグネチャ42を含む。シグネチャ42の必要なセット40は、トリガ信号36の識別子Pを使用する特定のクロストーク補償器エレメントによって選択され、セット内の特定のシグネチャ42は、クロストーク補償器エレメントの同一性、したがって隣接するピクセルの識別性に対応するように選択される。次いで、選択されたシグネチャは、トリガ信号の高さE、または立ち上がり時間などのトリガ信号の何らかの他の特性を利用してスケーリングされるか、または他の方法で変更され、トリガ信号のタイミングtを用いて正確に信号Siに適用され、信号Siとスケーリングされたシグネチャの重ね合わせとして、同期(すなわち時間ベースのアライメント)を提供する。
選択されたシグネチャ42は、例えば、トリガ信号の高さEと線形にスケーリングされるか、または他の何らかの方法、例えば、データに最良に適合する非線形関数を使用して、または理論的もしくは近似に従ってスケーリングすることができる。シグネチャ42が適用されるべき信号とのタイミングアライメントは、トリガのタイミングtのみを用いて実行され得、または、アライメントは、例えば、データに対する最良適合に基づいて、または理論または近似に従って、高さEの関数として変化させることができる。
したがって、各シグネチャ42は、エネルギー、高さ、または第1ピクセルのイベントフィーチャ35の他の態様に従ってスケーリングし得る別のピクセルからの信号上の第1ピクセルからの信号におけるイベントフィーチャ35の効果を表す時間ベースの関数またはデータ系列であってもよい。シグネチャ42は、解析関数のパラメータのセット、値の時系列、または他の方法で格納することができる。
2つ以上のイベントが、特定のピクセルからの信号において互いに近接して発生する場合、一般に「積み上げ(pile up)」と呼ばれるが、近接イベントを互いに区別することは困難であり、関連する光子または粒子のエネルギーをイベントフィーチャから測定することはより困難である。そのような場合、本書に記載された技術を使用するクロストーク補償はまた、あまり効果的ではない。しかし、イベントの近位の性質に基づく因子を用いたシグネチャ42のスケーリングは、そのような場合のクロストーク補償を改善するために使用することができる。例えば、スケーリングは、補償されるイベントの周りの特定の時間ウィンドウ内の同じ信号内に見出されるイベント、例えばそのようなイベントの数に基づいてもよく、あるいはそのような各イベントの補償されるべきイベントからの時間間隔によって重み付けされた近似イベント(好ましくはそのようなウィンドウ内で)の数をカウントするメトリックに基づいてもよいし、または他の方法でもよい。
シグネチャ42は、典型的には、クロストークが補償されるべきイベントの前後に、ある有限の時間を伸ばして時間制限される。例えば、図3に示すイベントフィーチャを先に参照すると、イベントから潜在的なクロストークを十分にカバーするために、各シグネチャは、1マイクロ秒のオーダーの時間範囲を有することができる。
それらシグネチャがどのように生成されるかに依存して、クロストークシグネチャの一部または全部はベースラインに戻らず、ピクセル信号に適用されたときにオフセットを追加することができ、それが行なわれなければ、複数のイベント補償動作にわたってピクセル信号におけるオフセットの累積が生じることがある。これを補償するために、ピクセル信号から減算されるか、別個の記録として保持されるかにかかわらず、ピクセル信号がリセットされるまで、シグネチャの各適用に対して(各イベントのシグネチャに適用されるエネルギーまたは他のスケーリングを可能にする)等価な反対のオフセットが蓄積され得る。
図2においては、クロストーク補償データ32が各ピクセルP1…Pnについて示されているが、あたかも別々のシグネチャ42を有する8つの隣接ピクセルがあるかのように示されているが、実際には、すべてのピクセルが8つの隣接するピクセルを有するわけではないので、例えば、エッジピクセル及びコーナーピクセルのように、すべての隣接ピクセルを補償するために必要とされるシグネチャがより少ないことに留意されたい。シグネチャ42は、ピクセル間で様々な方法で共有することもできる。例えば、特定の装置では、システムによって少数の別個のシグネチャ形状しか必要とされないこと、場合によっては特定のピクセルに適用可能なスケーリングまたは他の変換を用いて、クロストーク補償器エレメントに正しいシグネチャ形状を選択させるために、インデックス装置構成を使用することができること、について見出すことができる。
上述したように、クロストークシグネチャ42及び対応するクロストーク補償構成はまた、非隣接ピクセルに対して提供されるか、またはプリアンプ22と、検出器18とアナライザ30との間の他の電気導体との間のクロストークなどの信号経路内の他のソースに対して提供され得、隣接するピクセルについて論じたのと本質的に同じ方法で、クロストーク効果が最も近い隣接ピクセルを超えて広がる可能性がある。さらに、図2のクロストーク補償データ32は、ピクセルの直線的な配置構成を意味するが、線形、三角形、六角形及び他の構成も見出すことができ、それに応じてクロストーク補償データ32を構成することができる。
補償器38によって使用されるクロストーク補償データ32は、様々な方法で導き出すことができる。例えば、装置10の試験または較正の使用から生じる信号S1…Snは、イベントフィーチャを検出するために、かつ各信号S1…Snのための他の、例えば隣接するピクセルのそれぞれに生じる(時系列または関数としての)対応するクロストーク信号を測定するために、使用され得る。このような複数のクロストークイベントは、次に、平均化されるか、または他の方法で組み合わされて、他のまたは隣接するピクセルの各々について代表的なクロストークシグネチャを得ることができる。この解析は、典型的には、粒子または光子エネルギーの範囲に対応するイベントフィーチャについて、各信号、互いのまたは隣接するピクセルについてのクロストーク信号を測定すること、得られたシグネチャのためのエネルギーを有する特定の倍率を決定すること、または、特定の信号について他のピクセルまたは隣接ピクセルに共通の適切な倍率を得ること、または、検出器からのすべての信号に対して、他のすべてのピクセルまたは隣接するピクセルに共通の適切な倍率を得ること、を含むことができる。
クロストーク補償器エレメント38−1…38−nによる必要なクロストーク補償の後、図2の補償器38に続く信号スケッチによって示されるように、隣接するピクセルの信号に見られるクロストークの影響はかなり低減される。次に、これらの補償された信号39は、1つ以上の解析エレメント44に送られ、イベント検出、例えば各イベントのエネルギーを決定するイベント測定、及びイベントカウントなどの機能を実行する。これらの解析エレメント44は、イベント検出器エレメント34−1…34−nによって既に決定されたトリガ信号のような検出データを使用することができ、または、既に上述したプロセスとは独立して信号を処理することができる。解析エレメント44に送られたクロストーク補償された形式の信号S1…S2は、クロストークのレベルが低減され、解析エレメントの出力46が改善される。典型的には、解析エレメント44の出力は、検出器に入射した粒子または光子のエネルギーの分布を表すエネルギースペクトルを含み、例えば、ピクセルのすべての結合された領域に関するスペクトル、及び/または複数のピクセルの組み合わせの各ピクセルに関して分離するスペクトルを含むことができる。
解析エレメント44が補償された信号に対して作用するいくつかの方法は、R. Farrow らの「X-ray signal processing electronics for solid state detectors」(Scientific Instruments, Volume 66, Issue 2 (1995) pp. 2307-2309)、及びR. Farrow らの「XSPRESS: A new generation of detector systems for EXAFS studies」(Physica B, Volume 208 (1995), pp 256-258)で論じられている。
アナライザ内のイベント検出器34及びクロストーク補償器38の一方または両方は、例えば、アナライザ30に到達する前または後のいずれかに既にデジタル化された信号に作用する、必要に応じて1つ以上の汎用マイクロプロセッサ上で動作するソフトウェアを用いて実現(実装)することができる。しかしながら、検出器18が粒子または光子8の個々の衝撃のエネルギーを識別し決定することができる非常に高い速度は、回路内に直接的にまたはリアルタイムファームウェアを使用して実施されるハードウェアソリューションが好ましいことを意味する。このような高速実装は、現代のシンクロトロン光源のような設備のユーザにとって日常的に利用可能な安定した高フラックス環境を利用するのに役立つ。例えば、本発明は、レジスタ転送レベルVHDLまたはVerilogで設計されたフィールドプログラマブルゲートアレイ(FPGA)ファームウェアを使用して実装することができる。比較的短いクロストークシグネチャは、多くのシグネチャに同時にアクセスできるFPGAの内部メモリに容易に収まる。他の実施態様では、グラフィックス処理ユニット(GPU)及び/または他のデジタル信号プロセッサを使用することができる。
本発明は、高解像度X線及びガンマ線分光法のためのシンクロトロン光源、並びに本書の他の場所で述べたような他の用途に使用することができる。このような光源の輝度を高めるために、ビームラインは、X線フラックスが増加するにつれて、飽和しない、またはエネルギー分解能を失う検出器を必要とする。フラックスの増加に追随するために、検出器システムの性能に使用される2つの主なアプローチがある。それらは、与えられた立体角内のピクセル数を増加させ、性能低下なしに与えられたピクセル数がピクセルカウントレート当たりより高くなるようにすることである。
英国におけるダイヤモンドライトソース(DLS)は、上記で参照されたXSPRESSアナライザのアップデートされた機能拡張版であるXSPRESS2アナライザに結合された、社内設計の高速低ノイズプリアンプを装備した多数のピクセル(9,36及び64ピクセル)を備えたモノリシックセグメント化されたキャンベラインダストリーズ製HPGe検出器を使用して、主X線分光ビームラインでこれらのアプローチの両方を組み合わせる。
X線蛍光分光法の重要な仕様の少なくとも2つ(エネルギー分解能及びピーク/バックグラウンド比)は、ピクセル当たりのカウント速度が増加するにつれて著しく低下し、入射X線8のエネルギーが増加するにつれて顕著になることが発明者によって観察された。検出器ピクセルまたはプリアンプからの未処理信号の厳重な検査は、隣接するピクセル上で発生するX線相互作用から各ピクセルへのかなりのクロストークがあり、性能が低下することを示している。本発明者らは、そのようなクロストークアーチファクトの除去が、特により高い固有のカウントレートで、改善された性能をもたらすことを見出した。
様々な計数率での異なる化学元素(従って、異なる特性X線蛍光放出エネルギーを有する)の試料からのX線蛍光の測定を、ダイヤモンドライトソースのビームライン上で行った。検出器のピクセル20が同時に定期的にリセット(グローバルリセット)されるように動作するXspress2アナライザ30に接続されたキャンベラインダストリーズ製36ピクセルHPGe検出器18を用いて、図1に示す装置を用いて蛍光を測定した。
図1及び図2で示したような装置用に、クロストーク補償データ32を計算するために、比較的低い計数率を有する蛍光X線データを使用して、特定のピクセル及び隣接するピクセル間のコンカレント(並行)イベントフィーチャの数を最小限に抑えた。低計数率条件の下で多数の隣接ピクセルに対するプリアンプ22からのコンカレント電気信号は、デジタル化され記憶される。このデータは、XSPRESS2アナライザの解析エレメント44をエミュレートするように構成されたMatlab(RTM)ベースのソフトウェアモデルへの入力として使用され、グローバルリセット方式を利用してすべてのピクセルからのデータが同期された。
図3において、14ビットと80MHzでサンプリングされたピクセルの1つからのデジタル化信号の一部の例を示す。ステップの形態の一連のイベントフィーチャ35が見られ、各ステップは、X線の入射に起因するピクセルにおける電荷の蓄積を表し、ステップの高さはX線エネルギーに対応する。
次いで、Matlabモデルを用いて、図2のイベント検出プロセス34と同様または同等のイベント検出プロセスを適用して、各ピクセルの信号における各イベントフィーチャ35の高さE及びタイミングtの第1パス推定を決定した。特定の中心ピクセルの信号で検出されたイベントについては、隣接するピクセルで検出されたイベントと重複するイベントの使用を避けるよう注意しながら、8つの隣接ピクセルの各々における時系列としてのクロストーク特性を測定した。各隣接ピクセルについて、クロストークシグネチャは、イベントの高さEに対して正規化され、平均化されて、中心ピクセルにおけるイベントに関してその隣接ピクセルに対するクロストークシグネチャ42を形成する。このプロセスを検出器内のすべてのピクセルについて繰り返すことにより、各中心ピクセルに対して、その中心ピクセルに隣接する各ピクセルへの予想される時系列またはクロストークの関数を表すエネルギー正規化クロストークシグネチャ42のセット40からなるクロストーク較正データ32のセットをもたらす(用語「中心」は、厳密にはそれらの隣接ピクセルに対して「中心」ではないエッジまたはコーナーピクセルを排除することを意図していないことに留意されたい)。
図4には、このようなクロストークシグネチャ42のセットの例として、21,22,23,27,29,33,34及び35で示される直線アレイ内の隣接するピクセルを有する、ピクセル28と表示される「中心」ピクセルについて示されている。各シグネチャ42を示すグラフの水平軸は、数値範囲40〜110がほぼ1マイクロ秒を表すように、連続的なデジタル化試料の単位でマーキングされ、垂直軸は、補償されるべき検出されたイベントフィーチャのエネルギーに正規化される。したがって、例えばピクセル33に適用される補償信号は、ピクセル28に対して検出された対応するイベントフィーチャのサイズの約1%でピークになるのに対して、ピクセル34に適用される補償信号は、ピクセル28についての対応するイベントフィーチャの約5%でピークになる。
図4のシグネチャ42の全てが少なくとも微妙に異なる形状を有することが分かるので、クロストーク較正データ32は、隣接する各ピクセルに対して異なる別個のデータ系列または各中心ピクセルに対応する各セット40の各隣接ピクセルに対して異なるパラメータ化された関数シグネチャを含むことが適切であり得る。しかしながら、図4のシグネチャのいくつかの間の類似性も留意され、この類似性を利用して、様々な方法で格納および使用されるシグネチャの数を減らすことができる。
上記のように各ピクセルに関してクロストーク較正データ32を計算した後、このデータを用いて、イベントフィーチャが同時に頻発し、隣接するピクセルに重複するように、比較的高い計数率のX線検出を伴う実験から得られた補償されていない信号データに対してクロストーク補償を行った。このクロストーク補償は、先に図2に関連して上述した様々な方法で図2のアナライザ30に同じプロセスを実施できることに留意しながら、記載されたMatlabモデルを用いて最初に実行された。
図2のイベント検出器32に対応するクロストーク補償の第1段階の間に、各ピクセルの生(raw)信号のイベントが検出され、各イベントの高さE及び時間tを出力するために測定された。次に、図2の補償器38に対応する第2の段階の間に、高さEを用いて、イベントが検出された中心ピクセルに対する各隣接ピクセルのクロストークシグネチャを線形にスケーリングし、正確なアライメントを確実にするために検出されたイベントフィーチャ時間tを用いて、それらの隣接するピクセルに対する補償されていない信号からスケーリングされたシグネチャ42を差し引いて、すべてのピクセルについて補償された信号を生成した。
続いて、図2の解析エレメント44のMatlabモデルを用いてピクセルの補償信号を解析し、その結果を未補償の生信号に適用した同じ解析と比較した。この解析には、適応レートフィルタリングに続いて、さらにダウンストリームイベントフィーチャのシェイピングとバックグラウンドのリークランプ補正を使用し、各ピクセルに入射する各X線のエネルギーを導出するために使用される高さのイベントフィーチャを検出する。
このようなイベントを多数組み合わせることで、比較のために、クロストーク補償処理を用いた場合と用いなかった場合の放射線のエネルギースペクトルを提供する。図5は、KαおよびKβピークの周りの、モリブデン試料についてこのようにして導出されたエネルギースペクトルの一部を示す。円で表された点は、クロストーク補償なしで導出されたエネルギースペクトルを示し、実線は、検出器からの同じ電気信号の結果を示しているが、上述したようにクロストーク補償で処理されている。の比較はスペクトル応答の大幅な改善を示し、エネルギー分解能は通常35%向上し、ピークトゥバックグラウンド(低エネルギーテール)比は通常2倍になる。
上述のクロストーク補償の方法は、イベントのエネルギーまたは他の測定可能な特性に従ってスケーリングまたは変更されるにもかかわらず、イベントごとにほぼ同じままであるクロストークの特性にある程度依存する。しかしながら、本発明者らは、一般に、隣接するピクセルにおける所定のクロストークフィーチャと、導出されたクロストークシグネチャ42との間に良好な一致が存在することを見出した。
本発明の範囲から逸脱することなく、記載された実施形態に様々な変更を加えることができる。例えば、クロストーク補償シグネチャを使用して、隣接するピクセルまたは他のピクセルの組み合わせの信号を補償するいくつかの詳細な実施形態について説明したが、電離粒子または光子のエネルギーのようなそれらの信号の特性を導出するために、これらの信号のさらなる解析またはその後の解析の前に、信号の解析の一部としてクロストーク補償のためにクロストーク補償シグネチャを使用して、代わりにそのような特性を導出することもできる。
Claims (24)
- 電離放射線の粒子または光子のエネルギーを検出するための半導体検出器から信号を受信するように構成されたアナライザであって、
前記半導体検出器は、複数のピクセルを備え、各ピクセルは前記ピクセルに入射する電離粒子または光子に対するピクセルの応答を表す一連のイベントフィーチャを含む信号を出力するように構成され、前記アナライザによって受信された信号は、相互間にクロストークを呈し、それによって少なくとも一部の前記イベントフィーチャが、少なくとも1つの他のピクセルからの信号に対応するクロストークフィーチャを生じさせ、
タイム・ドメイン関数として提供される1つ以上のクロストーク補償シグネチャを含み、各ピクセルについて、そのピクセルからの信号におけるイベントフィーチャを検出するように構成されるとともに、該対応する検出されたイベントフィーチャとアライメントされた1つ以上のクロストーク補償シグネチャのうちの選択されたものを用いて、ピクセルのうちの1つ以上の他のピクセルからのコンカレント信号の解析を補償するように構成される、アナライザ。 - 前記イベントフィーチャの1つ以上のプロパティに関して前記クロストーク補償シグネチャを正規化するように、かつ、前記検出された各イベントフィーチャの前記1つ以上のプロパティを測定するように、かつ、前記コンカレント信号の補償解析に使用するために前記測定されたプロパティに従って前記1つ以上のクロストーク補償シグネチャの前記選択されたものをスケーリングするように、さらに構成される請求項1に記載のアナライザ。
- 前記1つ以上のプロパティが、前記光子または前記電離放射線の他の粒子のエネルギーに対応するプロパティを含み、該対応するイベントフィーチャを生じさせる、請求項1または2に記載のアナライザ。
- 前記ピクセルのうちの1つ以上の他のピクセルからのコンカレント信号の解析は、それぞれのアライメントされたクロストーク補償シグネチャとの前記コンカレント信号のそれぞれの重ね合わせを含む、請求項1〜3いずれかに記載のアナライザ。
- 各イベントフィーチャのタイミングを検出し、検出されたタイミングを用いて対応するクロストーク補償シグネチャを前記コンカレント信号とアライメントさせることを含む、請求項1〜4いずれかに記載のアナライザ。
- 前記1つ以上のクロストーク補償シグネチャは、各ピクセルに対して、前記ピクセルの1つ以上の他のピクセルからのコンカレント電気信号におけるクロストークを補償するために用いるように構成される少なくとも1つの補償シグネチャを含む、請求項1〜5いずれかに記載のアナライザ。
- 前記1つ以上のクロストーク補償シグネチャは、各ピクセルの補償シグネチャのセットを含み、前記ピクセルの他の1つのピクセルのコンカレント電気信号におけるクロストークを補償するために、セットの各構成要素を使用するように構成されている請求項1〜6いずれかに記載のアナライザ。
- 前記ピクセルのうちの他の1つは、1つ以上の隣接するピクセルを含む、請求項1〜7いずれかに記載のアナライザ。
- クロストークを引き起こすコンカレント(並行)検出されたイベントフィーチャの1つ以上のプロパティに従ってスケーリングされるクロストーク補償シグネチャのうちの1つ以上の選択されたものを用いて、前記ピクセルの1つ以上の他のピクセルからのコンカレント電気信号の解析におけるクロストークを補償するように構成されている、請求項1〜8いずれかに記載のアナライザ。
- 前記クロストークを引き起こすコンカレント検出されたイベントフィーチャの前記1つ以上のプロパティは、前記検出されたイベントフィーチャを生じさせる前記電離粒子または光子のエネルギーの測定値を含む、請求項9に記載のアナライザ。
- 各検出されたイベントフィーチャは、対応するピクセルによって出力された信号のステップを含み、ステップの高さは、検出されたイベントフィーチャを生じさせる電離粒子または光子のエネルギーを表す、請求項1〜10いずれかに記載のアナライザ。
- 請求項1〜11いずれかに記載の前記アナライザに前記信号を出力するように構成された前記半導体検出器を備える装置。
- 前記半導体検出器は、セグメント化HPGe検出器である、請求項12に記載の装置。
- 試料から蛍光的に散乱または放射されたX線を検出するように構成された請求項12または13に記載の装置を含むシンクロトロンビームライン。
- 複数のピクセルを含む電離放射線の半導体検出器から受信した信号を解析する方法であって、
各ピクセルは、粒子または光子のエネルギーに依存する前記ピクセルに入射する電離粒子または光子に対するピクセルの応答を表す一連のイベントフィーチャを含む前記信号を出力するように構成され、
前記方法は、
タイム・ドメイン関数として1つ以上のクロストーク補償シグネチャを提供するステップと、
各ピクセルから前記出力された信号を受信するステップと、
各ピクセルについて、そのピクセルからの信号におけるイベントフィーチャを検出するステップと、
各ピクセルに対し、対応するイベントフィーチャと同期された1つ以上のクロストーク補償シグネチャのうちの選択されたものを用いて、検出されたイベントフィーチャからのクロストークのための前記ピクセルの他のピクセルの信号の解析を補償するステップと、
を含む方法。 - 前記検出されたイベントフィーチャのタイミングを検出するステップと、前記検出されたタイミングを用いて、クロストーク補償シグネチャを対応するイベントフィーチャと同期させて、クロストークのための前記ピクセルの他のピクセルの信号の解析を補償するステップと、をさらに含む、請求項15に記載の方法。
- 前記検出されたイベントフィーチャからのクロストークのための前記ピクセルの他のピクセルの信号の補償解析が、前記1つ以上のクロストーク補償シグネチャの前記信号との重ね合わせを含む、請求項15または16に記載の方法。
- 各検出されたイベントについて、前記ピクセルの他のピクセルは、前記検出されたイベントの前記ピクセルに隣接するピクセルを含む、請求項15〜17いずれかに記載の方法。
- 各検出されたイベントフィーチャに対して、
前記イベントフィーチャの1つ以上のプロパティを測定するステップと、
測定された1つ以上のプロパティに従って、選択された1つ以上のクロストーク補償シグネチャをスケーリングするステップと、
スケーリングされた選択されたクロストークシグネチャを前記ピクセルの他のピクセルのコンカレント信号と組み合わせるステップと、
をさらに含む、請求項15〜18いずれかに記載の方法。 - 前記1つ以上のプロパティが、前記検出されたイベントフィーチャを生じさせる前記電離粒子または光子のエネルギーに対応するプロパティを含む、請求項19に記載の方法。
- イベントフィーチャが検出されるべき各ピクセルについて、前記1つ以上のクロストーク補償シグネチャは、前記ピクセルから、前記ピクセルの異なる他のピクセルへのクロストークを補償するために使用される一連の補償シグネチャを含む、請求項15〜20いずれかに記載の方法。
- 複数の他の、例えば隣接するピクセルの各々に対する各ピクセルの検出されたイベントフィーチャに起因するクロストーク信号の複数のインスタンスを測定することによって前記タイム・ドメイン関数のクロストークシグネチャを生成するステップと、
前記各ピクセルから他のピクセルのそれぞれの異なる1つへのクロストークのための異なる前記補償シグネチャを形成すべく前記クロストーク信号の前記複数のインスタンスを合成するステップと、
を含む請求項15〜21いずれかに記載の方法。 - 前記補償された信号から前記電離粒子または光子のエネルギースペクトルを導出するステップをさらに含む、請求項15〜22いずれかに記載の方法。
- 請求項23に記載の方法を用いてエネルギースペクトルを導出するステップを含む、X線蛍光分光法を用いて試料を特徴付ける方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB1408853.8 | 2014-05-19 | ||
GBGB1408853.8A GB201408853D0 (en) | 2014-05-19 | 2014-05-19 | Analysis of signals from pixellated detectors of ionizing radiation |
PCT/GB2015/051470 WO2015177534A1 (en) | 2014-05-19 | 2015-05-19 | Analysis of signals from pixellated detectors of ionizing radiation |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2017521676A true JP2017521676A (ja) | 2017-08-03 |
Family
ID=51135071
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017513366A Pending JP2017521676A (ja) | 2014-05-19 | 2015-05-19 | 電離放射線のピクセル化検出器からの信号の解析 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10168438B2 (ja) |
JP (1) | JP2017521676A (ja) |
GB (2) | GB201408853D0 (ja) |
WO (1) | WO2015177534A1 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201709298D0 (en) | 2017-06-12 | 2017-07-26 | Science And Tech Facilities Council | Depth Correction in Pixellated detectors of ionizing radiation |
US11676685B2 (en) | 2019-03-21 | 2023-06-13 | Illumina, Inc. | Artificial intelligence-based quality scoring |
US11210554B2 (en) | 2019-03-21 | 2021-12-28 | Illumina, Inc. | Artificial intelligence-based generation of sequencing metadata |
US11423306B2 (en) * | 2019-05-16 | 2022-08-23 | Illumina, Inc. | Systems and devices for characterization and performance analysis of pixel-based sequencing |
US11593649B2 (en) | 2019-05-16 | 2023-02-28 | Illumina, Inc. | Base calling using convolutions |
EP4107735A2 (en) | 2020-02-20 | 2022-12-28 | Illumina, Inc. | Artificial intelligence-based many-to-many base calling |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105995A (ja) * | 2004-10-07 | 2006-04-20 | Siemens Ag | ディジタルx線撮影装置およびディジタルx線撮影装置におけるx線像の撮影方法 |
JP2007175154A (ja) * | 2005-12-27 | 2007-07-12 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
US20070235657A1 (en) * | 2005-09-30 | 2007-10-11 | The Regents Of The University Of Michigan | Three-dimensional, position-sensitive radiation detection |
JP2007536549A (ja) * | 2004-05-03 | 2007-12-13 | コンパニ ジェネラル デ マティエール ニュークリエール | 2パラメータスペクトルの補正方法 |
US20130168562A1 (en) * | 2011-07-01 | 2013-07-04 | Commissariat A L'energie Atomique Et Aux Ene Alt | Ionizing radiation detection device with a semi-conductor detector having and improved spectrometric response |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3752986A (en) * | 1972-03-03 | 1973-08-14 | J Fletcher | Compton scatter attenuation gamma ray spectrometer |
DE2921780C3 (de) * | 1979-05-29 | 1982-02-11 | Siemens AG, 1000 Berlin und 8000 München | Verfahren und Schaltungsanordnung zur Übertragung von Daten über eine Zweidrahtleitung |
JPS6383687A (ja) * | 1986-09-29 | 1988-04-14 | Toshiba Corp | 核医学デ−タ処理装置 |
EP0583118A3 (en) * | 1992-07-30 | 1996-03-27 | Summit World Trade Corp | Gamma camera |
US5440322A (en) * | 1993-11-12 | 1995-08-08 | In Focus Systems, Inc. | Passive matrix display having reduced image-degrading crosstalk effects |
US5982847A (en) * | 1996-10-28 | 1999-11-09 | Utah State University | Compact X-ray fluorescence spectrometer for real-time wear metal analysis of lubrucating oils |
US6288399B1 (en) * | 1997-11-12 | 2001-09-11 | Cti Pet Systems, Inc. | Depth of interaction detector block for high resolution positron emission tomography |
US6140649A (en) * | 1997-12-23 | 2000-10-31 | General Electric Company | Imaging attenuation correction employing simultaneous transmission/emission scanning |
US6097785A (en) * | 1998-10-30 | 2000-08-01 | The United States Of America As Represented By The Secretary Of The Navy | Cone penetrometer utilizing an X-ray fluorescence metals sensor |
US6249003B1 (en) * | 1999-06-18 | 2001-06-19 | General Electric Company | Imaging attenuation correction method employing multiple energy scan masks and windows |
US6791090B2 (en) * | 2000-05-17 | 2004-09-14 | Koninklijke Philips Electronics N.V. | Compton deconvolution camera |
JP2002071816A (ja) * | 2000-08-29 | 2002-03-12 | Japan Atom Energy Res Inst | 2次元放射線および中性子イメージ検出器 |
US7456878B1 (en) | 2005-05-11 | 2008-11-25 | Magnachip Semiconductor Ltd. | Method of color filter design and color reproduction under the effect of pixel crosstalk in CMOS image sensors |
US7825370B2 (en) | 2005-06-23 | 2010-11-02 | General Electric Company | Method and system for calibrating a computed tomography system |
EP2266250A1 (en) * | 2008-03-13 | 2010-12-29 | University College Cork-National University of Ireland, Cork | An optical communication system and method for secure data communication using quantum key distribution |
CN101750621A (zh) * | 2008-12-08 | 2010-06-23 | 皇家飞利浦电子股份有限公司 | 对来自像素化探测器的测量的评估 |
US20120015825A1 (en) * | 2010-07-06 | 2012-01-19 | Pacific Biosciences Of California, Inc. | Analytical systems and methods with software mask |
FR2968879A1 (fr) * | 2010-12-10 | 2012-06-15 | St Microelectronics Sa | Procede de caracterisation de pixels d'un capteur d'image |
US20120166128A1 (en) | 2010-12-28 | 2012-06-28 | Ikhlef Abdelaziz | Method and apparatus for detector calibration |
US20120173175A1 (en) | 2010-12-30 | 2012-07-05 | Deepak Devicharan | Photodiode array with algorithm-based crosstalk reduction |
US8908831B2 (en) * | 2011-02-08 | 2014-12-09 | Rapiscan Systems, Inc. | Covert surveillance using multi-modality sensing |
US8374814B2 (en) * | 2011-02-16 | 2013-02-12 | Rigaku Corporation | X-ray detection signal processing apparatus and method therefor |
US9246106B2 (en) * | 2011-04-05 | 2016-01-26 | The Board Of Trustees Of The Leland Stanford Junior University | Electron deficient molecules and their use in organic electronic applications |
US8610081B2 (en) * | 2011-11-23 | 2013-12-17 | General Electric Company | Systems and methods for generating control signals in radiation detector systems |
-
2014
- 2014-05-19 GB GBGB1408853.8A patent/GB201408853D0/en not_active Ceased
-
2015
- 2015-05-19 WO PCT/GB2015/051470 patent/WO2015177534A1/en active Application Filing
- 2015-05-19 GB GB1508576.4A patent/GB2528168B/en active Active
- 2015-05-19 JP JP2017513366A patent/JP2017521676A/ja active Pending
- 2015-05-19 US US15/312,124 patent/US10168438B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007536549A (ja) * | 2004-05-03 | 2007-12-13 | コンパニ ジェネラル デ マティエール ニュークリエール | 2パラメータスペクトルの補正方法 |
JP2006105995A (ja) * | 2004-10-07 | 2006-04-20 | Siemens Ag | ディジタルx線撮影装置およびディジタルx線撮影装置におけるx線像の撮影方法 |
US20070235657A1 (en) * | 2005-09-30 | 2007-10-11 | The Regents Of The University Of Michigan | Three-dimensional, position-sensitive radiation detection |
JP2007175154A (ja) * | 2005-12-27 | 2007-07-12 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
US20130168562A1 (en) * | 2011-07-01 | 2013-07-04 | Commissariat A L'energie Atomique Et Aux Ene Alt | Ionizing radiation detection device with a semi-conductor detector having and improved spectrometric response |
Non-Patent Citations (2)
Title |
---|
A. BRAMBILLA, ET AL.: "CdTe Linear Pixel X-Ray Detector With Enhanced Spectrometric Performance for High Flux X-Ray Imaging", IEEE TRANSACTIONS ON NUCLEAR SCIENCE, vol. Volume: 59, Issue: 4, JPN7019000602, 18 July 2012 (2012-07-18), US, pages 1552 - 1558, ISSN: 0003988757 * |
SEZGIN AYDIN: "PERFORMANCE IMPROVEMENT OF AN AGATA SEGMENTEDGERMANIUM DETECTOR", PHD THESIS, JPN7019000601, 2009, TR, ISSN: 0003988756 * |
Also Published As
Publication number | Publication date |
---|---|
US20170102468A1 (en) | 2017-04-13 |
GB2528168B (en) | 2021-09-22 |
GB2528168A (en) | 2016-01-13 |
WO2015177534A1 (en) | 2015-11-26 |
US10168438B2 (en) | 2019-01-01 |
GB201508576D0 (en) | 2015-07-01 |
GB201408853D0 (en) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2017521676A (ja) | 電離放射線のピクセル化検出器からの信号の解析 | |
Zhu et al. | Sub-pixel position sensing for pixelated, 3-D position sensitive, wide band-gap, semiconductor, gamma-ray detectors | |
JP5914381B2 (ja) | X線データ処理装置、x線データ処理方法およびx線データ処理プログラム | |
Liao et al. | Pulse shape discrimination using EJ-299-33 plastic scintillator coupled with a Silicon Photomultiplier array | |
JP5425071B2 (ja) | パルスプロセッサのエネルギー測定フィルタの応答を調整する方法およびこの方法を実行するパルスプロセッサ、エネルギー分散型放射線分光分析システム | |
WO2014173812A1 (en) | Pulse processing circuit with correction means | |
WO2012029496A1 (ja) | 放射線計測装置 | |
Warburton et al. | Digital pulse processing: new possibilities in nuclear spectroscopy | |
US10641909B2 (en) | Method for processing a pulse generated by a detector of ionizing radiation | |
Abbene et al. | Performance of a digital CdTe X-ray spectrometer in low and high counting rate environment | |
Al-Adili et al. | On the Frisch–Grid signal in ionization chambers | |
US8927939B2 (en) | Radiation measurement apparatus and nuclear medicine diagnosis apparatus | |
CN110998368A (zh) | 电离辐射的像素化检测器中的深度校正 | |
JP2013512443A5 (ja) | ||
Lipschutz et al. | Digital data acquisition for the low energy neutron detector array (LENDA) | |
Stevanato et al. | High rate read-out of LaBr (Ce) scintillator with a fast digitizer | |
Bruyneel et al. | Crosstalk corrections for improved energy resolution with highly segmented HPGe-detectors | |
Alharbi | Energy resolution improvement of CdTe detectors by using the principal component analysis technique | |
Kim et al. | Impact of drift time variation on the Compton image from large-volume CdZnTe crystals | |
Papp et al. | A robust digital signal processor: Determining the true input rate | |
Guthoff et al. | Instrumentation for beam radiation and luminosity measurement in the CMS experiment using novel detector technologies | |
Nakhostin et al. | Determination of gas amplification factor by digital waveform analysis of avalanche counter signals | |
Sabbatucci et al. | First principles pulse pile-up balance equation and fast deterministic solution | |
Nakhostin et al. | Time walk correction of CdTe detectors using depth sensing technique | |
Wonders et al. | Characterization of a mixed-sinusoid multiplexing scheme with silicon photomultipliers and an inorganic scintillator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180517 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190305 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20191015 |