The influence of attenuation and scatter compensation on the apparent distribution of Tc-99m sestamibi in cardiac slices.

BACKGROUND: Our objective was to study the differences in relative count distributions in the left ventricular walls with attenuation compensation (AC) versus AC and triple-energy-window scatter compensation (SC), compared with standard filtered backprojection (FBP). METHODS AND RESULTS: Two hundred patients identified as having normal cardiac perfusion with FBP after undergoing either pharmacologically or physiologically induced stress were included in this study. Projection data were reconstructed with FBP, 10 iterations of ordered-subset expectation-maximization (OSEM) with AC, and OSEM with AC+SC. A comparison was made of average percentage of maximum counts within each of 9 regions of CEqual (Marconi Medical Systems, Inc, Cleveland, Ohio) polar maps (ie, the apex, 4 midventricular regions, and 4 basal regions). Compared with OSEM(AC), a slight decrease at the apex exists when SC is included. The elevated inferior-to-anterior count ratio in the midventricular and basal regions noted with OSEM(AC) decreased to close to 1.0 with OSEM(AC+SC). The anterior-to-lateral ratio for both regions was closest to 1.0 for OSEM(AC+SC). In the midventricular region, the lateral-to-septal ratio decreased further below 1.0 with OSEM(AC+SC) than it did with OSEM(AC). This was the only basal ratio not to improve to close to 1.0 with OSEM(AC+SC). In a subset of patients identified at the time of clinical reading as having a possible attenuation-caused decrease in the inferior region, AC elevated the inferior-to-anterior ratio to above 1.0 for the midventricular region. AC+SC resulted in a ratio of near 1.0 for this region. In another subset of patients identified as having anterior attenuation artifacts, compensation methods (either AC or AC+SC) failed to show an improvement compared with FBP. CONCLUSIONS: AC and SC improve the uniformity of the polar map, especially by bringing the inferior-to-anterior ratio closer to 1.0. Further investigation is necessary to determine the cause of the increased midventricular septal polar map count. In addition, the subset of patients identified as having breast-like attenuation artifacts causing a decreased polar map count in the anterior wall (relative to the inferior wall) also needs further attention.

An interior point iterative maximum-likelihood reconstruction algorithm incorporating upper and lower bounds with application to SPECT transmission imaging.

The algorithm we consider here is a block-iterative (or ordered subset) version of the interior point algorithm for transmission reconstruction. Our algorithm is an interior point method because each vector of the iterative sequence [x(k)], k = 0, 1, 2, ... satisfies the constraints a(j) < x(j)k < b(j), j = 1, ..., J. Because it is a block-iterative algorithm that reconstructs the transmission attenuation map and places constraints above and below the pixel values of the reconstructed image, we call it the BITAB method. Computer simulations using the three-dimensional mathematical cardiac and torso phantom, reveal that the BITAB algorithm in conjunction with reasonably selected prior upper and lower bounds has the potential to improve the accuracy of the reconstructed attenuation coefficients from truncated fan beam transmission projections. By suitably selecting the bounds, it is possible to restrict the over estimation of coefficients outside the fully sampled region, that results from reconstructing truncated fan beam projections with iterative transmission algorithms such as the maximum-likelihood gradient type algorithm.

Improved image quality and computation reduction in 4-D reconstruction of cardiac-gated SPECT images.

Spatiotemporal reconstruction of cardiac-gated SPECT images permits us to obtain valuable information related to cardiac function. However, the task of reconstructing this four-dimensional (4-D) data set is computation intensive. Typically, these studies are reconstructed frame-by-frame: a nonoptimal approach because temporal correlations in the signal are not accounted for. In this work, we show that the compression and signal decorrelation properties of the Karhunen-Loève (KL) transform may be used to greatly simplify the spatiotemporal reconstruction problem. The gated projections are first KL transformed in the temporal direction. This results in a sequence of KL-transformed projection images for which the signal components are uncorrelated along the time axis. As a result, the 4-D reconstruction task is simplified to a series of three-dimensional (3-D) reconstructions in the KL domain. The reconstructed KL components are subsequently inverse KL transformed to obtain the entire spatiotemporal reconstruction set. Our simulation and clinical results indicate that KL processing provides image sequences that are less noisy than are conventional frame-by-frame reconstructions. Additionally, by discarding high-order KL components that are dominated by noise, we can achieve savings in computation time because fewer reconstructions are needed in comparison to conventional frame-by-frame reconstructions.

LROC analysis of detector-response compensation in SPECT.

Localization ROC (LROC) observer studies examined whether detector response compensation (DRC) in ordered-subset, expectation-maximization (OSEM) reconstructions helps in the detection and localization of hot tumors. Simulated gallium (Ga-67) images of the thoracic region were used in the study. The projection data modeled the acquisition of attenuated 93- and 185-keV photons with a medium-energy parallel-hole collimator, but scatter was not modeled. Images were reconstructed with five strategies: 1) OSEM with no DRC; 2) OSEM preceded by restoration filtering; 3) OSEM with iterative DRC; 4) OSEM with an ideal DRC; and 5) filtered backprojection (FBP) with no DRC. All strategies included attenuation correction. There were four LROC studies conducted. In a study using a single tumor activity, the ideal DRC offered the best performance, followed by iterative DRC, restoration filtering, OSEM with no DRC, and FBP. Statistical significance at the 5% level was found between all pairs of strategies except for restoration filtering and OSEM with no DRC. A similar ranking was found for a more realistic study using multiple tumor activities. Additional studies considered the effects of OSEM iteration number and tumor activity on the detection improvement that iterative DRC offered with respect to OSEM with no DRC.