JP2022155147A - 核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム - Google Patents
核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム Download PDFInfo
- Publication number
- JP2022155147A JP2022155147A JP2021058499A JP2021058499A JP2022155147A JP 2022155147 A JP2022155147 A JP 2022155147A JP 2021058499 A JP2021058499 A JP 2021058499A JP 2021058499 A JP2021058499 A JP 2021058499A JP 2022155147 A JP2022155147 A JP 2022155147A
- Authority
- JP
- Japan
- Prior art keywords
- probability value
- scattering
- function
- electron density
- geometric information
- 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
Links
- 238000009206 nuclear medicine Methods 0.000 title claims abstract description 32
- 238000012545 processing Methods 0.000 title claims description 137
- 238000000034 method Methods 0.000 title claims description 27
- 238000003745 diagnosis Methods 0.000 title abstract description 14
- 230000005251 gamma ray Effects 0.000 claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims description 60
- 230000005855 radiation Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000002405 diagnostic procedure Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 253
- 238000012636 positron electron tomography Methods 0.000 description 58
- 238000004364 calculation method Methods 0.000 description 52
- 238000001514 detection method Methods 0.000 description 28
- 238000009826 distribution Methods 0.000 description 22
- 238000012937 correction Methods 0.000 description 13
- 238000003860 storage Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 238000013528 artificial neural network Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 238000004020 luminiscence type Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- XKUYOJZZLGFZTC-UHFFFAOYSA-K lanthanum(iii) bromide Chemical compound Br[La](Br)Br XKUYOJZZLGFZTC-UHFFFAOYSA-K 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 229910014323 Lanthanum(III) bromide Inorganic materials 0.000 description 1
- 229910052765 Lutetium Inorganic materials 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000012879 PET imaging Methods 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- OSGMVZPLTVJAFX-UHFFFAOYSA-N [Gd].[Lu] Chemical compound [Gd].[Lu] OSGMVZPLTVJAFX-UHFFFAOYSA-N 0.000 description 1
- ANDNPYOOQLLLIU-UHFFFAOYSA-N [Y].[Lu] Chemical compound [Y].[Lu] ANDNPYOOQLLLIU-UHFFFAOYSA-N 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- OHSVLFRHMCKCQY-UHFFFAOYSA-N lutetium atom Chemical compound [Lu] OHSVLFRHMCKCQY-UHFFFAOYSA-N 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 229940121896 radiopharmaceutical Drugs 0.000 description 1
- 239000012217 radiopharmaceutical Substances 0.000 description 1
- 230000002799 radiopharmaceutical effect Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- 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
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/04—Positioning of patients; Tiltable beds or the like
- A61B6/0407—Supports, e.g. tables or beds, for the body or parts of the body
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4258—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector for detecting non x-ray radiation, e.g. gamma radiation
-
- 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/452—Computed tomography involving suppression of scattered radiation or scatter correction
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Nuclear Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】画質を向上させること。【解決手段】実施形態に係る核医学診断装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。また、前記散乱イベントは多重散乱イベントを含む。【選択図】図1
Description
本明細書及び図面に開示の実施形態は、核医学診断装置、医用画像処理装置、核医学診断方法及びプログラムに関する。
現在のPET装置の多くは、被写体内での散乱を経ずに検出器に入射したガンマ線検出事象である直接入射事象を基に、被写体内の放射線分布を推測し画像再構成を行う。一般的に画像再構成のアルゴリズムは散乱線補正機能を有するが、多くの場合、散乱線補正の対象は、一回散乱までの散乱である。しかしながら、散乱線補正の対象が一回散乱までの場合、画質が低下する場合がある。
一方で、多重散乱を経た低エネルギーガンマ線検出事象を除去する散乱線補正を精度行うためには、高発光量のシンチレータと、当該高発光量シンチレータに対応する多数の受光セルが要求される場合がある。
加えて、例えばTOF(Time Of Flight)-PETにおいて高時間分解能を得ようとする場合には、シンチレータとしては、短時間発光シンチレータが要求され、また光検出器としては、個々のセルの素早い応答や短時間での回復が要求される場合がある。しかしながら、このような条件を満たすシンチレータは、高価なものになりうる。
本明細書及び図面の開示の実施形態が解決しようとする課題の一つは、画質を向上させることである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態に係る核医学診断装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。また、前記散乱イベントは多重散乱イベントを含む。
以下、図面を参照しながら、核医学診断装置、医用画像処理装置、核医学診断方法及びプログラムの実施形態について詳細に説明する。
(第1の実施形態)
図1は、実施形態に係る核医学診断装置としてのPET装置100の構成を示す図である。図1に示すように、実施形態に係るPET装置100は、架台装置10と、医用画像処理装置としてのコンソール装置20とを備える。架台装置10は、検出器1と、タイミング情報取得回路102と、天板103と、寝台104と、寝台駆動部106とを備える。
図1は、実施形態に係る核医学診断装置としてのPET装置100の構成を示す図である。図1に示すように、実施形態に係るPET装置100は、架台装置10と、医用画像処理装置としてのコンソール装置20とを備える。架台装置10は、検出器1と、タイミング情報取得回路102と、天板103と、寝台104と、寝台駆動部106とを備える。
検出器1は、被写体Pの陽電子から放出された消滅ガンマ線が発光体(シンチレータ)と相互作用することにより励起状態となった物質が再び基底状態に遷移する際に再放出される光であるシンチレーション光(蛍光)を検出することにより、放射線を検出する検出器である。検出器1は、被写体P内の陽電子から放出された消滅ガンマ線の放射線のエネルギー情報を検出する。検出器1は、被写体Pの周囲をリング状に取り囲むように複数配置され、例えば複数の検出器ブロックからなる。
検出器1の具体的な構成の一例としては、フォトンカウンティング方式、アンガー型の検出器であり、例えば、シンチレータと、光検出素子と、ライトガイドとを有する。すなわち、検出器1に含まれるピクセルのそれぞれは、シンチレータと、発生したシンチレーション光を検出する光検出素子とを有する。
シンチレータは、被写体P内の陽電子から放出されて入射した消滅ガンマ線をシンチレーション光(scintillation photons、optical photons)に変換し、出力する。シンチレータは、例えば、LaBr3(Lanthanum Bromide)、LYSO(Lutetium Yttrium Oxyorthosilicate)、LSO(Lutetium Oxyorthosilicate)、LGSO(Lutetium Gadolinium Oxyorthosilicate)等やBGO等の、TOF計測およびエネルギー計測に適するシンチレータ結晶によって形成され、例えば、2次元に配列される。
光検出素子としては、例えばSiPM(Silicon photomultiplier)や、光電子増倍管が用いられる。光電子増倍管は、シンチレーション光を受光して光電子を発生させる光電陰極、発生した光電子を加速する電場を与える多段のダイノード、及び、電子の流れ出し口である陽極を有し、シンチレータから出力されたシンチレーション光を増倍して電気信号に変換する。
また、架台装置10は、タイミング情報取得回路102により、検出器1の出力信号から計数情報を生成し、生成した計数情報を、コンソール装置20の記憶部130に格納する。なお、検出器1は、複数のブロックに区分けされ、タイミング情報取得回路102を備える。
タイミング情報取得回路102は、検出器1の出力信号をデジタルデータに変換し、計数情報を生成する。この計数情報には、消滅ガンマ線の検出位置、エネルギー値、及び検出時間が含まれる。例えば、タイミング情報取得回路102は、シンチレーション光を同じタイミングで電気信号に変換した複数の光検出素子を特定する。そして、タイミング情報取得回路102は、消滅ガンマ線が入射したシンチレータの位置を示すシンチレータ番号(P)を特定する。消滅ガンマ線が入射したシンチレータ位置を特定する手段は、各光検出素子の位置及び電気信号の強度に基づいて重心演算を行うことによって特定してもよい。また、シンチレータと光検出素子の各々の素子サイズが対応している場合には、出力が得られた光検出素子に対応するシンチレータを消滅ガンマ線が入射したシンチレータ位置として特定すればよい。
また、タイミング情報取得回路102は、各光検出素子から出力された電気信号の強度を積分計算することで、検出器1に入射した消滅ガンマ線のエネルギー値(E)を特定する。また、タイミング情報取得回路102は、検出器1によって消滅ガンマ線によるシンチレーション光が検出された検出時間(T)を特定する。なお、検出時間(T)は、絶対時刻であってもよいし、撮影開始時点からの経過時間であってもよい。このように、タイミング情報取得回路102は、シンチレータ番号(P)、エネルギー値(E)、及び検出時間(T)を含む計数情報を生成する。
なお、タイミング情報取得回路102は、例えば、CPU(Central Processing Unit)、GPU(Graphical Processing Unit)或いは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路により実現される。タイミング情報取得回路102は、タイミング情報取得部の一例である。
天板103は、被写体Pが載置されるベッドであり、寝台104の上に配置される。寝台駆動部106は、処理回路150の寝台制御機能105kによる制御の下、天板103を移動させる。例えば、寝台駆動部106は、天板103を移動させることで、被写体Pを架台装置10の撮影口内に移動させる。
コンソール装置20は、操作者によるPET装置100の操作を受け付け、PET画像の撮影を制御するとともに、架台装置10によって収集された計数情報を用いてPET画像を再構成する。図1に示すように、コンソール装置20は、処理回路150と、入力装置110と、ディスプレイ120と、記憶部130とを備える。なお、コンソール装置20が備える各部は、バスを介して接続される。処理回路150の詳細については後述する。
入力装置110は、PET装置100の操作者によって各種指示や各種設定の入力に用いられるマウスやキーボード等であり、入力された各種指示や各種設定を、処理回路150に転送する。例えば、入力装置110は、撮影開始指示の入力に用いられる。
ディスプレイ120は、操作者によって参照されるモニター等であり、処理回路150による制御の下、被写体の呼吸波形やPET画像を表示したり、操作者から各種指示や各種設定を受け付けるためのGUI(Graphical User Interface)を表示したりする。
記憶部130は、PET装置100において用いられる各種データを記憶する。記憶部130は、例えば、メモリで構成され、一例として、RAM(Random Access Memory)、フラッシュメモリ(flash memory)等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。記憶部130は、シンチレータ番号(P)、エネルギー値(E)、及び検出時間(T)が対応づけられた情報である計数情報、同時計数情報の通し番号であるコインシデンスNo.に計数情報の組が対応づけられた同時計数情報、再構成されたPET画像等を記憶する。
続いて、実施形態に係る背景について説明する。
現在のPET装置の多くは、基本的に被写体内での散乱を経ずに検出器に入射したガンマ線検出事象である直接入射事象を基に、被写体内の放射線分布を推測し画像再構成を行う。一般的に画像再構成のアルゴリズムは散乱線補正機能を有するが、多くの場合、散乱線補正の対象は、一回散乱までの散乱である。しかしながら、散乱線補正の対象が一回散乱までの場合、画質が低下する場合がある。
例えば、図2において、曲線30は、検出器に入射したガンマ線の信号強度をエネルギーの関数として示した図である。領域32は、直接入射事象及び単一散乱事象に対応するエネルギー領域である。これに対して、領域31は、多重散乱に対応するエネルギー領域である。多重散乱に対応するエネルギー領域である領域31の事象は、直接入射事象及び単一散乱事象に対応するエネルギー領域である領域32の事象と比較して、信号強度自体が小さく、また複数の種類の事象が重ね合わせられているので、散乱性補正の難易度が高くなる。
多重散乱を経た低エネルギーガンマ線検出事象を除去する散乱線補正を精度行うための一つの方法として、入射エネルギーを精度よく計測するための高発光量のシンチレータと、当該高発光量シンチレータに対応する多数の受光セルを用意する方法が考えられる。この場合、例えば、シンチレーション光子を電気信号に変換する光検出器としてSiPMを用いる場合、高発光量シンチレーションに対応するため、数千から数万の多数の受光セルが要求される。
加えて、例えばTOF(Time Of Flight)-PETにおいて高時間分解能を得ようとする場合、シンチレータとしては、短時間発光シンチレータが要求され、また光検出器としては、個々のセルの素早い応答や短時間での回復が要求される。加えて、光検出器の後段には信号に基づきエネルギーを算出するための電子回路が必要となり、多くの場合、当該装置に適したASICの開発を行う必要が生じる。
このように、TOF-PETにおいて、多重散乱を含む散乱線補正を行う場合、一つのPET装置100を構成するエネルギー計測系は、複雑なものとなり、また膨大なチャネル数からなるものになる。また、それぞれのチャネルの動作には固有のばらつきがあるため、これらの動作のばらつきを補正して計測精度を得るには、運用時に複雑な手順の構成を実施しなければならない。また、TOF-PETに要求される高発光量かつ短発光時間という条件を同時に満たすシンチレータは、高価なものになる。
従って、高発光量かつ短発光時間という条件を同時に満たすシンチレータを用いずとも、多重散乱事象を効果的に補正することができれば、エネルギー計測系によるPET装置の設計上の制約を大幅に減らすことができる。
実施形態に係る核医学診断装置は、かかる背景に鑑みて、新しい方法でPET装置における多重散乱事象を補正するものである。具体的には、実施形態に係る核医学診断装置は、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいて直接入射イベントに対応する第1の確率値及び多重散乱イベントを含む散乱イベントに対応する第2の確率値を推定し、推定したこれらの確率値と同時計数データとに基づいて、PET画像を再構成する。
これにより、画質を向上させることができ、または、PET装置のエネルギー計測系による設計上の制約を大幅に減らすことができる。例えば、現状より廉価なシンチレータを用いてPET装置を設計することが可能になり、またPET装置の性能を向上することができる。
続いて、図3及び図4を用いて、第1の実施形態に係るPET装置100について詳しく説明する。図3は、第1の実施形態に係るPET装置100の有する処理回路150の構成について説明した図である。また、図4は、実施形態に係るPET装置100が行う処理について説明したフローチャートである。
図3に示されるように、第1の実施形態において、処理回路150は、検出データ取得機能150a、同時計数情報データ生成機能150b、幾何情報取得機能150c、電子密度関数取得機能150d、直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150g、初期画像生成機能150h、再構成機能150i、システム制御機能150j、寝台制御機能150kを有する。なお、システム制御機能150j及び寝台制御機能150k以外の各機能については、図4を用いて後ほど詳しく説明する。
実施形態では、検出データ取得機能150a、同時計数情報データ生成機能150b、幾何情報取得機能150c、電子密度関数取得機能150d、直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150g、初期画像生成機能150h、再構成機能150i、システム制御機能150j、寝台制御機能150kにて行われる各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶部130へ記憶されている。処理回路150はプログラムを記憶部130から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路150は、図3の処理回路105内に示された各機能を有することになる。
なお、図3においては単一の処理回路150にて、検出データ取得機能150a、同時計数情報データ生成機能150b、幾何情報取得機能150c、電子密度関数取得機能150d、直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150g、初期画像生成機能150h、再構成機能150i、システム制御機能150j、寝台制御機能150kにて行われる処理機能が実現されるものとして説明するが、複数の独立したプロセッサを組み合わせて処理回路150を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。換言すると、上述のそれぞれの機能がプログラムとして構成され、1つの処理回路150が各プログラムを実行する場合であってもよい。別の例として、特定の機能が専用の独立したプログラム実行回路に実装される場合であってもよい。
上記説明において用いた「プロセッサ」という文言は、例えば、CPU(Central Processing Unit)、GPU(Graphical Processing Unit)或いは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサは記憶部130に保存されたプログラムを読み出し実行することで機能を実現する。
なお、図3の処理回路150に限られず、後述する図5、図7、図9の処理回路150についても、同様の実装が可能である。
なお、図3において、検出データ取得機能150a、同時計数情報データ生成機能150b及び初期画像生成機能150hは、第1の取得部の一例である。また、幾何情報取得機能150c、電子密度関数取得機能150dは、第2の取得部の一例である。直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150gは、推定部の一例である。再構成機能150iは、再構成部の一例である。システム制御機能150j及び寝台制御機能150kは、制御部の一例である。
処理回路105は、システム制御機能105cにより、架台装置10及びコンソール装置20を制御することによって、PET装置100の全体制御を行う。例えば、処理回路105は、システム制御部105cにより、PET装置100における撮影を制御する。
処理回路105は、寝台制御機能105kにより、寝台駆動部105を制御する。
続いて、図4を用いて、処理回路150が有する各機能を説明しながら、第1の実施形態に係るPET装置100が実行する処理の流れについて説明する。ステップS100において、処理回路150は、電子密度関数取得機能150dにより、被写体の電子密度関数を取得する。一例として、処理回路150は、電子密度関数取得機能150dにより、被写体のCT画像を、電子密度関数として取得する。一例として、処理回路150は、電子密度関数取得機能150dにより、図示しないCT装置により被写体を撮影した結果得られた被写体のCT画像を被写体の電子密度関数として取得する。また、別の例として、実施形態に係る核医学診断装置は、CT装置としての機能を有するPET-CT装置であり、当該PET-CT装置によるCT撮影により被写体を撮影した結果得られた被写体のCT画像を、処理回路150は、電子密度関数取得機能150dにより取得してもよい。
なお、被写体の電子密度関数を表すデータは、CT撮像により得られたデータに限られない。一例として、被写体の電子密度関数を表すデータは、磁気共鳴画像であってもよい。例えば、処理回路150は、電子密度関数取得機能150dにより、図示しないMRI装置により被写体を撮影し、例えばセグメンテーション処理を行った磁気共鳴画像を、被写体の電子密度関数として取得する。
このように、実施形態に係るPET装置100は、電子密度関数取得機能150dにより、ステップS100において、被写体の電子密度関数を、PET撮像により得られたデータとは別に取得する。これにより、処理回路150は、ガンマ線の多重散乱を推定するための基礎となる情報を取得することができ、既存の方法と比較して、以降のステップにおける多重散乱の推定の精度が向上する。
続いて、ステップS110において、処理回路150は、幾何情報取得機能150cにより、ガンマ線検出器である検出器1の幾何情報を取得する。ここで、検出器1の幾何情報とは、検出器1のそれぞれの相対的位置関係の情報及び、検出器1と被写体Pとの位置関係の情報を指す。一例として、処理回路150は、幾何情報取得機能150cにより、記憶部130にあらかじめ記憶されている検出器1のそれぞれの相対的位置関係の情報を、記憶部130から取得することにより、検出器1のそれぞれの相対的位置関係の情報を取得する。また、処理回路150は、幾何情報取得機能150cにより、寝台駆動部106を通じて被写体Pの位置を取得することにより、検出器1と被写体Pとの位置関係の情報を取得する。なお、例えばステップS110の処理と、ステップS100の処理とは、実行の順序は逆でもよい。また、後述の図10においても、同様に、ステップS110の処理と、ステップS100の処理とは、実行の順序は逆でもよい。
続いて、ステップS120において、処理回路150は、ステップS100で電子密度関数取得機能150dにより取得した被写体の電子密度と、ステップS110で幾何情報取得機能150cにより取得した検出器1の幾何情報とに基づいて、直接入射確率計算機能150eにより、ガンマ線が検出器1に直接入射する確率である第1の確率値を計算する。ここで、ガンマ線が検出器1に直接入射するとは、生成されたガンマ線の対が検出器1に到達するまでに途中で散乱された散乱回数が0回であることを意味する。
なお、ステップS120においては、処理回路150は、検出器視野内に設定した各ボクセルから、単位放射能が放射されたと仮定して計算を行う。なお、各ボクセルに対する実際の放射能の量は、ステップS150において推定される。
処理回路150は、直接入射確率計算機能150eにより、被写体の電子密度関数及びガンマ線検出器の幾何情報の少なくとも一方に基づいて、検出器視野内のボクセルに単位放射能があると仮定した場合の被写体における直接入射イベントに対応する第1の確率値を推定する。一例として、処理回路150は、直接入射確率計算機能150eにより、被写体の電子密度関数及びガンマ線検出器の幾何情報の少なくとも一方に基づいて、検出器視野内のボクセルに単位放射能があると仮定した場合の検出器1のそれぞれへの直接入射確率である第1の確率値を、例えば放射輸送方程式、モンテカルロシミュレーション、ニューラルネットワーク等により、算出する。
また、処理回路150は、ステップS100で電子密度関数取得機能150dにより取得した被写体の電子密度と、ステップS110で幾何情報取得機能150cにより取得した検出器1の幾何情報とに基づいて、散乱入射確率計算機能150fにより、ガンマ線が散乱を経て検出器1に入射する確率を計算する。ここで、散乱を経るとは、単一散乱及び多重散乱を経た入射の両方を含む。
換言すると、処理回路150は、散乱入射確率計算機能150fにより、被写体の電子密度関数及びガンマ線検出器の幾何情報の少なくとも一方に基づいて、検出器視野内のボクセルに単位放射能があると仮定した場合の、被写体での、多重散乱イベントを含む散乱イベントに対応する第2の確率値を推定する。具体的には、処理回路150は、散乱入射確率計算機能150fにより、検出器視野内のボクセルに単位放射能があると仮定した場合の検出器1の散乱入射確率を、例えば放射輸送方程式、モンテカルロシミュレーション、ニューラルネットワーク等により、算出する。
続いて、処理回路150は、システム行列生成機能150gにより、直接入射確率計算機能150eにより算出された直接入射確率と、散乱入射確率計算機能150fにより算出された散乱入射確率とに基づいて、システム行列Hを生成する。ここで、システム行列Hは、iをi番目のLOR(Line Of Response)とし、jをj番目のボクセルとしたときに、その要素がHijで表される行列である。
なお、第1の実施形態及び第2の実施形態では、散乱の効果は、散乱行列として、システム行列に取り込まれた形で表現されるが、第3の実施形態では、散乱の効果は、例えばShifted Poission Model等を用いて、システム行列とは別個の散乱項として表現されるという違いがある。
ステップS130において、処理回路150は、タイミング情報取得回路102から、計数情報を取得する。具体的には、処理回路150は、検出データ取得機能150aにより、タイミング情報取得回路102から、ガンマ線の検出データを取得する。ここで、処理回路150が検出データ取得機能150aによりタイミング情報取得回路102から取得するガンマ線の検出データは、例えば、シンチレータ番号(P)、エネルギー値(E)、及び検出時間(T)を含む計数情報である。かかる計数情報には、直接入射によるデータおよび散乱入射によるデータの両方が含まれる。
続いて、処理回路150は、同時計数情報データ生成機能150bにより、検出データ取得機能150aにより取得した計数情報に基づいて、同時計数情報データを生成する。
すなわち、ステップS130において、処理回路150は、検出データ取得機能150a及び同時計数情報データ生成機能150bにより、ガンマ線検出器1への直接入射イベントと被写体Pでの散乱イベントとを含む同時計数データを取得する。
続いて、ステップS140において、処理回路150は、初期画像生成機能150hにより、被写体放射能分布の初期分布を設定する。一例として、処理回路150は、初期画像生成機能150bにより、複数の被写体のPET画像に対して平均操作を行って平均的な放射能分布を作成し、当該平均的な放射能分布を、初期分布として設定する。また、別の例として、処理回路150は、初期画像生成機能150hより、ステップS100で取得した被写体のCT画像において、電子密度が一定の閾値以上である場所について、一定の放射能が分布し、その他の場所について、放射能が分布しないとする初期分布を設定する。また、別の例として、処理回路150は、初期画像生成機能150hより、ステップS100で取得した被写体のCT画像から算出した電子密度分布を、被写体放射能分布の初期分布とする。
続いて、ステップS150において、処理回路150は、画像再構成機能150iにより、ステップS140で設定した被写体放射能の初期分布と、ステップS120で生成したシステム行列と、ステップS130で取得した同時計数データとに基づいて、被写体の放射能分布を推定し、PET画像を再構成し再構成画像を生成する。換言すると、処理回路150は、画像再構成機能150iにより、ステップS140で設定した被写体放射能の初期分布と、ステップS120で推定された第1の確率値と第2の確率値と、ステップS130で取得した同時計数データとに基づいてPET画像を再構成する。一例として、処理回路150は、画像再構成機能150iにより、ステップS140で設定した被写体放射能の初期分布と、ステップS120で推定された第1の確率値と第2の確率値とに基づくシステム行列と、ステップS130で取得した同時計数データとに基づいてPET画像を再構成する。一例として、処理回路150は、画像再構成機能150iにより、ステップS140で設定した被写体放射能の初期分布と、ステップS120で生成したシステム行列とに基づいて、ステップS130で取得した同時計数情報を再現する被写体の放射能分布を、最尤法、ニューラルネットワーク等を用いて推定し、再構成画像を生成する。
続いて、ステップS160において、処理回路150は、システム制御機能150jより、ステップS150で生成した再構成画像を、ディスプレイ120に表示させる。
以上のように、第1の実施形態では、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいて直接入射イベントに対応する第1の確率値及び多重散乱イベントを含む散乱イベントに対応する第2の確率値を推定し、推定したこれらの確率値と同時計数データとに基づいて、PET画像を再構成する。第1の実施形態によると、散乱線が有する被写体情報を生かした撮像・再構成を行うPET装置を構成することができる。また、被写体の電子密度関数やガンマ線検出器の幾何情報があることで、検出ガンマ線のエネルギーを計測するためにPETに実装させる高価なエネルギー計測回路が必須でなくなり、設計上の自由度が広がる。
一例として、従来と比較して廉価なシンチレータを用いて高画質なPET画像を得ることができる核医学診断装置を製作することが可能となる。また、別の例として、発光量は少ないが、ガンマ線入射に対する応答が非常に速い、などのTOF機能に有利な特性を持つシンチレータを、実施形態に係る方法では採用することが可能となり、シンチレータの選択の幅が広がる。
加えて、第1の実施形態に係る方法では、生成されたシステム行列が、散乱過程と減弱過程を整合的に扱うため、例えば従来法のように、散乱過程と減弱過程との間の相対的な強度を手作業でチューニングする操作が不要になる。また、第1の実施形態に係る方法では、多重散乱を含めた全事象を対象とした画像再構成であるので、統計ノイズを抑制した画像を得ることができる。加えて、第1の実施形態に係る方法では、放射性薬剤の低投与量化による低被ばく撮像も期待できる。
(第2の実施形態)
続いて、図5及び図6を用いて、第2の実施形態について説明する。第1の実施形態では、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいてシステム行列を推定し、推定したシステム行列に基づいて、PET画像を再構成する場合を説明した。第2の実施形態も、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいてシステム行列を求める点では共通するが、詳細において、第1の実施形態と異なる手続きでシステム行列を求める。すなわち、システム行列のうち、直接入射の寄与部分は、被写体の電子密度関数の形状によらない部分であるので、システム行列のうち、直接入射の寄与部分に関しては、処理回路は、別途あらかじめ計算を行っておき、散乱入射の寄与部分を後ほど加算することでシステム行列を生成する。これにより、システム行列の計算を簡略化することができる。
続いて、図5及び図6を用いて、第2の実施形態について説明する。第1の実施形態では、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいてシステム行列を推定し、推定したシステム行列に基づいて、PET画像を再構成する場合を説明した。第2の実施形態も、被写体の電子密度関数及びガンマ線検出器の幾何情報を取得し、取得したこれらのデータに基づいてシステム行列を求める点では共通するが、詳細において、第1の実施形態と異なる手続きでシステム行列を求める。すなわち、システム行列のうち、直接入射の寄与部分は、被写体の電子密度関数の形状によらない部分であるので、システム行列のうち、直接入射の寄与部分に関しては、処理回路は、別途あらかじめ計算を行っておき、散乱入射の寄与部分を後ほど加算することでシステム行列を生成する。これにより、システム行列の計算を簡略化することができる。
図5に、第2の実施形態における処理回路150の構成が示されている。図5において、直接入射確率計算機能150e以外の機能については、第1の実施形態と同様であるので詳しい説明は省略する。直接入射確率計算機能150eについて、第1の実施形態においては、処理回路150は、検出器の幾何情報と被写体の電子密度関数の両方を用いて直接入射確率を推定したが、第2の実施形態では、処理回路150は、検出器の幾何情報のみを用いて直接入射確率を推定する。
なお、図5において、図3と同様に、検出データ取得機能150a、同時計数情報データ生成機能150b及び初期画像生成機能150hは、第1の取得部の一例である。また、幾何情報取得機能150c、電子密度関数取得機能150dは、第2の取得部の一例である。直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150gは、推定部の一例である。再構成機能150iは、再構成部の一例である。システム制御機能150j及び寝台制御機能150kは、制御部の一例である。
図6に、第2の実施形態に係るPET装置100が実行する処理の流れが示されている。ここで、ステップS70、ステップS90及びステップS120以外の処理については、処理の内容は第1の実施形態と共通であるから、繰り返しての説明は省略する。
はじめに、ステップS70において、処理回路150は、幾何情報取得機能150cにより、検出器幾何情報を取得する。ステップS70の処理は、第1の実施形態におけるステップS110と同様の処理であるが、第2の実施形態においては、このタイミングで処理が行われる。
続いて、ステップS90において、処理回路150は、直接入射確率計算機能150eにより、ステップS70において幾何情報取得機能150cにより取得された検出器幾何情報に基づいて、第1の確率値である直接入射確率を算出する。この際、処理回路150は、直接入射確率計算機能150eにより、被写体の電子密度関数を用いずに第1の確率値である直接入射確率を算出する。続いて、ステップS90において、処理回路150は、システム行列生成機能150gにより、算出した第1の確率値に基づいて、システム行列のうち、第1の確率値である直接入射確率に対応する寄与を算出する。すなわち、処理回路150は、システム行列生成機能150gにより、検出器幾何情報に基づいて、直接入射確率のシステム行列を生成する。
続いて、ステップS100及びステップS110において、第1の実施形態と同様の処理を行い、処理回路150は、電子密度関数取得機能150dにより、被写体の電子密度関数を取得する。
続いて、ステップS120において、処理回路150は、散乱入射確率計算機能150eにより、第1の実施形態と同様の処理を行い、第2の確率値である散乱入射確率を計算する。続いて、処理回路150は、システム行列生成機能150gにより、ステップS90であらかじめ算出された第1の確率値である直接入射確率と、第1の実施形態と同様にステップS120で算出された第2の確率値である散乱入射確率とに基づいて、システム行列を生成する。その後、ステップS130からステップS160において、第1の実施形態と同様の処理を行う。
このように、第2の実施形態では、処理回路150は、再構成機能150iにより、被写体の電子密度関数の取得に先立ってあらかじめ算出された第1の確率値に基づくシステム行列に基づいてPET画像を再構成する。これにより、被写体ごとにシステム行列のうち直接入射確率を算出する必要がなくなり、システム行列の計算を簡略化することができる。これにより、再構成処理における計算負荷を軽減することができる。
(第3の実施形態)
第1及び第2の実施形態では、散乱の効果が、散乱行列として、システム行列の値に組み込まれて表現される場合について説明した。しかしながら、実施形態はこれに限られない。第3の実施形態では、例えばShiftedPoissonModel等が採用され、散乱の効果が、システム行列に対する補正項である散乱項として表現される場合について説明する。
第1及び第2の実施形態では、散乱の効果が、散乱行列として、システム行列の値に組み込まれて表現される場合について説明した。しかしながら、実施形態はこれに限られない。第3の実施形態では、例えばShiftedPoissonModel等が採用され、散乱の効果が、システム行列に対する補正項である散乱項として表現される場合について説明する。
図7に、第3の実施形態における処理回路150の構成が示されている。図7において、散乱項計算機能150l以外の機能については、第2の実施形態と同様であるので詳しい説明は省略する。
ここで、jをj番目のLORを示す添え字として、処理回路150が散乱項計算機能150lにより計算する散乱項Sjは、例えば以下の式(1)で与えられる。
ここで、i及びmはそれぞれi番目、m番目のボクセルを表す添え字であり、kは、画像再構成において、k番目のイタレーションを表す添え字であり、λi
kはk番目のイタレーションにおけるi番目のボクセルにおける放射能の値の推定値を表し、Hijはi番目のボクセルとj番目のLORとの間のシステム行列であり、gjは、j番目のLORにおける検出器が検出するイベント数である。すなわち、処理回路150が散乱項計算機能150lにより計算する散乱項Sjは、散乱線に関するパラメータであって、PET画像の再構成時に、各LORごとに導入されるシステム行列に対する補正項である。
なお、図7において、検出データ取得機能150a、同時計数情報データ生成機能150b及び初期画像生成機能150hは、第1の取得部の一例である。また、幾何情報取得機能150c、電子密度関数取得機能150dは、第2の取得部の一例である。直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150gは、推定部の一例である。散乱項計算機能150iは、同定部の一例である。再構成機能150iは、再構成部の一例である。システム制御機能150j及び寝台制御機能150kは、制御部の一例である。
図8に、第3の実施形態に係るPET装置100が実行する処理の流れが示されている。ここで、ステップS90、ステップS125及びステップS150以外の処理については、処理の内容は第2の実施形態と共通であるから、繰り返しての説明は省略する。
第2の実施形態と同様、ステップS90において、処理回路150は、直接入射確率計算機能150eにより、ステップS70において幾何情報取得機能150cにより取得された検出器幾何情報に基づいて、第1の確率値である直接入射確率を推定し、システム行列生成機能150gにより、検出器幾何情報に基づいて、システム行列を生成する。
ここで、第3の実施形態では、散乱がシステム行列の成分に組み込まれた第2の実施形態とは異なり、散乱の効果は、システム行列に対する補正項である散乱項として表現されるので、システム行列の成分とはならない。従って、ステップS90において、処理回路150は、システム行列生成機能150gにより、すでにシステム行列を生成する。
ステップS125において、処理回路150は、ステップS100で電子密度関数取得機能150dにより取得した被写体の電子密度と、ステップS70で幾何情報取得機能150cにより取得した検出器1の幾何情報とに基づいて、散乱入射確率計算機能150fにより、ガンマ線が単一散乱または多重散乱を経て検出器に入射する確率である第2の確率を計算する。
換言すると、処理回路150は、散乱入射確率計算機能150fにより、被写体の電子密度関数とガンマ線検出器の幾何情報とに基づいて、被写体での、多重散乱イベントを含む散乱イベントに対応する第2の確率値を推定する。具体的には、処理回路150は、散乱入射確率計算機能150fにより、多重散乱イベントを含む散乱イベントによる散乱入射確率を、例えば放射輸送方程式、モンテカルロシミュレーション、ニューラルネットワーク等により、算出する。続いて、処理回路150は、散乱項計算機能150lにより、ステップS70で取得した検出器幾何情報と、ステップS90で推定された第1の確率値である直接入射確率と、ステップS100で取得した被写体の電子密度関数と、被写体の放射能分布の初期推定値とに基づいて散乱計算を行い、散乱線に関するパラメータである散乱項を同定する。ここで、散乱項に関するパラメータである散乱項は、多重散乱項を含む。
また、ステップS150では、処理回路150は、再構成機能150iにより、例えば式(1)の右辺を評価することにより、式(1)の左辺の値を推定し、ステップS140で設定された被写体の放射能分布の初期推定値と、ステップ90で生成されたシステム行列と、ステップS125で生成された散乱線に関するパラメータである散乱項とに基づいて、被写体の放射能分布を推定し、PET画像を再構成する。
なお、例えばステップS100において、処理回路150が電子密度関数150dにより取得するデータは、CT画像に限られず、例えば磁気共鳴画像であってもよい。一例として、処理回路150は、図示しないセグメンテーション機能により、磁気共鳴画像に対してセグメンテーション処理を行って部位の自動抽出/判定を行う。続いて、処理回路150は、当該磁気共鳴画像について、抽出/判定された部位ごとに、磁気共鳴信号値からCT値への換算を行い、CT画像に類似する信号値を有する磁気共鳴画像を生成する。処理回路150は、電子密度関数150dにより、当該磁気共鳴画像を、被写体の電子密度を示す画像として取得する。
また、処理回路150は、システム行列生成機能150gにより、ステップS100で取得された被写体のCT画像に対して、例えばセマンティックセグメーション処理を行って、散乱領域を特定し、特定した散乱領域に基づいて第2の確率値を推定してもよい。これにより、画像再構成に伴う計算負荷を軽減することができる。
以上述べたように、第3の実施形態においては、処理回路150は、散乱項計算機能150iにより、検出器幾何情報と電子密度関数とに基づいて散乱線に関するパラメータである散乱項を、多重散乱を含めて算出し、別途算出されたシステム行列と散乱項とを用いて、PET画像を再構成する。このように、多重散乱を含む項を散乱項としてシステム行列と別に算出し再構成を行うことで、多重散乱を含む散乱を含んだデータに対して、精度よくPET再構成を行うことが可能となる。
(第4の実施形態)
続いて、第4の実施形態について説明する。同時計数データを生成する際、ガンマ線検出器からの出力結果に対してローパスフィルタや波高値弁別器(ディスクリメータ)を用いてノイズ除去して同時計数データを生成する場合がある。ここで、例えば波高値弁別器が波高値カットオフに用いたパラメータの値は、検出器ごとに異なる場合がある。第4の実施形態では、検出器ごとに波高値カットオフに用いたパラメータの値を用いて画像再構成を行うことで、再構成されるPET画像のさらなる画質を向上させることができる。
続いて、第4の実施形態について説明する。同時計数データを生成する際、ガンマ線検出器からの出力結果に対してローパスフィルタや波高値弁別器(ディスクリメータ)を用いてノイズ除去して同時計数データを生成する場合がある。ここで、例えば波高値弁別器が波高値カットオフに用いたパラメータの値は、検出器ごとに異なる場合がある。第4の実施形態では、検出器ごとに波高値カットオフに用いたパラメータの値を用いて画像再構成を行うことで、再構成されるPET画像のさらなる画質を向上させることができる。
図9に、第4の実施形態における処理回路150の構成が示されている。図9において、カットオフ処理機能150n及びカットオフパラメータ設定機能150m以外の機能については、第1の実施形態と同様であるので詳しい説明は省略する。
カットオフ処理機能150nは、検出器データ取得機能150aにより得られた計数データに対して、ローパスフィルタまたは波高弁別器(ディスクリミネータ)によりノイズカット処理を行う機能である。一例として、処理回路150は、カットオフ処理機能150nにより、例えばガンマ線検出器データ取得機能150aにより得られた計数データに対して、定められた波高値カットオフパラメータ未満の計数データを除去する。例えば、検出器データ取得機能150aにより得られた計数データが定められた波高値カットオフパラメータ未満であった場合、出力される計数データは0となる。
また、処理回路150は、カットオフパラメータ設定機能150mにより、かかるカットオフパラメータの値を設定する。
なお、図9において、検出データ取得機能150a、同時計数情報データ生成機能150b、カットオフパラメータ設定機能150m、カットオフ処理機能150n及び初期画像生成機能150hは、第1の取得部の一例である。また、幾何情報取得機能150c、電子密度関数取得機能150dは、第2の取得部の一例である。直接入射確率計算機能150e、散乱入射確率計算機能150f、システム行列生成機能150gは、推定部の一例である。再構成機能150iは、再構成部の一例である。システム制御機能150j及び寝台制御機能150kは、制御部の一例である。
図10に、第4の実施形態に係るPET装置100が実行する処理の流れが示されている。ここで、ステップS80、ステップS85、ステップS120、ステップS130及びステップS150以外の処理については、処理の内容は第1の実施形態と共通であるから、繰り返しての説明は省略する。
ステップS80において、処理回路150は、カットオフパラメータ設定機能150mにより、波高値カットオフパラメータを設定する。また、ステップS130において、処理回路150は、第1の実施形態と同様に、検出データ取得機能150aにより、タイミング情報取得回路102から、ガンマ線の検出データを計数情報として取得する。続いて、処理回路150は、同時計数情報データ生成機能150bにより、ガンマ線検出データ取得機能150aにより取得した計数情報と、ステップS80で設定された波高値パラメータとに基づいて、同時計数情報データを生成する。すなわち、処理回路150は、同時計数情報データ生成機能150bにより、ガンマ線検出器からの出力結果に対して波高値カットオフを適用して同時計数データを生成する。
一方、ステップS85において、処理回路150は、ステップS80で設定された波高値カットオフパラメータを、直接入射確率計算機能150e及び散乱入射確率計算機能150fなどの散乱入射計算ユニットに通知する。
また、ステップS120において、処理回路150が直接入射確率計算機能150eにより直接入射確率を算出し、散乱入射確率計算機能150fに基づいて、散乱入射確率を算出するが、処理回路150は、直接入射確率計算機能150eにより、ステップS85で通知された波高値パラメータの値に基づいて、直接入射確率を第1の確率値として推定する。また、処理回路150は、散乱入射確率計算機能150fにより、ステップS85で通知された波高値パラメータの値に基づいて、散乱入射確率を第2の確率値として推定する。
また、処理回路150は、システム行列生成機能150gにより、直接入射確率計算機能150eにより算出された直接入射確率と、散乱入射確率計算機能150fにより算出された散乱入射確率とに基づいて、システム行列を生成するが、この際、ステップS85で通知された波高値パラメータの値に基づいてシステム行列を生成する。すなわち、処理回路150は、検出器毎のカットオフパラメータの差異を、ステップS120において生成するシステム行列に反映させる。
また、ステップS150において、処理回路150は、画像再構成機能150iにより、ステップS140で設定した被写体放射能の初期分布と、ステップS120で生成したシステム行列と、ステップS130で取得した同時計数データと、カットオフパラメータ設定機能により設定された波高値パラメータとに基づいて、被写体の放射能分布を推定し、PET画像を再構成し再構成画像を生成する。
以上のように、第4の実施形態では、検出器ごとに波高値カットオフに用いたパラメータの値を用いて散乱の評価及び画像再構成を行うことで、再構成されるPET画像のさらなる画質を向上させることができる。
以上説明した少なくとも1つの実施形態によれば、画質を向上させることができる。
以上の実施形態に関し、発明の一側面および選択的な特徴として以下の付記を開示する。
(付記1)
本発明の一つの側面において提供される核医学診断装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。前記散乱イベントは多重散乱イベントを含む。
本発明の一つの側面において提供される核医学診断装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。前記散乱イベントは多重散乱イベントを含む。
(付記2)
前記推定部は、放射輸送方程式、モンテカルロシミュレーションまたはニューラルネットワークを用いて、前記第1の確率値を推定してもよい。
前記推定部は、放射輸送方程式、モンテカルロシミュレーションまたはニューラルネットワークを用いて、前記第1の確率値を推定してもよい。
(付記3)
前記推定部は、放射輸送方程式、モンテカルロシミュレーションまたはニューラルネットワークを用いて、前記第2の確率値を推定してもよい。
前記推定部は、放射輸送方程式、モンテカルロシミュレーションまたはニューラルネットワークを用いて、前記第2の確率値を推定してもよい。
(付記4)
再構成部は、前記第1の確率値に基づくシステム行列に基づいて前記PET画像を再構成してもよい。
再構成部は、前記第1の確率値に基づくシステム行列に基づいて前記PET画像を再構成してもよい。
(付記5)
再構成部は、前記第1の確率値と前記第2の確率値とに基づくシステム行列に基づいて前記PET画像を再構成してもよい。
再構成部は、前記第1の確率値と前記第2の確率値とに基づくシステム行列に基づいて前記PET画像を再構成してもよい。
(付記6)
推定部は、前記幾何情報に基づいて前記第1の確率値を推定し、前記幾何情報と前記電子密度関数とに基づいて第2の確率値を推定してもよい。
推定部は、前記幾何情報に基づいて前記第1の確率値を推定し、前記幾何情報と前記電子密度関数とに基づいて第2の確率値を推定してもよい。
(付記7)
前記幾何情報と前記電子密度関数とに基づいて散乱線に関するパラメータを同定する同定部を更に備え、
再構成部は、前記パラメータに基づいて前記PET画像を再構成してもよい。
前記幾何情報と前記電子密度関数とに基づいて散乱線に関するパラメータを同定する同定部を更に備え、
再構成部は、前記パラメータに基づいて前記PET画像を再構成してもよい。
(付記8)
前記再構成部は、Shifted Poisson Modelを用いて前記PET画像を再構成してもよい。
前記再構成部は、Shifted Poisson Modelを用いて前記PET画像を再構成してもよい。
(付記9)
第1の取得部は、前記ガンマ線検出器からの出力結果に対して波高値カットオフを適用して前記同時計数データを生成し、
推定部は、前記波高値カットオフに用いたパラメータに基づいて、前記第1の確率値及び前記第2の確率値を推定してもよい。
第1の取得部は、前記ガンマ線検出器からの出力結果に対して波高値カットオフを適用して前記同時計数データを生成し、
推定部は、前記波高値カットオフに用いたパラメータに基づいて、前記第1の確率値及び前記第2の確率値を推定してもよい。
(付記10)
推定部は、セグメンテーション処理を行って散乱領域を特定し、特定した散乱領域に基づいて前記第2の確率値を推定してもよい。セグメンテーション処理は、セマンティックセグメンテーション処理であってもよい。
推定部は、セグメンテーション処理を行って散乱領域を特定し、特定した散乱領域に基づいて前記第2の確率値を推定してもよい。セグメンテーション処理は、セマンティックセグメンテーション処理であってもよい。
(付記11)
電子密度関数は、CT画像であってもよい。
電子密度関数は、CT画像であってもよい。
(付記12)
電子密度関数は、磁気共鳴画像をもとに生成したものであってもよい。
電子密度関数は、磁気共鳴画像をもとに生成したものであってもよい。
(付記13)
本発明の一つの側面において提供される医用画像処理装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。前記散乱イベントは多重散乱イベントを含む。
本発明の一つの側面において提供される医用画像処理装置は、第1の取得部と、第2の取得部と、推定部と、再構成部とを備える。第1の取得部は、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する。第2の取得部は、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する。推定部は、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定する。再構成部は、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する。前記散乱イベントは多重散乱イベントを含む。
(付記14)
本発明の一つの側面において提供される核医学診断方法は、核医学診断装置により実行される核医学診断方法であって、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得し、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定し、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成することを含み、前記散乱イベントは多重散乱イベントを含む。
本発明の一つの側面において提供される核医学診断方法は、核医学診断装置により実行される核医学診断方法であって、ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得し、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定し、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成することを含み、前記散乱イベントは多重散乱イベントを含む。
(付記15)
本発明の一つの側面において提供されるプログラムは、コンピュータに、ガンマ線検出器への直接入射イベントと被写体での多重散乱イベントを含む散乱イベントとを含む同時計数データを取得し、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定し、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する処理を実行させる。
本発明の一つの側面において提供されるプログラムは、コンピュータに、ガンマ線検出器への直接入射イベントと被写体での多重散乱イベントを含む散乱イベントとを含む同時計数データを取得し、前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントとに対応する第2の確率値とを推定し、前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する処理を実行させる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
150 処理回路
150a 検出データ取得機能
150b 同時計数情報生成機能
150c 幾何情報取得機能
150d 電子密度関数取得機能
150e 直接入射確率計算機能
150f 散乱入射確率計算機能
150g システム行列生成機能
150h 初期画像生成機能
150i 再構成機能
150k 寝台制御機能
150a 検出データ取得機能
150b 同時計数情報生成機能
150c 幾何情報取得機能
150d 電子密度関数取得機能
150e 直接入射確率計算機能
150f 散乱入射確率計算機能
150g システム行列生成機能
150h 初期画像生成機能
150i 再構成機能
150k 寝台制御機能
Claims (10)
- ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する第1の取得部と、
前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する第2の取得部と、
前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントに対応する第2の確率値とを推定する推定部と、
前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する再構成部とを備え
前記散乱イベントは多重散乱イベントを含む、核医学診断装置。 - 前記再構成部は、前記第1の確率値に基づくシステム行列に基づいて前記PET画像を再構成する、請求項1に記載の核医学診断装置。
- 前記再構成部は、前記第1の確率値と前記第2の確率値とに基づくシステム行列に基づいて前記PET画像を再構成する、請求項1に記載の核医学診断装置。
- 前記推定部は、前記幾何情報に基づいて前記第1の確率値を推定し、前記幾何情報と前記電子密度関数とに基づいて第2の確率値を推定する、請求項1乃至3のいずれか1つに記載の核医学診断装置。
- 前記幾何情報と前記電子密度関数とに基づいて散乱線に関するパラメータを同定する同定部を更に備え、
前記再構成部は、前記パラメータに基づいて前記PET画像を再構成する、請求項2に記載の核医学診断装置。 - 前記第1の取得部は、前記ガンマ線検出器からの出力結果に対して波高値カットオフを適用して前記同時計数データを生成し、
前記推定部は、前記波高値カットオフに用いたパラメータに基づいて、前記第1の確率値及び前記第2の確率値を推定する、請求項1に記載の核医学診断装置。 - 前記推定部は、セグメンテーション処理を行って散乱領域を特定し、特定した散乱領域に基づいて前記第2の確率値を推定する、請求項1に記載の核医学診断装置。
- ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得する第1の取得部と、
前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得する第2の取得部と、
前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントに対応する第2の確率値とを推定する推定部と、
前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する再構成部と、を備え
前記散乱イベントは多重散乱イベントを含む、
医用画像処理装置。 - 核医学診断装置により実行される核医学診断方法であって、
ガンマ線検出器への直接入射イベントと被写体での散乱イベントとを含む同時計数データを取得し、
前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、
前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントに対応する第2の確率値とを推定し、
前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成することを含み、
前記散乱イベントは多重散乱イベントを含む、核医学診断方法。 - コンピュータに、
ガンマ線検出器への直接入射イベントと被写体での多重散乱イベントを含む散乱イベントとを含む同時計数データを取得し、
前記被写体の電子密度関数及び前記ガンマ線検出器の幾何情報を取得し、
前記電子密度関数及び前記幾何情報の少なくとも一方に基づいて前記被写体における前記直接入射イベントに対応する第1の確率値と前記散乱イベントに対応する第2の確率値とを推定し、
前記第1の確率値と前記第2の確率値と、前記同時計数データとに基づいてPET画像を再構成する処理を実行させるプログラム。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021058499A JP2022155147A (ja) | 2021-03-30 | 2021-03-30 | 核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム |
CN202210319855.5A CN115137381A (zh) | 2021-03-30 | 2022-03-29 | 核医学诊断装置及方法、医用图像处理装置、存储介质 |
US17/656,957 US20220319068A1 (en) | 2021-03-30 | 2022-03-29 | Nuclear medicine diagnosis apparatus and nuclear medicine diagnosis method |
EP22165423.9A EP4092626A1 (en) | 2021-03-30 | 2022-03-30 | Nuclear medicine diagnosis apparatus, nuclear medicine diagnosis method, and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021058499A JP2022155147A (ja) | 2021-03-30 | 2021-03-30 | 核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022155147A true JP2022155147A (ja) | 2022-10-13 |
Family
ID=80999696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021058499A Pending JP2022155147A (ja) | 2021-03-30 | 2021-03-30 | 核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム |
Country Status (4)
Country | Link |
---|---|
US (1) | US20220319068A1 (ja) |
EP (1) | EP4092626A1 (ja) |
JP (1) | JP2022155147A (ja) |
CN (1) | CN115137381A (ja) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7310407B2 (en) * | 2004-09-03 | 2007-12-18 | Juni Jack E | Nuclear medical imaging device |
EP2748641B1 (en) * | 2011-10-06 | 2019-04-03 | Koninklijke Philips N.V. | Data-driven optimization of event acceptance/rejection logic |
CN109414238B (zh) * | 2016-06-12 | 2023-01-03 | 上海联影医疗科技股份有限公司 | 正电子发射断层成像***及其图像重建方法 |
US9916670B1 (en) * | 2016-09-13 | 2018-03-13 | Toshiba Medical Systems Corporation | Fast, efficient, and list-mode compatible tomographic image reconstruction using a novel quadratic surrogate |
US10973486B2 (en) * | 2018-01-08 | 2021-04-13 | Progenics Pharmaceuticals, Inc. | Systems and methods for rapid neural network-based image segmentation and radiopharmaceutical uptake determination |
US10768318B2 (en) * | 2018-10-17 | 2020-09-08 | Canon Medical Systems Corporation | Method and apparatus to determine energy correction arising from multi-channel detection in a gamma detector exhibiting non-linear energy measurement and/or cross-talk among channels |
-
2021
- 2021-03-30 JP JP2021058499A patent/JP2022155147A/ja active Pending
-
2022
- 2022-03-29 CN CN202210319855.5A patent/CN115137381A/zh active Pending
- 2022-03-29 US US17/656,957 patent/US20220319068A1/en active Pending
- 2022-03-30 EP EP22165423.9A patent/EP4092626A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
CN115137381A (zh) | 2022-10-04 |
US20220319068A1 (en) | 2022-10-06 |
EP4092626A1 (en) | 2022-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20230333269A1 (en) | System and method for image reconstruction in positron emission tomography | |
US9753146B2 (en) | Pixel identification for small pitch scintillation crystal arrays | |
US8558176B2 (en) | Medical image diagnosis apparatus and image reconstruction method | |
US10403010B2 (en) | Methods and systems for reconstructing images | |
WO2014180734A2 (en) | Apparatus and method for the evaluation of gamma radiation events | |
US20130251102A1 (en) | Image diagnosis apparatus and method | |
US20130136328A1 (en) | Methods and systems for enhanced tomographic imaging | |
JP2011252855A (ja) | 核医学イメージング装置 | |
US10219765B2 (en) | Nuclear medicine diagnostic apparatus and control method thereof | |
JP7317586B2 (ja) | 医用画像処理装置、方法及びプログラム | |
WO2015022606A1 (en) | A hybrid method based on simulation and experimental data to normalize pet data | |
CN108474862A (zh) | 具有lu谱减除的能量校准 | |
US10222490B2 (en) | PET scanner with emission and transmission structures in a checkerboard configuration | |
US20170119322A1 (en) | Medical image diagnosis system, structural image diagnosis apparatus, and nuclear medical image diagnosis apparatus | |
JP6502105B2 (ja) | 核医学画像撮影装置及び制御方法 | |
JPWO2016002084A1 (ja) | 画像再構成処理方法 | |
EP4080460A1 (en) | Medical image processing apparatus, medical image diagnosis apparatus, and program | |
US20220104781A1 (en) | Nuclear medicine diagnostic apparatus | |
JP2022155147A (ja) | 核医学診断装置、医用画像処理装置、核医学診断方法及びプログラム | |
JP2003294842A (ja) | 放射線投影データの補正方法 | |
JP2017512997A (ja) | 陽電子放射断層撮影(pet)における欠落画素の補償 | |
JP2017207482A (ja) | 医用画像診断装置 | |
EP4181069A1 (en) | Medical image processing device, medical image processing method, and program | |
CN112022189B (zh) | 图像重建方法、装置、设备及pet*** | |
JP2013007585A (ja) | 陽電子放出コンピュータ断層撮影装置及びX線CT(ComputedTomography)装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240208 |