Improvement of brain single photon emission tomography (SPET) using transmission data acquisition in a four-head SPET scanner

Evaluation of the number of SPECT projections in the ordered subsets-expectation maximization image reconstruction method

Filtered back projection (FBP) method, maximum likelihood-expectation maximization(ML-EM) method, and ordered subsets-expectation maximization (OS-EM) method are currently used for reconstruction of SPECT images in clinical studies. In the ML-EM method, images of good quality can be reconstructed even with a small sampling number of projection data, when compared with FBP. Shorter acquisition time and less radionuclide dose are preferable in the clinical setting if image quality is the same. In this study, we attempted to find optimal conditions for reconstruction of OS-EM images with commonly used sampling numbers of 30, 60 and 120 (step angles: 12 degrees, 6 degrees, and 3 degrees, respectively), with acquisition counts/projection of 30, 60, 120 and 240 each. We adjusted the pixel counts of reconstructed images to be constant, by setting combination of sampling number and counts/projection (120 sampling number for 30 counts/projection, 60 for 60, and 30 for 120). Among the 3 acquisition conditions, the small sampling number of 30 had large acquisition counts per direction, resulting in low signal to noise ratio. Under this condition, the resolution was slightly low, but the uniformity of images was high. The combination of OS-EM and smaller sampling projection number may be clinically useful with reduction of the examination time, which is also beneficial to reduce dead time for gamma-camera rotation.

A study on attenuation correction using Tc-99m external TCT source in Tc-99m GSA liver SPECT

PURPOSE

In attenuation correction of ECT images by transmission CT (TCT) with an external 99mTc gamma-ray source, simultaneous TCT/ECT data acquisition is difficult, when the same radionuclide such as 99mTc-tetrofosmin or 99mTc-GSA is used as the tracer. In this case, TCT is usually acquired before administration of the tracer, and ECT is acquired separately after the tracer injection. However, misregistration may occur between the TCT and ECT images, and the repetition of examinations add to the mental and physical stress of the patients. In this study, to eliminate this problem, we evaluated whether attenuation correction of ECT images can be achieved by acquiring TCT and ECT simultaneously, then acquiring ECT alone, and preparing an attenuation map by subtracting the latter from the former using 99mTc-GSA liver ECT.

METHOD

The ECT system used was a three-head gamma camera equipped with one cardiac fan beam collimator and two parallel beam collimators. External gamma-ray source for TCT of 99mTc was 740 MBq, and ECT of 99mTc-GSA was 185 MBq. First, pure TCT data were acquired for the original TCT-map, then, ECT/TCT data were acquired for the subtracted TCT-map, and finally, pure ECT data were acquired. The subtracted attenuation map was produced by subtracting the pure ECT image from the TCT/ECT image, and attenuation correction of the ECT image was done using both this subtracted TCT map and attenuation map from pure TCT. These two attenuation corrected images and non-corrected images were compared. Hot rods phantom, a liver phantom with a defect, and 10 patients were evaluated.

RESULTS

Attenuation corrected ECT values using the subtraction attenuation map showed an error of about 5% underestimation compared with ECT values of the images corrected by original attenuation map at the defect in the liver phantom. A good correlation of y = 22.65 + 1.06x, r = 0.958 was observed also in clinical evaluation.

CONCLUSION

By means of the method proposed in this study, it is possible to perform simultaneous TCT/ECT data acquisition for attenuation correction using Tc-99m external source in Tc-99m GSA liver SPECT. Moreover, it is thought that this method decreases the mental and physical stress of the patients.

Truncation correction of fan beam transmission data for attenuation correction using parallel beam emission data on a 3-detector SPECT system

BACKGROUND

When the simultaneous transmission computed tomography (TCT)/single photon emission CT (SPECT) acquisition protocol is applied to myocardial studies using a 3-detector SPECT, the narrow effective field of view of a fan beam collimator used for TCT acquisition may cause truncation artifacts on TCT images. In this paper, we propose a new method of correcting for the truncation of TCT.

METHODS

The truncated parts of the TCT projection data are corrected using quadratic functions, based on the properties that the integral of non-truncated TCT projection data is constant at any projection angle and the position of the centre of gravity is focused on a fixed point. The usefulness of our method was investigated in phantom and human studies using a 3-detector SPECT equipped with one cardiac fan beam collimator for TCT and two parallel beam collimators for SPECT. We used Tl as a tracer for SPECT and Tc as an external source for TCT.

RESULTS

The phantom and human studies showed that our method can adequately correct for the truncation of TCT data acquired using a fan beam collimator in a 3-detector SPECT, as long as there is no truncation in SPECT data.

CONCLUSION

Our method appears to be useful for improving the SPECT images obtained using simultaneous TCT/SPECT acquisition in a 3-detector SPECT. However, further studies will be necessary to establish the clinical usefulness of this method.

Segmented attenuation correction for myocardial SPECT

PURPOSE

One of the main factors contributing to the accuracy of attenuation correction for SPECT imaging using transmission computed tomography (TCT) with an external gamma-ray source is the radionuclide count. To reduce deterioration of TCT images due to inadequate radionuclide counts, a correction method, segmented attenuation correction (SAC), in which TCT data are transformed into several components (segments) such as water, lungs and spine, providing a satisfactory attenuation correction map with less counts, has been developed. The purpose of this study was to examine the usefulness of SAC for myocardial SPECT with attenuation correction.

METHODS

A myocardial phantom filled with Tc-99m was scanned with a triple headed SPECT system, equipped with one cardiac fan beam collimator for TCT and two parallel hole collimators for ECT. As an external gamma-ray source for TCT, 740 MBq of Tc-99m was also used. Since Tc-99m was also used for ECT, the TCT and ECT data were acquired separately. To make radionuclide counts, the TCT data were acquired in the sequential repetition mode, in which a 3-min-rotation was repeated 7 times followed by a 10-min-rotation 4 times (a total of 61 minutes). The TCT data were reconstructed by adding some of these rotations to make TCT maps with various radionuclide counts. Three types of SAC were used: (a) 1-segment SAC in which the body structure was regarded as water, (b) 2-segment SAC, in which the body structure was regarded as water and lungs, and (c) 3-segment SAC, in which the body structure was regarded as water, lungs and spine. We compared corrected images obtained with non-segmentation methods, and with 1- to 3-segment SACs. We also investigated the influence of radionuclide counts of TCT (3, 6, 9, 12, 15, 18, 21, 31, 41, 51, 61 min acquisition) on the accuracy of the attenuation correction.

RESULTS

Either 1-segment or 2-segment SAC was sufficient to correct the attenuation. When non-segmentation TCT attenuation methods were used, rotations of at least 31 minutes were required to obtain sufficiently large counts for TCT. When the 3-segment SAC was used, the minimal acquisition time for a satisfactory TCT map was 7 min.

CONCLUSION

The 3-segment SAC was effective for attenuation correction, requiring fewer counts (about 1/5 of the value for non-segmentation TCT), or less radiation for TCT.

Attenuation correction of myocardial SPECT images with X-ray CT: effects of registration errors between X-ray CT and SPECT

PURPOSE

Attenuation correction with an X-ray CT image is a new method to correct attenuation on SPECT imaging, but the effect of the registration errors between CT and SPECT images is unclear. In this study, we investigated the effects of the registration errors on myocardial SPECT, analyzing data from a phantom and a human volunteer.

METHODS

Registerion (fusion) of the X-ray CT and SPECT images was done with standard packaged software in three dimensional fashion, by using linked transaxial, coronal and sagittal images. In the phantom study, an X-ray CT image was shifted 1 to 3 pixels on the x, y and z axes, and rotated 6 degrees clockwise. Attenuation correction maps generated from each misaligned X-ray CT image were used to reconstruct misaligned SPECT images of the phantom filled with 201Tl. In a human volunteer, X-ray CT was acquired in different conditions (during inspiration vs. expiration). CT values were transferred to an attenuation constant by using straight lines; an attenuation constant of 0/cm in the air (CT value = -1,000 HU) and that of 0.150/cm in water (CT value = 0 HU). For comparison, attenuation correction with transmission CT (TCT) data and an external gamma-ray source (99mTc) was also applied to reconstruct SPECT images.

RESULTS

Simulated breast attenuation with a breast attachment, and inferior wall attenuation were properly corrected by means of the attenuation correction map generated from X-ray CT. As pixel shift increased, deviation of the SPECT images increased in misaligned images in the phantom study. In the human study, SPECT images were affected by the scan conditions of the X-ray CT.

CONCLUSION

Attenuation correction of myocardial SPECT with an X-ray CT image is a simple and potentially beneficial method for clinical use, but accurate registration of the X-ray CT to SPECT image is essential for satisfactory attenuation correction.

Performance evaluation of OSEM reconstruction algorithm incorporating three-dimensional distance-dependent resolution compensation for brain SPECT: a simulation study

Iterative reconstruction techniques such as an ordered subsets-expectation maximization (OSEM) algorithm can easily incorporated various physical models of attenuation or scatter. We implemented OSEM reconstruction algorithm incorporating compensation for distance-dependent blurring due to the collimator in SPECT. The algorithm was examined by computer simulation to estimate the accuracy for brain perfusion study.

METHODS

The detector response was assumed to be a two-dimensional Gauss function and the width of the function varied linearly with the source-to-detector distance. The attenuation compensation (AC) was also included. To investigate the properties of the algorithm, we performed computer simulations with the point source and digital brain phantoms. In the point source phantom, the uniformity of FWHM for the radial, tangential and longitudinal directions was evaluated on the reconstruction image. As for the brain phantom, quantitative accuracy was estimated by comparing the reconstructed images with the true image by the mean square error (MSE) and the ratio of gray and white matter counts (G/W). Both noise free and noisy simulations were examined.

RESULTS

In the point source simulation, FWHM in radial, tangential and longitudinal directions were 14.7, 14.7 and 15.0 mm at the image center and were 15.9, 9.83 and 10.6 mm at a distance of 15 cm from the center by using FBP, respectively. On the other hand, they were 8.12, 8.12 and 7.83 mm at the image center, and were 7.45, 7.44 and 7.01 mm at 15 cm from the center by OSEM with distance-dependent resolution compensation (DRC). An isotropic and stationary resolution was obtained at any location by OSEM with DRC. The spatial resolution was also improved about 6.5 mm by OSEM with DRC at the image center. In the brain phantom simulation, the blurring at the edge of the brain structure was eliminated by using OSEM with both DRC and AC. The G/W was 2.95 and 2.68 for noise free and noisy cases, respectively, when no compensation was performed. But the values for G/W without and with noise became 3.45 and 3.21 with AC only and were improved to 3.75 and 3.71 with both AC and DRC. The G/W approached the true value (4.00) by using OSEM with both AC and DRC even when there was statistical noise.

CONCLUSION

In conclusion, OSEM reconstruction including the distance-dependent resolution compensation algorithm was reasonably successful in achieving isotropic and stationary resolution and improving the quantitative accuracy for brain perfusion SPECT.

Receiver operating characteristic (ROC) analysis of images reconstructed with iterative expectation maximization algorithms

PURPOSE

The quality of images reconstructed by means of the maximum likelihood-expectation maximization (ML-EM) and ordered subset (OS)-EM algorithms, was examined with parameters such as the number of iterations and subsets, then compared with the quality of images reconstructed by the filtered back projection method.

METHODS

Phantoms showing signals inside signals, which mimicked single-photon emission computed tomography (SPECT) images of cerebral blood flow and myocardial perfusion, and phantoms showing signals around the signals obtained by SPECT of bone and tumor were used for experiments. To determine signals for recognition, SPECT images in which the signals could be appropriately recognized with a combination of fewer iterations and subsets of different sizes and densities were evaluated by receiver operating characteristic (ROC) analysis. The results of ROC analysis were applied to myocardial phantom experiments and scintigraphy of myocardial perfusion.

RESULTS

Taking the image processing time into consideration, good SPECT images were obtained by OS-EM at iteration No. 10 and subset 5. CONCLULSION: This study will be helpful for selection of parameters such as the number of iterations and subsets when using the ML-EM or OS-EM algorithms.

Attenuation correction in quantitative SPECT of cerebral blood flow: a Monte Carlo study

Monte Carlo simulation has been used to produce projections from a voxel-based brain phantom, simulating a 99mTc-HMPAO single photon emission computed tomography (SPECT) brain investigation. For comparison, projections free from the effects of attenuation and scattering were also simulated, giving ideal transaxial images after reconstruction. Three methods of attenuation correction were studied: (a) a pre-processing method, (b) a post-processing uniform method and (c) a post-processing non-uniform method using a density map. The accuracy of these methods was estimated by comparison of the reconstructed images with the ideal images using the normalized mean square error, NMSE, and quantitative values of the regional cerebral blood flow, rCBF. A minimum NMSE was achieved for the effective linear attenuation coefficient mu(eff) = 0.07 (0.09) cm(-1) for the uniform(pre) method, the effective mass attenuation coefficient mu(eff)/rho = 0.08 (0.10) cm2 g(-1) for the uniform(post) method and mu(eff)/rho = 0.12 (0.13) cm2 g(-1) for the non-uniform(post) method. Values in parentheses represent the case of dual-window scatter correction. The non-uniform(post) method performed better, as measured by the NMSE, both with and without scatter correction. Furthermore, the non-uniform(post) method gave, on average, more accurate rCBF values. Although the difference in rCBF accuracy was small between the various methods, the same method should be used for patient studies as for the reference material.

Experimental and numerical investigation of the 3D SPECT photon detection kernel for non-uniform attenuating media

Experimental tests for non-uniform attenuating media are performed to validate theoretical expressions for the photon detection kernel, obtained from a recently proposed analytical theory of photon propagation and detection for SPECT. The theoretical multi-dimensional integral expressions for the photon detection kernel, which are computed numerically, describe the probability that a photon emitted from a given source voxel will trigger detection of a photon at a particular projection pixel. The experiments were performed using a cylindrical water-filled phantom with large cylindrical air-filled inserts to simulate inhomogeneity of the medium. A point-like, a short thin cylindrical and a large cylindrical radiation source of 99Tcm were placed at various positions within the phantom. The values numerically calculated from the theoretical kernel expression are in very good agreement with the experimentally measured data. The significance of Compton-scattered photons in planar image formation is discussed and highlighted by these results. Using both experimental measurements and the calculated values obtained from the theory, the kernel's size is investigated. This is done by determining the square N x N pixel neighbourhood of the gamma camera that must be connected to a particular radiation source voxel to account for a specific fraction of all counts recorded at all camera pixels. It is shown that the kernel's size is primarily dependent upon the source position and the properties of the attenuating medium through Compton scattering events, with 3D depth-dependent collimator resolution playing an important but secondary role, at least for imaging situations involving parallel hole collimation. By considering small point-like sources within a non-uniform elliptical phantom, approximating the human thorax, it is demonstrated that on average a 12 cm x 12 cm area of the camera plane is required to collect 85% of the total count recorded. This is a significantly larger connectivity than the 3 cm x 3 cm area required if scattering contributions are ignored and only the 3D depth-dependent collimator resolution is considered.

Simultaneous emission and transmission measurements as an adjunct to dynamic planar gamma camera studies

Anatomical imaging provides useful information which complements functional imaging performed using a gamma camera. We have previously used transmission measurements in single-photon emission tomography acquired simultaneously with the emission scan using either a plane flood source or a moving line source for attenuation and scatter correction. This approach is equally applicable in planar imaging and provides useful information to assist in detecting patient motion and in defining regions of interest in dynamic studies. We have adapted a moving transmission line source to acquire dynamic geometric mean measurements in the study of the mucociliary clearance of inhaled technetium-99m labelled colloids with a single-headed rotating gamma camera. The line source makes a return pass for each emission acquisition frame (alternating anterior/posterior views), each pass being initiated by a signal from the gamma camera. The result is a dynamic sequence of emission and transmission measurements obtained from a single acquisition. In this application transmission measurements are used to define the lung outline for clearance determination and to check for subject movement throughout the duration of the study.

Segmentation of the body and lungs from Compton scatter and photopeak window data in SPECT: a Monte-Carlo investigation

In SPECT imaging of the chest, nonuniform attenuation correction requires use of a patient specific attenuation (mu) map. Such a map can be obtained by estimating the regions of (1) the lungs and (2) the soft tissues and bones, and then assigning an appropriate value of attenuation coefficient (mu) to each region. The authors proposed a method to segment such regions from the Compton scatter and photopeak window SPECT slices of Tc-99m Sestamibi studies. The Compton scatter slices are used to segment the body outline and to estimate the regions of the lungs. Locations of the back bone and sternum are estimated from the photopeak window slices to assist in the segmentation. To investigate the accuracy of using Compton scatter slices in estimating the regions of the body and the lungs, a Monte-Carlo SPECT simulation of an anthropomorphic phantom with an activity distribution and noise characteristics similar to patient data was conducted. Energy windows of various widths were simulated for use in locating a suitable Compton scatter window for imaging, The effects of attenuation correction using a mu map based on segmentation were also studied. The results demonstrated for the activity and mu maps studied herein that: (1) reasonable contrast could be obtained from Compton scatter data for the segmentation of the lung regions, (2) true positive rates of 99% and 89% for determining the body and lung regions, respectively, with total error rates of 4% and 29%, could be achieved, (3) usage of a mu map based on segmentation for attenuation correction improved relative quantification over filtered backprojection, (4) variations in the assigned mu value of 40% smaller or 40% larger in the lung regions had an insignificant impact on the results of relative quantification, (5) a wide energy window away from the photopeak window for recording scattered events could benefit both the segmentation of the lung regions and the attenuation correction of the activity in the myocardium region, and (6) usage of a smaller than true mu value in the lung regions of an assigned mu map might benefit attenuation correction for absolute quantification.

An evaluation of the accelerated expectation maximization algorithms for single-photon emission tomography image reconstruction

We previously reported that brain single-photon emission tomography (SPET) images could be improved by using an attenuation coefficient map constructed with transmission data and the iterative expectation maximization (EM) algorithm. However, the conventional EM algorithm (CEM) typically requires 30-80 iterations to provide acceptable results, limiting its clinical applicability. Several methods have been proposed to accelerate the EM algorithm. The purpose of this study was to search for a practical method for accelerating the EM algorithm. The methods investigated here include the accelerated EM algorithm (ACEM) using additive correction, ACEM using multiplicative correction, and Tanaka's filtered iterative reconstruction method (FIR). These methods were assessed by simulated SPET studies of a phantom incorporating nonuniform attenuation and by reference to clinical brain SPET data. In the simulation studies, the above methods were evaluated by using three parameters (root mean square error, log likelihood value, and contrast recovery coefficient); the results showed that FIR had an advantage over other methods in terms of all parameters. The results obtained using the clinical data demonstrated that FIR could reconstruct acceptable images in only five iterations. These results show that FIR offers significant advantages over CEM or other ACEMs, indicating that FIR can make the EM algorithm practical for clinical use in SPET.