JP7049709B2 - 信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置 - Google Patents

信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置 Download PDF

Info

Publication number
JP7049709B2
JP7049709B2 JP2020570150A JP2020570150A JP7049709B2 JP 7049709 B2 JP7049709 B2 JP 7049709B2 JP 2020570150 A JP2020570150 A JP 2020570150A JP 2020570150 A JP2020570150 A JP 2020570150A JP 7049709 B2 JP7049709 B2 JP 7049709B2
Authority
JP
Japan
Prior art keywords
signal
fitting
fitting model
final
model
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
JP2020570150A
Other languages
English (en)
Other versions
JP2021535454A (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.)
Raycan Technology Co Ltd
Original Assignee
Raycan Technology Co Ltd
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 Raycan Technology Co Ltd filed Critical Raycan Technology Co Ltd
Publication of JP2021535454A publication Critical patent/JP2021535454A/ja
Application granted granted Critical
Publication of JP7049709B2 publication Critical patent/JP7049709B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measurement Of Radiation (AREA)
  • Complex Calculations (AREA)

Description

本開示は、信号処理の技術分野、特に信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置に関する。
このセクションにおける説明は、本開示に関する背景情報を提供するのみであり、先行技術を構成するものではない。
陽電子放出断層撮影(Positron Emission Tomography、略称PET)は、臨床画像診断に放射性元素を用いる技術である。この技術のプロセスは、血流に添加される、又は生体組織の代謝プロセスに参加することができる化合物に陽電子放出放射性核種を標識すること、次いで放射性核種で標識された化合物を対象に注入することを含む。対象の体内で放射性核種によって放出された陽電子は約1mm移動し、次いで対象の体内にある負の電子と結合して電子対を消滅させ、このようにガンマ光子を生じさせるが、これをシンチレーション結晶によって受けて可視光に変換することができ、これが次に再構築のために光電変換器によって電気信号に変換され、これによって放射性核種の濃縮部位を決定し、活発な代謝の領域を特定して放射性核種の活性を評価する助けとなる。
PET又は他の関連分野において、信号を再構築するため、サンプリング後に信号のサンプリングポイントに対してフィッティング処理を実行することが必要である。先行技術において、サンプリングポイントをフィッティングするためのフィッティングモデルは、主に経験又は散布図に基づいて選択される。具体的には、各サンプリングポイントについて、図1(a)~図1(c)に示すように、信号の対象パラメータの測定値を横座標値とし、様々なフィッティングモデルの下でサンプリングポイントをフィッティングすることによって得られた対象パラメータのフィッティング値を縦座標値とすることによって、すべてのサンプリングポイントの様々な散布図を様々なフィッティングモデルの下で描くことができる。フィッティングモデルは、散布図における散布点間の分散度及び線形性を観察することによって選択される。
本開示を実現するプロセスにおいて、本発明者らは、先行技術において少なくとも以下の問題を発見した。
先行技術において、信号パラメータ及びフィッティング曲線パラメータは散布図を通して直感的に観察することができるが、散布点の分散度を測定することは困難であり、しばしば観察方法だけによって直感的な結論が得られる。結果として、後続の再構築結果が十分に正確でない可能性があり、これが信号再構築の精度に影響することになる。
本開示の目的は、信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置を提供して、信号再構築結果の正確性及び信号再構築の精度を改善することである。
上の目的を達成するため、本開示の実施形態において、以下の技術的解決策が提供される。
信号のためのフィッティングモデルを決定するための方法が、サンプリングされた信号をこの信号の収集された特性情報に基づいてセグメント化して複数の信号セグメントを得るステップと、複数のフィッティングモデルを用いることによって複数の信号セグメントにおいてサンプリングポイントをフィッティングするステップと、フィッティング結果に基づいて信号セグメントのそれぞれにおいて対象パラメータについてのフィッティング値を取得するステップと、対象パラメータについてのフィッティング値のそれぞれを、フィッティングモデルのそれぞれの下で、対象パラメータについての取得された測定値と順次比較し、比較結果に基づいて、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定するステップとを含む。
好ましくは、信号は、電気信号、光信号又は音声信号を含む。電気信号は、放射線検出器からの電気パルス信号を含む。
好ましくは、特性情報は、信号の振幅、立ち上がりスロープ、立ち上がり時間、立ち下がりスロープ及び/又は立ち下がり時間を含む。
好ましくは、複数のフィッティングモデルは、線形指数フィッティングモデル、二重指数フィッティングモデル及びグラフフィッティングモデルの少なくとも1つを含む。
好ましくは、グラフフィッティングモデルは、サンプリング閾値に対して上部で採用される三角フィッティングモデルと、サンプリング閾値に対して下部で採用される線形指数フィッティングモデル又は二重指数フィッティングモデルと、を含む。線形指数フィッティングモデルは、N点線形N点指数フィッティングモデル又はN点線形NプラスM点指数フィッティングモデルを含む。二重指数フィッティングモデルは、2N点二重指数フィッティングモデル又は2NプラスM点二重指数フィッティングモデルを含み、Nは1より大きい正の整数であり、Mは1以上の正の整数である。
好ましくは、対象パラメータは、信号のエネルギー、時間、エネルギー分解能及び/又は時間分解能を含む。
好ましくは、最終フィッティングモデルを決定するステップは、フィッティングモデルのそれぞれの下での信号セグメントのそれぞれにおける対象パラメータのフィッティング値のそれぞれと、対象パラメータについての取得された測定値と、の間の差値を順次計算し、差値に基づいて、複数のフィッティングモデルの中で、信号セグメントのそれぞれを再構築するための最終フィッティングモデルを決定するステップを含む。
好ましくは、最終フィッティングモデルを決定するステップは、複数のフィッティングモデルの中で、信号セグメントのそれぞれについての最終フィッティングモデルとして、信号セグメントのそれぞれにおける対象パラメータについての絶対差値が最小であるフィッティングモデルを決定するステップ、又は複数のフィッティングモデルの中で、信号セグメントのそれぞれについての最終フィッティングモデルとして、信号セグメントのそれぞれにおける対象パラメータについての差分布が最小であるフィッティングモデルを決定するステップを含む。
好ましくは、この方法は、最終フィッティングモデルと信号セグメントのそれぞれとの間の対応関係に従って最終フィッティングモデルを記憶するステップをさらに含む。
信号を再構築するための方法が次のステップ、すなわち、上の方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定するステップと、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するステップとを含む。
信号を再構築するための方法が次のステップ、すなわち、上の方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、実際にサンプリングされた信号のための最終フィッティングモデルを決定するステップと、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するステップとを含む。
好ましくは、参照信号と実際にサンプリングされた信号を比較するステップは、実際にサンプリングされた信号の特性情報を利用することによって、実際にサンプリングされた信号と同一又は類似の特性情報を有する参照信号を検索するステップと、参照信号が検索された後、参照信号と最終フィッティングモデルとの間の対応関係に従って参照信号のための最終フィッティングモデルを検索し、参照信号のための検索された最終フィッティングモデルを、実際にサンプリングされた信号のための最終フィッティングモデルとして決定するステップとを含む。
信号フィッティングモデルを決定するための装置が、サンプリングされた信号をこの信号の収集された特性情報に基づいてセグメント化して複数の信号セグメントを得るように構成されたセグメント化ユニットと、複数のフィッティングモデルを用いることによって複数の信号セグメントにおいてサンプリングポイントをフィッティングするように構成されたフィッティングユニットと、フィッティング結果に基づいて信号セグメントのそれぞれにおいて対象パラメータについてのフィッティング値を取得するように構成された取得ユニットと、対象パラメータについてのフィッティング値のそれぞれを、フィッティングモデルのそれぞれの下で、対象パラメータについての取得された測定値と比較し、比較結果に基づいて、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定するように構成された比較及び決定ユニットとを含む。
好ましくは、比較及び決定ユニットは、フィッティングモデルのそれぞれの下での信号セグメントのそれぞれにおける対象パラメータについてのフィッティング値のそれぞれと、対象パラメータについての取得された測定値と、の間の差値を順次計算し、差値に基づいて、複数のフィッティングモデルの中で、信号セグメントのそれぞれを再構築するための最終フィッティングモデルを決定するように具体的に構成することができる。
好ましくは、この装置は、最終フィッティングモデルと信号セグメントのそれぞれとの間の対応関係に従って最終フィッティングモデルを記憶するように構成された記憶ユニットをさらに含む。
信号を再構築するための装置が、上の方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成された決定ユニットと、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するように構成された再構築ユニットとを含む。
信号を再構築するための装置が、上の方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成された比較及び決定ユニットと、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するように構成された再構築ユニットとを含む。
好ましくは、比較及び決定ユニットは、実際にサンプリングされた信号の特性情報を利用することによって、実際にサンプリングされた信号と同一又は類似の特性情報を有する参照信号を検索して、参照信号が検索された後、参照信号と最終フィッティングモデルとの間の対応関係に従って参照信号のための最終フィッティングモデルを検索し、参照信号のための検索された最終フィッティングモデルを、実際にサンプリングされた信号のための最終フィッティングモデルとして決定するように具体的に構成されている。
本開示の上の実施形態によって提供される技術的解決策から、サンプリングされた信号が信号の特性情報に従ってセグメント化され、複数のフィッティングモデルを用いることによって得られる信号の対象パラメータのフィッティング値が測定値と比較されるため、最終フィッティングモデルが比較結果に従って決定されることが分かる。信号の異なる特性を考慮することによって異なるフィッティングモデルを採用することによって、信号の後続のフィッティング結果がより正確になる可能性があり、したがって信号再構築結果の正確性及び信号再構築の精度を改善するための目的を達成することができる。
本開示又は先行技術の実施形態における技術的解決策をより明確に説明するため、実施形態又は先行技術の説明に用いることが要求される添付の図面を以下に簡単に紹介する。示したような図面は、本開示に記載されたようないくつかの実施形態の単なる例示であるということが明らかである。創造的な作業を伴うことなく、示したような図面に対する様々な代替案を理解することができるということが当業者によって理解されるべきである。
先行技術において4点線形及び5点指数フィッティングモデルを用いてサンプリングポイントがフィッティングされた散布図である。 先行技術において8点二重指数フィッティングモデルを用いてサンプリングポイントがフィッティングされた散布図である。 先行技術においてグラフフィッティングモデルを用いてサンプリングポイントがフィッティングされた散布図である。 本開示の一実施形態による信号のためのフィッティングモデルを決定するための方法の概略フローチャートである。 N点線形及びN点指数フィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 N点線形及びNプラスM点指数フィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 2N点二重指数フィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 2NプラスM点二重指数フィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 下方曲線がN点線形及びNプラスM点指数フィッティングモデルによってフィッティングされた、グラフフィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 下方曲線が2N点二重指数フィッティングモデルによってフィッティングされた、グラフフィッティングモデルを用いてサンプリングポイントをフィッティングすることによって得られたフィッティング曲線である。 4点線形及び5点指数フィッティングモデルの下での電気パルス信号a及びbのエネルギー差分布のヒストグラムである。 8点二重指数フィッティングモデルの下での電気パルス信号a及びbのエネルギー差分布のヒストグラムである。 グラフフィッティングモデルの下での電気パルス信号a及びbのエネルギー差分布のヒストグラムである。 4点線形及び5点指数フィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図である。 8点二重指数フィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図である。 グラフフィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図である。 60~130Vの振幅範囲内の4点線形及び4点指数フィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図を示す。 60~130Vの振幅範囲内の4点線形及び5点指数フィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図を示す。 60~130Vの振幅範囲内の、下方曲線が9点二重指数フィッティングモデルによってフィッティングされているグラフフィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図を示す。 60~130Vの振幅範囲内の、下方曲線が8点二重指数フィッティングモデルによってフィッティングされた、グラフフィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図を示す。 本開示の一実施形態による信号を再構築するための方法の概略フローチャートである。 本開示におけるグラフフィッティングモデルを用いて電気パルス信号をフィッティングすることによって得られたフィッティング曲線、及び電気パルス信号の実際の測定曲線を示す。 既存の8点二重指数フィッティングモデルを用いて電気パルス信号をフィッティングすることによって得られたフィッティング曲線、及び電気パルス信号の実際の測定曲線を示す。 図7における信号を再構築するための方法を用いることによって得られたフィッティング曲線を示す。 既存の多閾値サンプリング方法を用いることによって得られたフィッティング曲線を示す。 本開示の一実施形態による信号を再構築するための他の方法の概略フローチャートである。 本開示の一実施形態による信号のためのフィッティングモデルを決定するための装置を示す。 本開示の一実施形態による信号を再構築するための装置を示す。 本開示の一実施形態による信号を再構築するための他の装置を示す。
本開示の実施形態における技術的解決策を、添付の図面を参照して以下の説明において明確且つ包括的に説明する。記載の実施形態は、これらのすべてを網羅的に詳述するのではなく、本開示のいくつかの例示的な実施形態を説明するためにのみ提供され、これらは本開示又は特許請求の範囲を限定するように解釈されるべきではないということが明らかである。本明細書に記載の実施形態に対する様々な代替案が、創造的な作業を伴うことなく、そして本開示の範囲から逸脱することなく、当業者によって使用され得るということが理解されるべきである。
特に、ある要素が他の要素「に配置されて」いると言うとき、これは他の要素に直接配置されていても、又は中間要素があってもよい。ある要素が他の要素に「結合又は接続されて」いると言うとき、これは他の要素に直接結合又は接続されていてもよく、又は中間要素がある。本明細書で用いる「結合又は接続」という用語は、電気的結合又は接続及び/又は機械的又は物理的結合又は接続を含むことができる。本明細書で用いる「含む(comprise or include)」という用語は、特徴、ステップ又は要素の存在を指すが、1又は2以上の他の特徴、ステップ又は要素の存在又は追加を排除するものではない。本明細書で用いる「及び/又は」という用語は、関連する列挙した項目の1又は2以上のありとあらゆる組み合わせを含む。「a」、「an」、「one」、「the」及び同様の用語は、文中に明確に記載されていない限り、単数及び複数を含む。
特に指定しない限り、本明細書で用いるすべての技術用語及び科学用語は、本開示に関連する技術分野における当業者によって通常理解されるような一般的な意味を有する。本明細書で用いる用語は、具体的な実施形態を説明することを目的とするが、本開示を限定するように意図されるものではない。
加えて、本明細書で用いる「第1」、「第2」、「第3」などの用語は、説明を目的とし、類似の物を互いに区別するためだけのものであり、これらはその順序を表さず、相対的な重要性の指標又は暗示として理解することもできない。加えて、特に指示しない限り、本明細書で用いる「複数」は2又は3以上を意味する。
以下、本開示の実施形態による、信号のためのフィッティングモデルを決定するための方法、信号を再構築するための方法、及びその装置を、添付の図面を参照して詳細に説明する。
図2に示すように、本開示の一実施形態において、信号のためのフィッティングモデルを決定するための方法が提供され、これは以下のステップを含む。
ステップS1:サンプリングされた信号をこの信号の収集された特性情報に基づいてセグメント化して複数の信号セグメントを得る。
信号をサンプリングした後、信号の収集された特性情報に基づいて信号をセグメント化することができる。具体的には、サンプリングされた信号は、信号の振幅(すなわち、振幅値)、立ち上がりスロープ、立ち上がり時間、立ち下がりスロープ、及び/又は立ち下がり時間などに基づいて、複数の信号セグメントにセグメント化することができる。例えば、信号は、振幅に基づいて、0~M1、M1~M2、M2~M3、…、Mk~Mn、及び>Mnなどのような、複数の信号セグメントにセグメント化することができる。他の一例では、信号は、時間に基づいて、0~T1、T1~T2、T2~T3、…、Ti~Tj、及び>Tjなどのような、複数の信号セグメントにセグメント化することができる。M1、M2、M3、Mk及びMnは、信号の異なる振幅(電圧又は電流などのような)を表し、T1、T2、T3、Ti及びTjは、信号の異なる時間値を表す。各信号セグメントの間隔は同じでも異なってもよいということが留意されるべきである。
信号は、電気信号、光信号、音声信号などであり得る。信号は、同じ装置からの同じ1つの信号でも、異なる装置からの様々な信号であってもよい。電気信号は、シンチレーション結晶検出器又はガス電離検出器のような放射線検出器からの電気パルス信号を含むことができる。
特性情報は、信号をサンプリングした結果から得られるデータとすることができる。特性情報は、信号の振幅、立ち上がりスロープ、立ち上がり時間、立ち下がりスロープ及び/又は立ち下がり時間などを含むことができる。
S2:複数のフィッティングモデルを用いることによって複数の信号セグメントにおいてサンプリングポイントをフィッティングする。
信号がセグメント化された後、複数のフィッティングモデルを用いて各信号セグメントにおいてサンプリングポイントをフィッティングして、複数のフィッティングモデルの下で完全な信号のフィッティング曲線を取得し、これによって信号の再構築波形を取得することができる。
フィッティングモデルは、線形及び指数フィッティングモデル、二重指数フィッティングモデル及びグラフフィッティングモデルなどの少なくとも1つのタイプを含むことができるが、これらに限定されない。グラフフィッティングモデルは、サンプリング閾値に対して上部で採用される三角フィッティングモデル、及びサンプリング閾値に対して下部で採用される線形及び指数フィッティングモデル又は二重指数フィッティングモデルを含むことができるが、これらに限定されない。線形及び指数フィッティングモデルは、N点線形及びN点指数フィッティングモデル又はN点線形及びNプラスM点指数フィッティングモデルなどを含むことができるが、これらに限定されない。二重指数フィッティングモデルは、2N点二重指数フィッティングモデル又は2NプラスM点二重指数フィッティングモデルなどを含むことができるが、これらに限定されない。Nは、多閾値サンプリング方法によって信号をサンプリングすることによって得られる立ち上がりエッジ部分又は立ち下がりエッジ部分上のサンプリングポイントの数を表し、これは1より大きい正の整数であり、Mは、多閾値方法によって信号をサンプリングすることに基づいて、時間間隔に従って信号をサンプリングすることによってさらに得られるサンプリングポイントの数を表し、これは1以上の正の整数である。多閾値サンプリング方法及び時間間隔サンプリング方法の具体的なプロセスに関しては先行技術を参照することができ、これはここでは詳述しない。
N点線形及びN点指数フィッティングモデル:最大振幅のサンプリングポイントを除いて、信号の立ち上がりエッジ部分上のN個のサンプリングポイントを基準点と見なし、これらをフィッティングして直線を生じさせ、信号の立ち下がりエッジ部分上のN個のサンプリングポイントを基準点と見なし、これらをフィッティングして指数曲線方程式に従って指数曲線を生じさせ、したがって生じさせられた直線及び指数曲線は、図3(a)に示すように、信号の再構築波形を構成し、N=4である。
N点線形及びNプラスM点指数フィッティングモデル:信号の立ち上がりエッジ部分上のN個のサンプリングポイントを基準点と見なし、これらをフィッティングして直線を生じさせ、信号の立ち下がりエッジ部分上のN個のサンプリングポイント及びM個のサンプリングポイントを基準点と見なし、これらをフィッティングして指数曲線方程式に従って指数曲線を生じさせ、したがって生じさせられた直線及び指数曲線は、図3(b)に示すように、信号の再構築波形を構成し、N=4及びM=1である。
2N点二重指数フィッティングモデル:信号に対して多閾値サンプリングを実行することによって得られた2N個のサンプリングポイントを基準点と見なし、これらをフィッティングして二重指数曲線方程式に従って二重指数曲線を生じさせ、したがって得られた二重指数曲線は、図3(c)に示すように、信号の再構築波形を構成し、N=4である。
2NプラスM点二重指数フィッティングモデル:信号に対して多閾値サンプリングを実行することによって得られた2N個のサンプリングポイント、及び時間間隔で信号をサンプリングすることによって得られたM個のサンプリングポイントを基準点と見なし、これらをフィッティングして二重指数曲線方程式に従って二重指数曲線を生じさせ、したがって得られた二重指数曲線は、図3(d)に示すように、信号の再構築波形を構成し、N=4及びM=1である。
グラフフィッティングモデル:サンプリングされた2N+M個のサンプリングポイントを基準点と見なし、信号の最大サンプリング閾値のような特定のサンプリング閾値に対して上部において、立ち上がりエッジ部分上の最後のサンプリングポイント、ピークでのサンプリングポイント、及び立ち下がりエッジ部分上の最初のサンプリングポイント(「最初」及び「最後」は時系列で示す)のような、3つのサンプリングポイントを直線で互いと結合して三角を形成し、信号のサンプリング閾値に対して下部において、サンプリングポイントを線形及び指数又は二重指数フィッティングモデルなどに従ってフィッティングして、図3(e)及び3(f)に示すように、完全な信号再構築波形を得る。図3(e)において、最大サンプリング閾値に対して下部における信号は、線形及び指数モードに従ってフィッティングされ、図3(f)において、最大サンプリング閾値に対して下部における信号は、二重指数フィッティングモデルに従ってフィッティングされている。
S3:フィッティング結果に基づいて信号セグメントのそれぞれにおいて対象パラメータについてのフィッティング値を取得する。
様々なフィッティングモデルを用いてすべての信号セグメントにおいてすべてのサンプリングポイントをフィッティングした後、各信号セグメントにおける各サンプリングポイントでの対象パラメータについてのフィッティング値は、得られたフィッティング結果から直接得ることができ、又はフィッティング結果から得られたデータを計算することによって得ることもできる。
フィッティング結果は、フィッティング曲線及び/又はフィッティングデータであり得る。
対象パラメータは、実際の用途に従って考慮する必要があるパラメータを指すことができ、これは、信号のエネルギー、時間、及び/又はエネルギー分解能などを含むことができるが、これらに限定されない。エネルギー分解能は、信号エネルギースペクトルの全電力ピークのエネルギーの期待値に対する全電力ピークの半値幅の比率を指し、信号エネルギー間隔を互いに区別する能力を表すことができ、エネルギー分解能が低いほど、エネルギースペクトルにおけるエネルギー間隔を互いに区別する能力が強くなる。
対象パラメータが時間であるとき、時間値、すなわち、フィッティング値は、得られたフィッティング曲線から直接抽出することができる。対象パラメータがエネルギーであるとき、フィッティング値は、電圧又は電流のような、抽出された振幅を処理することによって得ることができる。例えば、エネルギーは、電圧振幅を積分することによって得ることができ、その具体的な計算プロセスは先行技術を参照することができ、これはここでは詳述しない。
同様に、対象パラメータがエネルギー分解能などであるとき、フィッティング値は、抽出された時間又は振幅を処理することによって得ることができ、その具体的なプロセスはここでは限定されない。
S4:対象パラメータについてのフィッティング値のそれぞれを、フィッティングモデルのそれぞれの下で、対象パラメータについての取得された測定値と順次比較し、比較結果に従って、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定する。
各フィッティングモデルの下ですべての信号セグメントにおいてすべてのサンプリングポイントで対象パラメータについてのフィッティング値を得た後、各フィッティングモデルの下での各信号セグメントにおける各サンプリングポイントでの対象パラメータについてのフィッティング値のそれぞれを、サンプリングポイントでの対象パラメータの測定値(すなわち、真の値)と順次比較することができる。次いで比較結果に従って、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定することができる。具体的には、各信号のため、各フィッティングモデルの下での各信号セグメントにおける各サンプリングポイントでの対象パラメータについてのフィッティング値の大きさを、対応する測定値の大きさと直接比較することができ、又は各フィッティングモデルの下での各信号セグメントにおける各サンプリングポイントでの対象パラメータについてのフィッティング値と、対応する測定値と、の間の差値を計算することができ、次いで差値に従って、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定することができる。
具体的には、各信号のため、第1に各フィッティングモデルの下での各信号セグメントにおける各サンプリングポイントでの対象パラメータについてのフィッティング値と、対応する測定値と、の間の差値を、事前設定された順序で、例えば、サンプリングポイントの横座標又は縦座標の値の大きさの昇順で、順次計算することができる。例えば、対象パラメータがエネルギーである状況において、各フィッティングモデルの下での各信号セグメントにおける各サンプリングポイントでのエネルギーについてのフィッティング値と、対応する測定値と、の間の差値を順次計算することができる。第2に、対象パラメータについてのすべての差値を計算した後、差値に基づいて信号セグメントのそれぞれを再構築するための最終フィッティングモデルを決定することができる。具体的には、複数のフィッティングモデルの中で、各信号セグメントにおける対象パラメータについての絶対差値が最小となるフィッティングモデルを、それぞれの信号セグメントについての最終フィッティングモデルとして決定することができる。あるいは、複数のフィッティングモデルの中で、信号セグメントのそれぞれにおける対象パラメータについての差分布が最小となるフィッティングモデルを、次の方法に従ってそれぞれの信号セグメントについての最終フィッティングモデルとして決定して、すべての信号セグメントについてのすべての最終フィッティングモデルを決定することができる。
各フィッティングモデル及び各信号のため、n行掛ける(×)2列の行列を確立することができ、行列の第1の列は、信号の各サンプリングポイントでの対象パラメータの特性値、すなわち、上の特性情報における各パラメータの値を記録することができ、行列の第2の列は、信号の各サンプリングポイントでの対象パラメータのフィッティング値と測定値との間の差値を記録することができ、又は2行×n列の行列を確立することができ、行列の第1の行は、信号の各サンプリングポイントでの対象パラメータの特性値を記録することができ、行列の第2の行は、信号の各サンプリングポイントでの対象パラメータのフィッティング値と測定値との間の差値を記録することができ、nは特性値の数を表し、次に、行列における特性値及び差値を、それぞれ、横座標及び縦座標として用いて、各フィッティングモデルの下での各信号セグメントにおける対象パラメータの差分布図を生じさせることができ、又は特性値をそれらの大きさに従ってまずランク付けすることができ、次いで各フィッティングモデルの下での各信号セグメントにおける対象パラメータの差分布図をランキング順序に基づいて生じさせることができ、最後に、差分布図における最小差分布に対応するフィッティングモデルを、それぞれの信号セグメントについての最終フィッティングモデルとして決定することができる。
最終フィッティングモデルは、複数のフィッティングモデルの1又は2以上とすることができるということが留意されるべきである。それぞれの信号セグメントについての最終フィッティングモデルは、互いに同じでも異なってもよい。
例えば、図4(a)~4(c)は、3つの異なるフィッティングモデル、すなわち、それぞれ、4点線形及び5点指数フィッティングモデル、8点二重指数フィッティングモデル及びグラフフィッティングモデルの下での電気パルス信号a及びbのエネルギー差分布のヒストグラムを示す。E及びEはそれぞれ電気パルス信号a及びbのエネルギーの測定値を表し、
Figure 0007049709000001
Figure 0007049709000002
及び
Figure 0007049709000003
はそれぞれ上の3つのフィッティングモデルの下での電気パルス信号aのエネルギーのフィッティング値を表し、
Figure 0007049709000004
Figure 0007049709000005
及び
Figure 0007049709000006
はそれぞれ上の3つのフィッティングモデルの下での電気パルス信号bのエネルギーのフィッティング値を表し、V及びVはそれぞれフィッティング曲線から抽出された電気パルス信号a及びbの振幅(すなわち、電圧)を表す。3つの図に示すように、電気パルス信号aについて、8点二重指数フィッティングモデル及びグラフフィッティングモデルのいずれかの下でのエネルギー差値は比較的小さく、したがってこれら2つのフィッティングモデルのいずれかを選択することができるということが分かる。電気パルス信号bについて、グラフフィッティングモデルの下でエネルギー差が最小であり、したがってグラフフィッティングモデルが最適な選択である。
例えば、図5(a)~5(c)は、3つの異なるフィッティングモデル、すなわち、それぞれ、4点線形及び5点指数フィッティングモデル、8点二重指数フィッティングモデル、及び8点二重指数フィッティングモデルの下で下方曲線がフィッティングされたグラフフィッティングモデルの下での複数の電気パルス信号のエネルギー差分布の図を示す。3つの図から、電気パルス信号のため、62V~74Vの振幅範囲内ではグラフフィッティングモデルの下でエネルギー差が最小であり、したがってこの場合、フィッティングにグラフフィッティングモデルを選択することが好ましく、136Vより大きい(>)振幅範囲内では4点線形及び5点指数フィッティングモデルの下でエネルギー差が最小であり、したがってこの場合、フィッティングに4点線形及び5点指数フィッティングモデルを選択することが好ましいということが分かる。
例えば、図6(a)~6(d)は、4つの異なるフィッティングモデル、すなわち、それぞれ、4点線形及び4点指数フィッティングモデル、4点線形及び5点指数フィッティングモデル、9点二重指数フィッティングの下で下方曲線がフィッティングされたグラフフィッティングモデル、及び8点二重指数フィッティングの下で下方曲線がフィッティングされたグラフフィッティングモデルの下での60V~130Vの振幅範囲内の複数の電気パルス信号のエネルギー差分布の図を示し、横座標は電気パルス信号の振幅(V)を表し、縦座標はエネルギー差を表す。4つの図から、電気パルス信号のため、60V~80Vの振幅範囲内では8点二重指数フィッティングの下で下方曲線がフィッティングされたグラフフィッティングモデルの下でエネルギー差が最小であり、したがってこの場合、フィッティングにグラフフィッティングモデルを選択することが好ましく、80V~130Vの振幅範囲内では4点線形及び5点指数フィッティングモデルの下でエネルギー差が最小であり、したがってこの場合、フィッティングに4点線形及び5点指数フィッティングモデルを選択することが好ましいということが分かる。
上記から、本開示の実施形態において、信号の特性情報に従って、サンプリングされた信号をセグメント化し、各信号セグメントにおける対象パラメータのフィッティング値を複数のフィッティングモデルの下で測定値と比較して、比較結果に基づいて各信号セグメントについての最終フィッティングモデルを決定することができ、これは、信号の異なる特性を考慮することによって異なるフィッティングモデルを採用することができるということを反映しているということが分かる。したがって、信号の後続のフィッティング結果がより正確になり、これによって信号再構築結果の正確性及び信号再構築の精度を改善することができる。また、ヒストグラムの形式を用いてフィッティングモデルを比較することによって様々なフィッティングモデルの中の比較可能性が向上し、より直感的にすることになる。
本開示のさらなる一実施形態において、この方法は、S5:信号セグメントのそれぞれと最終フィッティングモデルとの間の対応関係に従って、決定された最終フィッティングモデルを記憶するステップをさらに含むことができる。
すべての信号セグメントについての最終フィッティングモデルが決定された後、決定された最終フィッティングモデルは、信号セグメントと最終フィッティングモデルとの間の対応関係に従って記憶することができ、後の使用のためにこれらを直接検索することができるようになり、信号再構築の効率が改善される。
本開示の実施形態において、信号を再構築するための方法が提供される。図7に示すように、この方法は以下のステップを含むことができる。
P1:上の決定方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定する。
取得された信号をサンプリングした後、上の決定方法を利用して、実際にサンプリングされた信号のための最終フィッティングモデルを決定することができる。具体的な決定プロセスに関しては図2~図6の説明を参照することができ、これはここでは詳述しない。
最終フィッティングモデルは、線形及び指数フィッティングモデル、二重指数フィッティングモデル及びグラフフィッティングモデルの少なくとも1つを含むことができる。
P2:決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築する。
実際にサンプリングされた信号のための最終フィッティングモデルを決定した後、決定された最終フィッティングモデルを用いて各信号セグメントにおいてサンプリングポイントをフィッティングして、信号のフィッティング曲線を取得し、これによって信号再構築を実現する。
例えば、各信号セグメントについての最終フィッティングモデルがN点線形及びNプラスM点指数フィッティングモデルであると決定した後、N点線形及びNプラスM点指数フィッティングモデルを用いて各信号セグメントにおいてサンプリングポイントをフィッティングして、信号のフィッティング曲線を得ることができる。例えば、信号のための最終フィッティングモデルが2NプラスM点二重指数フィッティングモデル及びグラフフィッティングモデルであると決定した後、2つのモデルをそれぞれ用いて、それぞれの信号セグメントにおいてサンプリングポイントをフィッティングすることができ、又はより速いフィッティング速度又はより低いエネルギー消費を伴う可能性がある2つのモデルの1つを用いてサンプリングポイントをフィッティングして、信号のフィッティング曲線を得ることができる。N点線形及びNプラスM点指数フィッティングモデル、2NプラスM点二重指数フィッティングモデル及びグラフフィッティングモデルの下で信号のサンプリングポイントをフィッティングするために、ステップS2におけるフィッティングモデルの説明を参照することができ、これはここでは説明しない。
フィッティング曲線から抽出された、振幅及び時間のような、対応するデータを計算することによって、信号のエネルギー又はエネルギー分解能を得ることができる。計算方法に関しては先行技術を参照することができ、これはここでは詳述しない。
本開示の実施形態による信号を再構築するための方法を用いることによって、信号再構築結果の正確性及び信号再構築の精度を改善することができる。
上の実施形態の利点を説明する例を与える。
図8は、8点二重指数フィッティングモデルによって下方曲線がフィッティングされた、本開示におけるグラフフィッティングモデルを用いて電気パルス信号をフィッティングすることによって得られたフィッティング曲線、及び電気パルス信号の実際の測定曲線を示す。図9は、既存の8点二重指数フィッティングモデルを用いて電気パルス信号をフィッティングすることによって得られたフィッティング曲線、及び電気パルス信号の実際の測定曲線を示す。(T9、V9)は9番目のサンプリングポイントの座標を表し、実線は測定曲線を表し、破線はフィッティング曲線を表す。2つの図に示すように、8点二重指数フィッティングモデルによって下方曲線がフィッティングされた本開示におけるグラフフィッティングモデルを用いることによって得られたフィッティング曲線は実際の測定曲線に非常に近いが、既存の8点二重指数フィッティングモデルを用いることによってフィッティングされた曲線は実際の測定曲線から遠いということが分かる。本開示の信号再構築法を用いることによって、信号再構築結果の正確性を改善することができるということが分かる。
図10及び図11は、それぞれ、本開示の実施形態による信号を再構築するための方法を用いることによって得られたフィッティング曲線、及び既存の多閾値サンプリング方法を用いることによって得られたフィッティング曲線を示す。2つの図を互いに比較することによって、本開示の実施形態による信号を再構築するための方法によって得られたエネルギー分解能は、既存の多閾値サンプリング方法によって得られたエネルギー分解能よりも約1.6%低く、これはエネルギー精度の向上を示すということが分かる。また、本開示の実施形態による信号を再構築するための方法を用いることによって、1274Kev付近に現れる光電ピークを識別することができるが、これは既存の多閾値サンプリング方法を用いることによっては全く識別不可能である。光電ピークを識別する能力は、異なる核種を区別するために重要である。本開示の実施形態による信号を再構築するための方法を用いることによって、光電ピークを識別する能力を改善することができるということが分かり、これは異なる核種を区別する助けとなる。
本開示の実施形態において、信号を再構築するための他の方法が提供される。図12に示すように、この方法は以下のステップを含むことができる。
P1':上の決定方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、実際にサンプリングされた信号のための最終フィッティングモデルを決定する。
取得された信号をサンプリングした後、上の決定方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較することができる。信号の特性情報、例えば、振幅、時間、立ち上がりエッジの立ち上がり時間、立ち上がりスロープなどを用いることによって、実際にサンプリングされた信号と同一又は類似の特性情報を有する参照信号を検索することができる。参照信号が検索された後、参照信号とその最終フィッティングモデルとの間の対応関係に従って参照信号のための最終フィッティングモデルを検索し、参照信号のための検索された最終フィッティングモデルを、実際にサンプリングされた信号のための最終フィッティングモデルとして決定することができる。参照信号のための最終フィッティングモデルの具体的な決定プロセスに関しては図2~図6の説明を参照することができ、これはここでは詳述しない。
最終フィッティングモデルは、線形及び指数フィッティングモデル、二重指数フィッティングモデル及びグラフフィッティングモデルの少なくとも1つを含むことができる。
P2':決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築する。
このステップの詳細については、上のステップP2を参照することができ、これはここでは詳述しない。
本開示の実施形態において、信号のためのフィッティングモデルを決定するための装置が提供される。図13に示すように、この装置は、サンプリングされた信号をこの信号の収集された特性情報に基づいてセグメント化して複数の信号セグメントを得るように構成することができるセグメント化ユニット121と、複数のフィッティングモデルを用いて信号セグメントのそれぞれにおいてサンプリングポイントをフィッティングするように構成することができるフィッティングユニット122と、フィッティング結果に基づいて複数の信号セグメントにおいて対象パラメータについてのフィッティング値を取得するように構成することができる取得ユニット123と、対象パラメータについてのフィッティング値のそれぞれを、フィッティングモデルのそれぞれの下で、対象パラメータについての取得された測定値と順次比較し、比較結果に基づいて、複数のフィッティングモデルの中で、信号を再構築するための最終フィッティングモデルを決定するように構成することができる比較及び決定ユニット124とを含むことができる。
上記のユニットの詳細について、図2に示すような方法の実施形態を参照することができ、これはここでは詳述しない。
本開示の実施形態による上記の装置を利用することによって、続いて信号を再構築する精度を改善するための目的を達成することができる。また、様々なフィッティングモデルの中の比較可能性が向上し、より直感的にすることができる。
本開示の実施形態において、信号を再構築するための装置が提供される。図14に示すように、この装置は、上の決定方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成することができる決定ユニット131と、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するように構成することができる再構築ユニット132とを含むことができる。
上記のユニットの詳細について、図7に示すような方法の実施形態を参照することができ、これはここでは詳述しない。
本開示の実施形態において、信号を再構築するための他の装置が提供される。図15に示すように、この装置は、上の決定方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成することができる比較及び決定ユニット141と、決定された最終フィッティングモデルを利用することによって、実際にサンプリングされた信号を再構築するように構成することができる再構築ユニット142とを含むことができる。
上記のユニットの詳細について、図12に示すような方法の実施形態を参照することができ、これはここでは詳述しない。
本開示の実施形態による上記の装置を利用することによって、信号を再構築する精度を改善するための目的及びしたがって再構築結果の正確性を達成することができる。
上記の装置におけるユニットによって実装される機能は、コンピュータにおけるプロセッサによってメモリに記憶された命令を実行することによって実装することもできるということが留意されるべきである。
本開示の実施形態による上記の方法及び装置は、PET、多電圧閾値サンプリング(Multi-Voltage Threshold、略称MVT)などの分野に応用されるだけに限定されず、信号再構築が要求される任意の分野に応用することもできるということが留意されるべきである。
上の実施形態に記載された装置、ユニットなどは、コンピュータチップ及び/又はエンティティによって具体的に実装、又は特定の機能を備えた製品によって実装することができる。便宜上、装置の異なる機能的ユニットをそれぞれ説明している。もちろん、本開示を実装するとき、個々のユニットの機能は、同じ1又は2以上のコンピュータチップに具現化することができる。
上記の実施形態又はフローチャートに記載の方法ステップが本開示に提供されているが、より多くの又はより少ないステップを、従来の又は日常的な作業と共にこれらの方法に含むことができる。論理的に必要な因果関係がないステップの実行順序は、本開示における実施形態において提供されるものに限定されない。
本明細書における様々な実施形態は、様々な実施形態にわたって参照されている同一又は類似の部分で漸進的な方法で説明されているが、各実施形態の説明は、他の実施形態との違いに焦点を合わせている。
上記の実施形態は、当業者が本開示を理解及び実践することを容易にするように説明されている。これらの実施形態に様々な修正を加え、本明細書に記載の一般原理を創造的な作業なしに他の実施形態に応用することも当業者には明らかである。したがって、本開示は、上記の実施形態に限定されるものではなく、本開示の範囲から逸脱することなく本開示に従って当業者によって行われる改善及び修正は、本開示の保護範囲に入るべきである。

Claims (18)

  1. 信号のためのフィッティングモデルを決定するための方法であって、
    サンプリングされた信号を前記信号の収集された特性情報に基づいてセグメント化して複数の信号セグメントを得るステップと、
    複数のフィッティングモデルを用いることによって前記複数の信号セグメントにおいてサンプリングポイントをフィッティングするステップと、
    フィッティング結果に基づいて前記信号セグメントのそれぞれにおいて対象パラメータについてのフィッティング値を取得するステップと、
    前記対象パラメータについての前記フィッティング値のそれぞれを、前記フィッティングモデルのそれぞれの下で、前記対象パラメータについての取得された測定値と順次比較し、比較結果に基づいて、前記複数のフィッティングモデルの中で、前記信号を再構築するための最終フィッティングモデルを決定するステップと、
    を含むことを特徴とする、前記方法。
  2. 前記信号が、電気信号、光信号又は音声信号を含み、
    前記電気信号が、放射線検出器からの電気パルス信号を含むことを特徴とする、請求項1に記載の方法。
  3. 前記特性情報が、信号の振幅、立ち上がりスロープ、立ち上がり時間、立ち下がりスロープ及び/又は立ち下がり時間を含むことを特徴とする、請求項1に記載の方法。
  4. 前記複数のフィッティングモデルが、線形及び指数フィッティングモデル、二重指数フィッティングモデル並びにグラフフィッティングモデルの少なくとも1つのタイプを含むことを特徴とする、請求項1に記載の方法。
  5. 前記グラフフィッティングモデルが、サンプリング閾値に対して上部で採用される三角フィッティングモデルと、前記サンプリング閾値に対して下部で採用される線形及び指数フィッティングモデル又は二重指数フィッティングモデルと、を含み、
    前記線形及び指数フィッティングモデルが、N点線形及びN点指数フィッティングモデル又はN点線形及びNプラスM点指数フィッティングモデルを含み、
    前記二重指数フィッティングモデルが、2N点二重指数フィッティングモデル又は2NプラスM点二重指数フィッティングモデルを含み、Nが1より大きい正の整数であり、Mが1以上の正の整数であることを特徴とする、請求項4に記載の方法。
  6. 前記対象パラメータが、信号のエネルギー、時間、エネルギー分解能及び/又は時間分解能を含むことを特徴とする、請求項1に記載の方法。
  7. 前記最終フィッティングモデルを決定するステップが、フィッティングモデルのそれぞれの下での信号セグメントのそれぞれにおける対象パラメータについてのフィッティング値のそれぞれと、前記対象パラメータについての測定値と、の間の差値を順次計算し、前記差値に基づいて、前記複数のフィッティングモデルの中で、前記信号セグメントのそれぞれを再構築するための最終フィッティングモデルを決定するステップを含むことを特徴とする、請求項1に記載の方法。
  8. 前記最終フィッティングモデルを決定するステップが、
    前記複数のフィッティングモデルの中で、信号セグメントのそれぞれにおける対象パラメータについての絶対差値が最小であるフィッティングモデルを、前記信号セグメントのそれぞれについての最終フィッティングモデルとして決定するステップ、又は
    前記複数のフィッティングモデルの中で、前記信号セグメントのそれぞれにおける対象パラメータについての差分布が最小であるフィッティングモデルを、前記信号セグメントのそれぞれについての最終フィッティングモデルとして決定するステップ
    を含むことを特徴とする、請求項7に記載の方法。
  9. 前記信号セグメントのそれぞれと最終フィッティングモデルとの間の対応関係に従って、決定された最終フィッティングモデルを記憶するステップをさらに含むことを特徴とする、請求項1に記載の方法。
  10. 信号を再構築するための方法であって、
    請求項1~9のいずれかに記載の方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定するステップと、
    決定された最終フィッティングモデルを利用することによって、前記実際にサンプリングされた信号を再構築するステップと
    を含むことを特徴とする、前記方法。
  11. 信号を再構築するための方法であって、
    請求項1~9のいずれかに記載の方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、前記実際にサンプリングされた信号のための最終フィッティングモデルを決定するステップと、
    決定された最終フィッティングモデルを利用することによって、前記実際にサンプリングされた信号を再構築するステップと
    を含むことを特徴とする、前記方法。
  12. 前記参照信号と前記実際にサンプリングされた信号を比較するステップが、
    前記実際にサンプリングされた信号の特性情報を利用することによって、前記実際にサンプリングされた信号と同一又は類似の特性情報を有する参照信号を検索するステップと、
    前記参照信号が検索された後、前記参照信号と最終フィッティングモデルとの間の対応関係に従って前記参照信号のための最終フィッティングモデルを検索し、前記参照信号のための検索された最終フィッティングモデルを、前記実際にサンプリングされた信号のための最終フィッティングモデルとして決定するステップと
    を含むことを特徴とする、請求項11に記載の方法。
  13. 信号フィッティングモデルを決定するための装置であって、
    サンプリングされた信号を前記サンプリングされた信号の特性情報に基づいてセグメント化して複数の信号セグメントを得るように構成されたセグメント化ユニットと、
    複数のフィッティングモデルを用いて前記複数の信号セグメントにおいてサンプリングポイントをフィッティングするように構成されたフィッティングユニットと、
    フィッティング結果に基づいて各信号セグメントにおいて対象パラメータについてのフィッティング値を取得するように構成された取得ユニットと、
    前記対象パラメータについての前記フィッティング値のそれぞれを、前記フィッティングモデルのそれぞれの下で、前記対象パラメータについての取得された測定値と順次比較し、比較結果に基づいて、前記複数のフィッティングモデルの中で、前記信号の再構築のための最終フィッティングモデルを決定するように構成された比較及び決定ユニットと、
    を含むことを特徴とする、前記装置。
  14. 前記比較及び決定ユニットが、前記フィッティングモデルのそれぞれの下での信号セグメントのそれぞれにおける対象パラメータについてのフィッティング値のそれぞれと、前記対象パラメータについての取得された測定値と、の間の差値を順次計算し、前記差値に基づいて、前記複数のフィッティングモデルの中で、前記信号セグメントのそれぞれを再構築するための最終フィッティングモデルを決定するように具体的に構成されていることを特徴とする、請求項13に記載の装置。
  15. 前記信号セグメントのそれぞれと前記最終フィッティングモデルとの間の対応関係に従って、決定された最終フィッティングモデルを記憶するように構成された記憶ユニットをさらに含むことを特徴とする、請求項13に記載の装置。
  16. 信号を再構築するための装置であって、
    請求項1~9のいずれかに記載の方法を利用することによって、実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成された決定ユニットと、
    決定された最終フィッティングモデルを利用することによって、前記実際にサンプリングされた信号を再構築するように構成された再構築ユニットと
    を含むことを特徴とする、前記装置。
  17. 信号を再構築するための装置であって、
    請求項1~9のいずれかに記載の方法を利用することによって最終フィッティングモデルが決定される参照信号と、実際にサンプリングされた信号を比較して、前記実際にサンプリングされた信号のための最終フィッティングモデルを決定するように構成された比較及び決定ユニットと、
    決定された最終フィッティングモデルを利用することによって、前記実際にサンプリングされた信号を再構築するように構成された再構築ユニットと、
    を含むことを特徴とする、前記装置。
  18. 前記比較及び決定ユニットが、実際にサンプリングされた信号の特性情報を利用することによって、前記実際にサンプリングされた信号と同一又は類似の特性情報を有する参照信号を検索して、前記参照信号が検索された後、前記参照信号と最終フィッティングモデルとの間の対応関係に従って前記参照信号のための最終フィッティングモデルを検索し、前記参照信号のための検索された最終フィッティングモデルを、前記実際にサンプリングされた信号のための最終フィッティングモデルとして決定するように具体的に構成されていることを特徴とする、請求項17に記載の装置。
JP2020570150A 2018-08-27 2019-05-09 信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置 Active JP7049709B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201810983248.2 2018-08-27
CN201810983248.2A CN109350098B (zh) 2018-08-27 2018-08-27 信号的拟合方式的确定方法、重建方法和装置
PCT/CN2019/086237 WO2020042663A1 (zh) 2018-08-27 2019-05-09 信号的拟合方式的确定方法、重建方法和装置

Publications (2)

Publication Number Publication Date
JP2021535454A JP2021535454A (ja) 2021-12-16
JP7049709B2 true JP7049709B2 (ja) 2022-04-07

Family

ID=65349929

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020570150A Active JP7049709B2 (ja) 2018-08-27 2019-05-09 信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置

Country Status (5)

Country Link
US (1) US11842426B2 (ja)
EP (1) EP3795085A4 (ja)
JP (1) JP7049709B2 (ja)
CN (1) CN109350098B (ja)
WO (1) WO2020042663A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109350098B (zh) * 2018-08-27 2021-02-26 苏州瑞派宁科技有限公司 信号的拟合方式的确定方法、重建方法和装置
CN110368020B (zh) * 2019-07-03 2022-03-08 北京必安必恒科技发展有限公司 一种心音信号预处理方法及装置
CN114325524B (zh) * 2020-09-29 2023-09-01 上海联影医疗科技股份有限公司 磁共振图像重建方法、装置、系统及存储介质
CN112965120B (zh) * 2021-02-01 2024-04-12 西双版纳景海高速公路建设投资有限公司 一种地质勘探用瞬变电磁信号处理方法、装置及存储介质
CN113032732B (zh) * 2021-03-24 2023-01-24 中国兵器装备集团自动化研究所有限公司 一种基于相对误差分段的剂量率拟合方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010131417A1 (ja) 2009-05-14 2010-11-18 日本電気株式会社 データ要約システム、データ要約方法、及びデータ要約プログラム
JP2014102779A (ja) 2012-11-22 2014-06-05 Fujitsu Ltd データ補間装置、データ補間プログラム及びデータ補間方法
JP2014529723A (ja) 2011-06-14 2014-11-13 サザン イノヴェーション インターナショナル プロプライアトリー リミテッド 検出器出力データ中のパルスを特定するための方法および装置
WO2016111240A1 (ja) 2015-01-07 2016-07-14 日本電気株式会社 情報処理システム、変化点検出方法、および記録媒体
US20160216246A1 (en) 2012-06-05 2016-07-28 The University Of Utah Research Foundation Reduced parameter space kinetic modeling systems and methods

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08221388A (ja) * 1995-02-09 1996-08-30 Nec Corp フィッティングパラメータ決定方法
CN102262238B (zh) * 2011-04-19 2014-07-23 苏州瑞派宁科技有限公司 一种提取闪烁脉冲信息的方法及装置
GB2509168B (en) * 2012-12-21 2020-04-01 James Shue Min Yeh Method of quantifying the characteristics of an object treated with a contrast agent
CN103969675B (zh) * 2013-02-05 2017-08-04 苏州瑞派宁科技有限公司 数字化闪烁脉冲的基线校正方法及系统
PL227660B1 (pl) * 2013-08-30 2018-01-31 Uniwersytet Jagiellonski Sposób wyznaczania parametrów reakcji kwantu gamma w scyntylatorze tomografu PET
CN104146726B (zh) * 2014-08-19 2016-09-14 中国科学院高能物理研究所 Pet系统符合探测响应的生成方法
CN105824817B (zh) * 2015-01-05 2019-05-21 苏州瑞派宁科技有限公司 一种闪烁脉冲的数字化方法
EP3098780A1 (en) 2015-05-28 2016-11-30 Galgo Medical, S.L. A method, a system and a computer program for determining data defining an estimate of the thickness and density of a tissue structure of interest from imaging data
CN105588667B (zh) 2015-12-18 2018-07-10 上海市计量测试技术研究院 一种高精密热敏电阻温度计校准装置
CN108021535A (zh) * 2016-11-03 2018-05-11 北京国双科技有限公司 数据拟合方法和装置
CN107024711B (zh) * 2017-04-17 2019-02-26 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化信号的拟合方法
CN109350098B (zh) * 2018-08-27 2021-02-26 苏州瑞派宁科技有限公司 信号的拟合方式的确定方法、重建方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010131417A1 (ja) 2009-05-14 2010-11-18 日本電気株式会社 データ要約システム、データ要約方法、及びデータ要約プログラム
JP2014529723A (ja) 2011-06-14 2014-11-13 サザン イノヴェーション インターナショナル プロプライアトリー リミテッド 検出器出力データ中のパルスを特定するための方法および装置
US20160216246A1 (en) 2012-06-05 2016-07-28 The University Of Utah Research Foundation Reduced parameter space kinetic modeling systems and methods
JP2014102779A (ja) 2012-11-22 2014-06-05 Fujitsu Ltd データ補間装置、データ補間プログラム及びデータ補間方法
WO2016111240A1 (ja) 2015-01-07 2016-07-14 日本電気株式会社 情報処理システム、変化点検出方法、および記録媒体

Also Published As

Publication number Publication date
CN109350098A (zh) 2019-02-19
EP3795085A4 (en) 2022-03-09
CN109350098B (zh) 2021-02-26
US11842426B2 (en) 2023-12-12
JP2021535454A (ja) 2021-12-16
WO2020042663A1 (zh) 2020-03-05
US20210275113A1 (en) 2021-09-09
EP3795085A1 (en) 2021-03-24

Similar Documents

Publication Publication Date Title
JP7049709B2 (ja) 信号のためのフィッティングモデルを決定し、信号を再構築するための方法及びその装置
Raczyński et al. Novel method for hit-position reconstruction using voltage signals in plastic scintillators and its application to Positron Emission Tomography
US11944470B2 (en) Method and device for sampling a pulse signal, and computer program medium
US20130334429A1 (en) Imaging device using gamma rays, image signal processor, and image processing method for gamma ray measurement data
US9734600B2 (en) Attenuation map with scattered coincidences in Positron Emission Tomography
Gaitanis et al. PET image reconstruction: A stopping rule for the MLEM algorithm based on properties of the updating coefficients
CN106963407B (zh) 正电子发射断层成像系统及其图像重建方法
BR112013026944B1 (pt) Método e dispositivo para extrair informação por impulso de cintilação
US20130136328A1 (en) Methods and systems for enhanced tomographic imaging
CN105824817A (zh) 一种闪烁脉冲的数字化方法
JP2014149308A (ja) 生理学的リストモード核医学イメージングにおけるカウントの最適時間サンプリングのための逆データ再構成
CN1969295B (zh) 数据处理系统以及使用该系统的检查设备
EP3427667A1 (en) Time calibration in pet device
WO2018207739A1 (ja) 陽電子寿命測定機能付きpet装置、及び、pet装置における陽電子寿命測定方法
Cooper et al. SmartPET: Applying HPGe and pulse shape analysis to small-animal PET
Brambilla et al. Comparative analysis of iterative reconstruction algorithms with resolution recovery and new solid state cameras dedicated to myocardial perfusion imaging
JP2023517792A (ja) 信号サンプリング方法及び再構成方法、並びに装置
Lyra et al. Quantification of Myocardial Perfusion in 3D SPECT images-Stress/Rest volume differences
CN110226943B (zh) 光子到达探测器的参数计算方法、装置和计算机设备
Le Meunier et al. Validation of a PET Monte-Carlo simulator with random events and dead time modeling
CN105212956A (zh) 一种基于ist的次晶体级pet系统时间修正方法
Jiang et al. Automatic detection of scintillation light splashes using conventional and deep learning methods
Tayefi Ardebili et al. Evaluation of Modular J-PET sensitivity
CN117909722A (zh) 模型训练方法、光子检测方法、终端设备以及存储介质
CN108648807B (zh) 一种图像重建方法和装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210309

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210309

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: 20220307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220318

R150 Certificate of patent or registration of utility model

Ref document number: 7049709

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150