Evaluation of penalty design in penalized maximum-likelihood image reconstruction for lesion detection

Published on Oct 1, 2014in Journal of medical imaging
· DOI :10.1117/1.JMI.1.3.035501
Li Yang4
Estimated H-index: 4
(UC Davis: University of California, Davis),
Andrea Ferrero8
Estimated H-index: 8
(UC Davis: University of California, Davis)
+ 2 AuthorsJinyi Qi52
Estimated H-index: 52
(UC Davis: University of California, Davis)
Sources
Abstract
Detecting cancerous lesions is a major clinical application in emission tomography. Previously, we developed a method to design a shift-variant quadratic penalty function in penalized maximum-likelihood (PML) image reconstruction to improve the lesion detectability. We used a multiview channelized Hotelling observer (mvCHO) to assess the lesion detectability in three-dimensional images and validated the penalty design using computer simulations. In this study, we evaluate the benefit of the proposed penalty function for lesion detection using real patient data and artificial lesions. A high-count real patient dataset with no identifiable tumor inside the field of view is used as the background data. A Na-22 point source is scanned in air at variable locations and the point source data are superimposed onto the patient data as artificial lesions after being attenuated by the patient body. Independent Poisson noise is introduced to the high-count sinograms to generate 200 pairs of lesion-present and lesion-absent datasets, each mimicking a 5-min scan. Lesion detectability is assessed using a mvCHO and a human observer two-alternative forced choice (2AFC) experiment. The results show improvements in lesion detection by the proposed method compared with the conventional first-order quadratic penalty function and a total variation (TV) edge-preserving penalty function.
📖 Papers frequently viewed together
30 Citations
2012ISBI: International Symposium on Biomedical Imaging
3 Authors (Li Yang, ..., Jinyi Qi)
5 Citations
2004
1 Author (Jinyi Qi)
3 Citations
References27
Newest
#1Li Yang (UC Davis: University of California, Davis)H-Index: 4
#2Jian Zhou (UC Davis: University of California, Davis)H-Index: 12
Last. Jinyi Qi (UC Davis: University of California, Davis)H-Index: 52
view all 5 authors...
Detecting cancerous lesions is a major clinical application in emission tomography. In previous work, we have studied penalized maximum-likelihood (PML) image reconstruction for the detection task and proposed a method to design a shift-invariant quadratic penalty function to maximize detectability of a lesion at a known location in a two dimensional image. Here we extend the regularization design to maximize detectability of lesions at unknown locations in fully 3D PET. We used a multiview chan...
24 CitationsSource
May 2, 2012 in ISBI (International Symposium on Biomedical Imaging)
#1Li Yang (UC Davis: University of California, Davis)H-Index: 4
#2Jian Zhou (UC Davis: University of California, Davis)H-Index: 140
Last. Jinyi Qi (UC Davis: University of California, Davis)H-Index: 52
view all 3 authors...
Detecting cancerous lesion is a major clinical application in emission tomography. In a previous work, we have shown that penalized maximum likelihood image reconstruction can improve lesion detection at a fixed location by designing a shift-invariant quadratic penalty function. Here we extend this work to detection of tumors at unknown positions. We present a method to design a shift-variant quadratic penalty function that maximizes the detectability of lesions at all possible locations. We con...
5 CitationsSource
#1H.R. Shi (UM: University of Michigan)H-Index: 1
#2Jeffrey A. Fessler (UM: University of Michigan)H-Index: 78
Statistical methods for tomographic image reconstruction have improved noise and spatial resolution properties that may improve image quality in X-ray computed tomography (CT). Penalized weighted least squares (PWLS) methods using conventional quadratic regularization lead to nonuniform and anisotropic spatial resolution due to interactions between the weighting, which is necessary for good noise properties, and the regularizer. Previously, we addressed this problem for parallel-beam emission to...
30 CitationsSource
#1Johan Nuyts (Katholieke Universiteit Leuven)H-Index: 61
#2Christian Michel (Siemens)H-Index: 20
Last. Raf Verscuren (Katholieke Universiteit Leuven)H-Index: 3
view all 10 authors...
For positron emission tomography (PET) imaging, different reconstruction methods can be applied, including maximum likelihood (ML ) and maximum a posteriori (MAP) reconstruction. Postsmoothed ML images have approximately position and object independent spatial resolution, which is advantageous for (semi-) quantitative analysis. However, the complex object dependent smoothing obtained with MAP might yield improved noise characteristics, beneficial for lesion detection. In this contribution, MAP a...
23 CitationsSource
Statistical image reconstruction methods based on maximum a posteriori (MAP) principle have been developed for emission tomography. The prior distribution of the unknown image plays an important role in MAP reconstruction. The most commonly used prior are Gaussian priors, whose logarithm has a quadratic form. Gaussian priors are relatively easy to analyze. It has been shown that the effect of a Gaussian prior can be approximated by linear filtering a maximum likelihood (ML) reconstruction. As a ...
12 CitationsSource
Projection and backprojection are operations that arise frequently in tomographic imaging. Recently, we proposed a new method for projection and backprojection, which we call distance-driven, and that offers low arithmetic cost and a highly sequential memory access pattern. Furthermore, distance-driven projection and backprojection avoid several artefact-inducing approximations characteristic of some other methods. We have previously demonstrated the application of this method to parallel and fa...
392 CitationsSource
#1J.W. Stayman (UM: University of Michigan)H-Index: 5
#2Jeffrey A. Fessler (UM: University of Michigan)H-Index: 78
Imaging systems that form estimates using a statistical approach generally yield images with nonuniform resolution properties. That is, the reconstructed images possess resolution properties marked by space-variant and/or anisotropic responses. We have previously developed a space-variant penalty for penalized-likelihood (PL) reconstruction that yields nearly uniform resolution properties . We demonstrated how to calculate this penalty efficiently and apply it to an idealized positron emission t...
102 CitationsSource
#1Jinyi Qi (LBNL: Lawrence Berkeley National Laboratory)H-Index: 52
Bayesian methods based on the maximum a posteriori principle (also called penalized maximum-likelihood methods) have been developed to improve image quality in emission tomography. To explore the full potential of Bayesian reconstruction for lesion detection, we derive simplified theoretical expressions that allow fast evaluation of the detectability of a lesion in Bayesian reconstruction. This work is built on the recent progress on the theoretical analysis of image properties of statistical re...
47 CitationsSource
#1Harrison H. BarrettH-Index: 64
#2Kyle J. MyersH-Index: 30
Last. Satyapal Rathee (U of A: University of Alberta)H-Index: 23
view all 3 authors...
769 Citations
#1Craig K. Abbey (UA: University of Arizona)H-Index: 35
#2Harrison H. Barrett (UA: University of Arizona)H-Index: 64
We consider detection of a nodule signal profile in noisy images meant to roughly simulate the statistical properties of tomographic image reconstructions in nuclear medicine. The images have two sources of variability arising from quantum noise from the imaging process and anatomical variability in the ensemble of objects being imaged. Both of these sources of variability are simulated by a stationary Gaussian random process. Sample images from this process are generated by filtering white-nois...
268 CitationsSource
Cited By4
Newest
#1Juan Gao (CAS: Chinese Academy of Sciences)H-Index: 20
#2Qiyang Zhang (CAS: Chinese Academy of Sciences)H-Index: 7
Last. Hu Zhanli (CAS: Chinese Academy of Sciences)H-Index: 17
view all 10 authors...
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 algorith...
Source
#1Elisabetta De Bernardi (University of Milano-Bicocca)H-Index: 7
#2Federico FallancaH-Index: 17
Last. Valentino BettinardiH-Index: 34
view all 5 authors...
Purpose The effects of regularizing priors on the maximum likelihood (ML) reconstruction of activity patterns in Positron Emission Tomography (PET) were assessed. Methods Two edge-preserving priors (one originally proposed by Nuyts et al. and nowadays implemented and commercialized by General Electric Medical Systems as Q.Clear software, and a second one originally proposed by Rapisarda et al. and our group) were assessed and compared to a standard Ordered Subset (OS)-ML reconstruction, assumed ...
2 CitationsSource
#1Li Yang (Toshiba)H-Index: 2
#2Jian Zhou (Toshiba)H-Index: 2
Last. Wenli Wang (Toshiba)H-Index: 7
view all 4 authors...
Region of interest (ROI) quantification is one of the important tasks in positron emission tomography (PET) imaging. Statistical reconstruction methods based on the penalized maximum-likelihood (PML) principle have been developed to improve image quality over early stopping with EM-based methods and analytic methods. The penalty function that controls spatial smoothness plays an important role in PML reconstruction. The commonly used quadratic penalty often oversmooths sharp edges, while the edg...
1 CitationsSource
#1Li Yang (UC Davis: University of California, Davis)H-Index: 4
#2Guobao Wang (UC Davis: University of California, Davis)H-Index: 17
Last. Jinyi Qi (UC Davis: University of California, Davis)H-Index: 52
view all 3 authors...
Detecting cancerous lesions is a major clinical application of emission tomography. In a previous work, we studied penalized maximum-likelihood (PML) image reconstruction for lesion detection in static PET. Here we extend our theoretical analysis of static PET reconstruction to dynamic PET. We study both the conventional indirect reconstruction and direct reconstruction for Patlak parametric image estimation. In indirect reconstruction, Patlak parametric images are generated by first reconstruct...
4 CitationsSource