Evaluation of scatter compensation methods by their effects on parameter estimation from SPECT projections

Correction of multiplexing artefacts in multi-pinhole SPECT through temporal shuttering, de-multiplexing of projections, and alternating reconstruction

Objective.Single-photon emission computed tomography (SPECT) with pinhole collimators can provide high-resolution imaging, but is often limited by low sensitivity. Acquiring projections simultaneously through multiple pinholes affords both high resolution and high sensitivity. However, the overlap of projections from different pinholes on detectors, known as multiplexing, has been shown to cause artefacts which degrade reconstructed images.Approach.Multiplexed projection sets were considered here using an analytic simulation model of AdaptiSPECT-C-a brain-dedicated multi-pinhole SPECT system. AdaptiSPECT-C has fully adaptable aperture shutters, so can acquire projections with a combination of multiplexed and non-multiplexed frames using temporal shuttering. Two strategies for reducing multiplex artefacts were considered: an algorithm to de-multiplex projections, and an alternating reconstruction strategy for projections acquired with a combination of multiplexed and non-multiplexed frames. Geometric and anthropomorphic digital phantoms were used to assess a number of metrics.Main results.Both de-multiplexing strategies showed a significant reduction in image artefacts and improved fidelity, image uniformity, contrast recovery and activity recovery (AR). In all cases, the two de-multiplexing strategies resulted in superior metrics to those from images acquired with only mux-free frames. The de-multiplexing algorithm provided reduced image noise and superior uniformity, whereas the alternating strategy improved contrast and AR.Significance.The use of these de-multiplexing algorithms means that multi-pinhole SPECT systems can acquire projections with more multiplexing without degradation of images.

Resolution Enhancement in PET Reconstruction Using Collimation

Collimation can improve both the spatial resolution and sampling properties compared to the same scanner without collimation. Spatial resolution improves because each original crystal can be conceptually split into two (i.e., doubling the number of in-plane crystals) by masking half the crystal with a high-density attenuator (e.g., tungsten); this reduces coincidence efficiency by 4× since both crystals comprising the line of response (LOR) are masked, but yields 4× as many resolution-enhanced (RE) LORs. All the new RE LORs can be measured by scanning with the collimator in different configurations. In this simulation study, the collimator was assumed to be ideal, neither allowing gamma penetration nor truncating the field of view. Comparisons were made in 2D between an uncollimated small-animal system with 2-mm crystals that were assumed to be perfectly absorbing and the same system with collimation that narrowed the effective crystal size to 1 mm. Digital phantoms included a hot-rod and a single-hot-spot, both in a uniform background with activity ratio of 4:1. In addition to the collimated and uncollimated configurations, angular and spatial wobbling acquisitions of the 2-mm case were also simulated. Similarly, configurations with different combinations of the RE LORs were considered including (i) all LORs, (ii) only those parallel to the 2-mm LORs; and (iii) only cross pairs that are not parallel to the 2-mm LORs. Lastly, quantitative studies were conducted for collimated and uncollimated data using contrast recovery coefficient and mean-squared error (MSE) as metrics. The reconstructions show that for most noise levels there is a substantial improvement in image quality (i.e., visual quality, resolution, and a reduction in artifacts) by using collimation even when there are 4× fewer counts or - in some cases - comparing with the noiseless uncollimated reconstruction. By comparing various configurations of sampling, the results show that it is the matched combination of both improved spatial resolution of each LOR and the increase in the number of LORs that yields improved reconstructions. Further, the quantitative studies show that for low-count scans, the collimated data give better MSE for small lesions and the uncollimated data give better MSE for larger lesions; for high-count studies, the collimated data yield better quantitative values for the entire range of lesion sizes that were evaluated.

Rapid Optimization of SPECT Scatter Correction Using Model LROC Observers

The problem we address is the optimization and comparison of window-based scatter correction (SC) methods in SPECT for maximum a posteriori reconstructions. While sophisticated reconstruction-based SC methods are available, the commonly used window-based SC methods are fast, easy to use, and perform reasonably well. Rather than subtracting a scatter estimate from the measured sinogram and then reconstructing, we use an ensemble approach and model the mean scatter sinogram in the likelihood function. This mean scatter sinogram estimate, computed from satellite window data, is itself inexact (noisy). Therefore two sources of noise, that due to Poisson noise of unscattered photons and that due to the model error in the scatter estimate, are propagated into the reconstruction. The optimization and comparison is driven by a figure of merit, the area under the LROC curve (ALROC) that gauges performance in a signal detection plus localization task. We use model observers to perform the task. This usually entails laborious generation of many sample reconstructions, but in this work, we instead develop a theoretical approach that allows one to rapidly compute ALROC given known information about the imaging system and the scatter correction scheme. A critical step in the theory approach is to predict additional (above that due to to the propagated Poisson noise of the primary photons) contributions to the reconstructed image covariance due to scatter (model error) noise. Simulations show that our theory method yields, for a range of search tolerances, LROC curves and ALROC values in close agreement to that obtained using model observer responses obtained from sample reconstruction methods. This opens the door to rapid comparison of different window-based SC methods and to optimizing the parameters (including window placement and size, scatter sinogram smoothing kernel) of the SC method.

The variation of intrinsic spatial resolution across the UFOV of scintillation cameras

The aim of the present study was to investigate in detail the variation of the intrinsic spatial resolution across the useful field of view (UFOV) of gamma-cameras and to explore whether this variation could lead to observable effects in clinical images. Two gamma-cameras were used, without their collimators, to acquire 560 (99m)Tc point source images at different points across their UFOVs, in order to measure the intrinsic spatial resolution at each point. Possible clinical effects of the resolution variation were examined on images of a thyroid phantom using a LEHR collimator, acquired at different locations on the UFOV and at various distances from the collimator. The gamma-camera intrinsic resolution varied significantly across the UFOV, being generally lower at the central region and deteriorating at the edges. Pronounced local maxima and minima were found at points corresponding to the centers of the photomultiplier tubes (PMTs) and halfway in between. Maximum differences of more than 50% were observed between the points presenting the best and worst intrinsic resolution. Differences between neighboring points reached 17%. The effects of resolution variation were clearly observable on the thyroid phantom images. It was concluded that an appropriate correction algorithm might be necessary in order to correct for the variation of the intrinsic spatial resolution across the UFOV of gamma-cameras.

Collimator optimization for detection and quantitation tasks: application to gallium-67 imaging

We describe a new approach to the problem of collimator optimization in nuclear medicine; our methodology is illustrated for the challenging case of gallium-67 imaging. Collimator-design methods based on empirical rules, such as specification of an allowable level of single-septal penetration (SSP) at a fixed energy, are especially inappropriate for radionuclides characterized by an abundance of high-energy contaminant photons that scatter in the patient, collimator, and/or detector before detection within one of a few photopeak energy windows. Lead X-rays produced in the collimator are an additional source of contamination. We designed optimal collimation for 67Ga based on relevant clinical imaging tasks and a realistic simulation of photon transport in a phantom, collimator, and detector. Collimator designs were compared on the basis of performance in lesion detection, as predicted by a three-channel Hotelling observer (CHO), as well as in tumor and background activity estimation (EST), quantified by task-specific signal-to-noise ratios (SNRs). The optimal values of collimator lead content were 22.0 and 23.8 g/cm2, respectively, for CHO and EST, while the optimal geometric resolution values were 1.8 and 1.6 cm full-width at half-maximum (FWHM), respectively, at a distance of 23.5 cm. The resolution of a commercially available medium-energy low-penetration collimator (MELP) is 1.9 cm FWHM at this distance. The optimal values for SSP at 300 keV were 7.3% and 5.8% based on CHO and EST, respectively, compared to 5.2% for the MELP collimator. Compared with the commercial MELP collimator, the 67Ga collimator optimized for tumor detection or activity estimation tasks provided improved geometric spatial resolution with reduced geometric efficiency and, surprisingly, allowed an increased level of single-septal penetration.

Scatter modelling and compensation in emission tomography

In nuclear medicine, clinical assessment and diagnosis are generally based on qualitative assessment of the distribution pattern of radiotracers used. In addition, emission tomography (SPECT and PET) imaging methods offer the possibility of quantitative assessment of tracer concentration in vivo to quantify relevant parameters in clinical and research settings, provided accurate correction for the physical degrading factors (e.g. attenuation, scatter, partial volume effects) hampering their quantitative accuracy are applied. This review addresses the problem of Compton scattering as the dominant photon interaction phenomenon in emission tomography and discusses its impact on both the quality of reconstructed clinical images and the accuracy of quantitative analysis. After a general introduction, there is a section in which scatter modelling in uniform and non-uniform media is described in detail. This is followed by an overview of scatter compensation techniques and evaluation strategies used for the assessment of these correction methods. In the process, emphasis is placed on the clinical impact of image degradation due to Compton scattering. This, in turn, stresses the need for implementation of more accurate algorithms in software supplied by scanner manufacturers, although the choice of a general-purpose algorithm or algorithms may be difficult.

A new design for a SPECT small-animal imager

We demonstrate, using computer models, the feasibility of a new SPECT system for imaging small animals such as mice. This system consists of four modular scintillation cameras, four multiple-pinhole apertures, electronics, and tomographic reconstruction software. All of these constituents have been designed in our laboratory. The cameras are 120mm×120mm with a resolution of approximately 2mm, the apertures can have either single or multiple pinholes, and reconstruction is performed using the OS-EM algorithm. One major advantage of this system is the design flexibility it offers, as the cameras are easy to move and the aperture s are simple to modify. We explored a number of possible configurations. One promising configuration had the four camera faces forming four sides of a cube with multiple-pinhole apertures employed to focus the incoming high-energy photons. This system is rotated three times, so that data are collected from a total of sixteen camera angles. It is shown that this hybrid system has some superior properties to single-aperture-type systems. We conclude that this proposed system offers advantages over current imaging systems in terms of flexibility, simplicity, and performance.