(57)【要約】本発明は、持続時間ＤとエネルギーＥからなるメインインパルスの連鎖を含んだ信号の処理を含む測定方法を提案するものであり、該メインインパルスが時間的に互いに間隔を開けられ、それぞれが、持続時間Ｄｉと、エネルギーが加算性の特性を有している変数Ｅｉによって評価される基本インパルスのパイルアップで構成され得るものであり、前記基本インパルスが強度λの均質なポワソン過程に従った出現時間Ｔｉを有している測定方法であって、前記信号をディジタル化することからなる過程と、各メインインパルスそれぞれの持続時間ＤとエネルギーＥを測定することで持続時間−エネルギーのペア（Ｄ、Ｅ）を作成するようにすることからなる過程と、作成されたペア（Ｄ、Ｅ）から基本インパルスのエネルギーのペア（Ｄｉ、Ｅｉ）を判定することからなる過程と、判定されたペア（Ｄｉ、Ｅｉ）から、各基本インパルスそれぞれのエネルギーＥｉを演繹することからなる過程を含んでいることを特徴とする測定方法である。【選択図】図６
2007538241.tif 000002
【特許請求の範囲】
【請求項１】
持続時間ＤおよびエネルギーＥからなるメインインパルス（１７）の連鎖を含んだ信号（１１）処理を含んだ測定方法であり、該メインインパルスが時間的に互いに間隔を開け、それぞれが、持続時間Ｄｉと、エネルギーが加算性の特性を有している変数Ｅｉによって評価される基本インパルス（１５）のパイルアップで構成され得るものであり、前記基本インパルス（１５）が強度λの均質なポワソン過程に従った出現時間Ｔｉを有している測定方法であって、
−前記信号をディジタル化することで、該信号を表すデータ（２４）を得るようにすることからなる過程と、
これらのデータ（２４）に基づき、
−各メインインパルス（１７）それぞれの持続時間ＤとエネルギーＥを測定することで持続時間−エネルギーのペア（Ｄ、Ｅ）を作成することからなる過程と、
−作成されたペア（Ｄ、Ｅ）から基本インパルス（１５）のエネルギーのペア（Ｄｉ、Ｅｉ）を判定することからなる過程と、
−判定されたペア（Ｄｉ、Ｅｉ）から、各基本インパルス（１５）それぞれのエネルギーＥｉを演繹的に導くことからなる過程、
を含んでいることを特徴とする、測定方法。
【請求項２】
持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）の判定過程が、メインインパルスの持続時間−エネルギーのペア（Ｄ、Ｅ）の関数を基本インパルスの持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）の関数に関係づけるパイルアップリジェクションの式を解くことで実施されることを特徴とする、請求項１に記載の測定方法。
【請求項３】
前記式がメインインパルスの持続時間−エネルギーのペア（Ｄ、Ｅ）に対する離散確率法則を、基本インパルスの持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）に対する離散確率法則に関係づけることを特徴とする、請求項２に記載の測定方法。
【請求項４】
持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）の判定過程において、持続時間−エネルギーのペア（Ｄ、Ｅ）に対する離散確率法則が、作成された持続時間−エネルギーのペア（Ｄ、Ｅ）と前記各ペア（Ｄ、Ｅ）に関連づけられた発生数を用いて算定され、前記発生が前記測定に基づいて事前に判定されていることを特徴とする、請求項３に記載の測定方法。
【請求項５】
さらに前記式がポワソン過程の強度λに依存していることを特徴とする、請求項２〜４のいずれか一つに記載の測定方法。
【請求項６】
前記パイルアップリジェクションの式が以下の形をとり、
【数１】
2007538241.tif 000017
ここに、変数ｚおよびｓが持続時間およびエネルギーのＺ変換を表し、ｔが持続時間、Ｂ（ｚ、ｓ）が確率密度ｂt、eのＺ変換を表し、ｂt、eが、メインインパルスの持続時間が値ｔと等しくなり、エネルギーが値ｅと等しくなる確率を表し、Ｋt（ｓ）が確率密度ｈt、eに依存した関数のＺ変換であり、ｈt、eが、基本インパルスの持続時間が値ｔに等しくなり、エネルギーが値ｅに等しくなる確率を表し、そして、αが分析される光子源の活性に依存したパラメータを表すことを特徴とする、請求項２〜５のいずれか一つに記載の測定方法。
【請求項７】
ポワソン過程の強度λがメインインパルス（１７）の分離時間を測定することで判定されることを特徴とする、請求項１〜６のいずれか一つに記載の測定方法。
【請求項８】
パラメータαが、以下の式：
α＝ｅｘｐ（−λＴｅ）
ここに、Ｔｅがディジタル化過程に関連したサンプリング周期である、
を用いて判定されることを特徴とする、請求項６または請求項７に記載の測定方法。
【請求項９】
各基本インパルス（１５）のエネルギーＥｉの演繹過程が、これらのエネルギーの離散確率法則の判定過程を含んでいることを特徴とする、請求項１〜８のいずれか一つに記載の測定方法。
【請求項１０】
エネルギーの離散確率法則が、エネルギーに応じている、持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）の離散確率法則の周辺確率法則であることを特徴とする、請求項１〜９のいずれか一つに記載の測定方法。
【請求項１１】
前記パイルアップリジェクションの式における等式の右辺の分数を持続時間に従った、冪級数で展開することで、この級数の係数を判定するようになっていることをさらに特徴とする、請求項６〜１０のいずれか一つに記載の測定方法。
【請求項１２】
前記冪級数の係数が前記パイルアップリジェクションの式における等式の左辺項に含まれる係数αt-Kt(s)と同一であると識別されることで、基本インパルス（１５）の持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）に関する前記確率密度ｈt、eを判定するようになっていることを特徴とする、請求項１１に記載の測定方法。
【請求項１３】
冪級数の係数の前記判定および係数の前記識別過程がそれぞれ、少なくとも一つの離散畳込みの計算を含んでいることを特徴とする、請求項１１または請求項１２に記載の測定方法。
【請求項１４】
係数の前記識別過程がさらに肯定性の制約テスト（２１７）を含んでいることを特徴とする、請求項１２または請求項１３に記載の測定方法。
【請求項１５】
メインインパルスがメイン流束を表し、該メイン流束のそれぞれが基本流束のパイルアップを含み得るものであり、各基本流束が基本インパルスで表されることを特徴とする、請求項１〜１４のいずれか一つに記載の測定方法。
【請求項１６】
流束が光子束（３）であることを特徴とする、請求項１５に記載の測定方法。
【請求項１７】
光子（３）がガンマ光子であることを特徴とする、請求項１６に記載の測定方法。
【請求項１８】
光子（３）の各流束のエネルギーが、対応するインパルスのエネルギーによって算定されることと、各インパルスそれぞれのエネルギーが、各対応するインパルスそれぞれの下のスペクトル面（１３）に応じて算定されることを特徴とする、請求項１６または請求項１７に記載の測定方法。
【請求項１９】
物理現象を表すメインインパルスの連鎖を含んだ前記信号を測定する初期過程を含んでいることを特徴とする、請求項１〜１８のいずれか一つに記載の測定方法。
【請求項２０】
前記信号に関する情報を利用者に提供することからなる少なくとも一つの過程を含んでいることを特徴とする、請求項１〜１９のいずれか一つに記載の測定方法。
【請求項２１】
前記情報が基本流束に関するものであることを特徴とする、請求項１〜２０のいずれか一つに記載の測定方法。
【請求項２２】
前記情報が基本流束を表す基本インパルスのエネルギーＥｉに関するものであることを特徴とする、請求項２１に記載の測定方法。
【請求項２３】
持続時間ＤとエネルギーＥからなるメインインパルス（１７）の連鎖を含んだ信号（１１）の分析装置であり、該メインインパルスが互いに時間的に間隔を開けられ、持続時間ＤｉとエネルギーＥｉからなる基本インパルス（１５）のパイルアップから構成され得るものであり、前記基本インパルス（１５）が強度λの均質なポワソン過程に従った出現時間Ｔｉを有している分析装置であって、前記装置が、請求項１〜２２のいずれか一つに記載の測定方法を実施することに適した手段を具備していることを特徴とする分析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、全体として、複数のメインインパルスを含んだ信号を分析することに関するものであり、該メインインパルスはそれぞれが識別しようとしている未知の基本インパルスのパイルアップで構成され得るものである。
【０００２】
より特定的には、本発明は、持続時間ＤとエネルギーＥからなるメインインパルスの連鎖を含んだ信号の処理を含む測定方法に関するものであり、該メインインパルスは時間的に互いに間隔を開けられ、それぞれが持続時間Ｄｉで、エネルギーが数学用語での加算性の特性を有する変数Ｅｉによって評価される基本インパルスのパイルアップで構成され得るものであり、前記基本インパルスは、強度λ（原語は、intensite 。日本では、「平均」と呼称することが多い、ポワソン分布のパラメータ。以下同様）の均質なポワソン過程に従った出現時間Ｔｉを有している。
【背景技術】
【０００３】
一般的に、このような信号は流束、特に光子束を表すことができる。
【０００４】
既に、このタイプの測定方法は、とりわけガンマ線スペクトロメトリーの領域で知られている。
【０００５】
ガンマ線スペクトロメトリーは、一方では、ガンマ光子の放出源に含まれる放射性核種を特徴づけることを目的としており、他方では、この放出源の活性（該活性は典型的には、壊変毎秒の数値またはベクレルと呼ばれる既知の単位によって定義される）を測定することを目的としていることを確認しておく。
【０００６】
このような処理を実施することのできる装置の例が、図１に示されている。
【０００７】
電源２を備えた検出器１は、光子信号３を電気信号４に変換する。
【０００８】
次に、この電気信号は、前置増幅器５によって処理されるのだが、該前置増幅器によって信号対雑音比を高めることが可能になっている。
【０００９】
次に、取得用の電子回路ユニット７の少なくとも一つの入力に適した信号を提示するために増幅器６が用いられる。
【００１０】
このような電子回路ユニットは、典型的にはアナログ−ディジタル変換器で構成される。
【００１１】
この電子回路ユニットの出力で、そのようにしてディジタル化された信号が処理用ユニット８に供給され、該処理用ユニットにおいて、特定の操作が実施される。
【００１２】
例えば、処理用ユニット８は、とりわけ、信号のディジタル・フィルタを実施するのに適した電子回路を含むことができる。
【００１３】
最後に、視覚化システム９によって、利用者がガンマ線放出源の構成を分析するのを助けて、前記の取得の流れが完了する。
【００１４】
図２は、典型的には、このような視覚化システムで観察できるものを示している。
【００１５】
この図２では、それぞれがエネルギーを数量化した値を示しているさまざまなスペクトル線を見ることができる（例えばスペクトル線１０）。
【００１６】
これらのスペクトル線の集合が、特定のガンマ線を放出する放射性核種を特徴づけている（たとえば、セシウム１３７は６６２ｋｅＶの単エネルギーのスペクトル線を有しており、また、コバルト６０は１１７３ＭｅＶと１３３３ＭｅＶの二つの単エネルギーのスペクトル線を有している）。
【００１７】
しかし、このようなスペクトロメトリー装置を用いた放射性核種の識別、より一般的にはガンマ線放出源の識別は、いくつかの問題を呈することがある。
【００１８】
とりわけ、放射性核種に関連する理想的なスペクトル線のスペクトルを劣化させうる三つのタイプの擾乱が知られているのだが、該タイプとは、
−コンプトン効果と、
−理想的な信号に測定ノイズが加わることと、
−それぞれが光子を表すインパルスの潜在的なパイルアップである。
【００１９】
コンプトン効果と測定ノイズ
光子が検出器との相互作用に入ると、電気インパルスが生じ、そして該電気インパルスは先に図１で説明したように増幅される。
【００２０】
このような検出器で生成された経時信号１１の例は、図３に示されている。
【００２１】
この経時信号が、異なった振幅と持続時間の電気インパルスの連続で構成され、各インパルスそれぞれが理想的には唯一の光子と検出器の相互作用に対応していることを観察することができる。
【００２２】
このような経時信号１１によって検出された光子のエネルギーの測定が可能になるのだが、それは、電気インパルスのスペクトル面（たとえばインパルス１２のスペクトル面１３）が、対象となるインパルスに関連した光子のエネルギーに比例するからである。
【００２３】
したがって、さまざまな電気インパルスの経時信号１１から、観察されている光子源のエネルギーのスペクトル線の集合を演繹的に導くことができる。
【００２４】
しかし、二つの現象によってこの測定原理は擾乱を受ける。
【００２５】
第一には、コンプトン効果（この効果は当業者には既知である）のために、光子エネルギーの一部だけが検出器にかけられることがあるのである。
【００２６】
よって、このような効果を受けた光子のスペクトロメトリー分析では、実際に観察されるはずのエネルギーより低いエネルギーのスペクトル線が生じる。
【００２７】
第二に、光子のエネルギーの全体がかけられた場合であっても、前記光子と検出器の間における相互作用と、測定電子機器によってもたらされるノイズが本質的にランダムであるという性質のため、創出された電気インパルスのスペクトル面は前記エネルギーに対してほぼ比例に近いというにすぎない。
【００２８】
図４、５は、とりわけ、この二つの擾乱現象を示している。
【００２９】
図４は、本質的に単エネルギーを備えた光子を放出する光子源に由来する信号のエネルギー分布を示しており、該単一のエネルギーはピーク１４に対応している。
【００３０】
ここでは、検出器から出た信号が単エネルギー性の単一のスペクトル線による理想的な形状では現れず、ピーク１４のエネルギーより小さいエネルギーで、コンプトン効果に起因する、重要な連続した低い部分１６も呈していることが確認できる。
【００３１】
図５は、図２の場合に該当するような、複数のエネルギーを有する光子を放出するガンマ線光子源の検出に由来する信号のエネルギー分布を示している。
【００３２】
ここでもまた、とりわけコンプトン効果と測定ノイズによって誘発された二つの現象に由来する、ノイズの上限の存在が確認できる。
【００３３】
このノイズの上限は、エネルギーが横座標上で高まるにつれ、顕著に低下している。
【００３４】
この図５からは、分析されたガンマ線源に存在する放射性核種の識別がここでは同様に難しいことが明らかであり、測定されたスペクトルが理想的に得られるはずのスペクトルとはかけ離れた形状を有している。
【００３５】
インパルスのパイルアップ
ガンマ線源をスペクトロメトリーで分析するときのもう一つの擾乱の原因は、時間インパルスのパイルアップの問題に関わる。
【００３６】
唯一の光子が検出器と相互作用に入ると、生成されるインパルスは短い持続時間ではあるが、ゼロではない。
【００３７】
たとえば、高速の検出器の場合、この持続時間はマイクロ秒をちょうど下回る。
【００３８】
二つの光子が近接した瞬間に放出されるとき、遅速の検出器によって生成されたそれぞれのインパルスが部分的に、さらには全体的に重なることも起こりうる。
【００３９】
この現象は、とりわけ、タイプ２のカウンタタイプの検出器を用いるときに知られている。
【００４０】
これらのカウンタでは、別のインパルスが検出中であるのに対し、到着する各インパルスそれぞれが、形成途中のインパルスをさらに延長する効果を有している。
【００４１】
この現象の例は、図６に示されている。
【００４２】
この図６は、それぞれＴｎ＋１およびＴｎ＋２という時間に到着し、Ｕｎ＋１とＵｎ＋２という持続時間をした二つのインパルス１５、１５’の重複を含んだ時間信号を表している。
【００４３】
インパルス１５’はインパルス１５が終わっていないにも関わらず現れている。
【００４４】
検出器は、ここではインパルスの持続時間に対して比較的遅い反応を示しており、二つのインパルスを互いに区別していないことが分かる。
【００４５】
実際、該検出器は二つのインパルス１５および１５’を融合して、持続時間Ｕｎ１、２のメインインパルス１７しか形成しておらず、このことは、唯一の光子が検出器と相互関係に入ったと錯覚させることになる。
【００４６】
本明細書の以下の部分では明確さのために、基本インパルスは唯一の光子に結びつけられたインパルス（たとえば基本インパルス１５または１５’）を指すこととする。
【００４７】
そして、メインインパルスは、基本インパルスのパイルアップを含み得るようなインパルスに対応することとする。
【００４８】
したがって、パイルアップ現象の第一の帰結は、分析中のガンマ光子源の全体的な活性を過小算定するおそれがある。
【００４９】
ここで注意すべきは、光子源の全体的な活性が強く、または、検出器の反応が遅い分だけ、それだけこの現象が強くなることである。
【００５０】
本明細書の以下の部分では、「全体的な活性」および「計数率」という用語の定義が同一であると見なすこととする。
【００５１】
このような光子のパイルアップの第二の帰結は、厳密な意味でエネルギーの算定に関わる。
【００５２】
実際、このような検出器の場合、単一の光子にパイルアップされた光子エネルギーの和が割り当てられるのだが、これは、定義が同一であることからして、メインインパルス１７が基本インパルス１５および１５’のスペクトル面より大きなスペクトル面を有しているからである。
【００５３】
一般的には、パイルアップ現象はスペクトルの一部を間違って高いエネルギーに向けて移動させる。
【００５４】
例えば、図７と図８ではこのような移動が確認できるわけだが、該図面はそれぞれ、パイルアップ現象が起こっていない、２２ｋｅＶと５５ｋｅＶという二つのメインのスペクトル線を含んだ正規化スペクトルと、パイルアップ現象が存在する、同一のガンマ線源に由来する正規化スペクトルを示している。
【００５５】
分析されたガンマ線源は本来、図７で完全に識別されている二つのエネルギーＥ１およびＥ２を備えた光子を放出することに注意すべきである。
【００５６】
この場合は計数率の低いスペクトロメトリー分析に対応する（たとえば光子源から離れて配置された小さな検出器）。
【００５７】
逆に、計数率が高くなるときパイルアップ現象が現れ、エネルギーＥ１およびＥ２の線形結合に対応しているエネルギーを備えた寄生スペクトル線が徐々に出現するのが観察される（たとえば、２Ｅ１、Ｅ１＋Ｅ２、２  Ｅ２、２Ｅ１＋２Ｅ２など。たとえばスペクトル線１８および１９）。
【００５８】
さらに、これらの寄生スペクトル線の出現は、高エネルギー方向へ信号が散逸することから、エネルギーＥ１とＥ２に関連する活性のひずみを引き起こす。
【００５９】
これらの寄生スペクトル線が観察されるガンマ線源の放射性核種の識別過程を擾乱し、さらには、いくつかのスペクトル線が光子源によって放出されるエネルギーのどれにも対応せず、実際に放出されたスペクトル線の活性が過小評価されるという、結果を誤らせる危険性があることが容易に理解できる。
【００６０】
上記に説明した問題全体を抑制し、特にパイルアップに関連する問題を抑制するために多くの解決策が提案されている。
【００６１】
たとえば、とりわけ線形フィルタリングの操作を含んだスペクトロメータで実装されている測定方法（非特許文献１参照）が知られている。
【００６２】
提案されているフィルタリングの目的は、基本インパルスのパイルアップで形成されているすべてのメインインパルスがノイズの一部をなすという了解の上で、有効な信号に加わった前記ノイズをできる限り除去することである。
【００６３】
ある程度の改善が可能にはなるのだが、この測定方法は特化したもので、その利用範囲を制限するという妥協を必要とする。
【００６４】
実際、前記信号に存在するノイズの性質のために、提案されているフィルタリングは長時間の操作であればそれだけ効果的になるのだが、これは必然的に、パイルアップ現象が含まれる危険性を高めることになる。
【００６５】
結果的に、この測定方法は分析されるガンマ線源の活性に強く依存することになる。
【００６６】
この活性が弱ければ、パイルアップの危険性は下がるので、効果的なフィルタリングを実施することができる。
【００６７】
逆に、この活性がある閾値を超えるときには、フィルタリングが実施される前記時間を短くし、処理すべき信号におけるパイルアップの数を抑えるようにしなければならないのだが、これは、ノイズ除去の効果を減じ、したがって、最終的に得られるスペクトルの分解能を下げることになる。
【００６８】
結果的に、このような測定方法を利用したスペクトロメータは、ガンマ線源の計数率が上がれば、分解能の劣化にさらされることになる。
【００６９】
ＣＡＭＢＥＲＲＡ／ＥＵＲＩＳＹＳ、ＯＲＴＥＣまたはＸＩＡのような、多くの製造者が、パイルアップの問題を克服することのできるディジタル・スペクトロメータも開発している。
【００７０】
しかし、インパルスの処理方法は概念的に従来の測定方法に近いままであり、該測定方法はやはり、上記で提示したフィルタリングのような、ＦＩＲ（Ｆｉｎｉｔｅ  Ｉｎｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅの頭字語）タイプの、非反復的線形フィルタリングに基づいている。
【００７１】
これはとりわけ、製造者ＸＩＡの「Ｄｉｇｉｔａｌ  Ｇａｍｍａ  Ｆｉｎｄｅｒ」というスペクトロメータに当てはまる。
【００７２】
この点について、詳細は、参考文献（特許文献１、２、３、４、５）を参照のこと。
【００７３】
パイルアップの問題を部分的に解決することを可能にするその他の測定方法もある。
【００７４】
この測定方法は検出された信号の形状分析に基づくもの（特許文献６）または、所定の閾値に対して累積したエネルギーを比較することに基づくもの（特許文献７，８）である。
【００７５】
しかし、これらの測定方法は用いられる検出器のタイプに強く依存しており、結局、計数率が高いときにはあまり安定していない。
【００７６】
また、特定の検出器（とりわけテルル化カドミウムをベースとする検出器）によってもたらされるひずみを修正することを可能にする最近のもう一つの測定方法（特許文献９）も知られており、該検出器は、スペクトル面が光子によってもたらされるエネルギーに比例するインパルスは生成しないというものである。光子によってもたらされるエネルギーは二つの変数に基づいて算定される。すなわち、一つはスペクトル面であり、もう一つは、前記光子に関連したインパルスの上昇時間である。
【００７７】
しかし、上記に挙げたスペクトロメトリーのあらゆる測定方法と同じように、この測定方法は実験結果にのみ基づいている。
【００７８】
さらに、これらの測定方法は一般的に、処理された信号から、検出したパイルアップインパルスを除去することで十分とされている。
【００７９】
最終的な信号は、パイルアップのないものではあるが、いくつかの基本インパルスおよび該基本インパルスが表す情報が前記除去に伴って消えてしまう。
【００８０】
このように、これらの測定方法の不都合は、メインインパルスが表す情報を分析するために、該メインインパルスに含まれたすべての基本インパルスを識別することができないということである。
【非特許文献１】Ｇ．Ｆ．Ｋｎｏｌｌ『Ｒａｄｉａｔｉｏｎ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｍｅａｓｕｒｅｍｅｎｔ』Ｗｉｌｅｙ、第二版、１９８９年
【特許文献１】米国特許第５，８７３，０５４号明細書
【特許文献２】米国特許第５，７７４，５２２号明細書
【特許文献３】米国特許第５，６８４，８５０号明細書
【特許文献４】米国特許第５，３３７，６０３号明細書
【特許文献５】米国特許第５，２０６，１７４号明細書
【特許文献６】米国特許第５，８８４，２３４号明細書
【特許文献７】米国特許第５，０６７，０９０号明細書
【特許文献８】米国特許第５，３４９，１９３号明細書
【特許文献９】米国特許第５，８２１，５３８号明細書
【発明の開示】
【発明が解決しようとする課題】
【００８１】
したがって、本発明の目的はこの不都合を克服することである。
【００８２】
この目的を達成するため、本発明は、とりわけガンマ放射能の分野において、実験的アプローチではなく、厳正なアプローチに基づいた効果的な識別方法を提案するものである。
【課題を解決するための手段】
【００８３】
よって、本発明は持続時間ＤとエネルギーＥからなるメインインパルスの連鎖を含んだ信号の処理を含む測定方法を提案するものであり、該メインインパルスは時間的に互いに間隔を開けられ、それぞれが、持続時間Ｄｉと、エネルギーが加算性の特性を有した変数Ｅｉによって評価される基本パルスのパイルアップで構成され得るものであり、前記基本インパルスが、強度λの均質なポワソン過程に従った出現時間Ｔｉを有している測定方法であって、
−前記信号をディジタル化して、該信号を表すデータを得られるようにすることからなる過程と、
これらのデータに基づき、
−各メインインパルスそれぞれの持続時間ＤおよびエネルギーＥを測定して、持続時間−エネルギー（Ｄ、Ｅ）のペアを作成することからなる過程と、
−作成されたペア（Ｄ、Ｅ）から基本インパルスのエネルギーのペア（Ｄｉ、Ｅｉ）を判定することからなる過程と、
−判定されたペア（Ｄｉ、Ｅｉ）から各基本インパルスそれぞれのエネルギーＥｉを演繹的に導くことからなる過程を含んでいることを特徴としている。
【００８４】
この測定方法の非制限的ないくつかの好ましい特徴は、以下のとおりである。
【００８５】
−持続時間−エネルギー（Ｄｉ、Ｅｉ）のペアを判定する過程が、メインインパルスの持続時間−エネルギー（Ｄ、Ｅ）のペアの関数を、基本インパルスの持続時間−エネルギー（Ｄｉ、Ｅｉ）のペアの関数に結びつけるパイルアップリジェクションの式を解くことで実施されること。
【００８６】
−前記式が、メインインパルスの持続時間−エネルギー（Ｄ、Ｅ）のペアの離散確率法則を、基本インパルスの持続時間−エネルギー（Ｄｉ、Ｅｉ）のペアの離散確率法則に結びつけること。
【００８７】
持続時間−エネルギー（Ｄｉ、Ｅｉ）のペアを判定する過程において、持続時間−エネルギー（Ｄ、Ｅ）の離散確率法則が、作成された持続時間−エネルギー（Ｄ、Ｅ）のペアと、前記各ペア（Ｄ、Ｅ）それぞれに関連づけられた発生数を用いて算定され、前記発生が予め前記測定から決定されていること。
【００８８】
−さらに、前記式がポワソン過程の強度λに依存すること。
【００８９】
−前記式が以下の形をしていること。
【００９０】
【数１】
2007538241.tif 000003
【００９１】
この式において、変数ｚおよびｓは持続時間とエネルギーのＺ変換を表し、ｔは持続時間、Ｂ（ｚ、ｓ）は確率密度ｂt、eのＺ変換であり、ｂt、eは、メインインパルスの持続時間が値ｔと同じになり、エネルギーが値ｅと同じになる確率を表し、Ｋｔ（ｓ）は確率密度ｈt、eに依存した関数のＺ変換であり、ｈt、eは、基本インパルスの持続時間が値ｔと等しくなり、エネルギーが値ｅと等しくなる確率を表し、そして、αは分析される光子源の活性に依存したパラメータである。
【００９２】
−ポワソン過程の強度λが、メインインパルスの分離時間を測定することで判定されること。
【００９３】
−パラメータαが、
α＝ｅｘｐ（−λＴe）、
という式を用いて判定され、該式において、Ｔｅがディジタル化過程に関連したサンプリング周期を表すこと。
【００９４】
−各基本インパルスそれぞれのエネルギーＥｉの演繹過程が、エネルギーの離散確率法則を判定する過程を含んでいること。
【００９５】
−エネルギーの離散確率法則が、エネルギーに応じた、持続時間−エネルギーペア（Ｄｉ、Ｅｉ）の離散確率法則の周辺確率法則であること。
【００９６】
−前記パイルアップリジェクションの式における等式の右辺の分数を、持続時間に従った、冪級数で展開することで、この級数の係数を決定すること。
【００９７】
−基本インパルスの持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）に関する前記確率密度ｈt、eを判定するために、前記冪級数の係数が、前記パイルアップリジェクションの式における等式の左辺項に含まれるαt-Kt(s)という係数と同一であると識別されること。
【００９８】
前記冪級数の係数決定と前記係数の識別過程がそれぞれ、少なくとも一回の離散畳込みの計算を含んでいること。
【００９９】
−前記係数の識別過程がさらに、肯定性に関わる制約のテストを含んでいること。
【０１００】
−メインインパルスが、それぞれ基本流束のパイルアップを含み得るようなメインの流束を表し、各基本流束それぞれが基本インパルスで表されること。
【０１０１】
−流束が光子束であること。
【０１０２】
−光子がガンマ光子であること。
【０１０３】
−各光子束それぞれのエネルギーが対応するインパルスのエネルギーで表され、各エネルギーそれぞれが対応する各インパルスそれぞれの下にあるスペクトル面から判定されることで表されること。
【０１０４】
本測定方法が、物理現象を表すメインインパルスの連鎖を含んだ前記信号を測定する初期過程を含んでいること。
【０１０５】
本測定方法が、前記信号に関連する情報を利用者に提供することからなる少なくとも一つの過程を含んでいること。
【０１０６】
前記情報が基本流束に関するものであること。
【０１０７】
前記情報が基本流束を表す基本インパルスのエネルギーＥｉに関するものであること。
【０１０８】
本発明のもう一つの目的は、持続時間ＤとエネルギーＥからなるメインインパルスの連鎖を含んだ信号を分析する装置を提案することであり、該メインインパルスは時間的に互いに間隔を開けられ、持続時間ＤｉとエネルギーＥｉからなる基本インパルスのパイルアップで構成され得るものであり、前記基本インパルスが強度λの均質なポワソン過程に従った出現時間Ｔｉを有するものであって、前記装置が、単一または組み合わせで上記の好ましい特徴にしたがって測定方法を実施することに適した手段を具備していることを特徴としている。
【発明を実施するための最良の形態】
【０１０９】
本発明のその他の特徴、目的および利点は、非制限的例示として示され、添付図面を参照して行われている、以下の本発明の好ましい実施態様の詳細な説明を読むことでより明確になるものであり、該添付図面において、
−図１はガンマスペクトロメトリーのディジタル取得の流れの各要素を概略的に示し、
−図２は、図１のタイプの装置によって観察可能なガンマ線源のスペクトル線のスペクトルを示しており、
−図３は、例示として、ガンマ光子の検出器によって生成された時間信号の例を示しており、
−図４は、例示として、セシウム１３７の単エネルギー光子源を測定したエネルギー分布を示しており、
−図５は、例示として、複数のエネルギーを備えた光子を放出する活性源のエネルギー分布を示しており、
−図６は、パイルアップ現象が起こったときの、遅速検出器に起因する時間信号を概略的に示しており、
−図７は、いかなるパイルアップ現象も生じていないときの、ガンマ線源の正規化されたスペクトルを示しており、
−図８は、活性の増加によってパイルアップ現象が起こったときの、図７で分析されたガンマ線源の正規化されたスペクトルを示しており、
−図９は、本発明によるシステムによって測定された、二つのメインインパルスを含んだ信号の一部を概略的に示しており、
−図１０は、例示として、本発明によるシステムによってディジタル化されたときの、図９の信号を示しており、
−図１１は、本発明に係る測定方法のユニット略図を示しており、
−図１２は、本発明に係る測定方法によって処理された信号の活動シーケンスのＤとＥのペアの二次元ヒストグラムの例を示しており、
−図１３は、前記処理された信号の基本インパルスのＤｉとＥｉのペアの二次元ヒストグラムの例を示しており、
−図１４は、適切な電子装置において本発明によるパイルアップリジェクション方法の第一部分を実装することを可能にするアルゴリズムのフローチャートを示しており、
−図１５は、適切な電子装置において本発明によるパイルアップリジェクション方法の第二部分を実装することを可能にするアルゴリズムのフローチャートを示しており、
−図１６は、本発明に係る測定方法を実施していないときの、セシウム１３７の観察されたガンマ線エネルギーのスペクトルを示しており、
−図１７は、本発明による測定方法を実施しているときの、セシウム１３７の観察されたガンマ線エネルギーのスペクトルを示している。
【０１１０】
前提として規定するのは、本明細書の以下の部分では、活動シーケンスとは、少なくとも一つの基本インパルスが常に存在しているとすること、つまり、各時間それぞれにおいて、信号の振幅が、たとえば検出器のホワイトノイズのレベルに対応する閾値を超えている、この信号の連続した部分を指すこととする。
【０１１１】
結局、一つの活動シーケンスは、一つのメインインパルスに対応する。
【０１１２】
ただし、この新しい概念によって規定できるのは、以下では非常に有用なのだが、非活動シーケンスを、信号が所定の閾値（閾値は、やはり検出器のノイズレベルとすることができる）以下となっているこの信号の連続した部分である。
【０１１３】
たとえば、図９は、本発明によるシステムによって測定された信号の部分を概略的に示しているが、前記信号は二つの活動シーケンス２０および２１を含んでいる。
【０１１４】
活動シーケンス２０は、パイルアップを含んでいないのに対し、シーケンス２１は、二つの基本インパルスのパイルアップを含んでいる。
【０１１５】
また、注意すべきことには、インパルスの持続時間ＤあるいはＤｉが、前記インパルスの全体の幅に対応する時間を意味していることである。
【０１１６】
この点に関し、本明細書の以下の部分では、
−持続時間：一時的な式ではｔ、または式がＺ変換されているときにはｚ、
−エネルギー：一時的な式ではｅ、または式がＺ変換されているときにはｓ、
という表記によって、式における持続時間とエネルギーを指すこととする。
【０１１７】
ここで図１０を参照すると、本発明に係るシステムは、測定過程において検出器によって生成された信号の集合をサンプリングし、ディジタル化することで、この信号を表すディジタルデータの流束２４、あるいは、より一般的には、関連しているガンマ光子束が得られるようになっている。
【０１１８】
なお、留意すべきは、前記サンプリングは、この説明においては式の記述ならびにその理解を容易にするために、値が１に固定されるサンプリング周期Ｔｅで行われていることである。
【０１１９】
このディジタルデータの流束は、図１における処理用ユニットによって処理される。
【０１２０】
利用される回路として典型的なものは、プロセッサ、ＤＳＰ（英語の表現「Ｄｅｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ」の頭字語）、ＦＰＧＡ（英語の表現「Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ」の頭字語）または他のあらゆる同等の回路といった、ディジタルデータに対する操作を行うことに適した回路である。
【０１２１】
この処理用ユニットでは、測定された信号に含まれるインパルスをパイルアップリジェクションするために、本発明による測定方法が実装される。
【０１２２】
ここで「パイルアップリジェクション」は、基本インパルスのパイルアップによって生じる問題を部分的または全体的に解決するための、メインインパルスに含まれる情報を処理すること全体を意味している。
【０１２３】
本発明で提案するパイルアップリジェクションの方法は、
−基本インパルスが強度λの均質なポワソン過程に従った出現時間を有していることと、
−エネルギーを表す変数が加算性の特性を有している、
という仮説に基づいている。
【０１２４】
本発明に係る測定方法のユニット略図が、図１１に示されている。
【０１２５】
まず、活動シーケンスを表すディジタルデータを考えることにする。
【０１２６】
とりわけ、持続時間ＤならびにエネルギーＥが算定され、それぞれが特徴的な活動シーケンスに関連づけられた持続時間−エネルギーのペア（Ｄ、Ｅ）の集合３０が作成される。
【０１２７】
活動シーケンスのエネルギーＥは、このシーケンス下に存在する領域に基づいて評価されることに注意すべきである。
【０１２８】
この領域はディジタルデータの振幅と得られた時間分解能に基づいて計算されるのだが、後者はとりわけ、図１における取得用ユニットの性能に依存している。
【０１２９】
領域の算定精度が時間分解能によって、したがって、オーバーサンプリングの係数によって高くなることは明らかである。
【０１３０】
次に、過程の第一段階（ユニット２６および２８）は、ポワソン過程の前記強度λを判定することから構成される（このパラメータλはユニット２７で必要になる）。
【０１３１】
このため、活動シーケンスのペア（Ｄ、Ｅ）（ユニット２６）とこれらの活動シーケンスの到着時間から、各活動シーケンスそれぞれの間の分離持続時間が計算される。
【０１３２】
そして、活動シーケンスを分離する持続時間に基づいて前記強度λが算定される（ユニット２８）。
【０１３３】
この算定が可能なのは、本出願人が、これらの時間がパラメータλの指数確率の法則に従うことを特定しているからである。
【０１３４】
したがって、このパラメータは、持続時間Ｄに基づいて算定することのできる分離時間と測定されたメインインパルスの発生時間に基づいて指数法則を特徴づけることで算定することができる。
【０１３５】
ユニット２５は、ペア（Ｄ、Ｅ）によって発生数を計算し、活動シーケンスに関わる存在についての情報を得られるようにすることからなる。
【０１３６】
たとえば、とりわけ１０ペア（０．４ｍｓ、１Ｍｅｖ）、２ペア（０．４ｍｓ、１．２Ｍｅｖ）および４０ペア（１ｍｓ、１Ｍｅｖ）を得ることができる。
【０１３７】
こうしてペア（Ｄ、Ｅ）による発生数を含んだデータの表が得られることで、統計分布の算定推定量を得ることができる。
【０１３８】
このような分布の可能な表示が図１２に示されている。
【０１３９】
この図面において、処理された信号における活動シーケンスのＤとＥのペアによる二次元ヒストグラムを確認することができる。
【０１４０】
このようなヒストグラムがここでは例示として示され、この例がまったく限定的でないことは明らかである。
【０１４１】
この図面では、指標時間１４と指標エネルギー２８による活動シーケンスに対応するピーク１００をはっきり観察することができる。
【０１４２】
測定によると、この活動シーケンスは検出器から出た信号においておよそ１００００回出現している。
【０１４３】
このように、ある観点からでは、ユニット２５は持続時間とエネルギーによる二次元ヒストグラムを表すデータの表を作成することから構成される。
【０１４４】
次に、もう一つのユニット２７はペア（Ｄ、Ｅ）を基本インパルスの持続時間−エネルギーのペア（Ｄｉ、Ｅｉ）に結びつけることからなる。
【０１４５】
このように処理することで、検出器から出た信号に含まれる基本インパルスに関わる情報が初めて得られるようになる。
【０１４６】
より的確には、メインインパルスから基本インパルスへの移行は、メインインパルスのペア（Ｄ、Ｅ）の離散確率法則を基本インパルスのペア（Ｄｉ、Ｅｉ）の離散確率法則に結びつける式を用いて行われる。
【０１４７】
パイルアップリジェクションの式として指示されるこのような式は、以下のように示される。
【０１４８】
【数２】
2007538241.tif 000004
【０１４９】
・ｔは持続時間を表し、変数ｚおよびｓは持続時間ｔとエネルギーｅのＺ変換を表す。
【０１５０】
・Ｂ（ｚ、ｓ）は、メインインパルスの持続時間が持続時間ｔの指数で表される値に等しくなり、エネルギーが指標ｅ（ｔおよびｅは、冪である）で表される値に等しくなる確率密度ｂt、eのＺ変換を表す。
【０１５１】
Ｂ（ｚ、ｓ）の方程式は以下の形をとる。
【０１５２】
【数３】
2007538241.tif 000005
【０１５３】
・Ｋｔ（ｓ）は、基本インパルスの持続時間が持続時間の指数ｔで表される値に等しくなり、エネルギーが指数ｅで表される値に等しくなる確率密度ｈt、eに依存する関数である。
【０１５４】
この関数は以下の形で表される。
【０１５５】
【数４】
2007538241.tif 000006
【０１５６】
最終的に、パラメータαは分析された光子源の活性に結びつけられ、とりわけ、ポワソン過程の前記強度λに依存している。
【０１５７】
こうして、ユニット２６および２８でこのパラメータλを決定したことの利点が理解される。
【０１５８】
前記パラメータαに関し、パラメータλに応じて該パラメータを決定することは、
α＝ｅｘｐ（−λＴｅ）            （４）
という式（４）のおかげで実施され、該式において、Ｔｅは検出器から出た信号をサンプリングする信号周期である。
【０１５９】
ユニット２７が終了すると、信号に含まれる基本インパルスのペア（Ｄｉ、Ｅｉ）に関わるデータの表３２が得られる。
【０１６０】
換言すれば、この段階で、前記基本インパルスの持続時間とエネルギーの新しい二次元ヒストグラムを作成することができるのである。
【０１６１】
このようなヒストグラムは非限定的な例として図１３に示されている。
【０１６２】
図１２と比較すると、インパルスの持続時間が減少していることが確認できる。
【０１６３】
基本インパルスのパイルアップを含んだメインインパルスは前記基本インパルスのために消失している。
【０１６４】
ユニット２９はペア（Ｄｉ、Ｅｉ）の離散確率法則のエネルギーに従った周辺確率法則の計算を行う。
【０１６５】
こうして、ペア（Ｄｉ、Ｅｉ）に基づき、選択されたエネルギーＥｉに関連づけられた持続時間の和を出すことで、各エネルギーＥｉそれぞれが決定される。
【０１６６】
たとえば、ユニット２７によって以下の表で表されるペア（Ｄｉ、Ｅｉ）３２を見つけることが可能になっていると考えられる。
【０１６７】
【表１】
2007538241.tif 000007
【０１６８】
あらゆる持続時間が混在する中で、一つのインパルスが信号における１ＭｅｖのエネルギーＥｉを有する確率は０．７２である（０．１＋０．２５＋０．２＋０．１７）。
【０１６９】
方程式の形で表現すると、本明細書の以下の部分においてＧeと表記される、エネルギーｅの発生に関する前記判定は、以下のようになる。
【０１７０】
【数５】
2007538241.tif 000008
【０１７１】
該式において、ｋは持続時間の指数であり、ＮＢ＿Ｄは持続時間におけるチャネル数である。
【０１７２】
本発明で提案される測定方法は、最終的に、信号に含まれる基本インパルスのエネルギーＥｉのスペクトルに関するデータ３１を得るために、適切な装置に上記のさまざまな式を実装することを必要としている。
【０１７３】
前述のように、このような装置はＤＳＰ、ＦＰＧＡ、ＡＳＩＣ（英語の表現「Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ」の頭字語）、または他のあらゆる同等の電子回路とすることができる。
【０１７４】
適切な選択は、コスト、サイズ、性能、モジュラー性といった、目標とされる用途に関わる制約に応じてなされることになる。
【０１７５】
あらゆる場合に、前記実装はとりわけ、これから説明する不可欠な二つの過程を含んでいる。
【０１７６】
もちろん、以下の説明を読めば、当業者であればその他の可能な変形例があり、ここでなされる提案がまったく限定的ではないことが理解されるものである。
【０１７７】
式（１）によると、等式の左辺は、係数が括弧に囲まれた項に対応する、冪級数の形、つまり、αt-Kt(s)であることが分かる。
【０１７８】
この利点の多い特性を考慮すると、本発明で提案される実装の第一過程は、等式の右辺項を冪級数で展開し、この冪級数の係数の集合を決定することからなる。
【０１７９】
換言すれば、係数ｙt（ｓ）は以下の式で求められる。
【０１８０】
【数６】
2007538241.tif 000009
【０１８１】
ここに、ｔは持続時間に対応する、冪数である。
【０１８２】
第一過程：冪級数の係数を求めること
このような過程は帰納関係に基づくアルゴリズムを実施する。
【０１８３】
このため、前記の等式（１）の右辺の項は、フィルタＩＩＲ（「Ｉｎｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｐｏｎｓｅ」の頭字語）のインパルス応答と考えられ、その係数ｙt（ｓ）はエネルギーに応じて決定される。
【０１８４】
このようなフィルタのＹ／Ｘ伝達関数のＺ変換によって、前記の等式の右辺の項により、
X(z,s)=Y(z,s)(1-(αz+(1-α)zB)(z,s)))、        (7)
または、
X(z,s)=Y(z,s)-αzY(z、s)-(1-α)zB(z,s)Y(z,s)、  (8)
と書くことができる。
【０１８５】
持続時間ｚにのみしたがって式（８）のＺ逆変換を計算すると、
xt(s)=yt(s)-αyt-1(s)-((1-α)b(s)*y(s))t-1             (9)
または、
yt(s)=αyt-1(s)+((1-α)b(s)*y(s))t-1+xt(s)         (10)
が得られる。
【０１８６】
【数７】
2007538241.tif 000010
【０１８７】
ここで、記号*が畳込みの演算子を意味することに注意すべきである。
【０１８８】
式（１０）によると、所与のエネルギーｓについて係数ｙt（ｓ）を反復的に決定することができる。
【０１８９】
そのため、前記フィルタのインパルス応答が対象であるので、ｘt（ｓ）は二次元のディラック関数であり、式（１１）における係数ｂt（ｓ）はメインインパルスのヒストグラムに関するデータから分かり、パラメータαは式（４）から分かる。
【０１９０】
上述した、主として式（１０）に基づく装置で実装することができるアルゴリズムのフローチャートが図１４に示されている。
【０１９１】
第一の過程２００は、フィルタの入力を表し、二次元のディラック関数のエネルギーに従ったＺ変換の値をとるｘt（ｓ）におけるいくつかの変数を初期化することから構成される。
【０１９２】
また、持続時間ｔは１に初期化されている。
【０１９３】
過程２０１では、フィルタの出力ｙt（ｓ）が初期化される。
【０１９４】
より的確には、第一の持続時間（ｔ＝０）が固定され、変数ｙ0（ｓ）のエネルギーに関するすべての値が１に初期化される。
【０１９５】
この段階で、方程式（１０）に基づいてｙt（ｓ）を計算し始めることができる。
【０１９６】
この目的を達成するために、持続時間ｔとエネルギーｓをネスティングする二つのループが利用される。
【０１９７】
持続時間に対する第一ループは過程２０２、２０３および２０４を含み、エネルギーに対する第二ループは過程２０２から２０６を含む。
【０１９８】
過程２０２では、所与のエネルギーｓについて式（１０）が決定される。
【０１９９】
たとえば、この過程２０２における最初の行程で、計算されるのはｙ1（ｓ＝０）であり、その次に、持続時間ｔのループにおける最初の繰り返しの後、決定されるのはｙ2（０）である。
【０２００】
この第一ループによって、ｙt（０）の値の集合が計算される。
【０２０１】
過程２０３は、持続時間にかかるループを続けるべきかどうか判定することを可能にするテストである。
【０２０２】
このテストは、持続時間ｔの変数を事前に設定した最大持続時間でのチャネル数（ＮＢ＿Ｄと表記）を比較することから構成される。
【０２０３】
テスト結果が肯定的であれば、過程２０４が実行され、持続時間ｔの変数がインクリメントされる。
【０２０４】
そうでなければ、アルゴリズムは前記第一のループから脱け、過程２０５で、第二ループを続けるべきかどうかテストする。
【０２０５】
したがって、この過程２０５において、変数ｓは事前に設定された最大エネルギーでのチャネル数（ＮＢ＿Ｅ）と比較される。
【０２０６】
テスト結果が肯定的であれば、アルゴリズムは過程２０６に移行し、該過程において、エネルギーｓがインクリメントされ、また、持続時間の変数ｔが過程２００と同じように値１に初期化される。
【０２０７】
次に、アルゴリズムは再び過程２０２に戻り、値ｙt（１）の集合が決定されることになる。
【０２０８】
たとえば、計算された第一の係数はｙ1（１）であり、次に、ｔの最初のインクリメントの後、過程２０２がｙ2（１）などを計算する。
【０２０９】
アルゴリズムは、エネルギーにかかる第二のループを脱けるとき、最後は過程２０７で終了する。
【０２１０】
最終的に持続時間とエネルギーのマトリクスが得られるのだが、該マトリクスの要素は、フィルタＩＩＲのインパルス応答の求められた係数ｙt（ｓ）、または換言すれば、式（６）における前記冪級数の求められた係数である。
【０２１１】
第二過程：係数の識別
第二過程は、このようにして決定された係数を、パイルアップリジェクションの式（１）における等式の左辺にある、冪級数の係数によって識別することから構成される。
【０２１２】
（１）および（６）より、以下のように書くことができる。
【０２１３】
【数８】
2007538241.tif 000011
【０２１４】
したがって、αt-Kt(s)＝ｙt（ｓ）である。
【０２１５】
活性αを（４）で置き換えることにより、
ｙt（ｓ）＝ｅｘｐ（−λｔ）ｅｘｐ（λＫt（ｓ））が得られる。
【０２１６】
ここで項λＫt（ｓ）を含んだ指数関数が、冪級数で展開される。
【０２１７】
【数９】
2007538241.tif 000012
【０２１８】
そして、選択された装置においてこのような式の実装を可能にするために、冪級数の切り捨てを行う。
【０２１９】
【数１０】
2007538241.tif 000013
【０２２０】
最終的に、等式（１５）の第二項に、エネルギーに従ったＺ逆変換を適用することで、以下の式にたどり着く。
【０２２１】
【数１１】
2007538241.tif 000014
【０２２２】
該式において、ｋt、eはエネルギーｓに従ったＫt（ｓ）のＺ逆変換である。
【０２２３】
ここでは、式（１６）で用いられている表記(e)*ｎがエネルギーｅに従ったｎ番目の畳込みを意味していることに注意すべきである。
【０２２４】
ｙt、eが分かると、係数ｋt、eは式（１６）に基づいて決定され、次に、係数ｈt、eが二重微分を行うことで関係式（３）から演繹され、そして最後に、係数ｈt、eのエネルギーにしたがって持続時間の和を出すことで、信号に含まれた基本インパルスの各エネルギーＥｉそれぞれを決定するようになっている。
【０２２５】
ここでは、選択された装置においてこの第二の識別過程を実装することを可能にする回帰的アルゴリズムが提案される。
【０２２６】
改めて、当業者であればここで提案されている実施態様には多くの可能な変形例があることが認識できるものである。
【０２２７】
前記アルゴリズムのメイン過程が例示として図１５に示されている。
【０２２８】
ここで確認するのは、これらの過程が式（１６）と（３）に主として基づいて計算できることである。
【０２２９】
第一の過程２１０では、等式（１６）の右辺項における指数ｎの和の判定を始めるために、数式ｙt、eｅｘｐ（λｔ）が決定される。
【０２３０】
したがって、この数式は、等式の左辺の指数関数を使ったときに、式（１６）からくるものである。
【０２３１】
図１５における過程２１０に示された矢印は値を変数に割り当てる操作を意味しており、この場合にはｙt、eである。
【０２３２】
次に、過程２１１は式（１６）に含まれる項ｋt*n（ｓ）を初期化することから構成される。
【０２３３】
この項の値は、図１５においてｇt（ｅ）と表記された変数に保存される。
【０２３４】
ｇt（ｅ）は、係数ｋt、eのエネルギー（ｋt、e）(e)*nに従ったｎ番目の畳込みを保存する変数である。
【０２３５】
エネルギー、持続時間、持続時間における最大チャネル数ＮＢ＿Ｄ、エネルギーにおける最大チャネル数ＮＢ＿Ｅといったその他の変数が初期化される。
【０２３６】
次に、過程２１２に基づく二つのネスティング・ループが続く。
【０２３７】
エネルギーに対する第一のループは過程２１２、２１３および２１４を含んでいる。
【０２３８】
該ループによって、ｎ番目の畳込みまたは、所与の持続時間ｔに対する係数ｋt、eの（ｋt、e）(e)*nを決定することが可能になる。
【０２３９】
第二ループは過程２１２から２２１を含んでいる。
【０２４０】
該ループは持続時間ｔを範囲とし、該ループによって、持続時間にしたがって係数ｋt、eを計算することが可能になっている。
【０２４１】
これから、これら二つのループのそれぞれを詳述することとする。
【０２４２】
上記で触れたように、エネルギーに対する第一のループは過程２１２で開始される。
【０２４３】
この過程では、エネルギーにしたがって、ｇt-1（ｅ）を用いてｙt（ｅ）の畳込みが計算され、この値がｇt（ｅ）に保存される。
【０２４４】
エネルギーの指数に対するテストが過程２１３で行われ、エネルギーに対するこの第一のループの予期しない出力を検出するようになっている。
【０２４５】
進行中のエネルギーがＮＢ＿Ｅを下回っている限り、過程２１４でこのエネルギーがインクリメントされ、過程２１２から再始動される。
【０２４６】
ひとたび、進行中の持続時間ｔに対するエネルギーにしたがってｎ番目の畳込み（ｋt、e）(e)*nが決定されると、「ｔｅｍｐ」と表記されている、式（１６）における等式の右辺にある和の一部を計算し、変数に保存する過程２１５に移行する。
【０２４７】
より的確には、変数「ｔｅｍｐ」は、以下の値をとる。
【０２４８】
【数１２】
2007538241.tif 000015
【０２４９】
ｔは、第二ループにおいて進行中の持続時間であり、ＭＡＸ＿Ｃは指数関数の冪級数展開の切り捨て指数である。
【０２５０】
この変数「ｔｅｍｐ」は次に、エネルギーにしたがい、進行中の持続時間ｔに対する係数ｋt、eの集合を計算するために過程２１６で用いられる。
【０２５１】
前記係数ｋt、eの計算は、
ｋt、e＝（ｙt、e−ｔｅｍｐ）／λである。
【０２５２】
したがって、この係数は、式（１６）の前記和と、同じ式の和の前記部分の差分を出すことで計算される。
【０２５３】
これは、以下の式を計算することになる。
【０２５４】
【数１３】
2007538241.tif 000016
【０２５５】
いったんエネルギーにしたがい、所与の持続時間ｔに対する係数ｋt、eが決定されると、過程２１７で肯定性の制約テストを行うことで、得られた値が確率密度と一貫していることを保証するようになっている。
【０２５６】
テストは以下の不等式が確認されるかどうか判定することからなるのだが、つまり、
ｋt、e−２ｋt-1、e＋ｋt-2、e＞０？
および、ｔ＞１？
ということである。
【０２５７】
ある意味では、一つ前の持続時間ｔ−１においてディジタル二重微分が正かどうか判定される。
【０２５８】
否定的の場合、アルゴリズムは値０に対応する係数ｈt、eを固定する過程２１９に移行する。
【０２５９】
次に、過程２１８で、これらの決定された係数ｋt、eを変数ｇt（ｅ）に保存する。
【０２６０】
逆に、テスト２１７が肯定的であれば、過程２１９は避けられ、直接過程２１８に移行する。
【０２６１】
過程２２０では、持続時間に対する第二ループのテストによって、係数ｋt、e（またはｇt（ｅ））の持続時間に従った計算が終了しているかどうか判定される。
【０２６２】
否定的の場合、過程２２１に移行し、該過程において、持続時間はインクリメントされ、エネルギーの指数は再度値０に初期化される。
【０２６３】
次に、この新しい持続時間によって過程２１２の最初から再開する。
【０２６４】
逆に、計算が終わっていれば、持続時間に対する第二ループから脱ける。
【０２６５】
実際、エネルギーおよび持続時間に応じた係数ｋt、eはこのとき決定されている。
【０２６６】
上記で触れたように、ここで関係式（３）によって係数ｈt、eが演繹される。
【０２６７】
この二重の誘導微分の操作が過程２２２に実装されている。
【０２６８】
この過程の最後に、処理された信号に含まれる基本インパルスの持続時間／エネルギーのペアの確率密度を表す係数ｈt、eの集合が得られる。
【０２６９】
最後の過程２２３では、基本インパルスのエネルギーが判定され、エネルギーに従った変数ｍに該エネルギーが保存される。
【０２７０】
この目的を達成するために、係数ｈt、eは持続時間にしたがって加算され、各エネルギーそれぞれについても同様である。
【０２７１】
この過程が、基本インパルスの持続時間／エネルギーのペアにおける離散確率密度のエネルギーに従った前記周辺確率の計算に対応することに注目すべきである。
【０２７２】
メモリに保存された基本インパルスのエネルギーは最終的に視覚化手段を用いて利用者に与えられることで、該手段が存在する放射性核種を識別し、診断を確立するようになっている（プリンタ、ディスプレイまたは機能が同等のその他のあらゆる装置）。
【０２７３】
ここで、当業者であれば明らかに一定数の変形例を見つけることができるものである。
【０２７４】
たとえば、本発明に係る測定方法によって与えられる情報は中間過程に関するものとすることができる（メインインパルスまたは基本インパルスの持続時間−エネルギーのペアのヒストグラムのグラフ表示など）。
【０２７５】
では、上記で提示したアルゴリズムを用いた、本発明に係る測定方法を実施することで得られた非制限的ないくつかの結果を提示することにする。
【０２７６】
図１６は、セシウム１３７について観察されるガンマ線源のスペクトル４９を示している。
【０２７７】
スペクトロメータ装置はパイルアップが生じるものである。
【０２７８】
さらに、この装置は本発明によるパイルアップリジェクション方法をまだ利用していない。
【０２７９】
このスペクトルが、エネルギーＥ３での基本スペクトル線５０、ならびに、エネルギーＥ３の二つの光子のパイルアップに対応するエネルギーＥ４での第二のスペクトル線５１を含んでいることを確認することができる。
【０２８０】
図１７は、本発明による測定方法が前記装置で用いられたときの、同一のガンマ線源（セシウム１３７）のスペクトル５２を示している。
【０２８１】
エネルギーＥ４のスペクトル線５１が消えている。つまり、エネルギーＥ３に対応するスペクトル線だけがスペクトルに残っており、その結果、該スペクトルがセシウム１３７のスペクトルに一致することになる。
【０２８２】
さらに、図１７では、ノイズ５３の上限が図１６で観察されるものに比べて大きく減少していることが分かる。
【０２８３】
情報としては、この上限はとりわけパイルアップとコンプトン効果の組み合わせに起因するものである。
【０２８４】
したがって、これらの得られた結果は満足できるものである。
【０２８５】
これらの結果から、本発明による測定方法の効果、並びに、適切な電子装置において提案されたアルゴリズムを実装することによる産業上の実用性も明らかになった。
【０２８６】
より一般的には、本発明は、さらに複数の利点を提供するものであり、すなわち、
−放出されたエネルギーおよびそれらの発生をより良く識別すること、
−導入部で触れた除去によるフィルタリング方法とは異なり、信号の有効部分を除去することなく、この測定された信号全体を対象とすること、
−インパルスの形状についていかなる仮説も必要ではないこと、
−光子源の活性が弱かろうと高かろうと、該活性とは独立して方法を強く信頼できること（しかし方法の良好な作動に必要な条件は基本インパルスの恒常的なパイルアップがあってはならないということである）、
−測定される信号に直接本方法を適用する可能性、あるいは該方法を後で実行する可能性、
−非常に接近した利用であり、このことによって本発明に係る測定方法が先行技術の測定方法よりもディジタル的になっていることである。
【０２８７】
もちろん、本発明は上記に説明し、図面に示された実施態様にまったく限定されるものではない。
【０２８８】
特に、本発明はその他のスペクトロメトリーの領域（たとえば、アルファ・スペクトロメトリー、中性子スペクトロメトリー、そしてより一般的にはあらゆるタイプの微粒子に対するあらゆるスペクトロメトリー）において、さまざまなタイプの検出器、とりわけ異なったサイズおよび形状の検出器に対して実施することができる。
【０２８９】
さらに一般的には、本発明は、各信号それぞれが用語の一般的な意味で流束を表し得るものであるとき、メインの信号にパイルアップによって含まれている基本信号の貢献が求められる応用などに用いることができる。
【０２９０】
非限定的には、光子束、蛍光元素の流束または通信システムの待ち行列における要素の流束などであってよい。
【０２９１】
たとえば待ち行列の場合、本発明において説明した実施態様に等価なものとして、持続時間Ｄの変数に接続時間を関連づけ、エネルギーＥの変数に通過帯域を関連づけることができる。
【０２９２】
また、エネルギーＥがインパルス下におけるスペクトル面の算定によって評価されることはまったく義務的ではない。
【０２９３】
本発明で説明された実施態様の場合、スペクトル面による判定が完全に適合していたのだが、それは、単純な関係式がこのスペクトル面をガンマ光子のエネルギーに結びつけるからである。
【０２９４】
しかし、異なった応用の場合にはこれは必ずしも当てはまらない。
【０２９５】
ここで、当業者であれば理解できるが、本発明の概念が非常に一般的な意味で、メインインパルスにおける基本インパルスの重みを表示するためのエネルギーの利用に基づいていることなのである。
【０２９６】
したがって、この重みの算定方法は応用に応じて大きく異なることがあることは明白である。
【０２９７】
非制限的な例としては、重みは基本インパルスの振幅、特定の持続時間などに対応することもあるだろう。
【０２９８】
ここで、本発明に係る測定方法を適用するために守られるべき以下の制約だけ確認することとするが、該制約とはすなわち、
−基本インパルスの出現時間が均質なポワソン過程に従っていること、
−メインインパルスが所定の持続時間Ｄを有していること、
−基本インパルスのエネルギーを表す変数Ｅｉが加算性の特性を有していることである。
【図面の簡単な説明】
【０２９９】
【図１】ガンマスペクトロメトリーのディジタル取得の流れの概略図
【図２】図１の装置によって観察可能なガンマ線源のスペクトル
【図３】ガンマ光子の検出器によって生成された経時信号の例
【図４】セシウム１３７の単エネルギー光子源を測定したエネルギー分布
【図５】複数のエネルギーを備えた光子を放出する活性光子源のエネルギー分布
【図６】パイルアップ現象が起こったときの、遅速検出器に起因する時間信号の概略図
【図７】パイルアップ現象がないときの、ガンマ線源の正規化スペクトル
【図８】活性の増加によってパイルアップ現象が起こったときの、図７で分析されたガンマ線源の正規化スペクトル
【図９】本発明によるシステムによって測定された、二つのメインインパルスを含んだ信号の一部の概略図
【図１０】本発明によるシステムによってディジタル化されたときの、図９の信号を示す図
【図１１】本発明に係る測定方法のブロック略図を示す図
【図１２】本発明に係る測定方法によって処理された信号の活動シーケンスのＤとＥのペアの二次元ヒストグラムの例
【図１３】前記処理された信号の基本インパルスのＤｉとＥｉのペアの二次元ヒストグラムの例
【図１４】適切な電子装置において本発明によるパイルアップリジェクション方法の第一部分を実装することを可能にするアルゴリズムのフローチャート
【図１５】適切な電子装置において本発明によるパイルアップリジェクション方法の第二部分を実装することを可能にするアルゴリズムのフローチャート
【図１６】本発明に係る測定方法を実施していないときの、セシウム１３７の観察されたガンマ線エネルギーのスペクトル
【図１７】本発明に係る測定方法を実施しているときの、セシウム１３７の観察されたガンマ線エネルギーのスペクトル
【図１】
2007538241.tif 000018
【図２】
2007538241.tif 000019
【図３】
2007538241.tif 000020
【図４】
2007538241.tif 000021
【図５】
2007538241.tif 000022
【図６】
2007538241.tif 000023
【図７】
2007538241.tif 000024
【図８】
2007538241.tif 000025
【図９】
2007538241.tif 000026
【図１０】
2007538241.tif 000027
【図１１】
2007538241.tif 000028
【図１２】
2007538241.tif 000029
【図１３】
2007538241.tif 000030
【図１４】
2007538241.tif 000031
【図１５】
2007538241.tif 000032
【図１６】
2007538241.tif 000033
【図１７】
2007538241.tif 000034
【国際調査報告】
2007538241.tif 000035
2007538241.tif 000036
2007538241.tif 000037
2007538241.tif 000038
2007538241.tif 000039
2007538241.tif 000040