DocumentCode :
1207485
Title :
Toward Realistic and Practical Ideal Observer (IO) Estimation for the Optimization of Medical Imaging Systems
Author :
He, Xin ; Caffo, Brian S. ; Frey, Eric C.
Author_Institution :
Dept. of Radiol., Johns Hopkins Sch. of Med., Baltimore, MD
Volume :
27
Issue :
10
fYear :
2008
Firstpage :
1535
Lastpage :
1543
Abstract :
The ideal observer (IO) employs complete knowledge of the available data statistics and sets an upper limit on observer performance on a binary classification task. However, the IO test statistic cannot be calculated analytically, except for cases where object statistics are extremely simple. Kupinski et al. have developed a Markov chain Monte Carlo (MCMC) based technique to compute the IO test statistic for, in principle, arbitrarily complex objects and imaging systems. In this work, we applied MCMC to estimate the IO test statistic in the context of myocardial perfusion SPECT (MPS). We modeled the imaging system using an analytic SPECT projector with attenuation, distant-dependent detector-response modeling and Poisson noise statistics. The object is a family of parameterized torso phantoms with variable geometric and organ uptake parameters. To accelerate the imaging simulation process and thus enable the MCMC IO estimation, we used discretized anatomic parameters and continuous uptake parameters in defining the objects. The imaging process simulation was modeled by precomputing projections for each organ for a finite number of discretely-parameterized anatomic parameters and taking linear combinations of the organ projections based on continuous sampling of the organ uptake parameters. The proposed method greatly reduces the computational burden and allows MCMC IO estimation for a realistic MPS imaging simulation. We validated the proposed IO estimation technique by estimating IO test statistics for a large number of input objects. The properties of the first- and second-order statistics of the IO test statistics estimated using the MCMC IO estimation technique agreed well with theoretical predictions. Further, as expected, the IO had better performance, as measured by the receiver operating characteristic (ROC) curve, than the Hotelling observer. This method is developed for SPECT imaging. However, it can be adapted to any linear imaging system.
Keywords :
Markov processes; Monte Carlo methods; Poisson distribution; biomedical imaging; cardiology; phantoms; single photon emission computed tomography; IO test statistics; MCMC IO estimation; Markov chain Monte Carlo technique; Poisson noise statistics; binary classification; detector-response modeling; ideal observer; medical imaging; myocardial perfusion SPECT; parameterized torso phantoms; receiver operating characteristic curve; Attenuation; Biomedical imaging; Computational modeling; Image analysis; Monte Carlo methods; Myocardium; Statistical analysis; Statistics; System testing; Torso; Ideal observer; Markov chain Monte Carlo (MCMC); Algorithms; Artificial Intelligence; Computer Simulation; Data Interpretation, Statistical; Image Enhancement; Image Interpretation, Computer-Assisted; Models, Biological; Models, Statistical; Numerical Analysis, Computer-Assisted; Pattern Recognition, Automated; Reproducibility of Results; Sensitivity and Specificity; Signal Processing, Computer-Assisted;
fLanguage :
English
Journal_Title :
Medical Imaging, IEEE Transactions on
Publisher :
ieee
ISSN :
0278-0062
Type :
jour
DOI :
10.1109/TMI.2008.924641
Filename :
4505799
Link To Document :
بازگشت