Simultaneous maximum a posteriori reconstruction of attenuation andactivity distributions from emission sinograms

In order to perform attenuation correction in emission tomography an attenuation map is required. We propose a new method to compute this map directly from the emission sinogram, eliminating the transmission scan from the acquisition protocol. The problem is formulated as an optimization task where the objective function is a combination of the likelihood and an a priori probability. The latter uses a Gibbs prior distribution to encourage local smoothness and a multimodal distribution for the attenuation coefficients. Since the attenuation process is different in positron emission tomography (PET) and single photon emission tomography (SPECT), a separate algorithm for each case is derived. The method has been tested on mathematical phantoms and on a few clinical studies. For PET, good agreement was found between the images obtained with transmission measurements and those produced by the new algorithm in an abdominal study. For SPECT, promising simulation results have been obtained for nonhomogeneous attenuation due to the presence of the lungs.     

An EM algorithm for estimating SPECT emission and transmission parameters from emissions data only

A maximum-likelihood (ML) expectation-maximization (EM) algorithm (called EM-IntraSPECT) is presented for simultaneously estimating single photon emission computed tomography (SPECT) emission and attenuation parameters from emission data alone. The algorithm uses the activity within the patient as transmission tomography sources, with which attenuation coefficients can be estimated. For this initial study, EM-IntraSPECT was tested on computer-simulated attenuation and emission maps representing a simplified human thorax as well as on SPECT data obtained from a physical phantom. Two evaluations were performed. First, to corroborate the idea of reconstructing attenuation parameters from emission data, attenuation parameters (mu) were estimated with the emission intensities (lambda) fixed at their true values. Accurate reconstructions of attenuation parameters were obtained. Second, emission parameters lambda and attenuation parameters mu were simultaneously estimated from the emission data alone. In this case there was crosstalk between estimates of lambda and mu and final estimates of lambda and mu depended on initial values. Estimates degraded significantly as the support extended out farther from the body, and an explanation for this is proposed. In the EM-IntraSPECT reconstructed attenuation images, the lungs, spine, and soft tissue were readily distinguished and had approximately correct shapes and sizes. As compared with standard EM reconstruction assuming a fix uniform attenuation map, EM-IntraSPECT provided more uniform estimates of cardiac activity in the physical phantom study and in the simulation study with tight support, but less uniform estimates with a broad support. The new EM algorithm derived here has additional applications, including reconstructing emission and transmission projection data under a unified statistical model.     

Modeling the axial extension of a transmission line source within iterative reconstruction via multiple transmission sources

Reconstruction algorithms for transmission tomography have generally assumed that the photons reaching a particular detector bin at a particular angle originate from a single point source. In this paper, we highlight several cases of extended transmission sources, in which it may be useful to approach the estimation of attenuation coefficients as a problem involving multiple transmission point sources. Examined in detail is the case of a fixed transmission line source with a fan-beam collimator. This geometry can result in attenuation images that have significant axial blur. Herein it is also shown, empirically, that extended transmission sources can result in biased estimates of the average attenuation, and an explanation is proposed. The finite axial resolution of the transmission line source configuration is modeled within iterative reconstruction using an expectation-maximization algorithm that was previously derived for estimating attenuation coefficients from single photon emission computed tomography (SPECT) emission data. The same algorithm is applicable to both problems because both can be thought of as involving multiple transmission sources. It is shown that modeling axial blur within reconstruction removes the bias in the average estimated attenuation and substantially improves the axial resolution of attenuation images.     

Monotonic penalized-likelihood image reconstruction for X-ray fluorescence computed tomography

In this paper, we derive a monotonic penalized-likelihood algorithm for image reconstruction in X-ray fluorescence computed tomography (XFCT) when the attenuation maps at the energies of the fluorescence X-rays are unknown. In XFCT, a sample is irradiated with pencil beams of monochromatic synchrotron radiation that stimulate the emission of fluorescence X-rays from atoms of elements whose K- or L-edges lie below the energy of the stimulating beam. Scanning and rotating the object through the beam allows for acquisition of a tomographic dataset that can be used to reconstruct images of the distribution of the elements in question. XFCT is a stimulated emission tomography modality, and it is thus necessary to correct for attenuation of the incident and fluorescence photons. The attenuation map is, however, generally known only at the stimulating beam energy and not at the energies of the various fluorescence X-rays of interest. We have developed a penalized-likelihood image reconstruction strategy for this problem. The approach alternates between updating the distribution of a given element and updating the attenuation map for that element's fluorescence X-rays. The approach is guaranteed to increase the penalized likelihood at each iteration. Because the joint objective function is not necessarily concave, the approach may drive the solution to a local maximum. To encourage the algorithm to seek out a reasonable local maximum, we include in the objective function a prior that encourages a relationship, based on physical considerations, between the fluorescence attenuation map and the distribution of the element being reconstructed.     

Incorporation of system resolution compensation (RC) in the ordered-subset transmission (OSTR) algorithm for transmission imaging in SPECT

In order to reconstruct attenuation maps with improved spatial resolution and quantitative accuracy, we developed an approximate method of incorporating system resolution compensation (RC) in the ordered-subset transmission (OSTR) algorithm for transmission reconstruction. Our method approximately models the blur caused by the finite intrinsic detector resolution, the nonideal source collimation and detector collimation. We derived the formulation using the optimization transfer principle as in the derivation of the OSTR algorithm. The formulation includes one forward-blur step and one back-blur step, which do not severely slow down reconstruction. The formulation could be applicable to various transmission geometries, such as point-source, line-source, and sheet-source systems. Through computer simulations of the MCAT phantom and transmission measurements of the air-filled Data Spectrum Deluxe single photo emission computed tomography (SPECT) Phantom on a system which employed a cone-beam geometry and a system which employed a scanning-line-source geometry, we showed that incorporation of RC increased spatial resolution and improved the quantitative accuracy of reconstruction. In simulation studies, attenuation maps reconstructed with RC correction improved the quantitative accuracy of emission reconstruction.     

A projector/backprojector with slice-to-slice blurring forefficient three-dimensional scatter modeling

Scatter correction is an important factor in single photon emission computed tomography (SPECT). Many scatter correction techniques, such as multiple-window subtraction and intrinsic modeling with iterative algorithms, have been under study for many years. Previously, we developed an efficient slice-to-slice blurring technique to model attenuation and system geometric response in a projector/backprojector pair, which was used in an ML-EM algorithm to reconstruct SPECT data. This paper proposes a projector/backprojector that models the three-dimensional (3-D) first-order scatter in SPECT, also using an efficient slice-to-slice blurring technique. The scatter response is estimated from a known nonuniform attenuation distribution map. It is assumed that the probability of detection of a first-order scattered photon from a photon that is emitted in a given source voxel and scattered in a given scatter voxel is proportional to the attenuation coefficient value at that voxel. Monte Carlo simulations of point sources and an MCAT torso phantom were used to verify the accuracy of the proposed projector/backprojector model. An experimental Jaszczak torso/cardiac phantom SPECT study was also performed. For a 64 x 64 x 64 image volume, it took 8.7 s to perform each iteration per slice on a Sun ULTRA Enterprise 3000 (167 MHz, 1 Gbyte RAM) computer, when modeling 3-D scatter, attenuation, and system geometric response functions. The main advantage of the proposed method is its easy implementation and the possibility of performing reconstruction in clinically acceptable time.     

Evaluation of an efficient compensation method for quantitative fan-beam brain SPECT reconstruction

Fan-beam collimators are designed to improve the system sensitivity and resolution for imaging small objects such as the human brain and breasts in single photon emission computed tomography (SPECT). Many reconstruction algorithms have been studied and applied to this geometry to deal with every kind of degradation factor. This paper presents a new reconstruction approach for SPECT with circular orbit, which demonstrated good performance in terms of both accuracy and efficiency. The new approach compensates for degradation factors including noise, scatter, attenuation, and spatially variant detector response. Its uniform attenuation approximation strategy avoids the additional transmission scan for the brain attenuation map, hence reducing the patient radiation dose and furthermore simplifying the imaging procedure. We evaluate and compare this new approach with the well-established ordered-subset expectation-maximization iterative method, using Monte Carlo simulations. We perform quantitative analysis with regional bias-variance, receiver operating characteristics, and Hotelling trace studies for both methods. The results demonstrate that our reconstruction strategy has comparable performance with a significant reduction of computing time.     

Simultaneous MR-Compatible Emission and Transmission Imaging for PET Using Time-of-Flight Information

Quantitative positron emission tomography (PET) imaging relies on accurate attenuation correction. Predicting attenuation values from magnetic resonance (MR) images is difficult because MR signals are related to proton density and relaxation properties of tissues. Here, we propose a method to derive the attenuation map from a transmission scan. An annulus transmission source is positioned inside the field-of-view of the PET scanner. First a blank scan is acquired. The patient is injected with FDG and placed inside the scanner. 511-keV photons coming from the patient and the transmission source are acquired simultaneously. Time-of-flight information is used to extract the coincident photons originating from the annulus. The blank and transmission data are compared in an iterative reconstruction method to derive the attenuation map. Simulations with a digital phantom were performed to validate the method. The reconstructed attenuation coefficients differ less than 5% in volumes of interest inside the lungs, bone, and soft tissue. When applying attenuation correction in the reconstruction of the emission data a standardized uptake value error smaller than 9% was obtained for all tissues. In conclusion, our method can reconstruct the attenuation map and the emission data from a simultaneous scan without prior knowledge about the anatomy or the attenuation coefficients of the tissues.     

Quantitative rotating multisegment slant-hole SPECT mammography with attenuation and collimator-detector response compensation

Rotating multisegment slant-hole (RMSSH) single photon emission computed tomography (SPECT) is suitable for detecting small and low-contrast breast lesions since it has much higher detection efficiency than conventional SPECT with a parallel-hole collimator and can image the breast at a closer distance. Our RMSSH SPECT reconstruction extends a previous rotation-shear transformation-based method to include nonuniform attenuation and collimator-detector response (CDR) compensation. To evaluate our reconstruction method, we performed two phantom simulation studies with 1) an isolated breast and 2) a breast phantom attached to the body torso. The reconstructed RMSSH SPECT images with attenuation and CDR compensation showed improved quantitative accuracy and less image artifacts than without. To evaluate the clinical efficacy of RMSSH SPECT mammography, we used a simulation study to compare with planar scintimammography in terms of the signal-to-noise ratio (SNR) value of a breast lesion. The RMSSH SPECT reconstruction images showed higher SNR value than the planar scintimammography images and even more so as we applied compensation for attenuation and collimator detector response. We conclude that attenuation and CDR compensation provide RMSSH SPECT mammography images with improved quality and quantitative accuracy.     

Dynamic imaging and tracer kinetic modeling for emission tomography using rotating detectors

When performing dynamic studies using emission tomography the tracer distribution changes during acquisition of a single set of projections. This is particularly true for some positron emission tomography (PET) systems which, like single photon emission computed tomography (SPECT), acquire data over a limited angle at any time, with full projections obtained by rotation of the detectors. In this paper, an approach is proposed for processing data from these systems, applicable to either PET or SPECT. A method of interpolation, based on overlapped parabolas, is used to obtain an estimate of the total counts in each pixel of the projections for each required frame-interval, which is the total time to acquire a single complete set of projections necessary for reconstruction. The resultant projections are reconstructed using traditional filtered backprojection (FBP) and tracer kinetic parameters are estimated using a method which relies on counts integrated over the frame-interval rather than instantaneous values. Simulated data were used to illustrate the technique's capabilities with noise levels typical of those encountered in either PET or SPECT. Dynamic datasets were constructed, based on kinetic parameters for fluoro-deoxy-glucose (FDG) and use of either a full ring detector or rotating detector acquisition. For the rotating detector, use of the interpolation scheme provided reconstructed dynamic images with reduced artefacts compared to unprocessed data or use of linear interpolation. Estimates for the metabolic rate of glucose had similar bias to those obtained from a full ring detector.     

A cone-beam filtered backprojection reconstruction algorithm for cardiac single photon emission computed tomography

A filtered backprojection reconstruction algorithm was developed for cardiac single photon emission computed tomography with cone-beam geometry. The algorithm reconstructs cone-beam projections collected from ;short scan' acquisitions of a detector traversing a noncircular planar orbit. Since the algorithm does not correct for photon attenuation, it is designed to reconstruct data collected over an angular range of slightly more than 180 degrees with the assumption that the range of angles is oriented so as not to acquire the highly attenuated posterior projections of emissions from cardiac radiopharmaceuticals. This sampling scheme is performed to minimize the attenuation artifacts that result from reconstructing posterior projections. From computer simulations, it is found that reconstruction of attenuated projections has a greater effect upon quantitation and image quality than any potential cone-beam reconstruction artifacts resulting from insufficient sampling of cone-beam projections. With nonattenuated projection data, cone beam reconstruction errors in the heart are shown to be small (errors of at most 2%).     

Reconstruction From Uniformly Attenuated SPECT Projection Data Using the DBH Method

An algorithm was developed for the 2-D reconstruction of truncated and nontruncated uniformly attenuated data acquired from single photon emission computed tomography (SPECT). The algorithm is able to reconstruct data from half-scan (180$^{circ}$ ) and short-scan (180$^{circ}$ +fan angle) acquisitions for parallel- and fan-beam geometries, respectively, as well as data from full-scan (360 $^{circ}$) acquisitions. The algorithm is a derivative, backprojection, and Hilbert transform (DBH) method, which involves the backprojection of differentiated projection data followed by an inversion of the finite weighted Hilbert transform. The kernel of the inverse weighted Hilbert transform is solved numerically using matrix inversion. Numerical simulations confirm that the DBH method provides accurate reconstructions from half-scan and short-scan data, even when there is truncation. However, as the attenuation increases, finer data sampling is required.      

Use of three-dimensional Gaussian interpolation in the projector/backprojector pair of iterative reconstruction for compensation of known rigid-body motion in SPECT

Due to the extended imaging times employed in single photon emission computed tomography (SPECT) and positron emission tomography (PET), patient motion during imaging is a common clinical occurrence. The fast and accurate correction of the three-dimensional (3-D) translational and rotational patient motion in iterative reconstruction is thus necessary to address this important cause of artifacts. We propose a method of incorporating 3-D Gaussian interpolation in the projector/backprojector pair to facilitate compensation for rigid-body motion in addition to attenuation and distance-dependent blurring. The method works as the interpolation step for moving the current emission voxel estimates and attenuation maps in the global coordinate system to the new patient location in the rotating coordinate system when calculating the expected projection. It also is employed for moving back the backprojection of the ratio of the measured projection to the expected projection and backprojection of the unit value (sensitivity factor) to the original location. MCAT simulations with known six-degree-of-freedom (6DOF) motion were employed to evaluate the accuracy of our method of motion compensation. We also tested the method with acquisitions of the data spectrum anthropomorphic phantom where motion during SPECT acquisition was measured using the Polaris IR motion tracking system. No motion artifacts were seen on the reconstructions with the motion compensation.     

Single-photon emission computed tomography

The subject of singl-photon emission computed tomography (SPECT) is generally reviewed. The basic interaction processes of gamma rays in matter are outlined, and the formation of conventional gamma-ray images is described. We then outline the extension of these concepts to the formation of three-dimensional or tomographic images. Of particular concern in emission tomography, the effects of gamma-ray attenuation and scattering are outlined. Several examples are given of practical SPECT systems, and representative results are given.     

Characterization of pinhole SPECT acquisition geometry

A method is presented to estimate the acquisition geometry of a pinhole single photon emission computed tomography (SPECT) camera with a circular detector orbit. This information is needed for the reconstruction of tomographic images. The calibration uses the point source projection locations of a tomographic acquisition of three point sources located at known distances from each other. It is shown that this simple phantom provides the necessary and sufficient information for the proposed calibration method. The knowledge of two of the distances between the point sources proves to be essential. The geometry is estimated by fitting analytically calculated projections to the measured ones, using a simple least squares Powell algorithm. Some mild a priori knowledge is used to constrain the solutions of the fit. Several of the geometrical parameters are however highly correlated. The effect of these correlations on the reconstructed images is evaluated in simulation studies and related to the estimation accuracy. The highly correlated detector tilt and electrical shift are shown to be the critical parameters for accurate image reconstruction. The performance of the algorithm is finally demonstrated by phantom measurements. The method is based on a single SPECT scan of a simple calibration phantom, executed immediately after the actual SPECT acquisition. The method is also applicable to cone-beam SPECT and X-ray CT.     

Attenuation Correction for SPECT: An Evaluation of Hybrid Approaches

Most methods that have been proposed for attenuation compensation in single-photon emission computed tomography (SPECT) either rely on simplifying assumptions, or use slow iteration to achieve accuracy. Recently, hybrid methods which combine iteration with simple multiplicative correction have been proposed by Chang and by Moore et al. In this study we evaluated these methods using both simulated and real phantom data from a rotating gamma camera. Of special concern were the effects of assuming constant attenuation distributions for correction and of using only 180 degrees of projection data in the reconstructions. Results were compared by means of image contrast, %RMS error, and a chi-square error statistic. Simulations showed the hybrid methods to converge after 1-2 iterations when 360 degrees data were used, less rapidly for 180 degrees data. The Moore method was more accurate than our modified Chang method for 180 degrees data. Phantom data indicated the importance of using an accurate attenuation map for both methods. The speed of convergence of the hybrid algorithms compared to traditional iterative techniques, and their accuracy in reconstructing photon activity, even with 180 degrees data, makes them attractive for use in quantitative analysis of SPECT reconstructions.     

Application of reconstruction-based scatter compensation tothallium-201 SPECT: implementations for reduced reconstructed imagenoise

Scatter compensation in Tl-201 single photon emission computed tomography (SPECT) presents an interesting challenge because of the multiple emission energies and relatively large proportion of scattered photons. In this paper, the authors present a simulation study investigating reconstructed image noise levels arising from various implementations of iterative reconstruction-based scatter compensation (RBSC) in Tl-201 SPECT. A two-stage analysis was used to study single and multiple energy window implementations of reconstruction-based scatter compensation, and RBSC was compared to the upper limits on performance for other approaches to handling scatter. In the first stage, singular value decomposition of the system transfer matrix was used to analyze noise levels in a manner independent of the choice of reconstruction algorithm, providing results valid across a wide range of regularizations. In the second stage, the data were reconstructed using maximum-likelihood expectation-maximization, and the noise properties of the resultant images were analyzed. The best RBSC performance was obtained using multiple energy windows, one for each emission photopeak, and RBSC outperformed the upper limit on subtraction-based compensation methods. Implementing RBSC with the correct choice of energy window acquisition scheme is promising method for performing scatter compensation for Tl-201 SPECT     

Resampling estimates of precision in emission tomography

Methods for estimating the regional variance in emission tomography images which arise from the Poisson nature of the raw data are discussed. The methods are based on the bootstrap and jackknife methods of statistical resampling theory. The bootstrap is implemented in time-of-flight PET (positron emission tomography); the same techniques can be applied to non-time-of-flight PET and SPECT (single-photon-emission computed tomography). The estimates are validated by comparing them to those obtained by repetition of emission scans, using data from a time-of-flight positron emission tomograph. Simple expressions for the accuracy of the estimates are given. The present approach is computationally feasible and can be applied to any reconstruction technique as long as the data are acquired in a raw, uncorrected form.