JP2024054856A - Ｐｅｔ装置、データ処理方法及びプログラム - Google Patents

Ｐｅｔ装置、データ処理方法及びプログラム Download PDF

Info

Publication number
JP2024054856A
JP2024054856A JP2023173213A JP2023173213A JP2024054856A JP 2024054856 A JP2024054856 A JP 2024054856A JP 2023173213 A JP2023173213 A JP 2023173213A JP 2023173213 A JP2023173213 A JP 2023173213A JP 2024054856 A JP2024054856 A JP 2024054856A
Authority
JP
Japan
Prior art keywords
loss correction
detector
coincidence
correction factor
count
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
Application number
JP2023173213A
Other languages
English (en)
Inventor
リ シャオリ
チャン イ
コルサマー ジェフリー
正樹 宮原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to EP23201886.1A priority Critical patent/EP4350396A1/en
Publication of JP2024054856A publication Critical patent/JP2024054856A/ja
Pending legal-status Critical Current

Links

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/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/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • 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/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647—Processing of scintigraphic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine (AREA)
  • Measurement Of Radiation (AREA)

Abstract

【課題】画質を向上させること。【解決手段】実施形態に係るＰＥＴ装置は、複数の検出器サブシステムと、処理回路とを備える。処理回路は、前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する。【選択図】図１

Description

本明細書及び図面に開示の実施形態は、ＰＥＴ装置、データ処理方法及びプログラムに関する。
ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置では、再構成画像のボクセル値は、放射能濃度の絶対単位で、適当な精度および正確さで較正される場合がある。ＰＥＴ装置は、複数のサブシステムおよび／または検出器ユニットを備え、そのそれぞれには、イベントの処理および／または一定量のデータの読取り、書込み、および転送のための時間が最小限となるのが望ましい。広範囲のファントムおよび線源の放射能全体にわたる定量的なＰＥＴ画像を生成するためには、計数率に依存し、かつサブシステムによって変化する様々な線源から起こる場合がある計数損失を、特性決定し、補正する必要がある。
米国特許出願公開第２０１４／０２３３６９３号明細書 米国特許出願公開第２０１４／０２３３６９４号明細書
本明細書及び図面の開示の実施形態が解決しようとする課題の一つは、画質を向上することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態に係るＰＥＴ装置は、複数の検出器サブシステムと、処理回路とを備える。処理回路は、前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する。
図１は、例示的ＰＥＴスキャナシステム内の検出器ユニットおよび様々なデータ処理ユニットを示す図である。 図２は、図１に示したスキャナシステムの同時計数損失補正係数の例示的計算を示す図である。 図３は、本開示による方法のフローチャートの一例を示す図である。 図４は、本開示の一実施形態によるポジトロン放射断層撮影（ＰＥＴ）スキャナの斜視図である。 図５は、本開示の一実施形態によるＰＥＴ装置および関連するハードウェアの概略図である。
以下、図面を参照しながら、ＰＥＴ装置、データ処理方法及びプログラムの実施形態について詳細に説明する。
従来の手法では、計数損失補正は、減衰線源測定から、またはシミュレーションから導出され得る、ＰＥＴスキャナシステム全体の計数損失補正（Count Loss Correction：ＣＬＣ）係数のルックアップテーブルを作成することによって実施される場合がある。しかし、この単純化された手法は、スキャナ内部の異なるサブシステムまたは検出器ユニットでの相対計数損失を変化させる可能性がある線源分布の空間的なばらつきは考慮されていない。
計数損失補正は、各サブシステムまたは各検出器ユニットの計数損失挙動をモデル化し、次いで、それらを一緒に組み合わせて、ＰＥＴスキャナシステム全体の計数損失挙動をモデル化することによって実施されることもある。しかし、ＰＥＴ装置における計数損失の複雑な性質が原因で、各サブシステムまたは各検出器ユニットの計数損失挙動を組み合わせる方法は、多くの課題がある。
特に、ＰＥＴ画像再構成には、同時計数検出が必要であり、同時計数損失補正係数は、通常、２つのシングル計数損失補正係数の積として計算される。しかし、この手法は、同時計数イベント時の２つのシングルイベントの計数損失が相関する場合、一方のシングルイベントが失われると、他方のシングルイベントも失われる可能性が高くなることを意味するため、正確さに欠ける。
さらに、ＰＥＴ装置では、異なる検出器ユニットが、異なる計数率を受信する可能性があり、したがって、異なる計数損失が発生する場合がある。したがって、システム全体の計数率に基づいてシステム計数損失補正係数を計算することは、正確でない。
加えて、82Rbスキャンでは数分といった、同位元素半減期と比べて長いスキャン時間間隔でのＰＥＴスキャンの場合、計数率は、スキャン中に著しく変化し、ゆえに、計数損失も、スキャン中に著しく変化する。したがって、平均計数率に基づいて総計数損失補正係数を計算することは、正確でない。
上述のように、再構成中に、補正係数が測定された計数率に適用されて、例えば、下記式（１）を使用して、実際の検出計数率が推定される。
Figure 2024054856000002
式中Csingle(fmeasured)は、シングル計数損失補正（ＣＬＣ）係数であり、また、通常、１より大きい値である。一般に、Csingle(fmeasured)の明確な解析形式は存在しない。そのため、実際には、近似的な解析形式またはルックアップテーブルのいずれかが代わりに使用される。
ＰＥＴ再構成では、光子イベントのペアが使用されるため、ペアの計数損失補正（ＣＬＣ）係数が計算される必要がある。ほとんどのシステムでは、シングルａおよびｂを含むイベントのペアのＣＬＣ係数は、下記式（２）によって得られる。
Figure 2024054856000003
式中、ｆａおよびｆbは、各検出器の計数率である。ペアのＣＬＣ係数は、再構成プロセスのイベントごとに推定されて、適用されるか、平均ペアＣＬＣ係数が、再構成画像にスケーリング係数として計算され、適用され得る。
全体的なペアのＣＬＣ係数は、システムシングル計数率の関数としてさらに下記式（３）のように単純化される。
Figure 2024054856000004
ここで、一実施形態では、１つの同時計数イベントでの２つのシングルイベントの計数損失が相関する場合、同時計数損失補正係数が、２つのシングル計数損失補正係数の最大値として計算される。
さらに、一実施形態では、異なる検出器ユニットで異なる計数損失が発生する場合、異なる検出器ユニット間の相対計数損失補正係数が、画像再構成に組み込まれて、より均一な画像が得られ、システム計数損失補正係数が、対応する検出器ユニットペアの同時計数率によって重み付けされた検出器ユニットペアの計数損失補正係数として計算される。
さらに、一実施形態では、計数率がスキャン中に著しく変化する場合、スキャナ全体の総計数損失補正係数が、対応する即時同時計数率によって重み付けされた即時計数損失補正係数として計算される。
図１に示すように、例示的ＰＥＴシステムでは、４つの象限データ処理ユニット（例えば、データ処理ユニット０からデータ処理ユニット３）のそれぞれは、いくつかの関連する検出器ユニットからデータを収集して、それらを処理する。さらに、データ処理ユニット５は、データ処理ユニット０から３のそれぞれからデータを収集して、それらを処理し、かつさらなる処理および制御が、データ処理ユニット５に接続された処理回路によって実施される。検出器ユニットと、データ処理ユニットと、処理回路との間のデータ転送は、例えば、光ファイバーケーブルまたはイーサネットを通して実施される。
（シングル計数損失補正係数）
サブシステムごとまたは検出器ユニットごとのシングル計数損失補正係数は、減衰線源測定から、またはシミュレーションを介して計算されるか、および／または各サブシステムの信号処理速度制限に基づく場合がある。
例えば、各検出器ユニットで、偶発的な不感時間が発生する場合がある。検出器ユニットごとのシングル計数損失補正係数は、検出器リングの中央に置かれた減衰ファントムから計算することができる。特に、減衰ファントムが中央に置かれている場合、各検出器ユニットは、類似したシングル計数率を受信し、ゆえに類似したシングル計数損失が発生する。さらに、同時計数損失補正係数は、システムシングル計数率または平均検出器ユニットシングル計数率から計算することができる。さらに、上記の式（３）に基づいて、検出器ユニットごとのシングル計数損失補正係数は、同時計数損失補正係数の平方根として計算することができる。また、システムペア率は、時間の関数として散乱フラクションが利用可能な場合、より正確な結果を得るために、（プロンプト率－偶発率）×（１－散乱フラクション）と計算することができる。
さらに、検出器サブシステムにイベント信号処理速度制限がある場合、サブシステムのシングル計数損失補正係数は、入力計数率をイベント信号処理速度閾値で除算することによって計算することができる。一実施形態では、異なる検出器ユニット間の相対計数損失補正係数が、画像再構成に組み込まれて、より均一な画像が得られる。
（同時計数損失補正係数）
サブシステムまたは検出器ユニットのペアごとの同時計数損失補正係数は、ペアのサブシステムごとまたはペアの検出器ユニットごとのシングル計数損失補正係数から計算することができる。
非相関計数損失の場合、例えば、偶発的な不感時間の場合、２つのシングルイベントの計数損失は、互いに独立している。同時計数損失補正係数は、同時計数ペアの２つのシングル計数損失係数の積である。
相関計数損失の場合、例えば、限られたバッファ深さによるオーバフローの場合、より高いシングル計数率でサブシステムに入射するシングルイベントが失われると、他のイベントも失われる。したがって、同時計数損失補正係数は、ペアの２つのシングル計数損失補正係数の最大値である。
非相関計数損失と相関計数損失との組み合わせの場合、非相関計数損失および相関計数損失は、組み合わせた同時計数損失補正係数を適切に計算するために、別々に決定される必要がある。例えば、計数損失の根本原因が分かっている場合、単純なモデルを使用して、非相関計数損失と、相関計数損失とを切り離すことができる。計数損失の根本原因が分かっていない場合、より高度な方法を使用して、非相関計数損失と、相関計数損失とを切り離す。例えば、以下に詳述するように、ニューラルネットワークを異なるファントムスキャンを用いて訓練して、特定のスキャンの非相関計数損失と相関計数損失とを切り離すことができる。
計数損失の根本原因が分かっている場合、一実施形態では、非相関計数損失と相関計数損失との組み合わせの同時計数損失補正係数は、以下のように計算される。まず、中央に置かれた減衰ファントムを使用して、非相関計数損失に基づくサブシステムごとのシングル計数損失補正係数を事前較正し、その結果を、シングル計数損失補正テーブルとして保存する。偶発的な不感時間に起因するような非相関計数損失は、通常、相関計数損失よりも低い計数率で発生することに留意されたい。したがって、アクティビティは、相関計数損失のみが発生するように制御することができる。
次に、相関計数損失に基づく検出器サブシステムごとのシングル計数損失補正係数は、イベント処理速度閾値を使用して計算することができる。ＰＥＴ装置がＮ個の検出器サブシステムを有し、かつ各サブシステムｉが、Ti（計数率の単位）のイベント処理速度制限値を有し、さらにサブシステムｉの測定された入力シングル計数率が、SRiであると仮定する。さらに、サブシステムペアｉ－ｊごとの同時発生率が、CRijであると仮定する。サブシステムｉの非相関シングル計数損失（Singles Count Loss：ＳＣＬ）補正係数であるSCLuncor_jは、入力シングル計数率SRiおよび事前較正されたシングル計数損失補正テーブルに基づいて計算することができる。さらに、サブシステムｉの相関シングル計数損失補正係数は、SCLcorr_j = SRi/Tiと計算することができる。
さらに、この実施形態では、総シングル計数損失補正係数サブシステムｉは、SCLi= SCLuncorr_i x SCLcorr_iのように計算することができる。すなわち、図５に示されるプロセッサ５０７等の処理回路は、検出器サブシステムｉごとに、検出器サブシステムｉで受信された測定済みのシングル計数率SRiに基づいて、シングル計数損失補正係数SCLiを決定するさらに、収集ごとのサブシステムペアｉ－ｊの同時計数損失（Coincidence Count Loss：ＣＣＬ）補正係数は、以下の式（４）のように計算することができる。
Figure 2024054856000005
次に、サブシステムｉごとの相対シングル計数損失補正係数は、SCLrel_i = SCLi / mean(SCLi)のように計算することができ、かつ画像再構成中に、該相対シングル計数損失補正係数に結晶効率テーブルからの値を乗じて、より均一な画像を得ることができる。スキャナの同時計数損失補正係数は、下記式（５）のように計算することができる。
Figure 2024054856000006
スキャナの同時計数損失補正係数を再構成画像に適用されるスケーリング係数として使用して、より良好な画像定量基準値、例えば基準更新値（Standard Update Value：ＳＵＶ）を得ることができることに留意されたい。
計数損失の根本原因が分かっていない場合、一実施形態では、非相関計数損失と相関計数損失との組み合わせの同時計数損失補正係数は、以下のように計算することができる。まず、ファントムスキャンが数回実施される。例えば、各ファントムのサイズが異なる場合、散乱イベントの量も異なる。さらに、ファントムが中央に置かれた場合と中央外に置かれた場合の両方のアクティビティが存在する場合、異なるサブシステム間で対称的な計数率と非対称的な計数率の両方が得られる。各ファントムスキャンは、様々なアクティビティをカバーするように一連の収集である必要がある。加えて、最も低いアクティビティは、計数損失がほとんど無いことに相当するはずである。ファントムスキャンごとに、以下の計数率およびアクティビティが、収集ごとに記録される必要がある。
収集ｋごとのアクティビティ：Activityk、収集ｋでの、サブシステムｉの入力シングル計数率：SRinput_i_k、収集ｋでの、サブシステムｉの出力シングル計数率：SRoutput_i_k及び収集ｋでの、サブシステムペアｉ－ｊの同時計数：CRij_k。
ファントムスキャンごとに、収集ごとの計数損失補正係数は、記録された計数率およびアクティビティから計算することができる。特に、一実施形態では、最も低いアクティビティでは計数損失が無い（Activity0）ことが想定される。収集ｋでサブシステムｉの計数損失が無い場合の予測シングル計数率は、記録されたアクティビティおよび計数率を用いて、以下の式（６）のように計算することができる。
Figure 2024054856000007
また、収集ｋでのサブシステムのシングル計数損失補正係数は、以下の式（７）のように計算することができる。
Figure 2024054856000008
収集ごとのサブシステムペアごとの計数損失が無い場合の予測同時計数率は、以下の式（８）のように計算することができる。収集ごとのサブシステムペアごとの同時計数損失補正係数は、以下の式（９）のように計算することができる。
Figure 2024054856000009
Figure 2024054856000010
ファントムスキャンごとに、収集ごとの計数損失補正係数は、非相関および相関シングル計数損失補正係数を使用して表すことができる。例えば、収集ｋでのサブシステムｉの非相関シングル計数損失係数補正は、SCLuncorr_j _k、かつ収集ｋでのサブシステムｉの相関シングル計数損失補正係数は、SCLcorr_j _kであると仮定する。次に、収集ｋでのサブシステムｉのシングル計数損失補正係数は、以下の式（１０）のように計算することができる。さらに、収集ｋでのサブシステムペアｉ－ｊの同時計数損失補正係数は、以下の式（１１）のように計算することができる。
Figure 2024054856000011
Figure 2024054856000012
さらに、一実施形態では、ニューラルネットワークが、異なるファントムスキャンを使用して、特定のスキャンの非相関および相関計数損失を推定するように訓練される。一例示的実装形態では、入力層は、（１）全てのファントムスキャンの記録された計数率およびアクティビティから推定される、収集ｋでのサブシステムｉのシングル計数損失補正係数（SCLi_k）、および／または（２）全てのファントムスキャンの記録された計数率およびアクティビティから推定される、収集ｋごとのサブシステムペアｉ－ｊの同時計数損失補正係数（CCLij_k）の入力を含む。この実施形態では、ニューラルネットワークの出力層は、（１）収集ｋでのサブシステムｉの非相関シングル計数損失補正係数（SCLuncorr_j _k）、および（２）収集ｋでのサブシステムｉの相関シングル計数損失補正係数（SCLcorr_j _k）の出力を含む。次に、曲線フィッティングまたは補間を使用して、検出器サブシステムｉごとのSCLuncorr_i_k対SRinput_i_kおよびSCLcorr_i_k対SRinput_i_kのプロットから、異なる入力シングル計数率でのSCLuncorr_iおよびSCLcorr_iを計算することができる。
したがって、一実施形態では、新しい患者またはファントムスキャンの場合、スキャナのサブシステムｉごとに、非相関および相関シングル計数損失補正係数（SCLuncor_jおよびSCLcorr_j）は、曲線フィッティングを用いて、前述したように、入力シングル計数率に基づいて計算することができる。さらに、サブシステムｉごとのシングル計数損失補正係数（SCLi= SCLuncorr_i x SCLcorr_i）およびサブシステムペアｉ－ｊごとの同時計数損失補正係数（CCLij= SCLuncorr_i x SCLuncorr_j x max(SCLcorr_i,SCLcorr_j)）は、上述したように計算することができる。最後に、サブシステムｉごとの相対シングル計数損失補正係数（SCLrel_i = SCLi / mean(SCLi)）およびスキャナの同時計数損失補正係数は、上述したように計算することができる。上述のように、画像再構成中に、相対シングル計数損失補正係数に結晶効率テーブルからの値を乗じて、より均一な画像を得ることができ、一方で、スキャナの同時計数損失補正係数を再構成画像に適用されるスケーリング係数として使用して、より良好な画像定量基準値、例えば基準更新値（ＳＵＶ）を得ることができる。
一実施形態では、スキャナシステム全体の同時計数損失補正係数は、重み付け係数としてサブシステムまたは検出器ユニットのペアごとのペア率を使用して、サブシステムまたは検出器ユニットのペアごとの同時計数損失補正係数から計算される。この実施形態では、サブシステムまたは検出器ユニットのペアごとのペア率は、例えば、プロンプト率から偶発率を減じて計算される。
一実施形態では、データフレーム全体の同時計数損失補正係数は、重み付け係数として即時ペア率を使用して、即時同時計数損失補正係数から計算される。一実施形態では、即時ペア率を計算する時間間隔は、放射性同位元素の半減期に基づいて決定される。この実施形態では、ペア率は、例えば、プロンプト率から偶発率を減じて計算される。
図２は、本開示の一実施形態による、図１に示したスキャナシステム全体の同時計数損失補正係数の例示的計算を示す。この例では、ＰＥＴ装置は、４つの象限に分割され、象限のそれぞれは１００Ｍｃｐｓのイベント処理速度制限を有する。したがって、この例では、象限のペアに対する入力シングル計数率が１００Ｍｃｐｓを超える場合に、相関計数損失が発生し始めることになる。
この例では、象限ごとの測定された入力シングル計数率は以下の通りである。第１象限：１６０Ｍｃｐｓ、第２象限：１２０Ｍｃｐｓ、第３象限：８０Ｍｃｐｓ、第４象限：１２０Ｍｃｐｓ。
さらに、象限ペアごとの測定された同時計数率は以下の通りである。象限１－２：４０Ｍｃｐｓ、象限１－３：８０Ｍｃｐｓ、象限１－４：４０Ｍｃｐｓ、象限２－４：８０Ｍｃｐｓ。
この例では、象限ごとのシングル計数損失補正係数は、上述したように決定される。この例では、少なくともイベント処理速度制限に基づいて、値は、象限１：１．６、象限２：１．２、象限３：１．０、象限４：１．２である。
さらに、この例では、画像再構成向けの象限ｉの相対計数損失補正係数は、全ての象限にわたって平均シングル計数損失補正係数によって分割された象限Ｉのシングル計数損失補正係数として計算することができる。この例では、値は、象限１：１．２８、象限２：０．９６、象限３：０．８、象限４：０．９６、と計算される。
さらに、この例では、相関計数損失の象限ペアごとの同時計数損失補正係数は、上述したように決定される。この例では、値は、象限１－２：１．６、象限１－３：１．６、象限１－４：１．６、象限２－４：１．２、と決定される。
最後に、スキャナ全体の同時計数損失補正は、重み付け係数として象限ペアごとの同時計数率を使用して、以下の式（１２）のように、加重平均として計算される。すなわち、処理回路は、スキャナ同時計数損失補正係数を、複数のペアの測定された同時計数率を使用して同時計数損失補正係数の重み付けされた組み合わせとして計算する。
Figure 2024054856000013
一般的な例として、ＰＥＴ装置はＮ個の検出器サブシステムを有し、かつ各システムがＴ（計数率の単位）のイベント処理速度制限値を有すると仮定する。したがって、この第２の例では、相関計数損失は、サブシステムへの入力シングル計数率が値Ｔを超過した場合に発生する。この例では、サブシステムｉごとの測定された入力シングル計数率がSRiであり、サブシステムペアｉ－ｊごとの測定された同時計数率がCRijであると想定されている。したがって、各サブシステムのイベント処理速度制限値に基づいて、サブシステムｉのシングル計数損失補正係数は、SCLcorr_i = SRi / Tである。
さらに、この例では、画像再構成向けのサブシステムｉの相対計数損失補正係数は、SCLrel_corr_i = SCLcorr_i / mean(SCLcorr_i)と計算することができ、ここで、mean(SCLcorr_i)は、サブシステム全体にわたる平均シングル計数損失補正係数である。
さらに、この例では、相関計数損失のサブシステムペアｉ－ｊごとの同時計数損失補正係数は、CCLij = max(SCLcorr_i, SCLcorr_j)のように計算される。
さらに、この例では、スキャナ全体の同時計数損失補正係数は、以下の式（１３）によって得られる加重平均である。
Figure 2024054856000014
本開示の発明は、既存の手法に優るいくつかの利点がある。特に、本明細書にて開示する本発明は、特に１つの同時計数イベントペアで２つのシングルイベントの計数損失が相関する場合の、実用的、かつより正確な計数損失補正の方法を提供する。
図３は、本開示の一実施形態による複数の検出器サブシステムを有するポジトロン放射断層撮影（ＰＥＴ）スキャナの計数損失補正の方法のフローチャートを示す。
ステップＳ３０１で、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信したシングル計数率が測定される。さらに、検出器サブシステムの複数のペアの検出器サブシステムペアごとに、検出器サブシステムペアの同時計数率が測定される。
ステップＳ３０２で、検出器サブシステムの複数のペアの検出器サブシステムペアごとに、相関シングル計数損失補正係数および非相関シングル計数損失補正係数は、上述したように、訓練済みニューラルネットワークを使用して生成された推定曲線を使用して、測定されたシングル計数率に基づいて予測される。
ステップＳ３０３で、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムのシングル計数損失補正係数は、上述したように、相関および非相関シングル計数損失補正係数に基づいて決定される。
次に、ステップＳ３０４で、検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数ＣＬ補正係数が、上述したように、相関および非相関シングル計数損失補正係数に基づいて、検出器サブシステムペアに対して決定される。
ステップＳ３０５で、複数の検出器サブシステムの検出器サブシステムごとに、相対シングル計数損失補正係数が、シングル計数損失補正係数に基づいて決定される。
ステップＳ３０６で、ＰＥＴ装置のスキャナ同時計数損失補正係数が、検出器サブシステムの複数のペアに対して決定された同時計数損失補正係数に基づいて計算される。
ステップＳ３０７で、ＰＥＴ画像は、計算されたスキャナ計数損失補正係数およびＰＥＴ装置を使用して実施された対象のスキャンから収集されたスキャンデータに基づいて再構成される。各サブシステムの相対シングル計数損失補正係数は、上述したように再構成で使用することができる。
本明細書にて開示する実施形態で使用可能なＰＥＴ装置を図４および５に示す。ＰＥＴ装置４００は、矩形の検出器モジュールとしてそれぞれが構成される、複数のガンマ線検出器（Gamma-Ray Detector：ＧＲＤ）４０１（例えば、ＧＲＤ１、ＧＲＤ２からＧＲＤＮ）を備える。一実装形態によれば、検出器は、リングに配置され、それにより、ガントリ４０４の周りに円形ボア４０２を形成する。この例では、リングは、４０個のＧＲＤ４０１を含んでいる。リングは、円形ボア４０２の所望のサイズなどの因子に応じて、異なる数のＧＲＤ４０１を有する場合がある。ＧＲＤ４０１は、ガンマ線を（例えば、可視の、赤外のおよび紫外の波長での）シンチレーション光子に変換するためのシンチレータ結晶アレイを含み、それらシンチレーション光子は、光検出器によって検出される。各ＧＲＤ４０１は、ガンマ線を吸収して、シンチレーション光子を放射する個々の検出器結晶の２次元アレイを備えるか、またはモノリシックアレイもしくはスラットアレイを備える場合がある。シンチレーション光子は、ＧＲＤ４０１に配設されているＳｉＰＭなどのデバイス（図示せず）の２次元アレイによっても検出される場合がある。ライトガイドが、検出器結晶アレイとＳｉＰＭとの間に配置可能である。本開示による結晶およびＳｉＰＭ配置については、下記にてより詳細に論じる。
図４は、対象ＯＢＪから放射されるガンマ線を検出するように配設されたＧＲＤを有するＰＥＴシステムの概略図を示す。ＧＲＤは、検出されたガンマ線のそれぞれに対応するタイミング、位置、およびエネルギーを測定することができる。図４の単一のＰＥＴ検出器リングは、ＰＥＴ装置の軸方向長さに沿って任意の数のＰＥＴ検出器リングを備えるように挿入される場合があることが理解できる。
図５は、撮像される対象ＯＢＪが台５０６上に置かれ、ＧＲＤ１からＧＲＤＮのＧＲＤモジュールが、対象ＯＢＪおよび台５０６の周りに周囲方向に配設される、ＰＥＴ装置５００の構造の例を示す。ＧＲＤは、ＰＥＴ検出器リングを備える場合があり、またガントリ５０４に固定して接続されている円形ボア５０２に固定して接続される場合がある。ガントリ５０４は、ＰＥＴ装置の多数の部品を収容する。ＰＥＴ装置のガントリ５０４はまた、対象ＯＢＪおよび台５０６がそれを通って通過できる、円形ボア５０２によって画定された開口を含み、また、消滅イベントにより対象ＯＢＪから反対方向に放射されるガンマ線は、ＧＲＤにより検出可能であり、かつ、タイミング情報およびエネルギー情報は、ガンマ線ペアの同時計数を決定するために使用される場合がある。
図５では、ガンマ線検出データを収集、記憶、処理、および分配する回路およびハードウェアも示されている。これらの回路およびハードウェアは、プロセッサ（処理回路）５０７、ネットワークコントローラ５０３、メモリ５０５、およびデータ収集システム（Data Acquisition System：ＤＡＳ）５０８、を備える。ＰＥＴ撮像装置はまた、ＧＲＤからの検出測定結果をデータ収集システム５０８、プロセッサ５０７、メモリ５０５、およびネットワークコントローラ５０３に送出するデータチャンネルを備える。データ収集システム５０８は、検出器からの検出データの収集、デジタル化、および送出を制御することができる。一実装形態では、データ収集システム５０８は、台５０６の動きを制御する。プロセッサ５０７は、配置誤差の特定、検出データの再構成前処理、画像再構成、および画像データの再構成後処理を含む機能を実施する。
一実施形態によれば、図５のＰＥＴ装置５００のプロセッサ５０７は、本明細書に記載の方法を実施するように構成される場合がある。プロセッサ５０７は、個別論理ゲート、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ）、または他の複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）として実装されることがあるＣＰＵを含む場合がある。ＦＰＧＡまたはＣＰＬＤの実装形態は、ＶＨＤＬ（登録商標）、Ｖｅｒｉｌｏｇ（登録商標）、または任意のその他のハードウェア記述言語で符号化されてもよく、かつ、その符号は、ＦＰＧＡまたはＣＰＬＤ内部の電子メモリに直接記憶されるか、または別個の電子メモリとして記憶されてよい。さらに、メモリ５０５は、ハードディスクドライブ、ＣＤ－ＲＯＭドライブ、ＤＶＤドライブ、ＦＬＡＳＨ（登録商標）ドライブ、ＲＡＭ、ＲＯＭまたは当該技術分野において既知の任意のその他の電子記憶装置とすることができる。メモリ５０５は、ＲＯＭ、ＥＰＲＯＭ（登録商標）、ＥＥＰＲＯＭ（登録商標）またはＦＬＡＳＨ（登録商標）メモリのような不揮発性のものであってもよい。メモリ５０５は、静的または動的ＲＡＭなどの揮発性メモリとすることができ、電子メモリおよびＦＰＧＡまたはＣＰＬＤとメモリとの間の相互動作を管理するために、マイクロコントローラまたはマイクロプロセッサなどのプロセッサが設けられてもよい。
あるいは、プロセッサ５０７のＣＰＵは、本明細書に記載の方法を実施するコンピュータ可読命令のセットを含むコンピュータプログラムを実行することができ、このプログラムは、上述の非一時的電子メモリおよび／またはハードディスクドライブ、ＣＤ、ＤＶＤ、ＦＬＡＳＨ（登録商標）ドライブまたは任意のその他の既知の記憶媒体のいずれかに記憶される。さらに、このコンピュータ可読命令は、実用アプリケーション、バックグラウンドデーモン、またはオペレーティングシステムのコンポーネント、あるいはそれらの組み合わせとして提供されてもよく、インテル社(登録商標)のＸｅｎｏｎプロセッサ（登録商標）、またはＡＭＤ社（登録商標）のＯｐｔｅｒｏｎプロセッサ（登録商標）などのプロセッサ、ならびに、Ｍｉｃｒｏｓｏｆｔ ＶＩＳＴＡ（登録商標）、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ（登録商標）、ＬＩＮＵＸ（登録商標）、Ａｐｐｌｅ（登録商標）、ＭＡＣ－ＯＳ（登録商標）および当業者に既知の他のオペレーティングシステムなどのオペレーティングシステムと連動して実行される。さらに、ＣＰＵは、命令を実施するために並行して協動する複数のプロセッサとして実装されることがある。命令は、メモリ５０５またはネットワークコントローラ５０３に設けられたメモリ（図示せず）内部に記憶されてよい。
一実装形態では、ＰＥＴ装置は、再構成画像を表示するためのディスプレイなどを含む場合がある。ディスプレイは、ＬＣＤディスプレイ、ＣＲＴディスプレイ、プラズマディスプレイ、ＯＬＥＤ、ＬＥＤまたは当該技術分野において既知の任意の他のディスプレイとすることができる。
インテル社（登録商標）のインテルイーサーネット（登録商標）ＰＲＯネットワークインターフェースカードなどのネットワークコントローラ５０３が、ＰＥＴ撮像装置の種々の部品間をインターフェースすることができる。追加として、ネットワークコントローラ５０３は、外部のネットワークともインターフェースすることができる。理解されるように、この外部ネットワークは、インターネットなどの公共ネットワーク、またはＬＡＮもしくはＷＡＮネットワークなどの私的ネットワーク、あるいはそれらの組み合わせとすることができ、かつ、ＰＳＴＮまたはＩＳＤＮサブネットワークを含むこともできる。この外部ネットワークは、イーサネット（登録商標）ネットワークのような有線方式、またはＥＤＧＥ、３Ｇおよび４Ｇ無線セルラーシステムを含むセルラーネットワークのような無線方式とすることもできる。無線ネットワークは、ＷｉＦｉ（登録商標）、Ｂｌｕｅｔｏｏｔｈ（登録商標）、または既知である通信の任意のその他の無線形態とすることができる。
以上の実施形態に関し、発明の一側面および選択的な特徴として以下の付記を開示する。
（付記１）複数の検出器サブシステムを有するポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）スキャナであって、
複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムのシングル計数損失補正係数を決定し、
検出器サブシステムの複数のペアの検出器サブシステムペアごとに、検出器サブシステムペアの同時計数損失補正係数を決定し、
検出器サブシステムの複数のペアの決定された同時計数損失補正係数に基づいて、ＰＥＴ装置のスキャナ同時計数損失補正係数を計算し、
計算されたスキャナ同時計数損失補正係数およびＰＥＴ装置を使用して実施された患者のスキャンから収集されたスキャンデータに基づいて画像を再構成する、ように構成された処理回路を備える、ＰＥＴ装置。
（付記２）処理回路は、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムの処理速度閾値に基づいて、検出器サブシステムのシングル計数損失補正係数を決定するようにさらに構成される、付記１に記載のＰＥＴ装置。
（付記３）処理回路は、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率に基づいて、検出器サブシステムのシングル計数損失補正係数を決定するようにさらに構成される、付記２に記載のＰＥＴ装置。
（付記４）処理回路は、検出器サブシステムの複数のペアの検出器サブシステムペアごとの測定された同時計数率に基づいて、ＰＥＴ装置のスキャナ同時計数損失補正係数を計算するようにさらに構成される、付記１に記載のＰＥＴ装置。
（付記５）処理回路は、ＰＥＴ装置のスキャナ同時計数損失補正係数を、重みの対応するセットを使用して検出器サブシステムの複数のペアに関して決定された同時計数損失補正係数の重み付けされた組み合わせとして計算するようにさらに構成され、該重みのセットは、検出器サブシステムの複数のペアの測定された同時計数率である、付記４に記載のＰＥＴ装置。
（付記６）処理回路は、
検出器サブシステムの複数のペアの特定の検出器サブシステムペアの計数損失が相関しているかどうかを判定し、
特定の検出器サブシステムペアに含まれている検出器サブシステムのシングル計数損失補正係数の最大値をとることによって、特定の検出器サブシステムペアの同時計数損失補正係数を決定するようにさらに構成される、
付記１に記載のＰＥＴ装置。
（付記７）処理回路は、
複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率に基づいて、相関シングル計数損失補正係数および非相関シングル計数損失補正係数を予測し、
複数の検出器サブシステムの検出器サブシステムごとに、相関シングル計数損失補正係数に非相関シングル計数損失補正係数を乗算することによって、検出器サブシステムのシングル計数損失補正係数を決定するようにさらに構成される、
付記１に記載のＰＥＴ装置。
（付記８）処理回路は、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率を検出器サブシステムでの相関シングル計数損失補正係数に関連付ける第１の関数、および検出器サブシステムで受信された測定済みのシングル計数率を検出器サブシステムでの非相関シングル計数損失補正係数に関連付ける第２の関数から、相関シングル計数損失補正係数を予測するようにさらに構成され、第１および第２の関数は、対応する特定の入力シングル計数率に対する、特定の相関シングル計数損失補正係数および特定の非相関シングル計数損失補正係数を出力するように訓練された訓練済みニューラルネットワークの出力から決定される、付記７に記載のＰＥＴ装置。
（付記９）処理回路は、複数の検出器サブシステムの検出器サブシステムごとに、相対シングル計数損失補正係数を決定し、計算されたスキャナ同時計数損失補正係数、決定された相対シングル計数損失補正係数、およびスキャンデータに基づいて、画像を再構成するようにさらに構成される、付記１に記載のＰＥＴ装置。
（付記１０）複数の検出器サブシステムを有するポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）スキャナ向けの計数損失補正の方法であって、
複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムのシングル計数損失補正係数を決定することと、
検出器サブシステムの複数のペアの検出器サブシステムペアごとに、検出器サブシステムペアの同時計数損失補正係数を決定することと、
検出器サブシステムの複数のペアの決定された同時計数損失補正係数に基づいて、ＰＥＴ装置のスキャナ同時計数損失補正係数を計算することと、
計算されたスキャナ同時計数損失補正係数およびＰＥＴ装置を使用して実施された患者のスキャンから収集されたスキャンデータに基づいて画像を再構成することと、
を含む、方法。
（付記１１）複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムの処理速度閾値に基づいて、検出器サブシステムのシングル計数損失補正係数を決定することをさらに含む、付記１０に記載の方法。
（付記１２）複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率に基づいて、検出器サブシステムのシングル計数損失補正係数を決定することをさらに含む、付記１１に記載の方法。
（付記１３）検出器サブシステムの複数のペアの検出器サブシステムペアごとの測定された同時計数率に基づいて、ＰＥＴ装置のスキャナ同時計数損失補正係数を計算することをさらに含む、付記１０に記載の方法。
（付記１４）ＰＥＴ装置のスキャナ同時計数損失補正係数を、重みの対応するセットを使用して検出器サブシステムの複数のペアに関して決定された同時計数損失補正係数の重み付けされた組み合わせとして計算することをさらに含み、該重みのセットは、検出器サブシステムの複数のペアの測定された同時計数率である、付記１３に記載の方法。
（付記１５）検出器サブシステムの複数のペアの特定の検出器サブシステムペアの計数損失が相関しているかどうかを判定することと、特定の検出器サブシステムペアに含まれている検出器サブシステムのシングル計数損失補正係数の最大値をとることによって、特定の検出器サブシステムペアの同時計数損失補正係数を決定することと、をさらに含む、付記１０に記載の方法。
（付記１６）複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率に基づいて、相関シングル計数損失補正係数および非相関シングル計数損失補正係数を予測することと、
複数の検出器サブシステムの検出器サブシステムごとに、相関シングル計数損失補正係数に非相関シングル計数損失補正係数を乗算することによって、検出器サブシステムのシングル計数損失補正係数を決定することと、
をさらに含む、付記１０に記載の方法。
（付記１７）予測するステップは、複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムで受信された測定済みのシングル計数率を検出器サブシステムでの相関シングル計数損失補正係数に関連付ける第１の関数、および検出器サブシステムで受信された測定済みのシングル計数率を検出器サブシステムでの非相関シングル計数損失補正係数に関連付ける第２の関数から、相関シングル計数損失補正係数を予測することをさらに含み、第１および第２の関数は、対応する特定の入力シングル計数率に対する、特定の相関シングル計数損失補正係数および特定の非相関シングル計数損失補正係数を出力するように訓練された訓練済みニューラルネットワークの出力から決定される、付記１６に記載の方法。
（付記１８）複数の検出器サブシステムの検出器サブシステムごとに、相対シングル計数損失補正係数を決定することと、計算されたスキャナ同時計数損失補正係数、決定された相対シングル計数損失補正係数、およびスキャンデータに基づいて、画像を再構成することと、をさらに含む、付記１７に記載の方法。
（付記１９）プログラムを記憶する非一時的コンピュータ可読媒体であって、
該プログラムは、複数の検出器サブシステムを有するポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）スキャナの処理回路によって実行されると、該処理回路に、計数損失補正の方法を実施させ、該方法は、
複数の検出器サブシステムの検出器サブシステムごとに、検出器サブシステムのシングル計数損失補正係数を決定することと、
検出器サブシステムの複数のペアの検出器サブシステムペアごとに、検出器サブシステムペアの同時計数損失補正係数を決定することと、
検出器サブシステムの複数のペアの決定された同時計数損失補正係数に基づいて、ＰＥＴ装置のスキャナ同時計数損失補正係数を計算することと、
計算されたスキャナ同時計数損失補正係数およびＰＥＴ装置を使用して実施された患者のスキャンから収集されたスキャンデータに基づいて画像を再構成することと、
を含む、コンピュータ可読媒体。
（付記２０）方法は、
検出器サブシステムの複数のペアの特定の検出器サブシステムペアの計数損失が相関しているかどうかを判定することと、
特定の検出器サブシステムペアに含まれている検出器サブシステムのシングル計数損失補正係数の最大値をとることによって、特定の検出器サブシステムペアの同時計数損失補正係数を決定することと、
をさらに含む、付記１９に記載のコンピュータ可読媒体。
（付記２１）
複数の検出器サブシステムと、
当該検出器サブシステムごとに、シングル計数損失補正係数を決定し、検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する処理回路と
を備える、ＰＥＴ装置。
（付記２２）
当該処理回路は、検出器サブシステムペアごとに測定された同時計数率に基づいて、スキャナ同時計数損失補正係数を計算する、付記２１に記載のＰＥＴ装置。
（付記２３）
当該処理回路は、前記スキャナ同時計数損失補正係数を、複数のペアの測定された同時計数率を使用して前記同時計数損失補正係数の重み付けされた組み合わせとして計算する、付記２２に記載のＰＥＴ装置。
（付記２４）
ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置が有する複数の検出器サブシステムの前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する、データ処理方法。
（付記２５）
ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置が有する複数の検出器サブシステムの前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する処理をコンピュータに実行させる、プログラム。
以上説明した少なくとも１つの実施形態によれば、画質を向上することができる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
４０１ ガンマ線検出器
５０２ 円形ボア
５０３ ネットワークコントローラ
５０４ ガントリ
５０５ メモリ
５０７ プロセッサ
５０８ データ収集システム

Claims (11)

  1. 複数の検出器サブシステムと、
    前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する処理回路と
    を備える、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置。
  2. 前記処理回路は、前記検出器サブシステムごとに前記検出器サブシステムの処理速度閾値に基づいて、前記シングル計数損失補正係数を決定する、請求項１に記載のＰＥＴ装置。
  3. 前記処理回路は、前記検出器サブシステムごとに、前記検出器サブシステムで受信された測定済みのシングル計数率に基づいて、前記シングル計数損失補正係数を決定する、請求項２に記載のＰＥＴ装置。
  4. 前記処理回路は、前記検出器サブシステムペアごとに測定された同時計数率に基づいて、前記スキャナ同時計数損失補正係数を計算する、請求項１に記載のＰＥＴ装置。
  5. 前記処理回路は、前記スキャナ同時計数損失補正係数を、前記複数のペアの前記測定された同時計数率を使用して前記同時計数損失補正係数の重み付けされた組み合わせとして計算する、請求項４に記載のＰＥＴ装置。
  6. 前記処理回路は、
    前記複数のペアの特定の検出器サブシステムペアの計数損失が相関しているかどうかを判定し、
    前記特定の検出器サブシステムペアに含まれている前記検出器サブシステムの前記シングル計数損失補正係数の最大値をとることによって、前記特定の検出器サブシステムペアの前記同時計数損失補正係数を決定する、請求項１に記載のＰＥＴ装置。
  7. 前記処理回路は、
    前記検出器サブシステムごとに、前記検出器サブシステムで受信された測定済みのシングル計数率に基づいて、相関シングル計数損失補正係数および非相関シングル計数損失補正係数を予測し、
    前記検出器サブシステムごとに、前記相関シングル計数損失補正係数に前記非相関シングル計数損失補正係数を乗算することによって、前記シングル計数損失補正係数を決定するようにさらに構成される、
    請求項１に記載のＰＥＴ装置。
  8. 前記処理回路は、前記検出器サブシステムごとに、前記測定済みのシングル計数率を前記相関シングル計数損失補正係数に関連付ける第１の関数、および前記測定済みのシングル計数率を前記非相関シングル計数損失補正係数に関連付ける第２の関数から、前記相関シングル計数損失補正係数を予測し、前記第１の関数および前記第２の関数は、対応する特定の入力シングル計数率に対する、特定の相関シングル計数損失補正係数および特定の非相関シングル計数損失補正係数を出力するように訓練された訓練済みニューラルネットワークの出力から決定される、請求項７に記載のＰＥＴ装置。
  9. 前記処理回路は、
    前記複数の検出器サブシステムの検出器サブシステムごとに、相対シングル計数損失補正係数を決定し、
    前記計算されたスキャナ同時計数損失補正係数、前記決定された相対シングル計数損失補正係数、および前記スキャンデータに基づいて、前記画像を再構成するようにさらに構成される、
    請求項１に記載のＰＥＴ装置。
  10. ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置が有する複数の検出器サブシステムの前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する、データ処理方法。
  11. ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）装置が有する複数の検出器サブシステムの前記検出器サブシステムごとに、シングル計数損失補正係数を決定し、前記検出器サブシステムの複数のペアの検出器サブシステムペアごとに、同時計数損失補正係数を決定し、前記同時計数損失補正係数に基づいて、スキャナ同時計数損失補正係数を計算し、被検体のスキャンから収集されたスキャンデータと前記スキャナ同時計数損失補正係数に基づいて画像を再構成する処理をコンピュータに実行させる、プログラム。
JP2023173213A 2022-10-05 2023-10-04 Ｐｅｔ装置、データ処理方法及びプログラム Pending JP2024054856A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
EP23201886.1A EP4350396A1 (en) 2022-10-05 2023-10-05 Pet apparatus, data processing method and program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US17/960,689 2022-10-05
US17/960,689 US20240125951A1 (en) 2022-10-05 2022-10-05 Method for correction of correlated count losses

Publications (1)

Publication Number Publication Date
JP2024054856A true JP2024054856A (ja) 2024-04-17

Family

ID=90627241

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2023173213A Pending JP2024054856A (ja) 2022-10-05 2023-10-04 Ｐｅｔ装置、データ処理方法及びプログラム

Country Status (2)

Country Link
US (1) US20240125951A1 (ja)
JP (1) JP2024054856A (ja)

Also Published As

Publication number Publication date
US20240125951A1 (en) 2024-04-18

Similar Documents

Publication Publication Date Title
RU2620862C2 (ru) Быстрая оценка рассеяния при реконструкции посредством позитронно-эмиссионной томографии
EP2984631B1 (en) Method for modeling and accounting for cascade gammas in images
CN103260521B (zh) 用于准确的输入函数估计的整合的工作流程
JP5920930B2 (ja) ガンマカメラ、ｓｐｅｃｔ装置、ｐｅｔ装置およびガンマ線計測画像生成方法
US9474501B2 (en) Hybrid method based on simulation and experimental data to normalize pet data
JP7286383B2 (ja) 陽電子放出撮像装置及び方法
US20200170605A1 (en) Method and apparatus for scatter correction in positron emission tomography (pet) imaging by performing a short pet scan in an extended region to estimate scatter coming from outside of the field of view (fov)
JP2021505907A (ja) 自動オンザフライ陽電子放出断層撮影（ｐｅｔ）スキャン計画及び最適化
US20210282732A1 (en) Method and apparatus for singles spectrum estimation and for dead-time correction in positron emission tomography (pet)
JP2021128143A (ja) ガンマ線放射イメージング装置、エネルギー校正方法及びプログラム
McIntosh et al. Validation of a GATE Model of $^{176} $ Lu Intrinsic Radioactivity in LSO PET Systems
JP2021173755A (ja) 医用画像処理装置、医用画像処理方法及びプログラム
WO2019145398A1 (en) A dead-time correction method in quantitative positron emission tomography (pet) reconstruction for various objects and radioactivity distributions
US20230281889A1 (en) Attenuation coefficient image generation method, nuclear medicine diagnostic apparatus, and trained model generation method
Bataille et al. Monte Carlo simulation for the ECAT HRRT using GATE
WO2009118843A1 (ja) 断層撮影装置、それを備えた撮影システム並びに撮影データ取得方法
JP2024054856A (ja) Ｐｅｔ装置、データ処理方法及びプログラム
EP4350396A1 (en) Pet apparatus, data processing method and program
US9750475B2 (en) Contour image generating device and nuclear medicine diagnosis apparatus
JP2024025712A (ja) 散乱フラクションの決定方法、放射線診断装置及び放射線診断システム
US20240023911A1 (en) Method to detect energy-dependent count loss
JP5073293B2 (ja) 断層撮影装置、それを備えた撮影システム並びに撮影データ取得方法
KR20200086814A (ko) 양전자 방출 단층 촬영기의 비행시간 정보를 이용하여 관측된 데이터에 포함된 반응깊이 정보를 보정하여 영상을 재구성하는 영상재구성방법
CN113253330B (zh) 伽玛射线放射成像装置及能量校准方法
JP7247782B2 (ja) 吸収係数画像推定方法、吸収係数画像推定プログラム、および、ポジトロンｃｔ装置