Positron emission tomography image reconstruction using feature extraction.

PURPOSE To reduce the cost of positron emission tomography (PET) scanning systems, image reconstruction algorithms for low-sampled data have been extensively studied. However, the current method based on total variation (TV) minimization regularization nested in the maximum likelihood-expectation maximization (MLEM) algorithm cannot distinguish true structures from noise resulting losing some fine features in the images. Thus, this work aims to recover fine features lost in the MLEM-TV algorithm from low-sampled data. METHOD A feature refinement (FR) approach previously developed for statistical interior computed tomography (CT) reconstruction is applied to PET imaging to recover fine features in this study. The proposed method starts with a constant initial image and the FR step is performed after each MLEM-TV iteration to extract the desired structural information lost during TV minimization. A feature descriptor is specifically designed to distinguish structure from noise and artifacts. A modified steepest descent method is adopted to minimize the objective function. After evaluating the impacts of different patch sizes on the outcome of the presented method, an optimal patch size of 7×7 is selected in this study to balance structure-detection ability and computational efficiency. RESULTS Applying MLEM-TV-FR algorithm to the simulated brain PET imaging using an emission activity phantom, a standard Shepp-Logan phantom, and mouse results in the increased peak signal-to-noise ratio (PSNR) and structural similarity (SSIM) as comparing to using the conventional MLEM-TV algorithm, as well as the substantial reduction of the used sampling numbers, which improves the computational efficiency. CONCLUSIONS The presented algorithm can achieve image quality superior to that of the MLEM and MLEM-TV approaches in terms of the preservation of fine structure and the suppression of undesired artifacts and noise, indicating its useful potential for low-sampled data in PET imaging.

[1]  Jinyi Qi,et al.  Iterative reconstruction techniques in emission computed tomography , 2006, Physics in medicine and biology.

[2]  Dong Liang,et al.  A feature refinement approach for statistical interior CT reconstruction , 2016, Physics in medicine and biology.

[3]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[4]  Yuanjun Wang,et al.  A new adaptive-weighted total variation sparse-view computed tomography image reconstruction with local improved gradient information. , 2018, Journal of X-ray science and technology.

[5]  Craig S. Levin,et al.  Compressed sensing for the multiplexing of large area silicon photomultiplier PET detectors: Acquisition and calibration , 2012 .

[6]  Otmar Scherzer,et al.  Inverse Total Variation Flow , 2007, Multiscale Model. Simul..

[7]  Joseph A. O'Sullivan,et al.  Alternating Minimization Algorithms for Transmission Tomography , 2007, IEEE Transactions on Medical Imaging.

[8]  Jungah Son,et al.  Gap compensation during PET image reconstruction by constrained, total variation minimization. , 2012, Medical physics.

[9]  Martin Burger,et al.  Reconstruction of short time PET scans using Bregman iterations , 2011, 2011 IEEE Nuclear Science Symposium Conference Record.

[10]  Martin Burger,et al.  Primal and Dual Bregman Methods with Application to Optical Nanoscopy , 2011, International Journal of Computer Vision.

[11]  O Mawlawi,et al.  MO-D-141-05: Dictionary Learning in Compressed Sensing Using Undersampled Data in PET Imaging. , 2013, Medical physics.

[12]  Guobao Wang,et al.  PET Image Reconstruction Using Kernel Method , 2015, IEEE Transactions on Medical Imaging.

[13]  Guobao Wang,et al.  Penalized Likelihood PET Image Reconstruction Using Patch-Based Edge-Preserving Regularization , 2012, IEEE Transactions on Medical Imaging.

[14]  Peter M Jakob,et al.  Compressed sensing for reduction of noise and artefacts in direct PET image reconstruction. , 2014, Zeitschrift fur medizinische Physik.

[15]  Richard M. Leahy,et al.  Analysis of Resolution and Noise Properties of Nonquadratically Regularized Image Reconstruction Methods for PET , 2008, IEEE Transactions on Medical Imaging.

[16]  Jing Tang,et al.  Sparsity-constrained PET image reconstruction with learned dictionaries , 2016, Physics in medicine and biology.

[17]  John W. Clark,et al.  Dictionary learning for data recovery in positron emission tomography , 2015, Physics in medicine and biology.

[18]  Jeffrey A. Fessler.  Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[19]  Dong Liang,et al.  PET Image Reconstruction from Under-sampled Data , 2017, 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC).

[20]  John W. Clark,et al.  Using compressive sensing to recover images from PET scanners with partial detector rings. , 2014, Medical physics.

[21]  M. Burger,et al.  Accurate EM-TV algorithm in PET with low SNR , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[22]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[23]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[24]  Li Yang,et al.  Evaluation of penalty design in penalized maximum-likelihood image reconstruction for lesion detection , 2014, Journal of medical imaging.

[25]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.