JP6188418B2 - 偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体 - Google Patents
偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体 Download PDFInfo
- Publication number
- JP6188418B2 JP6188418B2 JP2013103728A JP2013103728A JP6188418B2 JP 6188418 B2 JP6188418 B2 JP 6188418B2 JP 2013103728 A JP2013103728 A JP 2013103728A JP 2013103728 A JP2013103728 A JP 2013103728A JP 6188418 B2 JP6188418 B2 JP 6188418B2
- Authority
- JP
- Japan
- Prior art keywords
- tof
- entry
- count data
- fov
- list mode
- 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
Links
- 238000000034 method Methods 0.000 title claims description 39
- 230000009467 reduction Effects 0.000 title claims description 28
- 238000002600 positron emission tomography Methods 0.000 claims description 52
- 230000004044 response Effects 0.000 claims description 15
- 230000005855 radiation Effects 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000013078 crystal Substances 0.000 description 20
- 238000012545 processing Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 8
- 230000005251 gamma ray Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- AOYNUTHNTBLRMT-SLPGGIOYSA-N 2-deoxy-2-fluoro-aldehydo-D-glucose Chemical compound OC[C@@H](O)[C@@H](O)[C@H](O)[C@@H](F)C=O AOYNUTHNTBLRMT-SLPGGIOYSA-N 0.000 description 2
- 108090000248 Deleted entry Proteins 0.000 description 2
- 238000012879 PET imaging Methods 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 229940121896 radiopharmaceutical Drugs 0.000 description 2
- 239000012217 radiopharmaceutical Substances 0.000 description 2
- 230000002799 radiopharmaceutical effect Effects 0.000 description 2
- 230000004323 axial length Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012636 positron electron tomography Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/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
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Description
本発明の実施形態は、偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体に関する。
ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）は、ポジトロン（Positron）を放出する放射性医薬品が患者の身体に導入される核医学イメージングの１部門である。放射性医薬品が崩壊するにつれて、ポジトロンが発生する。具体的には、複数のポジトロンのそれぞれは、電子と反応し、対消滅イベント（annihilation event）として知られているイベントによって、ガンマ光子（gamma photon）のペア（pair）が同時に発生する。このガンマ光子は、同時計数線（line of coincidence）に沿って、略反対方向に進む。同時計数時間内に検出されたガンマ線のペアは、通常、ＰＥＴスキャナにより、対消滅イベントとして記録される。
飛行時間（Time-Of-Flight：ＴＯＦ）を用いたＴＯＦイメージングでは、同時計数対（coincident pair）の各ガンマ光子が検出される同時計数間隔（coincidence interval）の時間も測定される。飛行時間情報は、同時計数線に沿って検出されたイベントの対消滅位置を示す情報を提供する。複数の対消滅イベントから得られるデータは、一般的には、統計的（反復的）、又は、解析的な再構成アルゴリズムを用いて、スキャンされる被検体又は物体の画像の再構成や生成のために使用される。
図１は、ＰＥＴイメージング装置の幾何学的配置の一例を示す図である。図１は、放出されたポジトロンの横断座標（transaxial coordinate）及び軸座標（axial coordinate）と、３Ｄ検出器（3D detector）で測定された応答線（Line-Of-Response：ＬＯＲ）とを示した図である。なお、応答線（ＬＯＲ）は、上記の同時計数線と同義である。横断座標は、例えば、スキャナの軸方向や被検体の体軸方向に垂直な断面内で設定される座標であり、軸座標は、例えば、スキャナの軸方向や被検体の体軸方向に沿った断面内で設定される座標である。座標（ｘｅ，ｙｅ，ｚｅ）又は（ｓｅ，ｔｅ，ｚｅ）は、放出されたポジトロン「ｅ」の画像座標を示す。また、座標（ｘａ，ｙａ，ｚａ）は、ＬＯＲとして測定された一方のガンマ線を検出した検出器結晶「ａ」の位置を示し、座標（ｘｂ，ｙｂ，ｚｂ）は、ＬＯＲとして測定された他方のガンマ線を検出した検出器結晶「ｂ」の位置を示す。測定されたＬＯＲの投影座標は、非ＴＯＦでは、「（ｓ，φ，ｚ、θ）、ここでｚ＝（ｚａ＋ｚｂ）／２」で表すことができる。「φ」は、横断座標における線分ａｂの傾斜角に対応し、「θ」は、軸座標における線分ａｂの傾斜角に対応する。或いは、測定されたＬＯＲの投影座標は、ＴＯＦ−ＬＯＲ用のさらなる次元「ｔ」を含んでもよい。「ｔ」は、時間に対応する値であり、例えば、「ｔａ」は、同時計数時間内に「ａ」がガンマ光子を検出した時間に対応し、「ｔｂ」は、同時計数時間内に「ｂ」がガンマ光子を検出した時間に対応し、「ｔｅ」は、同時計数時間内に「ａ」がガンマ光子を検出した時間と「ｂ」がガンマ光子を検出した時間との時間差に対応する。「ｔ」の情報を用いることで、測定されたＬＯＲ上における「ｅ」の座標、すなわち、対消滅点の位置を推定することができる。
ＰＥＴでは、偶発同時計数（random coincidence）は、真の同時計数（true coincidence）の検出に使用される同時計数窓（coincidence window）の有限幅によって発生する。２つの無相関な単一イベント（single event）が、偶発的に同時計数窓内で検出された場合、それらは、誤って真の同時計数イベントとして識別されて、記録される可能性がある。偶発イベントの率は、以下の式（１）に示すように、各検出器の単一イベントの率と、同時計数窓のサイズとに比例する。なお、式（１）では、「Ｃｉｊ」は、ｉ番目の検出器とｊ番目の検出器とを結ぶＬＯＲにおける偶発同時計数率（random coincidences count rate）を示し、「ｒｉ」は、ｉ番目の検出器における単一イベントの計数率を示し、「ｒｊ」は、ｊ番目の検出器における単一イベントの計数率を示し、「τ」は、同時計数窓のサイズ（coincidence window size）を示す。
特に、３次元ＰＥＴ（３ＤＰＥＴ）スキャナを用いる撮影や、高放射能濃度での撮影では、偶発同時計数は、記録されたプロンプト同時計数イベント（prompt coincidence event）の大部分を含む可能性がある。プロンプト同時計数イベントは、真の同時計数イベント、散乱（scatter）の同時計数イベント及び偶発同時計数イベントを含む。適切に補正されない場合、偶発同時計数イベントは、再構成画像に相当な量の誤差を導入する可能性がある。
プロンプト同時計数イベント内の偶発イベントの量は、同時計数時間窓（coincidence timing window）と、同時計数を検出する２つの対向する結晶における単一イベントの計数率とによって決定される。同時計数時間窓は、収集される（再構成される）有効視野（Field-Of-View：ＦＯＶ）のサイズと、スキャナの時間分解能（timing resolution）とによって設定される。ＰＥＴスキャナの一般的なＦＯＶは、３次元円筒形である。この３次元円筒形は、スキャナの横断面方向に中心をおく円形領域を底面とし、軸方向における長さがスキャナの軸方向における長さと同じである。従って、デフォルトでは、収集されるＦＯＶは、横断面方向の円形領域における直径を引用する。ＰＥＴ全身イメージングの場合、ＦＯＶは、通常、５７６ｍｍ〜７００ｍｍである。時間分解能は、データ収集システムの結晶のタイプ、光学的光子検出器及びフロントエンド（front-end）のエレクトロニクス（electronics）に依存する。ＬＹＳＯ結晶の場合、典型的な時間分解能は、５００ｐｓ〜６５０ｐｓである。関連する同時計数時間窓は、この場合、全身イメージングの場合４ｎｓ〜６ｎｓである。ＦＤＧ（Ｆｌｕｏｒｏｄｅｏｘｙｇｌｕｃｏｓｅ：フルオロデオキシグルコース）を用いた全身撮影における偶発分画（偶発／プロンプト）は、被検体の軸方向部分の放射能濃度により、約３０％〜５０％となる。
しかし、固定した同時計数窓方法では、プロンプト同時計数データに必然的に大量の偶発イベントが導入される。例えば、単一イベントの計数率が全ての結晶で同じとなる固定したファントム／被検体の場合、偶発の総量は、収集されるＦＯＶの円形領域と、同時計数窓ＦＯＶのサイズとに正比例する。なお、同時計数窓ＦＯＶは、「２×１５０（ｍｍ／ｎｓ）×ｃｏｉｎＷｉｎ（ｎｓ）」として計算される。ただし、「ｃｏｉｎＷｉｎ」は、単位が「ｎｓ（ナノ秒）」の同時計数窓の長さである。
偶発イベント数を減らすためにさまざまな方法が、提案されている。
１つのアプローチでは、ＰＥＴスキャナの最大の同時計数時間窓より小さい短時間窓を用いて、利用可能な全てのＴＯＦビン（bin）より少ないＴＯＦビンの中央サブセット（subset）だけに投影データの収集を制限する。短時間窓ＦＯＶは、画像化される対象物の関心領域が、実際の対象物より小さくなることを表わす。換言すると、このアプローチは、ＴＯＦ情報から推定される対消滅点の位置が、再構成用のＦＯＶの中央部となる同時計数情報のみを採用する方法である。また、短時間窓ＦＯＶは、例えば、高計数率なＲｂ−８２（Rubidium：ルビジウム−８２）によるＰＥＴ心臓イメージングに適用される場合、リストモードデータの転送速度を下げる可能性がある。この短時間窓を用いる方法は、短時間窓を極端に小さく設定した場合、偶発イベントの数を減らすだけでなく、真のイベントの数も減らしてしまう。
もう一つのアプローチでは、個々の被検体のＣＴ（Computed Tomography：コンピュータ断層撮影）の形態学的情報を用いて同時計数時間窓を調整する。換言すると、このアプローチは、ＴＯＦ情報とＣＴ情報とを用いて、推定対消滅点が被検体の輪郭内に位置する同時計数情報のみを採用する方法である。ある応答線（Line-Of-Response：ＬＯＲ）について、ＴＯＦビンの下界および上界は、被検体の身体と交差するＬＯＲの入口点及び出口点を探索することで、決定される。ＴＯＦマスキング（masking）の効果は、主に偶発である被検体の外側のイベントを排除する、被検体の境界線に沿ったサイノグラムのトリミングにある。しかし、タイミングの統計学的不確実性（すなわち、完璧でない時間分解能）が、考慮に入れられていないので、真の同時計数の約１０〜１５％が、このＴＯＦマスキング手法で除去されてしまう。
第３のアプローチでは、異なるリング差の関数としての可変同時計数窓を用いて偶発イベントを減らす。この方法は、軸方向に長いＦＯＶとなるＰＥＴスキャナ、例えば、２０ｃｍを上回る軸方向ＦＯＶを有するＰＥＴスキャナに対してより適用可能である。
他の偶発削減方法は、異なる対になった検出器モジュール間で調整可能な同時計数窓を用いる方法、サイノグラム径方向のビン位置を用いる方法との両方又は一方を含み、ＦＯＶの外側に位置するＬＯＲを除去する。これらの方法は、時間分解能を考慮に入れていないので、やはり真の同時計数イベントを減らすことも避けられない。
本発明が解決しようとする課題は、偶発イベントを削除することができる偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体を提供することである。
実施形態の偶発イベント削減方法は、ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータ内の偶発イベントを削減する方法であり、取得ステップと、フィルタステップとを含む。取得ステップは、所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含むプロンプトリストモード計数データをＰＥＴスキャナが取得する。フィルタステップは、前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記取得ステップにより取得された前記プロンプトリストモード計数データから削除することで、当該プロンプトリストモード計数データをＰＥＴスキャナがフィルタ処理する。
以下、添付図面を参照して、偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体の実施形態を詳細に説明する。本明細書に記載する実施形態とこれに付随する多くの利点は、以下の詳細な説明を参照し添付図面と関連付けて考えれば、より完全に理解できる。
本明細書で開示する実施形態は、ＴＯＦ分解能及び収集されるＦＯＶの情報を利用するだけで、汎用的に、偶発イベントを削減する技術を対象としている。本実施形態では、収集されるＦＯＶ（再構成用のＦＯＶ）に貢献しない偶発同時計数イベントは、全て削除されるが、削除される偶発同時計数イベントは、真の同時計数イベントに影響を及ぼさない。
一実施形態では、ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータ内の偶発イベントを削減する偶発イベント削減方法が提供される。該方法は、（１）所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含むプロンプトリストモード計数データ（又は、遅延リストモード計数データ）を取得する取得ステップと、（２）再構成用のＦＯＶに対する接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、取得ステップにより取得されたプロンプトリストモード計数データ（又は、遅延リストモード計数データ）から削除することで、当該プロンプトリストモード計数データ（又は、遅延リストモード計数データ）をフィルタ処理するフィルタステップと、を含む。なお、プロンプトリストモード計数データは、同時計数用の時間窓を用いて同時計数された２つの計数情報（結晶位置、時間、エネルギー等）を１つのエントリとするリストである。プロンプトリストモード計数データには、真の同時計数情報とともに、散乱に由来する同時計数情報及び偶発同時計数情報が含まれる。また、遅延リストモード計数データは、一方のイベントの計数情報が含む時間に、同時計数用の時間窓より大きい時間を加算したうえで、同時計数用の時間窓を用いて同時計数とされた２つの計数情報を１つのエントリとするリストである。遅延リストモード計数データは、プロンプトリストモード計数データに含まれる偶発イベントを推定するために用いられる。
一実施形態では、接線方向のＴＯＦマスクは、ＰＥＴスキャナの時間分解能により定まる「推定放射点が存在する位置の確率分布」に基づいて設定される。一実施形態では、フィルタステップは、接線方向のＴＯＦマスクを表す以下の式（２）の不等式を満たさないエントリを削除する。
ただし、「ｔａ」及び「ｔｂ」それぞれは、エントリに含まれる「イベントａ」及び「イベントｂ」それぞれの到達時間である。「ｃ／２」は、時間を距離に変換するための変換係数（ｃは、光速）である。「θ」は、エントリにより表される応答線の横断面に対する軸方向の傾斜角である。「ｄ」は、再構成用のＦＯＶの直径である。「ｓ」は、エントリにより表される応答線の径方向の値である。「σＴＯＦ」は、空間領域に変換したＰＥＴスキャナの時間分解能の標準偏差である。「ｎ」は、予め設定された値であり、望ましくは、範囲「２、３」内の所定値であり、真のＬＯＲがマスクによって削除されない信頼区間を示す。
更に、一実施形態では、上記の方法は、フィルタ処理されたリストモード計数データで再構成を実行する再構成ステップを更に含む。
さて、図を参照すると、再構成ＦＯＶは、通常、横断面に一定の直径を有し、軸方向の長さが「スキャナの軸方向のＦＯＶを網羅する」長さを有する円筒として設定される。従って、横断面では、円形のＦＯＶは、２次元空間の座標系、例えば、図１に示すように、半径方向「ｓ」と接線方向「ｔ」との投影座標により定義することができる。
空間的分解能は、投影の半径方向及び接線方向で劇的に異なる。半径方向の分解能は、スキャナの横断面の結晶サイズによって、決定される。半径方向の分解能のＦＷＨＭ（Full Width at Half Maximum：半値全幅）は、通常、２〜４ｍｍである（例えば、４ｍｍの結晶ピッチが与えられる）。接線方向の分解能は、システムの時間分解能によって決定される。接線方向の分解能のＦＷＨＭは、空間領域に変換された場合、６０〜９０ｍｍである（例えば、４００〜６００ｐｓ時間分解能が与えられる）。
図２は、円筒における接線方向のＴＯＦマスクを横断面及び矢状面で示す図である。一実施形態では、図２の陰影領域として示された接線方向のＴＯＦマスクは、プロンプトデータにおける偶発イベントと、起こり得る散乱イベントの量を低減するために使用される。特に、上記の式（２）の不等式が満たされるとき、所与のＬＯＲｂ（図２に示す線ＡＢ）に対する真のデータは、フィルタによって保持されて削除されない。以下、上記の式（２）を再度記載する。
ただし、「ｔａ」及び「ｔｂ」それぞれは、エントリに含まれる「イベントａ」及び「イベントｂ」それぞれの到達時間（ｐｓ：picosecond：ピコ秒）である。「ｃ／２＝０．１５ｍｍ／ｐｓ」は、ｐｓをｍｍに変換するための変換係数である。「θ」は、応答線の横断面に対する軸方向の傾斜角である。「ｄ」は、再構成ＦＯＶの直径である。「ｓ」は、応答線と、再構成ＦＯＶの横断面における中心との径方向における距離である。「σＴＯＦ」は、空間領域に変換したＰＥＴスキャナの時間分解能の標準偏差である。「ｎ」は、予め設定された値であり、望ましくは、範囲「２、３」内の所定値である。ＴＯＦマスクは、本質的に、再構成ＦＯＶシリンダーを±ｔ次元に沿ってσＴＯＦだけ伸ばしている。ＴＯＦマスクは、再構成ＦＯＶシリンダーを、ＬＯＲに平行なＦＯＶの接線方向に沿って、「±ｎσＴＯＦ」だけ伸ばすこととなる。式（２）に示す項「｛（ｄ２／４）−ｓ２｝１／２」は、直径ｄの円ＦＯＶと交差する線Ａ’ Ｂ’の弦長の半分である。式（２）に示す項「（ｃ／２）×｜ｔａ−ｔｂ｜」は、放射点（対消滅点）からＬＯＲの中間点Ｏまでの距離、すなわち、セグメントＯＥの長さを表わす。
なお、放射点は、ＬＯＲ（線ＡＢ）上のどこでも良い。図２では、放射点が、ＴＯＦマスクの端にある場合を例示している。項「（ｃ／２）×｜ｔａ−ｔｂ｜×ｃｏｓθ」は、単に横断面内のＯＥの投影である。
ガウス分布の時間分解能の場合、「３σＴＯＦ＝３＊ＦＷＨＭＴＯＦ／２．３５５」であり、４５０ｐｓのＦＷＨＭＴＯＦに対して「３σＴＯＦ＝８６ｍｍ」となる。なお、「±３σＴＯＦ」は、９９．７％の信頼区間を達成するのに用いられる。信頼区間を９５．４％に減らした場合、「±２σＴＯＦ」が必要となり、より多くの偶発イベントと起こり得る散乱イベントとは、削除される。
偶発分布は、（ｓ、ｔ）次元で均一であると仮定すると、任意のＴＯＦマスクを使う前では、（ｓ、ｔ）次元の偶発イベントの数は、「ｄ×２×１５０（ｍｍ／ｎｓ）×ｃｏｉｎＷｉｎ（ｎｓ）」に比例している。接線方向のＴＯＦマスクを使った後では、（ｓ、ｔ）次元の偶発イベントの数は、「２ｎσＴＯＦｄ＋πｄ２／４」に比例する。
ＩＥＣ（International Electrotechnical Commission：国際電気標準会議）のファントムを用いた場合、偶発イベントの数は、接線方向のＴＯＦマスクを用いることで、５０％まで減らすことができ、同様に、ＴＯＦ再構成を用いない場合でも、雑音等価計数率（noise-equivalent-count rate）は、１２％まで増えたことが分かった。
なお、ＴＯＦマスクは、更に、別の形状の再構成ＦＯＶ、例えば、被検体の横断面により共形な再構成ＦＯＶに対しても、使用できる。
図３は、楕円筒における接線方向のＴＯＦマスクを横断面及び矢状面で示す図である。図３は、軸方向に均一に分布する楕円形状の横断面ＦＯＶ、すなわち、３次元の楕円形円筒を示した図である。図３の（ａ）は、横断図であり、図３の（ｂ）は、矢状図である。陰影のついたエリア（area）は、ｘ軸に対して横断面回転角φを有するＬＯＲに対する楕円ＦＯＶのためのＴＯＦマスクである。マスクは、楕円を±ｔ次元に沿って３σＴＯＦだけ伸ばすことによって生成される。なお、楕円等の非円形の再構成ＦＯＶの場合、ＴＯＦマスクの形状は、ＬＯＲの横断面角φに依存する。
図４は、一実施形態に係るＰＥＴにおいて、偶発イベントをフィルタ処理する方法のフローチャートである。図４は、本実施形態に係る方法を示したフローチャートである。ステップ４００では、所与の再構成用の有効視野（ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、ＴＯＦ情報（飛行時間情報）を含むＴＯＦプロンプトリストモード計数データ（プロンプトリストモード計数データ）を取得する。或いは、ステップ４００では、所与の再構成用のＦＯＶを有するＰＥＴスキャナについて、複数のエントリを含み、ＴＯＦ情報を含むＴＯＦ遅延リストモード計数データ（遅延リストモード計数データ）を取得する。
Ｎｏｎ−ＴＯＦリストモードデータには、一般的に、同時計数イベントごとに以下のフォーマット、｛ｘａ、ｚａ、ｘｂ、ｚｂ、ｅａ、ｅｂ｝がある。ここで、ｘａ及びｘｂは、それぞれイベントａ及びｂが入射した結晶の横断面での番号であり、ｚａ及びｚｂは、それぞれイベントａ及びｂが入射した結晶の軸方向での番号であり、ｅａ及びｅｂは、それぞれイベントａ及びｂのエネルギレベルである。同様に、ＴＯＦリストモードデータには、フォーマット｛ｘａ、ｚａ、ｘｂ、ｚｂ、ｅａ、ｅｂ、ｔｏｆ｝がある。ここで、ｔｏｆは、イベントａ及びｂの到着時間差、すなわち、ｔｏｆ＝ｔａ−ｔｂである。
なお、遅延同時計数窓方法は、ＴＯＦ遅延リストモードデータの生成に使用でき、また、本明細書で開示する偶発イベントフィルタ処理方法にも使用できる。このように、ステップ４００では、遅延リストモードデータを得ることができる。
ステップ４１０では、リストモードデータのエントリごとに径方向の距離ｓ及び軸方向の傾斜角θを算出する。
ステップ４２０では、フィルタ処理されたリストモード計数データを得るために、接線方向のＴＯＦマスク（例えば、図２に示されたマスク）の外側に存在する放射点を表すエントリを、得られたプロンプトリストモード計数データから削除することで、得られたプロンプトリストモード計数データをフィルタ処理する。或いは、ステップ４２０では、フィルタ処理されたリストモード計数データを得るために、接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、得られた遅延リストモード計数データから削除することで、得られた遅延リストモード計数データをフィルタ処理する。
すなわち、ステップ４２０では、フィルタ処理は、図４に示す不等式が満たされないエントリを削除することによって実行される。
なお、上述したように、図４に示す不等式においても、「ｔａ」及び「ｔｂ」それぞれは、エントリに含まれる「イベントａ」及び「イベントｂ」それぞれの到達時間である。「ｃ／２」は、時間を距離に変換するための変換係数である。「θ」は、エントリにより表される応答線の横断面に対する軸方向の傾斜角である。「ｄ」は、再構成用のＦＯＶの直径である。「ｓ」は、エントリにより表される応答線の径方向の値である。「σＴＯＦ」は、空間領域に変換したＰＥＴスキャナの時間分解能の標準偏差である。「ｎ」は、望ましくは、範囲「２、３」内の所定値である。
ステップ４３０では、ＰＥＴ画像再構成が、フィルタ処理されたリストモード計数データで実行される。具体的には、プロンプトリストモード計数データを接線方向のＴＯＦマスクによりフィルタ処理した場合、フィルタ処理後のプロンプトリストモード計数データを用いてＰＥＴ画像が再構成される。フィルタ処理後のプロンプトリストモード計数データには、偶発イベントの同時計数情報が削除されていることから、かかるＰＥＴ画像は、偶発イベントに起因するアーチファクトが低減された画像となる。
或いは、遅延リストモード計数データを接線方向のＴＯＦマスクによりフィルタ処理した場合、例えば、プロンプトリストモード計数データを用いて再構成された画像から、フィルタ処理後の遅延リストモード計数データを用いて再構成された画像を差分することで、最終的なＰＥＴ画像を得る。すなわち、フィルタ処理前の遅延リストモード計数データは、フィルタ処理前の遅延リストモード計数データから、接線方向のＴＯＦマスク以外の領域に推定対消滅点が位置するＬＯＲのエントリを削除したデータである。遅延リストモード計数データから削除されるエントリを「削除エントリ」と呼ぶと、削除エントリ以外のエントリは、偶発イベントに起因するアーチファクトを低減するために有効である。しかし、削除エントリを用いた場合、補正する必要の無い領域においても補正処理が行なわれるため、ＰＥＴ画像に他のアーチファクトが発生する場合がある。このように、遅延リストモード計数データを接線方向のＴＯＦマスクによりフィルタ処理することで、偶発イベントを低減することができる。
なお、本実施形態は、プロンプトリストモード計数データ及び遅延リストモード計数データの双方に対して、接線方向のＴＯＦマスクを用いたフィルタ処理が行なわれても良い。フィルタ処理後のプロンプトリストモード計数データには、偶発イベントの同時計数情報が完全に削除されていない場合がある。そこで、例えば、フィルタ処理後のプロンプトリストモード計数データを用いて再構成された画像から、フィルタ処理後の遅延リストモード計数データを用いて再構成された画像を差分することで、偶発イベントに起因するアーチファクトが更に低減されたＰＥＴ画像を得ることができる。
図５は、一実施形態に係るＰＥＴシステムのハードウェア構成を示す図である。図５は、本実施形態で使用できるＰＥＴのハードウェア構成例を示す図である。図５では、光電子増倍管１３５及び１４０が、ライトガイド１３０の上に配列され、シンチレーション結晶１０５のアレイが、ライトガイド１３０の下に配列される。シンチレーション結晶１２５の第２のアレイは、シンチレーション結晶１０５に対向して配置され、その上にライトガイド１１５と光電子増倍管１９５及び１１０とが配列される。
図５では、２つのガンマ線が、検査下の身体から放出されると、２つのガンマ線は、互いに略１８０°の反対方向に進む。個々のガンマ線の検出は、シンチレーション結晶１００及び１２０において同時に起こる。ガンマ線が、所定の制限時間内にシンチレーション結晶１００及び１２０において検出されたとき、シンチレーションイベントは、確定される。このように、ガンマ線タイミング検出システムは、２つのガンマ線をシンチレーション結晶１００及び１２０において同時に検出する。しかし、以下では、簡単のために、シンチレーション結晶１００におけるガンマ線検出について記述する。しかし、当業者には自明であるが、本明細書においてシンチレーション結晶１００に対して与えられる説明は、シンチレーション結晶１２０におけるガンマ線検出にも同様に適用できる。
各光電子増倍管１１０、１３５、１４０、１９５は、データ取得ユニット１５０にそれぞれ接続されている。データ取得ユニット１５０は、光電子増倍管からの信号を処理するように構成されたハードウェアを含む。データ取得ユニット１５０は、ガンマ線の到達時間を測定する。データ取得ユニット１５０は、システムクロック（図示せず）に対する識別パルスの時間を符号化する２つの出力（ＰＭＴ（Photomultiplier Tube：光電子増倍管）１３５／１４０の組合せのための１つ、及び、ＰＭＴ１１０／１９５の組合せのための１つ）を生成する。ＴＯＦ−ＰＥＴシステムの場合、データ取得ユニット１５０は、一般的に１５〜２５ｐｓの精度でタイムスタンプを生成する。データ取得ユニット１５０は、ＰＭＴごとの信号（データ取得ユニット１５０からの出力のうちの４つ）の振幅を測定する。
データ取得ユニット１５０の出力は、中央処理装置（Central Processing Unit：ＣＰＵ）１７０に供給され処理される。処理工程には、データ取得ユニット１５０の出力からのエネルギ及び位置と、各イベントに対するタイムスタンプ出力からの到着時間との推定工程が含まれる。処理工程は、エネルギ、位置及び時間の推定の精度を改善するために、従来技術の較正に基づく多くの較正ステップを適用して含んでも良い。
さらに、ＣＰＵ１７０は、上述の図４に示したフローチャートに従って偶発イベントを推定するための方法を実行するように構成される。
当業者には自明であるが、ＣＰＵ１７０は、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（Field-Programmable Gate Array：ＦＰＧＡ）、又は、他の複合プログラマブル論理回路（Complex Programmable Logic Device：ＣＰＬＤ）のような個別論理ゲートとして実装できる。ＦＰＧＡ又はＣＰＬＤの実装は、ＶＨＤＬ（VHSIC Hardware Description Language）、Ｖｅｒｉｌｏｇ、又は、他の任意のハードウェア記述言語でコード化されると良い。コードは、直接、ＦＰＧＡ又はＣＰＬＤ内に、或いは、個別の電子メモリとして電子メモリに格納されるとよい。更に、電子メモリは、ＲＯＭ（Read Only Memory）、ＥＰＲＯＭ（Electrically Programmable Read Only Memory）、ＥＥＰＲＯＭ（Electrically erasable Programmable Read only Memory）、又は、フラッシュ（ＦＬＡＳＨ）メモリ等の不揮発性でも良い。電子メモリは、また、静的または動的のＲＡＭなどの揮発性であっても良い。また、ＦＰＧＡ又はＣＰＬＤと電子メモリの間の対話のみならず電子メモリを管理するためにマイクロコントローラやマイクロプロセッサー等の処理装置を備えても良い。
或いは、ＣＰＵ１７０は、上記の電子メモリとハードディスク、ＣＤ、ＤＶＤ、フラッシュドライブ等の既知の記憶媒体の両方または一方のいずれかに格納されている一連のコンピュータ可読命令として実装してもよい。更に、コンピュータ可読命令は、ユーティリティアプリケーション、バックグラウンドデーモン、またはオペレーティングシステムの構成要素、あるいはこれらの組み合わせとして提供され、米国Ｉｎｔｅｌ社製のＸｅｏｎプロセッサ（登録商標）、又は、米国ＡＭＤ社製のＯｐｔｅｒｏｎプロセッサ（登録商標）等の処理装置、及び、Ｍｉｃｒｏｓｏｆｔ ＶＩＳＴＡ（登録商標）、ＵＮＩＸ（登録商標）、Ｓｏｌａｒｉｓ（登録商標）、ＬＩＮＵＸ（登録商標）、及びＡｐｐｌｅ ＭＡＣ−ＯＳ（登録商標）等の当業者に知られているオペレーティングシステムと連動して実行される。
一旦、ＣＰＵ１７０で処理されると、処理された信号は、電子格納部１８０への格納、及び、表示部１４５での表示の両方または一方が行われる。当業者には自明であるが、電子格納部１８０は、ハードディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、フラッシュドライブ、ＲＡＭ、ＲＯＭ、又は、技術的に周知の任意の他の電子格納部（electronic storage）でも良い。表示部１４５は、液晶ディスプレイ（liquid crystal display：ＬＣＤ）、ブラウン管ディスプレイ（cathod ray tube：ＣＲＴ）、プラズマディスプレイ、有機発光ダイオード（organic light emitting diode：ＯＬＥＤ）、発光ダイオード（ＬＥＤ）等の当技術分野で知られているディスプレイとして実装してもよい。このように、ここで説明した電子格納部１８０および表示部１４５の記載は例示にすぎず、決して本実施形態の進歩の範囲を制限するものではない。
また、図５は、ガンマ線検出システムが、他の外部デバイス及びユーザーの両方又は一方と接続するためのインターフェース１７５も含む。例えば、インターフェース１７５は、ユニバーサルシリアルバス（Ｕｎｉｖｅｒｓａｌ Ｓｅｒｉａｌ Ｂｕｓ：ＵＳＢ）インターフェース、パーソナルコンピュータメモリカード国際協会（Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ Ｍｅｍｏｒｙ Ｃａｒｄ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ａｓｓｏｃｉａｔｉｏｎ：ＰＣＭＣＩＡ）インターフェース、イーサネット（登録商標）インターフェース等の当技術分野で知られているインターフェースであって良い。また、インターフェース１７５は、有線又は無線であっても良く、キーボードとマウスの両方又は一方、等のユーザーと対話するためのヒューマンインターフェースデバイスを含んでよい。
以上、説明したとおり本実施形態によれば、偶発イベントを削除することができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
１１０、１３５、１４０、１９５ 光電子増倍管
１４５ 表示部
１５０ データ取得ユニット
１７０ 中央処理装置（Central Processing Unit：CPU）
１７５ インターフェース
１８０ 電子格納部
１４５ 表示部
１５０ データ取得ユニット
１７０ 中央処理装置（Central Processing Unit：CPU）
１７５ インターフェース
１８０ 電子格納部
Claims (13)
- ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータ内の偶発イベントを削減する偶発イベント削減方法であって、
所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含むプロンプトリストモード計数データをＰＥＴスキャナが取得する取得ステップと、
前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記取得ステップにより取得された前記プロンプトリストモード計数データから削除することで、当該プロンプトリストモード計数データをＰＥＴスキャナがフィルタ処理するフィルタステップと、
を含む、偶発イベント削減方法。 - 前記接線方向のＴＯＦマスクは、前記エントリにより表される応答線上で放射点が存在する位置の確率分布であり、前記ＰＥＴスキャナの時間分解能により定まる確率分布と、に基づいて前記再構成用のＦＯＶを接線方向に伸ばすことで設定される、請求項１に記載の偶発イベント削減方法。
- 前記フィルタステップは、「ｔａ」及び「ｔｂ」それぞれを、前記エントリに含まれる「イベントａ」及び「イベントｂ」それぞれの到達時間とし、「ｃ／２」を、時間を距離に変換するための変換係数とし、「θ」を、前記エントリにより表される応答線の横断面に対する軸方向の傾斜角とし、「ｄ」を、前記再構成用のＦＯＶの直径とし、「ｓ」を、前記エントリにより表される前記応答線の径方向の値とし、「σＴＯＦ」を、空間領域に変換した前記ＰＥＴスキャナの時間分解能の標準偏差とし、「ｎ」を、予め設定された値とした場合、前記接線方向のＴＯＦマスクを表す以下の式（１）の不等式を満たさないエントリを削除する、請求項１に記載の偶発イベント削減方法。
- ＰＥＴスキャナが前記フィルタ処理されたリストモード計数データで再構成を実行する再構成ステップ
を更に含む、請求項１に記載の偶発イベント削減方法。 - 前記取得ステップは、更に、複数のエントリを含み、各エントリが飛行時間の情報を含む遅延リストモード計数データを取得し、
前記フィルタステップは、更に、前記接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記取得ステップにより取得された前記遅延リストモード計数データから削除することで、当該遅延リストモード計数データをフィルタ処理する、請求項１に記載の偶発イベント削減方法。 - ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータの偶発イベントを削減する偶発イベント削減方法であって、
所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含む遅延リストモード計数データをＰＥＴスキャナが取得する取得ステップと、
前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記取得ステップにより取得された前記遅延リストモード計数データから削除することで、当該遅延リストモード計数データをＰＥＴスキャナがフィルタ処理するフィルタステップと、
を含む、偶発イベント削減方法。 - ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータ内の偶発イベントを削減する偶発イベント削減装置であって、
所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含むプロンプトリストモード計数データを記憶するメモリと、
前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記プロンプトリストモード計数データから削除することで、当該プロンプトリストモード計数データをフィルタ処理するプロセッサと、
を備える、偶発イベント削減装置。 - 前記接線方向のＴＯＦマスクは、前記エントリにより表される応答線上で放射点が存在する位置の確率分布であり、前記ＰＥＴスキャナの時間分解能により定まる確率分布と、に基づいて設定される、請求項７に記載の偶発イベント削減装置。
- 前記プロセッサは、「ｔａ」及び「ｔｂ」それぞれを、前記エントリに含まれる「イベントａ」及び「イベントｂ」それぞれの到達時間とし、「ｃ／２」を、時間を距離に変換するための変換係数とし、「θ」を、前記エントリにより表される応答線の横断面に対する軸方向の傾斜角とし、「ｄ」を、前記再構成用のＦＯＶの直径とし、「ｓ」を、前記エントリにより表される前記応答線の径方向の値とし、「σＴＯＦ」を、空間領域に変換した前記ＰＥＴスキャナの時間分解能の標準偏差とし、「ｎ」を、予め設定された値とした場合、前記接線方向のＴＯＦマスクを表す以下の式（２）の不等式を満たさないエントリを削除する、請求項７に記載の偶発イベント削減装置。
- 前記プロセッサは、前記フィルタ処理されたリストモード計数データで再構成を実行する、請求項７に記載の偶発イベント削減装置。
- 前記メモリは、更に、複数のエントリを含み、各エントリが飛行時間の情報を含む遅延リストモード計数データを記憶し、
前記プロセッサは、更に、前記接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、取得された前記遅延リストモード計数データから削除することで、当該遅延リストモード計数データをフィルタ処理する、請求項７に記載の偶発イベント削減装置。 - ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）で得られるリストモードデータ内の偶発イベントを削減する偶発イベント削減装置であって、
所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含む遅延リストモード計数データを記憶するメモリと、
前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記遅延リストモード計数データから削除することで、当該遅延リストモード計数データをフィルタ処理するプロセッサと、
を備える、偶発イベント削減装置。 - プログラムを格納する非一時的コンピュータ可読記憶媒体であって、
前記プログラムがプロセッサにより実行された場合、前記プロセッサは、
所与の再構成用の有効視野（Field-Of-View：ＦＯＶ）を有するＰＥＴスキャナについて、複数のエントリを含み、各エントリが飛行時間（Time-Of-Flight：ＴＯＦ）の情報を含むプロンプトリストモード計数データを取得する取得ステップと、
前記再構成用のＦＯＶに対するＴＯＦマスクであって、前記再構成用のＦＯＶを接線方向に伸ばすことで設定される接線方向のＴＯＦマスクの外側に存在する放射点を表すエントリを、前記取得ステップにより取得された前記プロンプトリストモード計数データから削除することで、当該プロンプトリストモード計数データをフィルタ処理するフィルタステップと、
を含む方法を実行する、非一時的コンピュータ可読記憶媒体。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US13/472,533 | 2012-05-16 | ||
| US13/472,533 US9291725B2 (en) | 2012-05-16 | 2012-05-16 | Random coincidence reduction in positron emission tomography using tangential time-of-flight mask |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013238603A JP2013238603A (ja) | 2013-11-28 |
| JP6188418B2 true JP6188418B2 (ja) | 2017-08-30 |
Family
ID=49582012
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013103728A Active JP6188418B2 (ja) | 2012-05-16 | 2013-05-16 | 偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US9291725B2 (ja) |
| JP (1) | JP6188418B2 (ja) |
| CN (1) | CN103547942B (ja) |
| WO (1) | WO2013172402A1 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104036137B (zh) * | 2014-06-11 | 2017-03-08 | 沈阳东软医疗系统有限公司 | 一种去除随机符合事件的方法及装置 |
| JP6535112B2 (ja) * | 2016-02-16 | 2019-06-26 | 日本電信電話株式会社 | マスク推定装置、マスク推定方法及びマスク推定プログラム |
| CN108109182B (zh) * | 2016-11-24 | 2021-08-24 | 上海东软医疗科技有限公司 | 一种pet图像重建方法和装置 |
| CN107464270B (zh) * | 2017-07-17 | 2020-08-11 | 东软医疗系统股份有限公司 | 一种图像重建方法和装置 |
| CN109188489B (zh) * | 2018-06-20 | 2020-04-28 | 江苏赛诺格兰医疗科技有限公司 | 随机事例的估计方法、装置及计算机可读存储介质 |
| WO2020059135A1 (ja) * | 2018-09-21 | 2020-03-26 | 株式会社島津製作所 | データ処理方法、プログラム、データ処理装置および陽電子放出断層撮像装置 |
| CN109613593B (zh) * | 2018-11-26 | 2020-10-20 | 北京永新医疗设备有限公司 | 伽马光子作用位置抽样方法、装置、电子设备及存储介质 |
| CN110215227B (zh) * | 2019-06-05 | 2022-10-14 | 上海联影医疗科技股份有限公司 | 时间窗设置方法、装置、计算机设备和存储介质 |
| US11576629B1 (en) | 2021-08-12 | 2023-02-14 | GE Precision Healthcare LLC | System and method for adaptive coincidence processing for high count rates |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6198104B1 (en) * | 1998-10-23 | 2001-03-06 | Adac Laboratories | Randoms correction using artificial trigger pulses in a gamma camera system |
| US6294788B1 (en) | 1999-08-18 | 2001-09-25 | Marconi Medical Systems, Inc. | Randoms correction in positron imaging |
| JP4695640B2 (ja) * | 2004-03-19 | 2011-06-08 | アドバンスト アプライド フィジックス ソリューションズ，インコーポレイテッド | 陽電子放出断面撮影におけるランダム同時発生を否認する方法および機器 |
| CN101163989B (zh) * | 2005-04-22 | 2013-04-10 | 皇家飞利浦电子股份有限公司 | 具有飞行时间性能的pet/mr扫描仪 |
| US7381959B2 (en) * | 2005-08-17 | 2008-06-03 | General Electric Company | Technique for reconstructing PET scan images |
| US7402807B2 (en) * | 2005-11-02 | 2008-07-22 | Siemens Medical Solutions Usa, Inc. | Method for reducing an electronic time coincidence window in positron emission tomography |
| WO2007054843A1 (en) | 2005-11-10 | 2007-05-18 | Koninklijke Philips Electronics, N.V. | Pet imaging using anatomic list mode mask |
| EP1977270B1 (en) * | 2006-01-09 | 2017-03-15 | Koninklijke Philips N.V. | Method of reconstructing time-of-flight pet images |
| JP4588042B2 (ja) * | 2007-03-16 | 2010-11-24 | 株式会社日立メディコ | 核医学診断装置 |
| US7737405B2 (en) | 2007-09-17 | 2010-06-15 | Siemens Medical Solutions Usa, Inc. | Time-of-flight (TOF) positron emission tomography (PET) reconstruction from time-truncated projection data |
-
2012
- 2012-05-16 US US13/472,533 patent/US9291725B2/en active Active
-
2013
- 2013-05-16 CN CN201380000774.8A patent/CN103547942B/zh active Active
- 2013-05-16 JP JP2013103728A patent/JP6188418B2/ja active Active
- 2013-05-16 WO PCT/JP2013/063626 patent/WO2013172402A1/ja active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| US9291725B2 (en) | 2016-03-22 |
| CN103547942A (zh) | 2014-01-29 |
| CN103547942B (zh) | 2017-09-26 |
| WO2013172402A1 (ja) | 2013-11-21 |
| JP2013238603A (ja) | 2013-11-28 |
| US20130311142A1 (en) | 2013-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6188418B2 (ja) | 偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体 | |
| Zhang et al. | Performance evaluation of the next generation solid-state digital photon counting PET/CT system | |
| JP5984246B2 (ja) | 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置に実行させるためのプログラム、及び陽電子放出コンピュータ断層撮影装置によって実行される方法 | |
| CN108474862B (zh) | 具有lu谱减除的能量校准 | |
| US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
| JP6125309B2 (ja) | 偶発同時計数推定方法及び偶発同時計数推定装置 | |
| JP7286383B2 (ja) | 陽電子放出撮像装置及び方法 | |
| US20090074281A1 (en) | Method for Energy Calculation and Pileup Determination for Continuously Sampled Nuclear Pulse Processing | |
| JP5925498B2 (ja) | 陽電子放出コンピュータ断層撮影装置、陽電子放出コンピュータ断層撮影装置で実行される方法及びプログラム | |
| Piccinelli et al. | Advances in single-photon emission computed tomography hardware and software | |
| US8520797B2 (en) | Medical imaging apparatus, control method, and computer program product | |
| JP6054050B2 (ja) | 核医学イメージング方法、核医学イメージング装置及び記憶媒体 | |
| JP7317586B2 (ja) | 医用画像処理装置、方法及びプログラム | |
| EP2902806B1 (en) | Nuclear medicine diagnostic device and medical data processing device | |
| US8519342B2 (en) | Nuclear medicine imaging apparatus, control method, and computer program product | |
| US20200085397A1 (en) | Attenuation coefficient image estimation method, attenuation coefficient image estimation program, and positron ct apparatus equipped with the same | |
| US11982779B2 (en) | Method and apparatus for guided pairing of multi-coincidences for time of flight positron emission tomography | |
| JP7247782B2 (ja) | 吸収係数画像推定方法、吸収係数画像推定プログラム、および、ポジトロンｃｔ装置 | |
| EP4075168A1 (en) | Nuclear medicine diagnosis device and nuclear medicine image data generation method | |
| WO2013168778A1 (ja) | 偶発同時計数推定方法及び偶発同時計数推定装置 | |
| US20220343566A1 (en) | Methods and systems for reconstructing a positron emission tomography image | |
| JP2024025712A (ja) | 散乱フラクションの決定方法、放射線診断装置及び放射線診断システム | |
| CN113253330A (zh) | 伽玛射线放射成像装置及能量校准方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20151102 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160222 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20160513 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20160929 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20161021 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20161206 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170202 |
|
| 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: 20170704 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170801 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6188418 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |