Biomedical subjects

Eric C Frey

Publications and source records attributed to Eric C Frey.

13 recordsLinked to original sources

PubMed2006

Comparison of conventional, model-based quantitative planar, and quantitative SPECT image processing methods for organ activity estimation using In-111 agents.

Accurate quantification of organ radionuclide uptake is important for patient-specific dosimetry. The quantitative accuracy from conventional conjugate view methods is limited by overlap of projections from different organs and background activity, and attenuation and scatter. In this work, we propose and validate a quantitative planar (QPlanar) processing method based on maximum likelihood (ML) estimation of organ activities using 3D organ VOIs and a projector that models the image degrading effects. Both a physical phantom experiment and Monte Carlo simulation (MCS) studies were used to evaluate the new method. In these studies, the accuracies and precisions of organ activity estimates for the QPlanar method were compared with those from conventional planar (CPlanar) processing methods with various corrections for scatter, attenuation and organ overlap, and a quantitative SPECT (QSPECT) processing method. Experimental planar and SPECT projections and registered CT data from an RSD Torso phantom were obtained using a GE Millenium VH/Hawkeye system. The MCS data were obtained from the 3D NCAT phantom with organ activity distributions that modelled the uptake of (111)In ibritumomab tiuxetan. The simulations were performed using parameters appropriate for the same system used in the RSD torso phantom experiment. The organ activity estimates obtained from the CPlanar, QPlanar and QSPECT methods from both experiments were compared. From the results of the MCS experiment, even with ideal organ overlap correction and background subtraction, CPlanar methods provided limited quantitative accuracy. The QPlanar method with accurate modelling of the physical factors increased the quantitative accuracy at the cost of requiring estimates of the organ VOIs in 3D. The accuracy of QPlanar approached that of QSPECT, but required much less acquisition and computation time. Similar results were obtained from the physical phantom experiment. We conclude that the QPlanar method, based on 3D organ VOIs and accurate models of the projection process, provided a substantial increase in accuracy of organ activity estimates from planar images compared to CPlanar processing and had accuracy approaching that of QSPECT.

Algorithms↗

PubMed2006

Model-based compensation for quantitative 123I brain SPECT imaging.

Previously we have developed a model-based method that can accurately estimate downscatter contamination from high-energy photons in 123I imaging. In this work we combined the model-based method with iterative reconstruction-based compensations for other image-degrading factors such as attenuation, scatter, the collimator-detector response function (CDRF) and partial volume effects to form a comprehensive method for performing quantitative 123I SPECT image reconstruction. In the model-based downscatter estimation method, photon scatter inside the object was modelled using the effective source scatter estimation (ESSE) technique, including contributions from all the photon emissions. The CDRFs, including the penetration and scatter components due to the high-energy 123I photons, were estimated using Monte Carlo (MC) simulations of point sources in air at various distances from the face of the collimator. The downscatter contamination was then compensated for during the iterative reconstruction by adding the estimated results to the projection steps. The model-based downscatter compensation (MBDC) was evaluated using MC simulated and experimentally acquired projection data. From the MC simulation, we found about 39% of the total counts in the energy window of 123I were attributed to the downscatter contamination, which reduced image contrast and caused a 1.5% to 10% overestimation of activities in various brain structures. Model-based estimates of the downscatter contamination were in good agreement with the simulated data. Compensation using MBDC removed the contamination and improved the image contrast and quantitative accuracy to that of the images obtained from 159 keV photons. The errors in absolute quantitation were reduced to within +/-3.5%. The striatal specific binding potential calculated based on the activity ratio to the background was also improved after MBDC. The errors were reduced from -4.5% to -10.93% without compensation to -0.55% to 4.87% after compensation. The model-based method provided accurate downscatter estimation and, when combined with iterative reconstruction-based compensations, accurate quantitation was obtained with minimal loss of precision.

Brain↗

PubMed2006

Three-class ROC analysis--a decision theoretic approach under the ideal observer framework.

Receiver operating characteristic (ROC) analysis is well established in the evaluation of systems involving binary classification tasks. However, medical tests often require distinguishing among more than two diagnostic alternatives. The goal of this work was to develop an ROC analysis method for three-class classification tasks. Based on decision theory, we developed a method for three-class ROC analysis. In this method, the objects were classified by making the decision that provided the maximal utility relative to the other two. By making assumptions about the magnitudes of the relative utilities of incorrect decisions, we found a decision model that maximized the expected utility of the decisions when using log-likelihood ratios as decision variables. This decision model consists of a two-dimensional decision plane with log likelihood ratios as the axes and a decision structure that separates the plane into three regions. Moving the decision structure over the decision plane, which corresponds to moving the decision threshold in two-class ROC analysis, and computing the true class 1, 2, and 3 fractions defined a three-class ROC surface. We have shown that the resulting three-class ROC surface shares many features with the two-class ROC curve; i.e., using the log likelihood ratios as the decision variables results in maximal expected utility of the decisions, and the optimal operating point for a given diagnostic setting (set of relative utilities and disease prevalences) lies on the surface. The volume under the three-class surface (VUS) serves as a figure-of-merit to evaluate different data acquisition systems or image processing and reconstruction methods when the assumed utility constraints are relevant.

Algorithms↗

PubMed2006

Three-class ROC analysis--the equal error utility assumption and the optimality of three-class ROC surface using the ideal observer.

Previously, we have developed a decision model for three-class receiver operating characteristic (ROC) analysis based on decision theory. The proposed decision model maximizes the expected decision utility under the assumption that incorrect decisions have equal utilities under the same hypothesis (equal error utility assumption). This assumption reduced the dimensionality of the "general" three-class ROC analysis and provided a practical figure-of-merit to evaluate the three-class task performance. However, it also limits the generality of the resulting model because the equal error utility assumption will not apply for all clinical three-class decision tasks. The goal of this study was to investigate the optimality of the proposed three-class decision model with respect to several other decision criteria. In particular, besides the maximum expected utility (MEU) criterion used in the previous study, we investigated the maximum-correctness (MC) (or minimum-error), maximum likelihood (ML), and Nyman-Pearson (N-P) criteria. We found that by making assumptions for both MEU and N-P criteria, all decision criteria lead to the previously-proposed three-class decision model. As a result, this model maximizes the expected utility under the equal error utility assumption, maximizes the probability of making correct decisions, satisfies the N-P criterion in the sense that it maximizes the sensitivity of one class given the sensitivities of the other two classes, and the resulting ROC surface contains the maximum likelihood decision operating point. While the proposed three-class ROC analysis model is not optimal in the general sense due to the use of the equal error utility assumption, the range of criteria for which it is optimal increases its applicability for evaluating and comparing a range of diagnostic systems.

Algorithms↗

PubMed2005

A Monte Carlo and physical phantom evaluation of quantitative In-111 SPECT.

Accurate estimation of the 3D in vivo activity distribution is important for dose estimation in targeted radionuclide therapy (TRT). Although SPECT can potentially provide such estimates, SPECT without compensation for image degrading factors is not quantitatively accurate. In this work, we evaluated quantitative SPECT (QSPECT) reconstruction methods that include compensation for various physical effects. Experimental projection data were obtained using a GE VH/Hawkeye system and an RSD torso phantom. Known activities of In-111 chloride were placed in the lungs, liver, heart, background and two spherical compartments with inner diameters of 22 mm and 34 mm. The 3D NCAT phantom with organ activities based on clinically derived In-111 ibritumomab tiuxetan data was used for the Monte Carlo (MC) simulation studies. Low-noise projection data were simulated using previously validated MC simulation methods. Fifty sets of noisy projections with realistic count levels were generated. Reconstructions were performed using the OS-EM algorithm with various combinations of attenuation (A), scatter (S), geometric response (G), collimator-detector response (D) and partial volume compensation (PVC). The QSPECT images from the various combinations of compensations were evaluated in terms of the accuracy and precision of the estimates of the total activity in each organ. For experimental data, the errors in organ activities for ADS and PVC compensation were less than 6.5% except the smaller sphere (-11.9%). For the noisy simulated data, the errors in organ activity for ADS compensation were less than 5.5% except the lungs (20.9%) and blood vessels (15.2%). Errors for other combinations of compensations were significantly (A, AS) or somewhat (AGS) larger. With added PVC, the error in the organ activities improved slightly except for the lungs (11.5%) and blood vessels (3.6%) where the improvement was more substantial. The standard deviation/mean ratios were all less than 1.5%. We conclude that QSPECT methods with appropriate compensations provided accurate In-111 organ activity estimates. For the collimator used, AGS was almost as good as ADS and may be preferable due to the reduced reconstruction time. PVC was important for small structures such as tumours or for organs in close proximity to regions with high activity. The improved quantitative accuracy from QSPECT methods has the potential for improving organ dose estimations in TRT.

Algorithms↗

PubMed2005

Partial volume effect compensation for quantitative brain SPECT imaging.

Partial volume (PV) effects degrade the quantitative accuracy of SPECT brain images. In this paper, we extended a PV compensation (PVC) method originally developed for brain PET, the geometric transfer matrix (GTM) method, to brain SPECT using iterative reconstruction-based compensations. In the GTM method a linear transform between the true regional activities and the measured results was assumed. Elements of the GTM were calculated by projecting and reconstructing maps with uniform regions representing different structures. However, with iterative reconstruction methods, especially when reconstruction-based compensation for detector response was applied, we found that it was important to treat the region maps as a perturbation to the reconstructed image in the estimation of the GTM. This modified method, termed perturbation-based GTM (pGTM) was evaluated using Monte Carlo (MC) simulated and experimentally acquired data. Results showed great improvement of the quantitative accuracy in brain SPECT imaging. For MC simulated data, PVC using pGTM reduced the underestimation of striatal activities from 30% to less than 1.2%. For experimental data, PVC using pGTM reduced the underestimation of striatal activities from 36% to less than 7.8%. The underestimation of the striatum to background activity ratio was also improved from 31% to 2.7%.

Algorithms↗

PubMed2005

Effect of attenuation correction on lesion detection using a hybrid PET system.

OBJECTIVE: The purpose of this study was to investigate the effect of attenuation correction (AC) on lesion detection for a hybrid PET system. MATERIAL AND METHOD: Experimental list-mode data were acquired from hot spheres inside a uniform cylindrical phantom with an elliptical cross-section using a Siemens E. CAM+ dual-camera hybrid PET system. Spheres with inner diameters of 0.8- and 1-cm and the cylindrical phantom were filled with F-18 to simulate lesions with lesion-to-background (L/B) ratios of 14:1 and 8:1, respectively, found in clinical PET studies. The list-mode data of each sphere size were regrouped into sinograms with peak-to-peak energy window settings at 30% and 20% for the 0.8- and 1-cm diameter lesion, respectively. They were then rebinned using the single slice rebinning method. Attenuation correction was applied assuming uniform attenuation. The sinograms with and without AC were reconstructed using 5 iterations of OS-EM algorithm with 8 angles/ subset and postfiltered with a Butterworth filter with n = 5 and fc = 0.52 cycles/cm. Human observer performance study and localization receiver operating characteristic (LROC) analysis were used to evaluate the reconstructed images for maximum lesion detection. Average areas under the LROC curves (A(LROC)) across 8 observers obtained with and without AC were determined. The null hypothesis that there was no difference between with AC and without AC was tested using a two-tailed t-test with 95% confidence interval. RESULTS: The results indicated that for the 0. 8-cm lesion with 14:1 L/B ratio, the A(LROC) decreases from 0.66 to 0.62 when AC is applied as compared to without AC andfrom 0.69 to 0.63 for the 1.0-cm lesion with 8:1 L/ B ratio, but no statistical significant difference (p > 0. 05). CONCLUSION: The authors conclude that for a phantom with hot lesions embedded in a uniform background, AC decreases lesion detectability compared to without AC using a hybrid PET system for small lesion sizes.

Humans↗

PubMed2004

Assessment of scatter compensation strategies for (67)Ga SPECT using numerical observers and human LROC studies.

UNLABELLED: 67Ga citrate is an oncologic SPECT imaging agent often used to diagnose or stage patients with non-Hodgkin's lymphoma. As (67)Ga decay involves the emission of multiple-energy gamma-rays, significant photon downscatter will be present within each photopeak energy window. We have previously shown that the inclusion of these scattered photons significantly degrades lesion detectability. The goal of this study was to investigate the extent to which this decrease in detectability can be reversed by applying scatter compensation strategies. METHODS: We have compared 5 different scatter compensation methods to the case of no scatter compensation in iterative SPECT image reconstruction. The strategies consisted of (a). perfect scatter rejection, (b). ideal scatter compensation, (c). triple-energy window (TEW) scatter estimation, (d). effective scatter source estimation (ESSE), and (e). postreconstruction scatter subtraction. Reconstruction parameters used for each method were first optimized using a channelized Hotelling numerical observer. Strategies were then ranked in terms of lesion detectability using a human observer localization receiver operating characteristic (LROC) study. An additional comparison was made comparing the human LROC rankings with a recently developed channelized nonprewhitening (CNPW) LROC numerical observer. RESULTS: Using the area-under-the-LROC-curve (A(LROC)) as the assessment criterion, our results indicate that the TEW and ESSE scatter compensation methods are able to significantly improve lesion detectability over no compensation (A(LROC) = 0.75 and 0.73 vs. 0.67, respectively). However, these compensations failed to achieve the same detectability as perfect scatter rejection (A(LROC) = 0.84). Both ideal scatter compensation and postreconstruction scatter subtraction resulted in numerical increases in detection accuracy that were not statistically significant from no scatter compensation. Good agreement is seen between the CNPW observer and human LROC studies (Spearman rank order coefficient, r(s) = 0.74), thus indicating that the LROC observer may be a good predictor of human observer performance in (67)Ga SPECT. CONCLUSION: Scatter compensation in (67)Ga SPECT imaging using techniques such as TEW or ESSE is able to improve lesion detectability compared with no scatter compensation. A recently developed numerical observer model appears to be a good predictor of human observer performance and may be used to perform imaging optimizations, thereby reducing the need for human LROC studies.

Citrates↗

PubMed2002

Application of task-based measures of image quality to optimization and evaluation of three-dimensional reconstruction-based compensation methods in myocardial perfusion SPECT.

In this paper, we apply the channelized Hotelling observer (CHO) using a defect detection task to the optimization and evaluation of three-dimensional iterative reconstruction-based compensation methods for myocardial perfusion single-photon emission computed tomography (SPECT). We used a population of 24 mathematical cardiac-torso phantoms that realistically model the activity and attenuation distribution in three classes of patients: females, and males with flat diaphragms and raised diaphragms. Projection data were generated and subsequently reconstructed using methods based on the ordered subsets-expectation maximization (OSEM) algorithm. The methods evaluated included compensation for attenuation, detector response blurring, and scatter in various combinations. We applied the CHO to optimize the number of iterations for OSEM and the cutoff frequency and order of a three-dimensional postreconstruction Butterworth filter. Using the optimal parameters, we then compared the compensation methods. The index of comparison in these studies was the area under the receiver operating characteristics curve (AUC) for the CHO. We found that attenuation compensation with either detector response or scatter compensation gave statistically significant increases in the AUC compared to attenuation compensation alone. The greatest increase in the AUC occurred when all three compensations were applied. These results indicate that compensation for detector response and scatter, in addition to attenuation compensation, will improve defect detectability in myocardial SPECT images.

Algorithms↗

PubMed2002

Efficient simulation of SPECT down-scatter including photon interactions with crystal and lead.

A major image degrading factor in simultaneous Dual Isotope (DI) SPECT or simultaneous Emission-Transmission (ECT-TCT) imaging, is the detection of photons emitted by the higher energy isotope in the energy window used for imaging the lower energy isotope. In Tc-99m/Tl-201 DI-SPECT typically tens of percents of the total detected down-scatter is caused by lead x rays. In Tc-99m/Gd-153 ECT-TCT, a comparable fraction of the down-scatter originates from Tc-99m photons which only partly deposit their energy in the detector crystal (i.e., due to crystal interactions). Efficient simulation methods which model down-scatter can be used to optimize DI-SPECT or ECT-TCT imaging acquisition or reconstruction protocols. In this paper we adapt a previously proposed efficient down-scatter simulation method, to include the interactions of photons with the detector crystal and collimator lead. To this end, point spread function tables including crystal and lead interactions are precalculated. Subsequently, photons are traced through the patient body until their last scatter position, and the precalculated responses are used to project the photons onto the detector plane, while photon attenuation in the patient is taken into account. The approach is evaluated by comparing simulated Tc-99m down-scatter projections with measured projections. Incorporation of photon interaction with crystal and lead leads to significantly improved accuracy of the shape of down-scatter responses, while differences in total counts between simulated and measured projections typically decrease from tens of percents to a couple of percents. Calculating 60 down-scatter projections of an extended distribution on a 64 x 64 x 64 grid takes about three minutes on a PC with two 1.2 GHz processors. We conclude that accurate and efficient simulation of down-scatter is now possible including the major effects of the nonuniform mass density of the patient as well as photon interactions with the crystal and collimator lead.

Biophysical Phenomena↗

PubMed2002

Optimum compensation method and filter cutoff frequency in myocardial SPECT: a human observer study.

UNLABELLED: Attenuation, photon scatter, and distance-dependent collimator-detector response are major degrading factors in myocardial SPECT images. The current study investigated whether compensation for these factors improves perfusion defect detectability, and compared the results for human observers with a previous study using a mathematical observer. METHODS: Four methods were investigated: attenuation compensation (AC); attenuation and detector response compensation; attenuation and scatter compensation; and attenuation, detector response, and scatter compensation (ADSC). For ADSC, 4 three-dimensional postreconstruction Butterworth filter cutoff frequencies were investigated for a pixel size of 0.62 cm: 0.12, 0.14, 0.16, and 0.22 pixel(-1). Five observers read images reconstructed using the 4 compensation methods. Receiver operating characteristics (ROC) analysis was used to determine the area under the ROC curve in each treatment studied. RESULTS: Reconstruction methods that incorporated scatter and detector response compensation had higher indices of detectability than AC alone. Over the range studied, a filter cutoff frequency of 0.14 pixel(-1) was optimal. A comparison of human observer results with an earlier channelized Hotelling observer study performed with the same images showed excellent agreement in trend and ranking of defect detectability. CONCLUSION: Compensation for detector response and scatter improves defect detectability compared with AC alone, although detectability may depend on phantom population choice and noise level. An optimal filter cutoff was found that is lower than what is typically used in a clinical setting. The channelized Hotelling observer is a good predictor of human observer performance and may reduce the need for tedious, time-consuming studies with human observers.

Area Under Curve↗

PubMedDate not supplied

Development of a 4-D digital mouse phantom for molecular imaging research.

PURPOSE: We develop a realistic and flexible 4-D digital mouse phantom and investigate its usefulness in molecular imaging research. METHODS: Organ shapes were modeled with non-uniform rational B-spline (NURBS) surfaces based on high-resolution 3-D magnetic resonance microscopy (MRM) data. Cardiac and respiratory motions were modeled based on gated magnetic resonance imaging (MRI) data obtained from normal mice. Pilot simulation studies in single-photon emission computed tomography (SPECT) and X-ray computed tomography (CT) were performed to demonstrate the utility of the phantom. RESULTS: NURBS are an efficient and flexible way to accurately model the anatomy and cardiac and respiratory motions for a realistic 4-D digital mouse phantom. The phantom is capable of producing realistic molecular imaging data from which imaging devices and techniques can be evaluated. CONCLUSION: The phantom provides a unique and useful tool in molecular imaging research. It can be used in the development of new imaging instrumentation, image acquisition strategies, and image processing and reconstruction methods.

Animals↗

PubMedDate not supplied

Comparison of 180 degrees and 360 degrees acquisition for myocardial perfusion SPECT with compensation for attenuation, detector response, and scatter: Monte Carlo and mathematical observer results.

BACKGROUND: The optimal projection data acquisition strategy for myocardial perfusion (MP) single photon emission computed tomography (SPECT) remains controversial. METHODS: We compared MP SPECT using 180 degrees and 360 degrees projection data obtained with the same acquisition time, reconstructed either with filtered back projection (FBP) or the iterative ordered-subsets expectation maximization (OS-EM) algorithm with various combinations of attenuation, detector response, and scatter compensation using mathematical observers and a myocardial defect detection task. We used Monte Carlo-simulated projection data from a population of 3-dimensional nurbs-based cardiac-torso (NCAT) phantoms with ranges of variability in patient anatomy, organ uptake, defect location, defect size, and noise level based on clinical data. Projection data from 180 degrees and 360 degrees acquisitions were generated by assuming the same acquisition time. After iterative or FBP reconstruction, standard postprocessing methods were applied. For each acquisition and reconstruction method, we optimized the number of iterations and cut-off frequency of the Butterworth filter using the Channelized Hotelling Observer methodology. The optimum set of parameters was that which gave the maximum area under the curve. RESULTS: For both acquisition protocols, OS-EM with compensations provided better performance than FBP or OS-EM without compensation. For FBP, the optimized 180 degrees acquisition provided a statistically significant increase in AUC as compared with optimized 360 degrees acquisition. For OS-EM, the AUCs for 180 degrees were slightly larger than for 360 degrees acquisitions when comparing images reconstructed with the same compensations. However, the differences were smaller and not statistically significant. CONCLUSION: With optimized reconstruction and filtering parameters, 180 degrees acquisition provided a statistically significant improvement over 360 degrees acquisition for FBP reconstruction. However, for OS-EM the differences were small and not statistically significant.

Algorithms↗