Resolution modeling in projection space using a factorized multi-block detector response function for PET image reconstruction.

Positron emission tomography (PET) images usually suffer from limited resolution and statistical uncertainties. However, a technique known as resolution modeling (RM) can be used to improve image quality by accurately modeling the system's detection process within the iterative reconstruction. In this study, we present an accurate RM method in projection space based on a simulated multi-block detector response function (DRF) and evaluate it on the Siemens hybrid MR-BrainPET system. The DRF is obtained using GATE simulations that consider nearly all the possible annihilation photons from the field-of-view (FOV). Intrinsically, the multi-block DRF allows the block crosstalk to be modeled. The RM blurring kernel is further generated by factorizing the blurring matrix of one line-of-response (LOR) into two independent detector responses, which can then be addressed with the DRF. Such a kernel is shift-variant in 4D projection space without any distance or angle compression, and is integrated into the image reconstruction for the BrainPET insert with single instruction multiple data (SIMD) and multi-thread support. Evaluation of simulations and measured data demonstrate that the reconstruction with RM yields significantly improved resolutions and reduced mean squared error (MSE) values at different locations of the FOV, compared with reconstruction without RM. Furthermore, the shift-variant RM kernel models the varying blurring intensity for different LORs due to the depth-of-interaction (DOI) dependencies, thus avoiding severe edge artifacts in the images. Additionally, compared to RM in single-block mode, the multi-block mode shows significantly improved resolution and edge recovery at locations beyond 10 cm from the center of BrainPET insert in the transverse plane. However, the differences have been observed to be low for patient data between single-block and multi-block mode RM, due to the brain size and location as well as the geometry of the BrainPET insert. In conclusion, the RM method proposed in this study can yield better reconstructed images in terms of resolution and MSE value, compared to conventional reconstruction without RM.

[1]  Andrew J. Reader,et al.  Impact of Image-Space Resolution Modeling for Studies with the High-Resolution Research Tomograph , 2008, Journal of Nuclear Medicine.

[2]  C. C. Watson.  New, faster, image-based scatter correction for 3D PET , 1999 .

[3]  Steven G. Ross,et al.  Application and Evaluation of a Measured Spatially Variant System Model for PET Image Reconstruction , 2010, IEEE Transactions on Medical Imaging.

[4]  W P Segars,et al.  Fast modelling of the collimator–detector response in Monte Carlo simulation of SPECT imaging using the angular response function , 2005, Physics in medicine and biology.

[5]  R. Laforest,et al.  Positron range modeling for statistical PET image reconstruction , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[6]  V. Bettinardi,et al.  Evaluation of a New Regularization Prior for 3-D PET Reconstruction Including PSF Modeling , 2012, IEEE Transactions on Nuclear Science.

[7]  U Pietrzyk,et al.  High performance volume-of-intersection projectors for 3D-PET image reconstruction based on polar symmetries and SIMD vectorisation. , 2015, Physics in medicine and biology.

[8]  Keishi Kitamura,et al.  Transaxial system models for jPET-D4 image reconstruction , 2005, Physics in medicine and biology.

[9]  A. Martineau,et al.  A method for accurate modelling of the crystal response function at a crystal sub-level applied to PET reconstruction , 2011, Physics in medicine and biology.

[10]  R. Fontaine,et al.  Fast, accurate and versatile Monte Carlo method for computing system matrix , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[11]  P K Marsden,et al.  Developments in component-based normalization for 3D PET. , 1999, Physics in medicine and biology.

[12]  F Hofheinz,et al.  A volume of intersection approach for on-the-fly system matrix calculation in 3D PET image reconstruction. , 2014, Physics in medicine and biology.

[13]  G. Delso,et al.  Performance Measurements of the Siemens mMR Integrated Whole-Body PET/MR Scanner , 2011, The Journal of Nuclear Medicine.

[14]  A. Dell'Acqua,et al.  Geant4 - A simulation toolkit , 2003 .

[15]  Jürgen Scheins,et al.  The Jülich Experience With Simultaneous 3T MR-BrainPET: Methods and Technology , 2019, IEEE Transactions on Radiation and Plasma Medical Sciences.

[16]  S Stute,et al.  GATE V6: a major enhancement of the GATE simulation platform enabling modelling of CT and radiotherapy , 2011, Physics in medicine and biology.

[17]  R. Leahy,et al.  High-resolution 3D Bayesian image reconstruction using the microPET small-animal scanner. , 1998, Physics in medicine and biology.

[18]  Jeih-San Liow,et al.  Variance reduction on randoms from coincidence histograms for the HRRT , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[19]  Thomas K. Lewellen,et al.  Modeling and incorporation of system response functions in 3-D whole body PET , 2006, IEEE Transactions on Medical Imaging.

[20]  Dan J Kadrmas,et al.  LOR-OSEM: statistical PET reconstruction from raw line-of-response histograms , 2004, Physics in medicine and biology.

[21]  E. Hoffman,et al.  Quantitation in Positron Emission Computed Tomography: 4. Effect of Accidental Coincidences , 1981, Journal of computer assisted tomography.

[22]  Michael E Casey,et al.  Generalized PSF modeling for optimized quantitation in PET imaging , 2017, Physics in medicine and biology.

[23]  Nadim Joni Shah,et al.  Analysis and Correction of Count Rate Reduction During Simultaneous MR-PET Measurements With the BrainPET Scanner , 2012, IEEE Transactions on Medical Imaging.

[24]  C Lartizien,et al.  GATE: a simulation toolkit for PET and SPECT. , 2004, Physics in medicine and biology.

[25]  R Marabini,et al.  System models for PET statistical iterative reconstruction: A review , 2016, Comput. Medical Imaging Graph..

[26]  Pablo Aguiar,et al.  Geometrical and Monte Carlo projectors in 3D PET reconstruction. , 2010, Medical physics.

[27]  Alfred O. Hero,et al.  Ieee Transactions on Image Processing: to Appear Penalized Maximum-likelihood Image Reconstruction Using Space-alternating Generalized Em Algorithms , 2022 .

[28]  Gaspar Delso,et al.  Design Features and Mutual Compatibility Studies of the Time-of-Flight PET Capable GE SIGNA PET/MR System , 2016, IEEE Transactions on Medical Imaging.

[29]  Long Zhang,et al.  Fast and memory-efficient Monte Carlo-based image reconstruction for whole-body PET. , 2010, Medical physics.

[30]  A. Chatziioannou,et al.  Fully 3D system model estimation of OPET by Monte Carlo simulation , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[31]  Vladimir Y. Panin,et al.  Fully 3-D PET reconstruction with system matrix derived from point source measurements , 2006, IEEE Transactions on Medical Imaging.

[32]  M. Rafecas,et al.  Use of a Monte Carlo-based probability matrix for 3-D iterative reconstruction of MADPET-II data , 2004, IEEE Transactions on Nuclear Science.

[33]  Jian Zhou,et al.  Quantitative image reconstruction for total-body PET imaging using the 2-meter long EXPLORER scanner , 2017, Physics in medicine and biology.

[34]  Jing Tang,et al.  Analytic system matrix resolution modeling in PET: an application to Rb-82 cardiac imaging , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[35]  F Hofheinz,et al.  Effects of cold sphere walls in PET phantom measurements on the volume reproducing threshold , 2010, Physics in medicine and biology.

[36]  Andrew J. Reader,et al.  EM algorithm system modeling by image-space techniques for PET reconstruction , 2003 .

[37]  Arman Rahmim,et al.  Resolution modeling in PET imaging: Theory, practice, benefits, and pitfalls. , 2013, Medical physics.

[38]  K-J Langen,et al.  High resolution BrainPET combined with simultaneous MRI , 2011, Nuklearmedizin.

[39]  H. Herzog,et al.  Reconstruction of PET data acquired with the BrainPET using STIR , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

[40]  A. Rahmim,et al.  Monte Carlo-based evaluation of inter-crystal scatter and penetration in the PET subsystem of three GE Discovery PET/CT scanners , 2011 .

[41]  Jeih-San Liow,et al.  Variance Reduction on Randoms from delayed coincidence histograms for the HRRT , 2005 .

[42]  Elsa D. Angelini,et al.  Locally weighted total variation denoising for ringing artifact suppression in pet reconstruction using PSF modeling , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[43]  Nadim Joni Shah,et al.  Fully-3D PET Image Reconstruction Using Scanner-Independent, Adaptive Projection Data and Highly Rotation-Symmetric Voxel Assemblies , 2011, IEEE Transactions on Medical Imaging.

[44]  Jian Zhou,et al.  Sinogram Blurring Matrix Estimation From Point Sources Measurements With Rank-One Approximation for Fully 3-D PET , 2017, IEEE Transactions on Medical Imaging.

[45]  Hans Herzog,et al.  Comparison of Template-Based Versus CT-Based Attenuation Correction for Hybrid MR/PET Scanners , 2015, IEEE Transactions on Nuclear Science.

[46]  R. Leahy,et al.  Accurate geometric and physical response modelling for statistical image reconstruction in high resolution PET , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[47]  D. Louis Collins,et al.  Twenty New Digital Brain Phantoms for Creation of Validation Image Data Bases , 2006, IEEE Transactions on Medical Imaging.

[48]  Nicola Belcari,et al.  Accurate and efficient modeling of the detector response in small animal multi-head PET systems , 2013, Physics in medicine and biology.

[49]  Martin A Lodge,et al.  Simultaneous measurement of noise and spatial resolution in PET phantom images. , 2010, Physics in medicine and biology.

[50]  S. Ross,et al.  Properties and Mitigation of Edge Artifacts in PSF-Based PET Reconstruction , 2011, IEEE Transactions on Nuclear Science.

[51]  A. Rahmim,et al.  Space-variant and anisotropic resolution modeling in list-mode EM reconstruction , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[52]  Paul Kinahan,et al.  Noise and signal properties in PSF-based fully 3D PET image reconstruction: an experimental evaluation , 2010, Physics in medicine and biology.