Regularization design in penalized maximum-likelihood image reconstruction for lesion detection in 3D PET.

Published on Jan 20, 2014in Physics in Medicine and Biology2.883
· DOI :10.1088/0031-9155/59/2/403
Li Yang4
Estimated H-index: 4
(UC Davis: University of California, Davis),
Jian Zhou12
Estimated H-index: 12
(UC Davis: University of California, Davis)
+ 2 AuthorsJinyi Qi52
Estimated H-index: 52
(UC Davis: University of California, Davis)
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 channelized Hotelling observer (mvCHO) to assess the lesion detectability in 3D images to mimic the condition where a human observer examines three orthogonal views of a 3D image for lesion detection. We derived simplified theoretical expressions that allow fast prediction of the detectability of a 3D lesion. The theoretical results were used to design the regularization in PML reconstruction to improve lesion detectability. We conducted computer-based Monte Carlo simulations to compare the optimized penalty with the conventional penalty for detecting lesions of various sizes. Only true coincidence events were simulated. Lesion detectability was also assessed by two human observers, whose performances agree well with that of the mvCHO. Both the numerical observer and human observer results showed a statistically significant improvement in lesion detection by using the proposed penalty function compared to using the conventional penalty function.
📖 Papers frequently viewed together
5 Authors (Li Yang, ..., Jinyi Qi)
348 Citations
2015ISBI: International Symposium on Biomedical Imaging
3 Authors (Li Yang, ..., Jinyi Qi)
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
#1Ljiljana PlatisaH-Index: 9
#2Asli KumcuH-Index: 7
Last. Wilfried PhilipsH-Index: 47
view all 7 authors...
This study explores detection performance trends of human observers with respect to two parameters: task complexity determined by the frequency content of background and signal, and image viewing mode: singleslice (ss) versus multi-slice (ms) stack-browsing image presentation. The images are 3D correlated Gaussian noise with a 3D Gaussian signal centered in the image volume. To vary task complexity, we consider three different noise kernels while keeping the signal spread constant across all ima...
5 CitationsSource
#1Ljiljana Platisa (UGent: Ghent University)H-Index: 9
#2Bart Goossens (UGent: Ghent University)H-Index: 17
Last. Wilfried Philips (UGent: Ghent University)H-Index: 47
view all 7 authors...
Current clinical practice is rapidly moving in the direction of volumetric imaging. For two-dimensional (2D) images, task-based medical image quality is often assessed using numerical model observers. For three-dimensional (3D) images, however, these models have been little explored so far. In this work, first, two novel designs of a multislice channelized Hotelling observer (CHO) are proposed for the task of detecting 3D signals in 3D images. The novel designs are then compared and evaluated in...
80 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
#1Bastien Guerin (UPMC: Pierre-and-Marie-Curie University)H-Index: 18
#2G. El FakhriH-Index: 13
We have developed and validated a realistic simulation of random coincidences, pixelated block detectors, light sharing among crystal elements and dead-time in 2D and 3D positron emission tomography (PET) imaging based on the SimSET Monte Carlo simulation software. Our simulation was validated by comparison to a Monte Carlo transport code widely used for PET modeling, GATE, and to measurements made on a PET scanner. Methods: We have modified the SimSET software to allow independent tracking of s...
24 CitationsSource
#1Hongye Liang (CDRH: Center for Devices and Radiological Health)H-Index: 4
#2Subok Park (CDRH: Center for Devices and Radiological Health)H-Index: 15
Last. Aldo Badano (CDRH: Center for Devices and Radiological Health)H-Index: 26
view all 5 authors...
Rationale and Objectives Statistics show that radiologists are reading more studies than ever before, creating the challenge of interpreting an increasing number of images without compromising diagnostic performance. Stack-mode image display has the potential to allow radiologists to browse large three-dimensional (3D) datasets at refresh rates as high as 30 images/second. In this framework, the slow temporal response of liquid crystal displays (LCDs) can compromise the image quality when the im...
20 CitationsSource
#1Norbert Avril (QMUL: Queen Mary University of London)H-Index: 35
#2Lee P. AdlerH-Index: 1
Breast cancer is the most common female malignancy in Western countries. The limitations of mammography, ultrasound and MRI do not allow reliable identification of primary breast cancer at early stages. Functional breast imaging with positron emission tomography (PET) and F-18 fluorodeoxyglucose (FDG) enables the visualization of increased glucose metabolism of breast cancer. However, despite the successful identification of primary breast cancer, FDG-PET provides a low sensitivity to detect sma...
51 CitationsSource
#1Jinyi Qi (UC Davis: University of California, Davis)H-Index: 52
#2Ronald H. Huesman (LBNL: Lawrence Berkeley National Laboratory)H-Index: 24
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 simp...
29 CitationsSource
#1Howard C. Gifford (UMass: University of Massachusetts Amherst)H-Index: 22
#2A. Lehovich (UMass: University of Massachusetts Amherst)H-Index: 3
Last. Michael A. King (UMMS: University of Massachusetts Medical School)H-Index: 55
view all 3 authors...
A human-model observer for tumor detection-localization studies featuring multislice-multiview (or volumetric) image displays has been introduced. This volumetric observer, an extension of multiclass linear observers previously tested with single-slice and multislice displays, produces rating and localization data by integrating perception measurements from the different image views. A channelized NPW (CNPW) version of the observer was evaluated against humans for a background-known-exactly (BKE...
4 CitationsSource
Cited By23
#1Xinhui Wang (Oakland University)H-Index: 2
#2Bao Yang (UR: University of Rochester)H-Index: 4
Last. Jing Tang (Oakland University)H-Index: 14
view all 4 authors...
Myocardial perfusion (MP) PET imaging plays a key role in risk assessment and stratification of patients with coronary artery disease. In this work, we proposed a patch-based artificial neural network (ANN) fusion approach that integrates information from the maximum-likelihood (ML) and the post-smoothed ML reconstruction to improve MP PET imaging. To enhance quantification and tasked-based MP defect detection, the proposed method fused features from patches of the ML and the post-smoothed ML re...
2 CitationsSource
#1Yu-Jung Tsai (UCL: University College London)H-Index: 4
#2Georg Schramm (Katholieke Universiteit Leuven)H-Index: 17
Last. Kris Thielemans (UCL: University College London)H-Index: 36
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...
3 CitationsSource
Recently, patch-based regularisation methods have been a topic for tomographic reconstruction, since they are expected to preserve features within patches. However, conventional patch-based methods tend to ignore important similarities when the two patches being compared have a large intensity difference, while they share similar features in a scaled domain, which often results in the loss of fine details. In this work, to overcome this problem, the authors measure the similarity using normalise...
1 CitationsSource
#1Xinhui Wang (UR: University of Rochester)H-Index: 2
#2Bao Yang (UR: University of Rochester)H-Index: 4
Last. Jing Tang (UR: University of Rochester)H-Index: 14
view all 6 authors...
: Myocardial perfusion (MP) PET imaging plays an important role in risk assessment and stratification of patients with coronary artery disease. In this work, we developed an anatomy-assisted maximum a posteriori (MAP) reconstruction method incorporating a wavelet-based joint entropy (WJE) prior for MP PET imaging. Using the XCAT phantom, we first simulated three MP PET datasets, one with normal perfusion and the other two with non-transmural and transmural regionally reduced perfusion of the lef...
1 CitationsSource
: Filtered backprojection (FBP) algorithms reduce image noise by smoothing the image. Iterative algorithms reduce image noise by noise weighting and regularization. It is believed that iterative algorithms are able to reduce noise without sacrificing image resolution, and thus iterative algorithms, especially maximum-likelihood expectation maximization (MLEM), are used in nuclear medicine to replace FBP algorithms. Methods: This short paper uses counter examples to show that this belief is not t...
1 CitationsSource
#1Xuezhu ZhangH-Index: 13
#2Qiyu PengH-Index: 9
Last. Jinyi QiH-Index: 52
view all 6 authors...
: The first generation Tachyon PET (Tachyon-I) is a demonstration single-ring PET scanner that reaches a coincidence timing resolution of 314 ps using LSO scintillator crystals coupled to conventional photomultiplier tubes. The objective of this study was to quantify the improvement in both lesion detection and quantification performance resulting from the improved time-of-flight (TOF) capability of the Tachyon-I scanner. We developed a quantitative TOF image reconstruction method for the Tachyo...
4 CitationsSource
#1Brian Hutton (UCL: University College London)H-Index: 66
#2Kjell Erlandsson (UCL: University College London)H-Index: 27
Last. Kris Thielemans (UCL: University College London)H-Index: 4
view all 3 authors...
In this article, we describe recent developments in the design of both single-photon emission computed tomography (SPECT) and positron emission tomography (PET) instrumentation that have led to the current range of superior performance instruments. The adoption of solid-state technology for either complete detectors [e.g., cadmium zinc telluride (CZT)] or read-out systems that replace photomultiplier tubes [avalanche photodiodes (APD) or silicon photomultipliers (SiPM)] provide the advantage of ...
27 CitationsSource
#1Yihuan Lu (Yale University)H-Index: 7
#2Kathryn Fontaine (Yale University)H-Index: 4
Last. Chi Liu (Yale University)H-Index: 17
view all 7 authors...
With the recent development of high resolution time-of-flight (TOF) positron emission tomography (PET) and resolution recovery incorporated in the reconstruction algorithm, accurate detection, and quantification of sub-centimeter nodules might become feasible. In this paper, we performed a comprehensive simulation and mini-Derenzo phantom study to explore the quantitative accuracy of sub-centimeter nodules using the Siemens Biograph mCT scanner. We simulated nodules ranging from 4 to 10 mm in di...
6 CitationsSource
#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...
12 CitationsSource
#1Grace J. Gang (Johns Hopkins University)H-Index: 13
#2Jeffrey H. Siewerdsen (Johns Hopkins University)H-Index: 67
Last. J. Webster Stayman (Johns Hopkins University)H-Index: 15
view all 3 authors...
This paper presents a joint optimization of dynamic fluence field modulation (FFM) and regularization in quadratic penalized-likelihood reconstruction that maximizes a task-based imaging performance metric. We adopted a task-driven imaging framework for prospective designs of the imaging parameters. A maxi-min objective function was adopted to maximize the minimum detectability index ( {d}\,') throughout the image. The optimization algorithm alternates between FFM (represented by low-dimensio...
14 CitationsSource