(57)【特許請求の範囲】
【請求項１】
臨床ポジトロン放出断層撮影（positron emission tomography（ＰＥＴ））デバイスに含まれるＬ個の検出器の信号をＭ個のグループにグループ分けするステップであって、各グループの検出器の数は、少なくとも２個、多くともＮ個であり、ａ番目のグループ内の検出器の信号数をＰ（ａ）とし、
【数１】
0006679715.tif 000005
であり、ａ番目のグループのｂ番目の検出器の信号をＳｉｇｎａｌ（ａ，ｂ）と示され、Ｍ≧２、Ｎ≧２、１≦ａ≦Ｍ、１≦ｂ≦ＮであるステップＳ１と、
ステップＳ１におけるＬ個の検出器の信号を第１のソース信号及び第２のソース信号に分割するステップＳ２と、
２つの読出チャンネルＡ及びＢを含む第１の信号伝送ラインを提供し、前記検出器の信号を前記第１の信号伝送ラインに入力し、前記第１の信号伝送ライン上の前記検出器の信号の間に少なくとも１つの第１の信号遅延ユニットを設けるステップＳ３と、
２つの読出チャンネルＣ及びＤを含む第２の信号伝送ラインを提供し、前記検出器の信号を前記第２の信号伝送ラインに入力し、前記第２の信号伝送ライン上の前記検出器の信号の間に少なくとも１つの第２の信号遅延ユニットを設けるステップＳ４と、
前記４つの読出チャンネルＡ、Ｂ、Ｃ及びＤからのパルスに基づいて、
前記読出チャンネルＡ及びＢにおける信号の到着時間の時間差を算出し、前記読出チャンネルＣ及びＤにおける信号の到着時間の時間差を算出し、信号ソースの行番号及び列番号を判定し、前記時間差に基づいて前記信号を生成するソース検出器を特定し、最終的なパルス情報を取得するステップＳ５と
を有する検出器信号読出チャンネル多重化方法。
【請求項２】
ステップＳ１において、前記Ｍ個のグループは、それぞれＮ個である同数の検出器を含み、Ｍ×Ｎ＝Ｌ、Ｍ≧２、Ｎ≧２である請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項３】
前記第１の信号伝送ライン及び前記第２の信号伝送ラインは、互いに直交しており、前記伝送ライン上の位置情報は、前記遅延ユニットによって分割される請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項４】
ステップＳ３及びステップＳ４において、前記第１の信号伝送ライン及び前記第２の信号伝送ラインには、空間位置が対応する２つの検出器グループがそれぞれ接続される請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項５】
ステップＳ３及びステップＳ４において、前記第１の信号伝送ライン及び前記第２の信号伝送ラインには、前記各検出器の２つの読出チャンネルがそれぞれ接続される請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項６】
ステップＳ３において、前記ａ番目のグループの検出器の第１のソース信号を第１の信号として機能し、前記Ｍ個のグループにおいてＭ個の第１の信号が生成され、前記Ｍ個の第１の信号が第１の信号伝送ラインに入力され、前記少なくとも１つの第１の信号遅延ユニットは、前記Ｍ個の第１の信号のうちの隣接する２つの第１の信号の間に設けられる請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項７】
前記第１の信号遅延ユニットが複数存在し、前記第１の信号遅延ユニットの少なくとも１つは、前記Ｍ個の第１の信号のうちの任意の２つの隣接する第１の信号の間に設けられる請求項６記載の検出器信号読出チャンネル多重化方法。
【請求項８】
ステップＳ４において、Ｍ個のグループのそれぞれの第ｂ番目の検出器の第２のソース信号を第２の信号として機能し、Ｎ個の第２の信号が生成され、前記Ｎ個の第２の信号が第２の信号伝送ラインに入力され、前記少なくとも１つの遅延ユニットは、前記Ｎ個の第２の信号の２つの隣接する第２の信号の間に設けられる請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項９】
前記第２の信号遅延ユニットが複数存在し、前記第２の信号遅延ユニットの少なくとも１つは、前記Ｎ個の第２の信号のうちの任意の２つの隣接する第２の信号の間に設けられる請求項８記載の検出器信号読出チャンネル多重化方法。
【請求項１０】
ステップＳ５は、減算回路によって、前記読出チャンネルＡ及びＢのパルス間の差分パルス、並びに前記読出チャンネルＣ及びＤのパルス間の差分パルスを取得することと、前記差分パルスのパルス幅によって、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表すこととを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項１１】
ステップＳ５において、減算回路によって、前記読出チャンネルＡ及びＢのパルス間の差分パルス、並びに前記読出チャンネルＣ及びＤのパルス間の差分パルスを取得することと、前記差分パルスの振幅によって、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表すこととを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項１２】
ステップＳ５の前に、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出することと、前記第１の信号伝送ライン、前記第２の信号伝送ライン、及び前記第１及び第２信号伝送ラインに与えられる前記第１及び第２の信号遅延ユニットの配置を含む回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出することと、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別することとを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項１３】
ステップＳ５の前に、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出することと、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出することと、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別することとを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項１４】
ステップＳ５の前に、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出することと、前記第１の信号伝送ライン、前記第２の信号伝送ライン、及び前記第１及び第２信号伝送ラインに与えられる前記第１及び第２の信号遅延ユニットの配置を含む回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出することと、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、前記各位置の周囲の一定の領域を設定することとを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【請求項１５】
ステップＳ５の前に、前記読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、前記読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出することと、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出することと、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、前記各位置の周囲の一定の領域を設定することとを含む請求項１記載の検出器信号読出チャンネル多重化方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本出願は、２０１５年６月２９日に中国国家知的所有権局に提出された「検出器信号読出チャンネル多重化方法」と題された中国特許出願第２０１５１０３６７１０５．５号に対する優先権を主張し、この文献は、参照により全体が本明細書に援用される。
【０００２】
本開示は、核検出及び核医学イメージングの技術分野に関し、特に、検出器信号読出チャンネル多重化方法に関する。
【背景技術】
【０００３】
臨床ポジトロン放出断層撮影（positron emission tomography：以下、ＰＥＴ）デバイスは、３０，０００個の検出器を含み、したがって、各検出器によって生成される電気信号を直接読み出すためには３０，０００個の電子処理チャンネルが必要である。多数の電子チャンネルを設けることは高コストであり、ＰＥＴ電子システム全体の工程の実現が困難になる。
【０００４】
このために、本分野では、検出器信号読出チャンネル多重化技術の研究開発に注目が集まっている。既存の主流の読出チャンネル多重化方法では、抵抗ネットワークを利用して、行及び列に基づいて出力された検出器信号を重み付けし、次いで出力信号を読み出す。
【０００５】
抵抗ネットワークによってチャンネルを多重化することにより、検出器の読出チャンネル数が減少し、チャンネル数が過剰になる問題が解決される。しかしながら、抵抗ネットワークでは、信号入力点が異なるために、抵抗ネットワークの同等の抵抗が検出器毎に異なり、異なる検出器の出力信号の振幅の間に大きな差が生じる。この場合、バックエンド読出回路は、大きなダイナミックレンジを必要とする。既存の処理回路は、通常、ダイナミックレンジが制限されており、すなわち、入力信号が小さいと信号対雑音比が低くなり、入力信号が大きいと飽和問題が生じる。
【発明の概要】
【発明が解決しようとする課題】
【０００６】
したがって、上述の技術的課題に鑑み、上述の欠点を克服し、チャンネル多重化された検出器によって出力される信号のダイナミックレンジ要求が過大になる問題を効果的に解決するように改良された、検出器信号を読み出すための構造を提供するチャンネル多重化方法を提供する必要がある。
【０００７】
そこで、本発明の目的は、使用される電子チャンネルの数を減らし、バックエンド読出回路のダイナミックレンジに対する要求を効果的に低減し、及びバックエンド読出回路の時間分解能に対する要求を低減する検出器信号読出チャンネル多重化方法を提供することである。
【課題を解決するための手段】
【０００８】
上述の目的を達成するために、本開示は、以下の技術的ソリューションを提供する。
【０００９】
すなわち、検出器信号読出チャンネル多重化方法を提供し、この方法は、
ステップＳ１：Ｌ個の検出器の信号をＭ個のグループにグループ分けするステップであって、各グループの検出器の数は、少なくとも２個、多くともＮ個であり、ａ番目のグループ内の検出器の信号数をＰ（ａ）とし、
【数１】
0006679715.tif 000002
であり、ａ番目のグループのｂ番目の検出器の信号をＳｉｇｎａｌ（ａ，ｂ）と示し、Ｍ≧２、Ｎ≧２、１≦ａ≦Ｍ、１≦ｂ≦Ｎであるステップと、
ステップＳ２：ステップＳ１におけるＬ個の検出器の信号を第１のソース信号及び第２のソース信号に分割するステップと、
ステップＳ３：２つの読出チャンネルＡ及びＢを含む第１の信号伝送ラインを提供し、検出器の信号を第１の信号伝送ラインに入力し、第１の信号伝送ライン上の検出器の信号の間に少なくとも１つの第１の信号遅延ユニットを設けるステップと、
ステップＳ４：２つの読出チャンネルＣ及びＤを含む第２の信号伝送ラインを提供し、検出器の信号を第２の信号伝送ラインに入力し、第２の信号伝送ライン上の検出器の信号の間に少なくとも１つの第２の信号遅延ユニットを設けるステップと、
ステップＳ５：４つの読出チャンネルＡ、Ｂ、Ｃ及びＤからのパルスに基づいて、信号を生成するソース検出器をシンボル化し、最終的なパルス情報を取得するステップとを有する。
【００１０】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ１において、Ｍ個のグループは、それぞれＮ個である同数の検出器を含み、Ｍ×Ｎ＝Ｌ、Ｍ≧２、Ｎ≧２である。
【００１１】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、第１の信号伝送ライン及び第２の信号伝送ラインは、互いに直交しており、伝送ライン上の位置情報は、遅延ユニットによって分割される。
【００１２】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ３及びステップＳ４において、第１の信号伝送ライン及び第２の信号伝送ラインには、空間位置が対応する２つの検出器グループがそれぞれ接続される。
【００１３】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ３及びステップＳ４において、第１の信号伝送ライン及び第２の信号伝送ラインには、各検出器の２つの読出チャンネルがそれぞれ接続される。
【００１４】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ３において、ａ番目のグループの検出器の第１のソース信号を第１の信号として機能し、Ｍ個のグループにおいてＭ個の第１の信号が生成され、Ｍ個の第１の信号が第１の信号伝送ラインに入力され、少なくとも１つの第１の信号遅延ユニットは、Ｍ個の第１の信号のうちの隣接する２つの第１の信号の間に設けられる。
【００１５】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、第１の信号遅延ユニットが複数存在し、第１の信号遅延ユニットの少なくとも１つは、Ｍ個の第１の信号のうちの任意の２つの隣接する第１の信号の間に設けられる。
【００１６】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ４において、Ｍ個のグループのそれぞれの第ｂ番目の検出器の第２のソース信号を第２の信号として機能し、Ｎ個の第２の信号が生成され、Ｎ個の第２の信号が第２の信号伝送ラインに入力され、少なくとも１つの遅延ユニットは、Ｎ個の第２の信号の２つの隣接する第２の信号の間に設けられる。
【００１７】
上述した検出器信号読出チャンネル多重化方法では、第２の信号遅延ユニットが複数存在し、第２の信号遅延ユニットの少なくとも１つは、Ｎ個の第２の信号のうちの任意の２つの隣接する第２の信号の間に設けられる。
【００１８】
上述した検出器信号読出チャンネル多重化方法では、ステップ５において、４つの読出チャンネルＡ、Ｂ、Ｃ、Ｄにおける各信号の到着時間を測定し、読出チャンネルＡ及びＢにおける信号の到着時間差と、読出チャンネルＣ及びＤにおける信号の到着時間差とを算出し、時間差に基づいて信号ソースの行番号及び列番号を判定し、２つの時間差を使用して信号を生成するソース検出器をシンボル化する。
【００１９】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップＳ５において、減算回路によって、読出チャンネルＡ及びＢのパルス間の差分パルス、並びに読出チャンネルＣ及びＤのパルス間の差分パルスを取得し、この差分パルスのパルス幅によって、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表す。
【００２０】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップ５において、減算回路によって、読出チャンネルＡ及びＢのパルス間の差分パルス、並びに読出チャンネルＣ及びＤのパルス間の差分パルスを取得し、差分パルスの振幅によって、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表す。
【００２１】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。
【００２２】
上述の検出器信号読出チャンネル多重化方法において、好ましくは、ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。
【００２３】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、各位置の周囲の一定の領域を設定する。
【００２４】
上述した検出器信号読出チャンネル多重化方法では、好ましくは、ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、各位置の周囲の一定の領域を設定する。
【００２５】
上述の技術的ソリューションは、従来技術に比べて以下の利点を有する。
【００２６】
１．検出器信号を読み取るためのチャンネル多重化方法は、
【００２７】
ステップＳ１：Ｌ個の検出器の信号をＭ個のグループにグループ分けするステップであって、各グループの検出器の数は、少なくとも２個、多くともＮ個であり、ａ番目のグループ内の検出器の信号数をＰ（ａ）とし、
【数２】
0006679715.tif 000003
であり、ａ番目のグループのｂ番目の検出器の信号をＳｉｇｎａｌ（ａ，ｂ）とシンボル化し、Ｍ≧２、Ｎ≧２、１≦ａ≦Ｍ、１≦ｂ≦Ｎであるステップと、ステップＳ２：ステップＳ１におけるＬ個の検出器の信号を第１のソース信号及び第２のソース信号に分割するステップと、ステップＳ３：２つの読出チャンネルＡ及びＢを含む第１の信号伝送ラインを提供し、検出器の信号を第１の信号伝送ラインに入力し、第１の信号伝送ライン上の検出器の信号の間に少なくとも１つの第１の信号遅延ユニットを設けるステップと、ステップＳ４：２つの読出チャンネルＣ及びＤを含む第２の信号伝送ラインを提供し、検出器の信号を第２の信号伝送ラインに入力し、第２の信号伝送ライン上の検出器の信号の間に少なくとも１つの第２の信号遅延ユニットを設けるステップと、ステップＳ５：４つの読出チャンネルＡ、Ｂ、Ｃ及びＤからのパルスに基づいて、信号を生成するソース検出器をシンボル化し、最終的なパルス情報を取得するステップとを有する。本開示では、使用される電子チャンネルの数が少なくなり、ＰＥＴシステム全体の工程の実現のコスト及び困難性が低減される。また、伝送ラインの時分割多重設計により、異なる位置の検出器によって出力される信号の波形と、読み出される波形とが１対１の形式で類似することが保証され、これによって、バックエンド読出回路のダイナミックレンジに対する要求が緩和される。更に、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００２８】
２．ステップＳ１において、Ｍ個のグループは、それぞれＮ個である同数の検出器を含み、Ｍ×Ｎ＝Ｌ、Ｍ≧２、Ｎ≧２である。この技術的ソリューションにおいて、各行の検出器の数は、実際の状況に応じて、同じであってもよく、異なっていてもよく、このことは、本方法に影響を与えない。
【００２９】
３．第１の信号伝送ライン及び第２の信号伝送ラインは、互いに直交しており、伝送ライン上の位置情報は、遅延ユニットによって分割される。この技術的ソリューションにおいて、信号は、直交する２本の伝送ラインに入力され、信号を出力する検出器の位置は、２本の伝送ライン上で判定された位置範囲に基づいて判定される。平面アレイ検出器の出力は、４つの読出チャンネルを介して取得される。
【００３０】
４．第１の信号伝送ライン及び第２の信号伝送ラインには、空間位置が対応する２つの検出器グループがそれぞれ接続され、又は第１の信号伝送ライン及び第２の信号伝送ラインには、各検出器の２つの読出チャンネルがそれぞれ接続される。この技術的ソリューションは、検出器を伝送線路に接続する２つの手法を提供し、検出器の実際の状況に基づいて適切な手法を選択できる。例えば、検出器が２つの読出チャンネルを含む場合、２つの読出チャンネルを伝送ラインに接続する手法を選択できる。
【００３１】
５．ステップＳ３では、ａ番目のグループの検出器の第１のソース信号を第１の信号として機能し、Ｍ個のグループからＭ個の第１の信号を生成し、Ｍ個の第１の信号を第１の信号伝送ラインに入力し、Ｍ個の第１の信号のうちの隣接する２つの第１の信号の間に、少なくとも１つの第１の信号遅延ユニットを設ける。この技術的ソリューションにより、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００３２】
６．第１の信号遅延ユニットが複数存在し、第１の信号遅延ユニットの少なくとも１つは、Ｍ個の第１の信号のうちの任意の２つの隣接する第１の信号の間に設けられる。この技術的ソリューションでは、各第１の信号に対して遅延ユニットを追加し、これにより、時間差の計算が容易になる。
【００３３】
７．ステップＳ４において、Ｍ個のグループのそれぞれの第ｂ番目の検出器の第２のソース信号を第２の信号として機能し、Ｎ個の第２の信号が生成され、Ｎ個の第２の信号が第２の信号伝送ラインに入力され、少なくとも１つの遅延ユニットは、Ｎ個の第２の信号の２つの隣接する第２の信号の間に設けられる。この技術的ソリューションにより、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００３４】
８．第２の信号遅延ユニットが複数存在し、第２の信号遅延ユニットの少なくとも１つは、Ｎ個の第２の信号のうちの任意の２つの隣接する第２の信号の間に設けられる。この技術的ソリューションでは、各第１の信号に対して遅延ユニットを追加し、これにより、時間差の計算が容易になる。
【００３５】
９．ステップ５において、４つの読出チャンネルＡ、Ｂ、Ｃ、Ｄにおける各信号の到着時間を測定し、読出チャンネルＡ及びＢにおける信号の到着時間差と、読出チャンネルＣ及びＤにおける信号の到着時間差とを算出し、時間差に基づいて信号ソースの行番号及び列番号を判定し、２つの時間差を使用して信号を生成するソース検出器をシンボル化する。この技術的ソリューションでは、時間差を計算する手法が複数ある。４つのパルスに関する情報が最初にデジタル化される場合、４つの読出チャンネルにおける到着時間に基づいて時間差を計算できる。パルスの情報を可能な限り早くデジタル化することによって後続するアナログ回路からの干渉を低減できる。
【００３６】
１０．ステップＳ５において、減算回路によって、読出チャンネルＡ及びＢのパルス間の差分パルス、並びに読出チャンネルＣ及びＤのパルス間の差分パルスを取得し、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表す差分パルスのパルス幅を取得する。この技術的ソリューションでは、時間差を表すために差分パルスのパルス幅を使用できる。ある位相差でパルスが出力されると、そのパルスの差分パルスの幅を用いて位相差を表すことができる。この方法により、バックエンド読出回路のチャンネル数を減らすことができる。
【００３７】
１１．ステップ５において、減算回路によって、読出チャンネルＡ及びＢのパルス間の差分パルス、並びに読出チャンネルＣ及びＤのパルス間の差分パルスを取得し、取得した差分パルスの振幅によって、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差、並びに読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差を表す。この技術的ソリューションでは、パルスの位相及び立ち上がり時間が適切である場合、差分パルスの振幅を用いて時間差を表すことができる。到着時間を測定する代わりに、ピーク検出を採用することによって、回路の複雑さを低減することができる。
【００３８】
１２．ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。この技術的ソリューションでは、回路構造のパラメータに基づいて確率分布関数上の検出器に対応する位置を計算することにより、試験プロセスにおける困難性が低減され、量産時の生産性が向上する。
【００３９】
１３．ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出して、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。この技術的ソリューションでは、異なる検出器の位置の間の分離線を用いて確率分布関数上の検出器に対応する分布領域を区別するため、システムの検出効率が向上する。各イベントは、位置によって判定される。
【００４０】
１４．ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、各位置の周囲の一定の領域を設定する。この技術的ソリューションでは、確率分布関数上の位置に対応する検出器の分布領域として、ある位置の周囲の一定の領域を設定することにより、誤認識されやすいイベントが排除され、信号対雑音比が改善される。
【００４１】
１５．ステップ５において、読出チャンネルＡ及びＢにおける信号の到着時刻間の時間差をｘとし、読出チャンネルＣ及びＤにおける信号の到着時刻間の時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいて、ｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の位置に対応する検出器の分布領域として、各位置の周囲の一定の領域を設定する。この技術的ソリューションでは、既存のデータに基づいて確率分布関数を取得する。この方法では、検出器アレイ内の様々な分布パラメータを反映させることができ、したがって、より正確な位置を得ることができ、これは検出器アレイにより適している。
【００４２】
以下、本開示の実施形態又は従来技術に基づく技術的ソリューションをより明確に示すために、本開示の実施形態又は従来技術の説明において必要とされる図面について簡単に説明する。明らかに、以下に説明する図面は、本開示の実施形態の一部のみを示すものであり、当業者であれば、これらの図面に基づいて、創造的な努力なしに他の図面を想到できることは明らかである。
【図面の簡単な説明】
【００４３】
【図１】両端に読出回路が設けられた伝送ラインへのパルス信号の入力を示す概略図である。
【図２】図１の伝送ライン上に遅延ユニットを追加した概略図である。
【図３】両端に読出回路を設けた伝送ラインにパルス信号を入力し、伝送ライン上に複数の遅延器を設けた概略図である。
【図４】検出器の積層構造を示す概略図である。
【図５】４×４の平面アレイを有する結晶アレイを示す概略図である。
【図６】本開示に基づく検出器信号読出チャンネル多重化方法における検出器回路基板上の第１の信号伝送ラインの概略構成図である。
【図７】本開示に基づく検出器信号読出チャンネル多重化方法における検出器回路基板上の第２の信号伝送ラインの概略構成図である。
【図８】本開示に基づく検出器信号読出チャンネル多重化方法において、第１の信号伝送ラインが第２の信号伝送ラインに直交することを示す概略図である。
【発明を実施するための形態】
【００４４】
本開示は、使用される電子チャンネルの数を減らし、バックエンド読出回路のダイナミックレンジに対する要求を効果的に低減し、及びバックエンド読出回路の時間分解能に対する要求を低減する検出器信号読出チャンネル多重化方法を提供する。
【００４５】
図１に示すように、両端に読出回路を有する伝送ラインにパルス信号を入力する場合、パルス信号は、入力点から両端に向かって伝送される。図２では、伝送ライン上に１つの遅延ユニットを追加している。伝送ライン上の信号の伝搬速度は、遅延器を通過する信号の速度よりも遥かに高いので、２つの端部の読出チャンネルにおける信号の到着時間の時間差は、伝送ライン上の遅延ユニットに依存する。更に、図３では、伝送ラインに複数の遅延ユニットを追加している。遅延ユニットの時間遅延は、入力点間で信号を送信する時間よりも遥かに大きいので、伝送ラインに入力される信号は、両端の各信号の到着時間の時間差に基づいて複数のグループにグループ分けされ、信号の時間差を測定することによって、すなわち、伝送ラインに信号が入力された位置区間（隣接する２つの遅延器の間の位置区間）を判定できる。図３において、参照符号１００は、クリスタルアレイを示し、参照符号２００は、シリコン光電子倍増管（ｓｉｌｉｃｏｎ  ｐｈｏｔｏｍｕｌｔｉｐｌｉｅｒ：ＳｉＰＭ）検出回路基板、参照符号３００は、ＳＩＰＭ検出回路基板２を示している。
【００４６】
上述のプロトタイプは、ＰＥＴ検出器に適用される。検出器には多くの種類があり、例えば、シリコン光電子増倍管（ＳｉＰＭ）検出器を用いることができる。ＳｉＰＭ検出器は、シリコン光電子増倍管（ＳｉＰＭ）回路基板と、ＳｉＰＭ回路基板上に結合されたクリスタルアレイとを含む。検出器の積層構造は、図４に示すとおりである。
【００４７】
図５は、４×４の平面アレイを有する結晶アレイを示す概略図である。本開示による技術的ソリューションは、結晶アレイを例として説明される。
【００４８】
図６〜図８において、参照符号４００は、読出チャンネル、参照符号５００は、伝送ライン、参照符号６００は、遅延ユニット、参照符号７００は、ＳＩＰＭを示している。本開示による検出器信号読出チャンネル多重化方法は、ステップＳ１からステップＳ５を含む。
【００４９】
ステップＳ１では、Ｌ個の検出器の信号をＭ個のグループにグループ分けする。各グループにおける検出器の数は、少なくとも２個、多くともＮ個であり、ａ番目のグループの検出器の信号数をＰ（ａ）とし、
【数３】
0006679715.tif 000004
であり、ａ番目のグループのｂ番目の検出器の信号をＳｉｇｎａｌ（ａ，ｂ）と示し、ここで、Ｍ≧２、Ｎ≧２、１≦ａ≦Ｍ、１≦ｂ≦Ｎである。
【００５０】
ステップＳ２では、ステップＳ１におけるＬ個の検出器の信号を第１のソース信号及び第２のソース信号に分割する。
【００５１】
ステップＳ３では、２つの読出チャンネルＡ、Ｂを含む第１の信号伝送ラインを提供し、検出器の信号を第１の信号伝送ラインに入力し、第１の信号伝送ライン上の各検出器の信号の間に少なくとも１つの第１の信号遅延ユニットを設ける。
【００５２】
ステップＳ４では、２つの読出チャンネルＣ及びＤを含む第２の信号伝送ラインを提供し、検出器の信号を第２の信号伝送ラインに入力し、第２の信号伝送ライン上の検出器の信号の間に少なくとも１つの第２の信号遅延ユニットを設ける。
【００５３】
ステップ５では、信号を生成するソース検出器をシンボル化し、４つの読出チャンネルＡ、Ｂ、Ｃ及びＤからのパルスに基づいて最終的なパルス情報を取得する。
【００５４】
本開示では、使用される電子チャンネルの数が低減され、これにより、ＰＥＴシステム全体の工程実施のコスト及び困難性が低減される。また、伝送ラインの時分割多重設計により、異なる位置の検出器によって出力される信号の波形と、読み出される波形とが１対１の形式で類似することが保証され、これによって、バックエンド読出回路のダイナミックレンジに対する要求が緩和される。更に、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００５５】
ステップＳ１において、Ｍ個のグループに含まれる検出器の数は、同じであってもよく、異なっていてもよい。Ｍ個のグループが同数の検出器を含む場合、この同数をＮと表すことができ、Ｍ×Ｎ＝Ｌ、Ｍ≧２、Ｎ≧２である。各行の検出器の数は、実際の状況に依存し、これらが同じであるか異なるかは、本方法に影響を与えない。
【００５６】
第１の信号伝送ライン及び第２の信号伝送ラインは、互いに直交しており、伝送ライン上の位置情報は、遅延ユニットによって分割される。ステップＳ３及びステップＳ４では、第１の信号伝送ライン及び第２の信号伝送ラインには、空間位置が対応する２つの検出器グループがそれぞれ接続される。あるいは、各検出器の２つの読出チャンネルは、第１の信号伝送ライン及び第２の信号伝送ラインにそれぞれ接続される。このように、直交する２本の伝送ライン（２本の伝送ラインが直交する場合、それぞれを行伝送ライン及び列伝送ラインと呼ぶことがある）を設け、伝送ラインを分割する遅延器に基づいて、伝送ライン上の位置（検出器グループ）を決定する。空間位置が対応する２つの検出器グループは、それぞれ２つの伝送ラインに接続してもよく、各検出器の２つの読出チャンネルを２つの伝送ラインにそれぞれ接続してもよい（この場合、２つの読出信号は、空間位置が対応する信号とみなすこともできる）。この技術的ソリューションでは、各信号が２つの直交する伝送ラインに入力され、信号を出力する検出器の位置は、２つの伝送ライン上の決定された位置範囲に基づいて決定され、４つの読出チャンネルを介して平面アレイ検出器の出力が取得される。
【００５７】
検出器は、ＳｉＰＭ検出器、フォトダイオード検出器、又は１つの出力チャンネルのみを有する他の検出器であってもよい。光電子増倍管（ＰＭＴ）のような１つの出力チャンネルのみを有する検出器には、２つのプロセスを適用できる。プロセスａでは、空間位置が対応する２つの検出器グループを２つの伝送ラインにそれぞれ接続する。プロセスｂでは、アナログ回路を介して各出力信号を２つの信号に分割し、２つの信号を２本の伝送ラインにそれぞれ入力する。シリコン光電子増倍管（略称ＳｉＰＭ）検出器のような２つの出力チャンネルを有する検出器の場合、各信号は、アノード及びカソードの２つのチャンネルから同時に読み出すことができる。ＳｉＰＭを行及び列に配置すると、各信号をアノード及びカソードから出力して２つの信号を生成できる。すなわち、第１のソース信号及び第２のソース信号は、同じ検出器から出力してもよく、空間位置が対応する２つの検出器から出力してもよい。
【００５８】
ステップＳ３では、ａ番目のグループの検出器の第１のソース信号を第１の信号として機能し、Ｍ個のグループからＭ個の第１の信号を生成する。Ｍ個の第１の信号は、第１の信号伝送ラインに入力される。Ｍ個の第１の信号のうちの隣接する２つの第１の信号の間に、少なくとも１つの第１の信号遅延ユニットを設ける。この技術的ソリューションにより、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００５９】
ステップＳ４では、Ｍ個のグループのそれぞれのｂ番目の検出器の第２のソース信号を第２の信号として機能し、Ｎ個の第２の信号を生成する。Ｎ個の第２の信号は、第２の信号伝送ラインに入力される。Ｎ個の第２の信号の２つの隣接する第２の信号の間に、少なくとも１つの遅延ユニットを設ける。この技術的ソリューションでは、遅延ユニットを追加することによって、バックエンド読出回路の時間分解能に対する要求が大幅に緩和される。
【００６０】
複数の第１の信号遅延ユニットと、複数の第２の信号遅延ユニットとを設けてもよく、Ｍ個の第１の信号のうちの任意の２つの隣接する第１の信号の間に少なくとも１つの第１の信号遅延ユニットを設け、Ｎ個の第２の信号のうちの任意の２つの隣接する第２の信号の間に少なくとも１つの第２の信号遅延ユニットを設ける。この技術的ソリューションでは、各第１の信号に対して遅延ユニットを追加し、これにより、時間差の計算が容易になる。
【００６１】
図４に示すように、γ光子がシンチレーション結晶を介して可視光光子に変換されるとき、２つの端部の検出器によって可視光光子が検出される。両端の検出回路に接続された伝送ラインは、互いに直交しており、２つの伝送ラインの位置情報を介してγ光子を変換する結晶ストリップを決定でき、すなわち、γ光子の入射位置が検出される。更に、２つの検出器によって検出された波形を比較することによって、エネルギ位置深度（ＤＯＩ）に関する情報を得ることもできる。
【００６２】
ステップＳ５は、４つの読出チャンネルＡ、Ｂ、Ｃ、Ｄにおける信号の到着時間を測定し、読出チャンネルＡ及びＢにおける信号の到着時間の時間差を算出し、読出チャンネルＣ及びＤにおける信号の到着時間の時間差を算出し、時間差に基づいて信号ソースの行番号及び列番号を判定し、２つの時間差を使用して信号を生成するソース検出器をシンボル化することによって実現される。この技術的ソリューションでは、時間差は、複数の手法で算出できる。４つのパルス情報が最初にデジタル化される場合、４つの読出チャンネルにおける到着時間に基づいて時間差を算出できる。可能な限り早くパルスをデジタル化することによって、後続するアナログ回路からの干渉を低減できる。
【００６３】
時間差を取得する手法は、以下に限定されるわけではないが、以下の２つを含む。
【００６４】
第１の手法では、減算回路が、読出チャンネルＡ及びＢのパルス間の差分パルスと、読出チャンネルＣ及びＤのパルス間の差分パルスとを取得し、差分パルスのパルス幅を用いて、読出チャンネルＡ及びＢの間の信号の到達時間の時間差と、読出チャンネルＣ及びＤの間の信号の到達時間の時間差とを表現する。差分パルスのパルス幅は、時間差を表すために使用される。同じパルスが一定の位相差で出力される場合、パルスの差分パルスの幅を用いて位相差を表すことができる。バックエンド読出回路のチャンネル数は、第１の手法によって削減できる。
【００６５】
第２の手法では、ステップＳ５において、減算回路が、読出チャンネルＡ及びＢのパルス間の差分パルスと、読出チャンネルＣ及びＤのパルス間の差分パルスとを取得し、差分パルスの振幅を用いて、読出チャンネルＡ及びＢの間の信号の到達時間の時間差と、読出チャンネルＣ及びＤの間の信号の到達時間の時間差とを表現する。パルスの位相及び立ち上がりエッジ時間が適切である場合、差分パルスの振幅を用いて時間差を表すことができる。到着時間を測定する代わりに、ピーク検出を採用して回路の複雑さを低減できる。
【００６６】
時間差に基づいて信号ソースの行番号及び列番号を決定し、２つの時間差に基づいて信号を生成するソース検出器をシンボル化することは、以下のような手法で実施される。
【００６７】
第１の手法では、ステップＳ５において、読出チャンネルＡ及び読出チャンネルＢの時間差をｘとし、読出チャンネルＣ及び読出チャンネルＤの時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいてｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出し、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。この確率分布関数は、後に、検出器の位置を認識するための基準となる。確率分布関数に従って検出器の位置を算出でき、検出器に対応する位置は、何らかの計算法によって確率分布関数上で「認識」される（例えば、白い斑点、白い斑点の周囲の領域、例えば、四角形内の領域を検出器に対応する位置として設定する）。後にパルスが検出されると、得られた（ｘ，ｙ）の値を確率密度グラフと照合し、このパルスのソースとなる検出器が判定又は区別される。異なる検出器の位置間の分離線は、確率分布関数上の検出器に対応する分布領域を区別するために使用され、これによりシステムの検出効率が改善される。各イベントは、位置によって決定される。
【００６８】
第２の手法では、読出チャンネルＡ及び読出チャンネルＢの時間差をｘとし、読出チャンネルＣ及び読出チャンネルＤの時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいてｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、ｇ（ｘ，ｙ）上の隣接する検出器間の分離線を算出し、ｇ（ｘ，ｙ）上の異なる検出器に対応する分布領域を区別する。検出器アレイ内の特定の検出器がパルス信号を出力すると、この信号は、４つの読出チャンネルで検出され、２つの時間差ｘ、ｙが算出される。確率密度関数を使用して、検出器の位置を認識できる。アレイ内の特定の検出器に対応する値（ｘ，ｙ）は、回路構造に基づいて取得することもできる。伝送ライン内のパルスの伝播速度は、演算によって求めてもよく、一定としてもよい。したがって、伝送ラインにおける検出器の信号の２つの読出チャンネルへの伝播距離に基づいて、その検出器に対応する時間差（ｘ，ｙ）を算出できる。この時間差（ｘ，ｙ）は、後のデータ収集におけるパルスに対応する検出器の位置を取得するための基準とすることもできる。回路構成に基づいて検出器の位置情報の較正を行うことができる。確率密度関数上の検出器に対応する位置を回路構造のパラメータに基づいて計算することにより、試験プロセスの困難性が低減され、量産時の生産性が向上する。
【００６９】
第３の手法では、ステップＳ５において、読出チャンネルＡ及び読出チャンネルＢの時間差をｘとし、読出チャンネルＣ及び読出チャンネルＤの時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、回路構造に基づいてｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、その位置の周囲の一定の領域をｇ（ｘ，ｙ）上の対応する検出器の分布領域とする。その位置の周囲の一定の領域を確率分布関数上の対応する検出器の分布領域として設定することにより、誤認識されやすいイベントが排除され、信号対雑音比が改善される。
【００７０】
第４の手法では、ステップＳ５において、読出チャンネルＡ及び読出チャンネルＢの時間差をｘとし、読出チャンネルＣ及び読出チャンネルＤの時間差をｙとし、確率分布関数ｇ（ｘ，ｙ）を算出し、確率分布関数ｇ（ｘ，ｙ）に基づいてｇ（ｘ，ｙ）上のＬ個の検出器に対応する位置を算出し、その位置の周囲の一定の領域をｇ（ｘ，ｙ）上の対応する検出器の分布領域とする。既存データに基づいて確率分布関数を求めることにより、検出器アレイ内の様々な分布パラメータを反映させることができ、したがって、より正確な位置を得ることができ、これは検出器アレイにより適している。
【００７１】
上述した第３の手法及び第４の手法において、「一定の領域」は、以下のように定義される。
【００７２】
１．確率分布関数において、検出器の決定された位置を中心とする半径Ｒの円を「一定の領域」とする。隣接する２つの検出器に対応する領域は、互いに重ならない。Ｒは、検出器に対応する斑点上に２次元ガウス曲面フィッティングを実行し、ガウス曲面フィッティングの関数の最大値の半分に相当する半径をＲとすることによって求めてもよい。Ｒは、他の手法で判定してもよく、具体的な値は、実際の状況に応じて合理的に判定される。
【００７３】
２．確率分布関数において、検出器の位置に対応する斑点が、ある値（例えば、最大値の１／３）より大きい領域を「一定の領域」とし、特定の値は、実際の状況に応じて合理的に決定される。
【００７４】
第１の手法及び第２の手法が採用される場合、すなわち、分離線を用いて、異なる検出器に対応する位置を区別する場合、この後に取得される全ての情報は、対応する検出器、すなわち対応する位置に分類される。２つの手法では、検出器は、より高い感度を有し、すなわち、検出されるパルスは、いずれも「無駄」にならない。
【００７５】
この多重化方法では、新たな検出パルスの時間差が、較正された位置情報の複数の点（ｘ，ｙ）に高度に一致する場合、どの点が時間差に最も適合しているかを判断することは困難であることがある。例えば、パルス情報がどの検出器に属するかを判定する際に、分離線に近いパルス情報は、正方形の真ん中のパルス情報に比べて、誤判定の可能性が高い。この場合、「この位置を中心とした一定の領域をｇ（ｘ，ｙ）上の対応する検出器の分布領域として設定する」第３、第４の手法では、位置の誤判定の可能性が高いパルスを直接的に除外することによって位置情報の誤りを低減できる。
【００７６】
確率分布関数ｇ（ｘ，ｙ）は、直接的には得られないが、統計的演算で得ることができる。まず一定量のデータを取得する必要がある。１０，０００個のパルスを検出すると、１０，０００個の（ｘ，ｙ）を取得することができ、表現に要求される精度の確率又は数で、１０，０００対の（ｘ，ｙ）をカウントする必要がある。例えば、ある特定の（ｘ，ｙ）における時間差の数が１００個であり、別の（ｘ，ｙ）における時間差の数が５０個である場合、得られた３次元グラフは、確率密度関数を表し、幾つかのデータを取得することによって確率密度関数を得て、検出器の位置情報を較正することができる。
【００７７】
上述したように、上記の実施形態は、本開示の技術的ソリューションを説明することのみを意図しており、本開示を限定するものではない。本開示は、上述の実施形態を参照して詳細に説明されているが、当業者であれば、上記実施形態の技術的ソリューションに変更を加えることができ、技術的ソリューションの幾つかの技術的特徴を同等に置き換えることができ、これらの変形又は置換は、本発明の実施形態の技術的ソリューションの精神及び範囲に含まれる。
【図１】
0006679715.tif 000006
【図２】
0006679715.tif 000007
【図３】
0006679715.tif 000008
【図４】
0006679715.tif 000009
【図５】
0006679715.tif 000010
【図６】
0006679715.tif 000011
【図７】
0006679715.tif 000012
【図８】
0006679715.tif 000013