Penalized maximum-likelihood image reconstruction for lesion detection

Published on Aug 21, 2006in Physics in Medicine and Biology3.609
路 DOI :10.1088/0031-9155/51/16/009
Jinyi Qi53
Estimated H-index: 53
(UC Davis: University of California, Davis),
Ronald H. Huesman33
Estimated H-index: 33
(LBNL: Lawrence Berkeley National Laboratory)
Detecting cancerous lesions is one major application in emission tomography. In this paper, we study penalized maximum-likelihood image reconstruction for this important clinical task. Compared to analytical reconstruction methods, statistical approaches can improve the image quality by accurately modelling the photon detection process and measurement noise in imaging systems. To explore the full potential of penalized maximum-likelihood image reconstruction for lesion detection, we derived simplified theoretical expressions that allow fast evaluation of the detectability of a random lesion. The theoretical results are used to design the regularization parameters to improve lesion detectability. We conducted computer-based Monte Carlo simulations to compare the proposed penalty function, conventional penalty function, and a penalty function for isotropic point spread function. The lesion detectability is measured by a channelized Hotelling observer. The results show that the proposed penalty function outperforms the other penalty functions for lesion detection. The relative improvement is dependent on the size of the lesion. However, we found that the penalty function optimized for a 5 mm lesion still outperforms the other two penalty functions for detecting a 14 mm lesion. Therefore, it is feasible to use the penalty function designed for small lesions in image reconstruction, because detection of large lesions is relatively easy.
馃摉 Papers frequently viewed together
#1Jinyi Qi (University of California, Berkeley)H-Index: 53
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 is the Gaussian prior, 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 ...
#1J. Oldan (SUNY: State University of New York System)H-Index: 1
#2S. Kulkarni (SUNY: State University of New York System)H-Index: 6
Last. Gene Gindi (SUNY: State University of New York System)H-Index: 26
view all 5 authors...
We compared the performance of a channelized Hotelling observer (CHO) to that of human observers to determine an optimal smoothing parameter /spl beta/ for an SKE/BKE detection task in a SPECT MAP (maximum a posteriori) reconstruction. The study is motivated in part by the recent development of theoretical methods that can rapidly predict CHO signal-to-noise ratios (SNRs) for MAP reconstructions. We found that a CHO not adjusted for internal noise effects was less predictive of the optimal smoot...
#1Jinyi Qi (LBNL: Lawrence Berkeley National Laboratory)H-Index: 53
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...
#1J.W. Stayman (UM: University of Michigan)H-Index: 3
#2Jeffrey A. Fessler (UM: University of Michigan)H-Index: 79
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...
#1Yuxiang Xing (SBU: Stony Brook University)H-Index: 18
#2Ing-Tsung Hsiao (CGU: Chang Gung University)H-Index: 26
Last. Gene Gindi (SBU: Stony Brook University)H-Index: 26
view all 3 authors...
We consider the calculation of lesion detectability using a mathematical model observer, the channelized Hotelling observer (CHO), in a signal-known-exactly/background-known-exactly detection task for single photon emission computed tomography (SPECT). We focus on SPECT images reconstructed with Bayesian maximum a posteriori methods. While model observers are designed to replace time-consuming studies using human observers, the calculation of CHO detectability is usually accomplished using a lar...
#1Jinyi Qi (LBNL: Lawrence Berkeley National Laboratory)H-Index: 53
Iterative image estimation methods have been widely used in emission tomography. Accurate estimate of the uncertainty of the reconstructed images is essential for quantitative applications. While theoretical approach has been developed to analyze the noise propagation from iteration to iteration, the current results are limited to only a few iterative algorithms that have an explicit multiplicative update equation. This paper presents a theoretical noise analysis that is applicable to a wide ran...
#1Jeffrey A. Fessler (UM: University of Michigan)H-Index: 79
What type of regularization method is optimal for penalized-likelihood image reconstruction when the imaging task is signal detection based on a channelized Hotelling (CHO) observer? To answer such questions, one would like to have a simple analytical expression (even if approximate) for the performance (SNR) of the CHO observer given different reconstruction methods. Bonetto, Qi, and Leahy (IEEE T-NS, Aug. 2000) derived and validated one such expression for penalized-likelihood (aka MAP) recons...
#1Jinyi Qi (LBNL: Lawrence Berkeley National Laboratory)H-Index: 53
#2Ronald H. Huesman (LBNL: Lawrence Berkeley National Laboratory)H-Index: 33
The low signal-to-noise ratio (SNR) in emission data has stimulated the development of statistical image reconstruction methods based on the maximum a posteriori (MAP) principle. Experimental examples have shown that statistical methods improve image quality compared to the conventional filtered backprojection (FBP) method. However, these results depend on isolated data sets. Here, the authors study the lesion detectability of MAP reconstruction theoretically, using computer observers. These the...
The maximum-likelihood (ML) approach in emission tomography provides images with superior noise characteristics compared to conventional filtered backprojection (FBP) algorithms. The expectation-maximization (EM) algorithm is an iterative algorithm for maximizing the Poisson likelihood in emission computed tomography that became very popular for solving the ML problem because of its attractive theoretical and practical properties. Recently, (Browne and DePierro, 1996 and Hudson and Larkin, 1991)...
#1J.S. Huber (LBNL: Lawrence Berkeley National Laboratory)H-Index: 13
#2Stephen E. Derenzo (LBNL: Lawrence Berkeley National Laboratory)H-Index: 33
Last. Thomas F. Budinger (LBNL: Lawrence Berkeley National Laboratory)H-Index: 62
view all 6 authors...
We present a conceptual design of a compact positron tomograph for prostate imaging using a pair of external curved detector banks, one placed above and one below the patient. The lower detector bank is fixed below the patient bed, and the top bank adjusts vertically for maximum sensitivity and patient access. Each bank is composed of 40 conventional block detectors, forming two arcs (44 cm minor, 60 cm major axis) that are tilted to minimize attenuation and positioned as close as possible to th...
Cited By30
#1Caroline CrockettH-Index: 3
#2Jeffrey A. Fessler (UM: University of Michigan)H-Index: 79
This review discusses methods for learning parameters for image reconstruction problems using bilevel formulations. Image reconstruction typically involves optimizing a cost function to recover a vector of unknown variables that agrees with collected measurements and prior assumptions. State-of-the-art image reconstruction methods learn these prior assumptions from training data using various machine learning techniques, such as bilevel methods. One can view the bilevel problem as formalizing hy...
#1Yu-Jung Tsai (UCL: University College London)H-Index: 7
#2Georg Schramm (Katholieke Universiteit Leuven)H-Index: 18
Last. Kris Thielemans (UCL: University College London)H-Index: 37
view all 9 authors...
In this study, we explore the use of a spatiallyvariant penalty strength in penalized image reconstruction using anatomical priors to reduce the dependence of lesion contrast on surrounding activity and lesion location. This work builds on a previous method to make the local perturbation response (LPR) approximately spatially invariant. While the dependence of lesion contrast on the local properties introduced by the anatomical penalty is intentional, the method aims to reduce the influence from...
#1Hossein Arabi (Geneva College)H-Index: 20
#2Habib Zaidi (University of Geneva)H-Index: 66
PET images commonly suffer from the high noise level and poor signal-to-noise ratio (SNR), thus adversely impacting lesion detectability and quantitative accuracy. In this work, a novel hybrid dual-domain PET denoising approach is proposed, which combines the advantages of both spatial and transform domain filtering to preserve image textures while minimizing quantification uncertainty. Spatial domain denoising techniques excel at preserving high-contrast patterns compared to transform domain fi...
#1Niccol貌 Camarlinghi (INFN: Istituto Nazionale di Fisica Nucleare)H-Index: 17
#2Giancarlo Sportelli (INFN: Istituto Nazionale di Fisica Nucleare)H-Index: 16
Last. Alberto Del GuerraH-Index: 19
view all 4 authors...
: Positron emission tomography (PET) iterative 3D reconstruction is a very computational demanding task. One of the main issues of the iterative reconstruction concerns the management of the system response matrix (SRM). The SRM models the relationship between the projection and the voxel space and its memory footprint can easily exceed hundreds of GB. Moreover, in order to make the reconstruction fast enough not to hinder its practical application, the SRM must be stored in the random access me...
#1Jorge Cabello (TUM: Technische Universit盲t M眉nchen)H-Index: 16
#2Sibylle Ziegler (TUM: Technische Universit盲t M眉nchen)H-Index: 50
The combination of positron emission tomography (PET) and MRI has attracted the attention of researchers in the past approximately 20 years in small-animal imaging and more recently in clinical research. The combination of PET/MRI allows researchers to explore clinical and research questions in a wide number of fields, some of which are briefly mentioned here. An important number of groups have developed different concepts to tackle the problems that PET instrumentation poses to the exposition o...
#1Hao Dang (Johns Hopkins University)H-Index: 13
#2J. Webster Stayman (Johns Hopkins University)H-Index: 15
Last. Jeffrey H. Siewerdsen (Johns Hopkins University)H-Index: 67
view all 11 authors...
: Task-based analysis of medical imaging performance underlies many ongoing efforts in the development of new imaging systems. In statistical image reconstruction, regularization is often formulated in terms to encourage smoothness and/or sharpness (e.g. a linear, quadratic, or Huber penalty) but without explicit formulation of the task. We propose an alternative regularization approach in which a spatially varying penalty is determined that maximizes task-based imaging performance at every loca...
Tube current modulation (TCM) is routinely adopted on diagnostic CT scanners for dose reduction. Conventional TCM strategies are generally designed for filtered-backprojection (FBP) reconstruction to satisfy simple image quality requirements based on noise. This work investigates TCM designs for model-based iterative reconstruction (MBIR) to achieve optimal imaging performance as determined by a task-based image quality metric. Additionally, regularization is an important aspect of MBIR that is ...
#1Bert-Ram Sah (UZH: University of Zurich)H-Index: 16
#2Paul Stolzmann (UZH: University of Zurich)H-Index: 57
Last. Patrick Veit-HaibachH-Index: 39
view all 11 authors...
PURPOSE To investigate the clinical performance of a block sequential regularized expectation maximization (BSREM) penalized likelihood reconstruction algorithm in oncologic PET/computed tomography (CT) studies. METHODS A total of 410 reconstructions of 41 fluorine-18 fluorodeoxyglucose-PET/CT studies of 41 patients with a total of 2010 lesions were analyzed by two experienced nuclear medicine physicians. Images were reconstructed with BSREM (with four different 尾 values) or ordered subset expec...
#1Kristen A. Wangerin (UW: University of Washington)H-Index: 9
#2Sangtae Ahn (GE: General Electric)H-Index: 23
Last. Ravindra Mohan Manjeshwar (GE: General Electric)H-Index: 16
view all 6 authors...
: We have previously developed a convergent penalized likelihood (PL) image reconstruction algorithm using the relative difference prior (RDP) and showed that it achieves more accurate lesion quantitation compared to ordered subsets expectation maximization (OSEM). We evaluated the detectability of low-contrast liver and lung lesions using the PL-RDP algorithm compared to OSEM. We performed a two-alternative forced choice study using a channelized Hotelling observer model that was previously val...
#1Li Yang (UC Davis: University of California, Davis)H-Index: 4
#2Guobao Wang (UC Davis: University of California, Davis)H-Index: 19
Last. Jinyi Qi (UC Davis: University of California, Davis)H-Index: 53
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...
This website uses cookies.
We use cookies to improve your online experience. By continuing to use our website we assume you agree to the placement of these cookies.
To learn more, you can find in our Privacy Policy.