A mathematical observer study for the evaluation and optimization of compensation methods for myocardial SPECT using a phantom population that realistically models patient variability

The goal of this study was to develop and apply a population of phantoms that realistically models patient variability and use it to optimize and evaluate different compensation methods used during reconstruction process with respect to defect detection in myocardial SPECT images. Various combinations of attenuation, detector response and scatter compensation were used in this study. A major difference between this and previous studies was that the level of realism was significantly increased by inclusion of variability in heart and organ uptakes, in the heart size and orientation, and in the defect size and contrast. In this study we used a population of 24 4-D NCAT phantoms (half male, half female) recently developed with statistical models for organ uptake and organ size based on clinical data. Almost noise-free projection data of the torso, heart, liver, lungs, and other organs were simulated for each phantom using the SIMSET MC simulation code. They were then combined to form 72 sets of projections for each phantom using randomly sampled activity ratios from a clinically realistic distribution. Poisson noise was then added to the projection data. We applied the channelized hotelling observer (CHO) and receiver operating characteristic (ROC) analysis to optimize iteration number for OSEM and cutoff frequency of a 3-D post-reconstruction Butterworth filter. We found that the area under the ROC curve (AUC) values were reduced compared to a previous study that included significantly less phantom variability, even though the defect contrast was higher and noise level was lower. The resulting AUC values were similar to those obtained using patient data. We found, in agreement with the previous study, that including compensation for more effects resulted in improved defect detectability. However, the optimal filter cutoff frequency was increased compared to the previous study. These studies demonstrate the importance of including realistic levels of phantom variability in myocardial perfusion studies using simulated data.

[1]  Benjamin M. W. Tsui,et al.  Modeling respiratory mechanics in the MCAT and spline-based MCAT phantoms , 1999 .

[2]  D. Lalush,et al.  A realistic spline-based dynamic heart phantom , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[3]  Eric C. Frey,et al.  Application of task-based measures of image quality to optimization and evaluation of three-dimensional reconstruction-based compensation methods in myocardial perfusion SPECT , 2002, IEEE Transactions on Medical Imaging.

[4]  E C Frey,et al.  Receiver operating characteristic evaluation of iterative reconstruction with attenuation correction in 99mTc-sestamibi myocardial SPECT images. , 2000, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[5]  Eric C Frey,et al.  Optimum compensation method and filter cutoff frequency in myocardial SPECT: a human observer study. , 2002, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[6]  C. Metz,et al.  Maximum likelihood estimation of receiver operating characteristic (ROC) curves from continuously-distributed data. , 1998, Statistics in medicine.

[7]  B M Tsui,et al.  Evaluation of cardiac cone-beam single photon emission computed tomography using observer performance experiments and receiver operating characteristic analysis. , 1993, Investigative radiology.

[8]  H.H. Barrett,et al.  Model observers for assessment of image quality , 1993, 2002 IEEE Nuclear Science Symposium Conference Record.

[9]  J. Hanley The Robustness of the "Binormal" Assumptions Used in Fitting ROC Curves , 1988, Medical decision making : an international journal of the Society for Medical Decision Making.

[10]  B. Tsui,et al.  A new method for modeling the spatially-variant, object-dependent scatter response function in SPECT , 1996, 1996 IEEE Nuclear Science Symposium. Conference Record.

[11]  E C Frey,et al.  Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction. , 1998, Physics in medicine and biology.

[12]  David Gur,et al.  Statistical issues in ROC curve analysis , 1990, Medical Imaging.

[13]  Matt A. King,et al.  Attenuation compensation for cardiac single-photon emission computed tomographic imaging: Part 2. Attenuation compensation algorithms , 1996, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[14]  C. Metz,et al.  Statistical significance tests for binormal ROC curves , 1980 .

[15]  Harrison H. Barrett,et al.  Hotelling trace criterion as a figure of merit for the optimization of imaging systems , 1986 .

[16]  William Paul Segars,et al.  Development of a new dynamic NURBS-based cardiac-torso (NCAT) phantom , 2001 .