3D RBI-EM reconstruction with spherically-symmetric basis function for SPECT rotating slat collimator

A single photon emission computed tomography (SPECT) rotating slat collimator with strip detector acquires distance-weighted plane integral data, along with the attenuation factor and distance-dependent detector response. In order to image a 3D object, the slat collimator device has first to spin around its axis and then rotate around the object to produce 3D projection measurements. Compared to the slice-by-slice 2D reconstruction for the parallel-hole collimator and line integral data, a more complex 3D reconstruction is needed for the slat collimator and plane integral data. In this paper, we propose a 3D RBI-EM reconstruction algorithm with spherically-symmetric basis function, also called 'blobs', for the slat collimator. It has a closed and spherically symmetric analytical expression for the 3D Radon transform, which makes it easier to compute the plane integral than the voxel. It is completely localized in the spatial domain and nearly band-limited in the frequency domain. Its size and shape can be controlled by several parameters to have desired reconstructed image quality. A mathematical lesion phantom study has demonstrated that the blob reconstruction can achieve better contrast-noise trade-offs than the voxel reconstruction without greatly degrading the image resolution. A real lesion phantom study further confirmed this and showed that a slat collimator with CZT detector has better image quality than the conventional parallel-hole collimator with NaI detector. The improvement might be due to both the slat collimation and the better energy resolution of the CZT detector.     

The cardiofocal collimator: a variable-focus collimator for cardiac SPECT

Investigators in nuclear medicine have long been in search of a practical method to increase the number of detected events in cardiac SPECT. A clinically practical method requires a simple data acquisition protocol, clinically acceptable reconstruction times, artifact levels near or below visual threshold, and the use of currently available cameras and computers. Towards this end, we have developed the Cardiofocal collimator, a variable-focus collimator for cardiac SPECT that increases the number of detected events from the heart by more than a factor of two compared to that of a parallel-hole collimator with equivalent resolution. In both the transverse and axial dimensions, the focusing is strongest at the centre of the collimator, and gradually relaxes to nearly parallel-hole collimation at the edge of the collimator. The variable-focus concept provides an increase in the number of counts from organs imaged near the centre of the collimator, where the heart will spend most of the time during a cardiac SPECT study, while adequately sampling enough of the background activity distribution to prevent truncation artifacts in the reconstructed images. Images are reconstructed in clinically acceptable times using a filtered backprojection reconstruction algorithm. The algorithm supports both full-scan (360 degrees) and short-scan (180 degrees plus the fan angle) acquisitions. The results of simulations and phantom studies are included to demonstrate the performance of the Cardiofocal collimator.     

Design and performance of a multi-pinhole collimation device for small animal imaging with clinical SPECT and SPECT-CT scanners

A multi-pinhole collimation device is developed that uses the gamma camera detectors of a clinical SPECT or SPECT-CT scanner to produce high-resolution SPECT images. The device consists of a rotating cylindrical collimator having 22 tungsten pinholes with 0.9 mm diameter apertures and an animal bed inside the collimator that moves linearly to provide helical or ordered-subsets axial sampling. CT images also may be acquired on a SPECT-CT scanner for purposes of image co-registration and SPECT attenuation correction. The device is placed on the patient table of the scanner without attaching to the detectors or scanner gantry. The system geometry is calibrated in-place from point source data and is then used during image reconstruction. The SPECT imaging performance of the device is evaluated with test phantom scans. Spatial resolution from reconstructed point source images is measured to be 0.6 mm full width at half maximum or better. Micro-Derenzo phantom images demonstrate the ability to resolve 0.7 mm diameter rod patterns. The axial slabs of a Micro-Defrise phantom are visualized well. Collimator efficiency exceeds 0.05% at the center of the field of view, and images of a uniform phantom show acceptable uniformity and minimal artifact. The overall simplicity and relatively good imaging performance of the device make it an interesting low-cost alternative to dedicated small animal scanners.     

Resolution recovery for list-mode reconstruction in SPECT

The purpose of the study was to evaluate the resolution recovery in the list-mode iterative reconstruction algorithm (LMIRA) for SPECT. In this study we compare the performance of the proposed method with other iterative resolution recovery methods for different noise levels. We developed an iterative reconstruction method which uses list-mode data instead of binned data. The new algorithm makes use of a more accurate model of the collimator structure. We compared the SPECT list-mode reconstruction with MLEM, OSEM and RBI, all including resolution recovery. For the evaluation we used Gaussian shaped sources with different FWHM at three different locations and three noise levels. For these distributions we calculated the reconstructed images for a different number of iterations. The absolute error for the reconstructed images was used to evaluate the performance. The performance of all four methods is comparable for the sources located in the centre of the field of view. For the sources located out of the centre, the error of the list-mode method is significantly lower than that of the other methods. Splitting the system model into a separate object-dependent and detector-dependent module gives us a flexible reconstruction method. With this we can very easily adapt the resolution recovery to different collimator types.     

An analytical reconstruction algorithm for multifocal converging-beam SPECT

Multifocal converging-beam collimation has been suggested for cardiac SPECT imaging to increase sensitivity over the heart without truncation of the activity distribution in the chest. In this study, an analytical reconstruction algorithm is derived for multifocal fan-beam and multifocal cone-beam tomography. In the algorithm, the projection data are differently weighted and filtered, depending on the distance from the detector. For a given image point, the set of filtered data corresponding to the distance between this point and detector is backprojected to determine the pixel value of the point. Thus, the backprojection is done only once at each projection view. To evaluate the algorithm, simulation studies are performed using a 3D Defrise slab phantom without considering photon attenuation and scatter, detector response, and statistical noise. Reconstructed images demonstrate that reasonable quality can be achieved with a modest focal-length variation rate and with a small radius of rotation. A collimator with a focal length increasing quickly near its centre provides better quality in the image region distant from the central plane of the cone geometry, but produces more severe artifacts at the centre of the reconstructed image, compared to a collimator with an initially slowly varying focal length.     

Brain SPECT with short focal-length cone-beam collimation

Single-photon emission-computed tomography (SPECT) imaging of deep brain structures is compromised by loss of photons due to attenuation. We have previously shown that a centrally peaked collimator sensitivity function can compensate for this phenomenon, increasing sensitivity over most of the brain. For dual-head instruments, parallel-hole collimators cannot provide variable sensitivity without simultaneously degrading spatial resolution near the center of the brain; this suggests the use of converging collimators. We have designed collimator pairs for dual-head SPECT systems to increase sensitivity, particularly in the center of the brain, and compared the new collimation approach to existing approaches on the basis of performance in estimating activity concentration of small structures at various locations in the brain. The collimator pairs we evaluated included a cone-beam collimator, for increased sensitivity, and a fan-beam collimator, for data sufficiency. We calculated projections of an ellipsoidal uniform background, with 0.9-cm-radius spherical lesions at several locations in the background. From these, we determined ideal signal-to-noise ratios (SNRCRB) for estimation of activity concentration within the spheres, based on the Cramer-Rao lower bound on variance. We also reconstructed, by an ordered-subset expectation-maximization (OS-EM) procedure, images of this phantom, as well as of the Zubal brain phantom, to allow visual assessment and to ensure that they were free of artifacts. The best of the collimator pairs evaluated comprised a cone-beam collimator with 20 cm focal length, for which the focal point is inside the brain, and a fan-beam collimator with 40 cm focal length. This pair yielded increased SNRCRB, compared to the parallel-parallel pair, throughout the imaging volume. The factor by which SNRCRB increased ranged from 1.1 at the most axially extreme location to 3.5 at the center. The gains in SNRCRB were relatively robust to mismatches between the center of the brain and the center of the imaging volume. Artifact-free reconstructions of simulated data acquired using this pair were obtained. Combining fan-beam and short-focusing cone-beam collimation should greatly improve dual-head brain SPECT imaging, especially for centrally located structures.     

Comparing planar image quality of rotating slat and parallel hole collimation: influence of system modeling

The main remaining challenge for a gamma camera is to overcome the existing trade-off between collimator spatial resolution and system sensitivity. This problem, strongly limiting the performance of parallel hole collimated gamma cameras, can be overcome by applying new collimator designs such as rotating slat (RS) collimators which have a much higher photon collection efficiency. The drawback of a RS collimated gamma camera is that, even for obtaining planar images, image reconstruction is needed, resulting in noise accumulation. However, nowadays iterative reconstruction techniques with accurate system modeling can provide better image quality. Because the impact of this modeling on image quality differs from one system to another, an objective assessment of the image quality obtained with a RS collimator is needed in comparison to classical projection images obtained using a parallel hole (PH) collimator. In this paper, a comparative study of image quality, achieved with system modeling, is presented. RS data are reconstructed to planar images using maximum likelihood expectation maximization (MLEM) with an accurate Monte Carlo derived system matrix while PH projections are deconvolved using a Monte Carlo derived point-spread function. Contrast-to-noise characteristics are used to show image quality for cold and hot spots of varying size. Influence of the object size and contrast is investigated using the optimal contrast-to-noise ratio (CNR(o)). For a typical phantom setup, results show that cold spot imaging is slightly better for a PH collimator. For hot spot imaging, the CNR(o) of the RS images is found to increase with increasing lesion diameter and lesion contrast while it decreases when background dimensions become larger. Only for very large background dimensions in combination with low contrast lesions, the use of a PH collimator could be beneficial for hot spot imaging. In all other cases, the RS collimator scores better. Finally, the simulation of a planar bone scan on a RS collimator revealed a hot spot contrast improvement up to 54% compared to a classical PH bone scan.     

Monte Carlo modelling of the performance of a rotating slit-collimator for improved planar gamma-camera imaging.

Planar imaging with a gamma camera is currently limited by the performance of the collimator. Spatial resolution and sensitivity trade off against each other; it is not possible with conventional parallel-hole collimation to have high geometric sensitivity and at the same time excellent spatial resolution unless field-of-view is sacrificed by using fan- or cone-beam collimators. We propose a rotating slit-collimator which collects one-dimensional projections from which the planar image may be reconstructed by the theory of computed tomography. The performance of such a collimator is modelled by Monte Carlo methods and images are reconstructed by a convolution and backprojection technique. The performance is compared with that of a conventional parallel-hole collimator and it is shown that higher spatial resolution with increased sensitivity is possible with the slit-collimator. For a point source a spatial resolution of some 6 mm at a distance of 100 mm from the collimator with a x7 sensitivity compared with a parallel-hole collimator was achieved. Applications to bone scintigraphy are modelled and an improved performance in hot-spot imaging is demonstrated. The expected performance in cold-spot imaging is analytically investigated. The slit-collimator is not expected to improve cold-spot imaging. Practical design considerations are discussed.     

Monte Carlo simulations of clinical PET and SPECT scans: impact of the input data on the simulated images

Monte Carlo simulations of emission tomography have proven useful to assist detector design and optimize acquisition and processing protocols. The more realistic the simulations, the more straightforward the extrapolation of conclusions to clinical situations. In emission tomography, accurate numerical models of tomographs have been described and well validated under specific operating conditions (collimator, radionuclide, acquisition parameters, count rates, etc). When using these models under these operating conditions, the realism of simulations mostly depends on the activity distribution used as an input for the simulations. It has been proposed to derive the input activity distribution directly from reconstructed clinical images, so as to properly model the heterogeneity of the activity distribution between and within organs. However, reconstructed patient images include noise and have limited spatial resolution. In this study, we analyse the properties of the simulated images as a function of the properties of the reconstructed images used to define the input activity distributions in (18)F-FDG PET and (131)I SPECT simulations. The propagation through the simulation/reconstruction process of the noise and spatial resolution in the input activity distribution was studied using simulations. We found that the noise properties of the images reconstructed from the simulated data were almost independent of the noise in the input activity distribution. The spatial resolution in the images reconstructed from the simulations was slightly poorer than that in the input activity distribution. However, using high-noise but high-resolution patient images as an input activity distribution yielded reconstructed images that could not be distinguished from clinical images. These findings were confirmed by simulated highly realistic (131)I SPECT and (18)F-FDG PET images from patient data. In conclusion, we demonstrated that (131)I SPECT and (18)F-FDG PET images indistinguishable from real scans can be simulated using activity maps with spatial resolution higher than that used in routine clinical applications.     

MR-based keyhole SPECT for small animal imaging

The rationale for multi-modality imaging is to integrate the strengths of different imaging technologies while reducing the shortcomings of an individual modality. The work presented here proposes a limited-field-of-view (LFOV) SPECT reconstruction technique that can be implemented on a multi-modality MR/SPECT system that can be used to obtain simultaneous MRI and SPECT images for small animal imaging. The reason for using a combined MR/SPECT system in this work is to eliminate any possible misregistration between the two sets of images when MR images are used as a priori information for SPECT. In nuclear imaging the target area is usually smaller than the entire object; thus, focusing the detector on the LFOV results in various advantages including the use of a smaller nuclear detector (less cost), smaller reconstruction region (faster reconstruction) and higher spatial resolution when used in conjunction with pinhole collimators with magnification. The MR/SPECT system can be used to choose a region of interest (ROI) for SPECT. A priori information obtained by the full field-of-view (FOV) MRI combined with the preliminary SPECT image can be used to reduce the dimensions of the SPECT reconstruction by limiting the computation to the smaller FOV while reducing artifacts resulting from the truncated data. Since the technique is based on SPECT imaging within the LFOV it will be called the keyhole SPECT (K-SPECT) method. At first MRI images of the entire object using a larger FOV are obtained to determine the location of the ROI covering the target organ. Once the ROI is determined, the animal is moved inside the radiofrequency (rf) coil to bring the target area inside the LFOV and then simultaneous MRI and SPECT are performed. The spatial resolution of the SPECT image is improved by employing a pinhole collimator with magnification >1 by having carefully calculated acceptance angles for each pinhole to avoid multiplexing. In our design all the pinholes are focused to the center of the LFOV. K-SPECT reconstruction is accomplished by generating an adaptive weighting matrix using a priori information obtained by simultaneously acquired MR images and the radioactivity distribution obtained from the ROI region of the SPECT image that is reconstructed without any a priori input. Preliminary results using simulations with numerical phantoms show that the image resolution of the SPECT image within the LFOV is improved while minimizing artifacts arising from parts of the object outside the LFOV due to the chosen magnification and the new reconstruction technique. The root-mean-square-error (RMSE) in the out-of-field artifacts was reduced by 60% for spherical phantoms using the K-SPECT reconstruction technique and by 48.5-52.6% for the heart in the case with the MOBY phantom. The K-SPECT reconstruction technique significantly improved the spatial resolution and quantification while reducing artifacts from the contributions outside the LFOV as well as reducing the dimension of the reconstruction matrix.     

Investigation of dynamic SPECT measurements of the arterial input function in human subjects using simulation, phantom and human studies

Computer simulations, a phantom study and a human study were performed to determine whether a slowly rotating single-photon computed emission tomography (SPECT) system could provide accurate arterial input functions for quantification of myocardial perfusion imaging using kinetic models. The errors induced by data inconsistency associated with imaging with slow camera rotation during tracer injection were evaluated with an approach called SPECT/P (dynamic SPECT from positron emission tomography (PET)) and SPECT/D (dynamic SPECT from database of SPECT phantom projections). SPECT/P simulated SPECT-like dynamic projections using reprojections of reconstructed dynamic (94)Tc-methoxyisobutylisonitrile ((94)Tc-MIBI) PET images acquired in three human subjects (1 min infusion). This approach was used to evaluate the accuracy of estimating myocardial wash-in rate parameters K(1) for rotation speeds providing 180° of projection data every 27 or 54 s. Blood input and myocardium tissue time-activity curves (TACs) were estimated using spatiotemporal splines. These were fit to a one-compartment perfusion model to obtain wash-in rate parameters K(1). For the second method (SPECT/D), an anthropomorphic cardiac torso phantom was used to create real SPECT dynamic projection data of a tracer distribution derived from (94)Tc-MIBI PET scans in the blood pool, myocardium, liver and background. This method introduced attenuation, collimation and scatter into the modeling of dynamic SPECT projections. Both approaches were used to evaluate the accuracy of estimating myocardial wash-in parameters for rotation speeds providing 180° of projection data every 27 and 54 s. Dynamic cardiac SPECT was also performed in a human subject at rest using a hybrid SPECT/CT scanner. Dynamic measurements of (99m)Tc-tetrofosmin in the myocardium were obtained using an infusion time of 2 min. Blood input, myocardium tissue and liver TACs were estimated using the same spatiotemporal splines. The spatiotemporal maximum-likelihood expectation-maximization (4D ML-EM) reconstructions gave more accurate reconstructions than did standard frame-by-frame static 3D ML-EM reconstructions. The SPECT/P results showed that 4D ML-EM reconstruction gave higher and more accurate estimates of K(1) than did 3D ML-EM, yielding anywhere from a 44% underestimation to 24% overestimation for the three patients. The SPECT/D results showed that 4D ML-EM reconstruction gave an overestimation of 28% and 3D ML-EM gave an underestimation of 1% for K(1). For the patient study the 4D ML-EM reconstruction provided continuous images as a function of time of the concentration in both ventricular cavities and myocardium during the 2 min infusion. It is demonstrated that a 2 min infusion with a two-headed SPECT system rotating 180° every 54 s can produce measurements of blood pool and myocardial TACs, though the SPECT simulation studies showed that one must sample at least every 30 s to capture a 1 min infusion input function.     

A dual modality approach to quantitative quality control in emission tomography

Routine quality control (QC) and optimization of image quality of reconstructed images in single photon emission computed tomography (SPECT) and positron emission tomography (PET) remains a relatively qualitative exercise. With the advent of combined SPECT/CT and PET/CT devices, and accurate post hoc co-registration algorithms, the potential exists to utilize high resolution structural information for QC evaluation in addition to their use for anatomical correlation in clinical studies. The aim of this work was to explore, in principle, the uses of x-ray CT data of QC phantoms used in SPECT and PET to develop more objective assessments of performance of the emission tomographic (ET) devices and reconstructed data. A CT reconstruction of a novel ET QC phantom was segmented into the various compartments it contained. Using software, the voxel values in the different compartments were then altered to correspond to the concentration of the radioactivity in the actual scan of the same phantom on the SPECT system. This produces a high resolution version of a 'perfect' ET scan. Image co-registration techniques were then used to spatially align the synthetic high resolution SPECT scan to the measured SPECT scan. Various parameters can then be objectively derived from the registered data, for example, image contrast, spatial resolution, spatial non-uniformity, etc. In this study, we have used this approach to estimate spatial resolution (full width at half maximum, FWHM) and recovered contrast in reconstructed images of a SPECT phantom. Two independent methods were used to measure spatial resolution, obtaining excellent agreement. In conclusion, the ability to produce high resolution synthetic phantoms in emission tomography QC affords an objective approach to assessing system performance and optimizing protocols which is readily automated and quantifiable.     

Cone beam collimation for single photon emission computed tomography: analysis, simulation, and image reconstruction using filtered backprojection

This paper presents an analysis of two cone beam configurations (having focal lengths of 40 and 60 cm) for the acquisition of single photon emission computed tomography (SPECT) projection data. A three-dimensional filtered backprojection algorithm is used to reconstruct SPECT images of cone beam projection data obtained using Monte Carlo simulations. The mathematical analysis resulted in on-axis point source sensitivities (calculated for a distance of 15 cm from the collimator surface) for cone beam configurations that were 1.4-3 times the sensitivities of parallel-hole and fan beam geometries having similar geometric resolutions. Cone beam collimation offers the potential for improved sensitivity for SPECT devices using large-field-of-view scintillation cameras.     

Noise, resolution, and sensitivity considerations in the design of a single-slice emission-transmission computed tomographic system

A prototype Emission-Transmission Computed Tomography (ETCT) system is being developed that will acquire single-slice x-ray transmission CT images simultaneously with single photon emission computed tomography (SPECT) images. This system will permit the correlation of anatomical information from x-ray CT with functional information from SPECT images. The patient-specific attenuation map derived from the x-ray CT images can be used to perform attenuation correction of the SPECT images, so that accurate quantitative information can be obtained. The fan-beam scanning geometry and the use of a segmented HPGe detector array impose special constraints on the design of the collimator for the system. Based on a signal detection model, an efficiency-resolution figure of merit (ERFM) as a function of the collimator geometric efficiency, system resolution width, and object diameter is defined. The ERFM is proportional to the square of the detection signal-to-noise ratio. The collimator design parameters can then be optimized by optimizing the ERFM for an anticipated object diameter. The collimator point-spread function, geometric efficiency, and resolution are calculated. The collimator optimized for the detection of a 1-cm object will have a single-slice point source efficiency of 1.2 X 10(-4), and a FWHM of 6.5 mm at the center of the reconstruction circle, at 12 cm from the collimator face. The minimum object contrast which will give a detection SNR of 5 is 74%, for a total accumulated count per slice of 2 X 10(6).     

Computerized mass detection for digital breast tomosynthesis directly from the projection images

Digital breast tomosynthesis (DBT) has recently emerged as a new and promising three-dimensional modality in breast imaging. In DBT, the breast volume is reconstructed from 11 projection images, taken at source angles equally spaced over an arc of 50 degrees. Reconstruction algorithms for this modality are not fully optimized yet. Because computerized lesion detection in the reconstructed breast volume will be affected by the reconstruction technique, we are developing a novel mass detection algorithm that operates instead on the set of raw projection images. Mass detection is done in three stages. First, lesion candidates are obtained for each projection image separately, using a mass detection algorithm that was initially developed for screen-film mammography. Second, the locations of a lesion candidate are backprojected into the breast volume. In this feature volume, voxel intensities are a combined measure of detection frequency (e.g., the number of projections in which a given lesion candidate was detected), and a measure of the angular range over which a given lesion was detected. Third, features are extracted after reprojecting the three-dimensional (3-D) locations of lesion candidates into projection images. Features are combined using linear discriminant analysis. The database used to test the algorithm consisted of 21 mass cases (13 malignant, 8 benign) and 15 cases without mass lesions. Based on this database, the algorithm yielded a sensitivity of 90% at 1.5 false positives per breast volume. Algorithm performance is positively biased because this dataset was used for development, training, and testing, and because the number of algorithm parameters was approximately the same as the number.of patient cases. Our results indicate that computerized mass detection in the sequence of projection images for DBT may be effective despite the higher noise level in those images.     

A multipinhole small animal SPECT system with submillimeter spatial resolution

Single photon emission computed tomography (SPECT) is an important technology for molecular imaging studies of small animals. In this arena, there is an increasing demand for high performance imaging systems that offer improved spatial resolution and detection efficiency. We have designed a multipinhole small animal imaging system based on position sensitive avalanche photodiode (PSAPD) detectors with the goal of submillimeter spatial resolution and high detection efficiency, which will allow us to minimize the radiation dose to the animal and to shorten the time needed for the imaging study. Our design will use 8 x 24 mm2 PSAPD detector modules coupled to thallium-doped cesium iodide [CsI(Tl)] scintillators, which can achieve an intrinsic spatial resolution of 0.5 mm at 140 keV. These detectors will be arranged in rings of 24 modules each; the animal is positioned in the center of the 9 stationary detector rings which capture projection data from the animal with a cylindrical tungsten multipinhole collimator. The animal is supported on a bed which can be rocked about the central axis to increase angular sampling of the object. In contrast to conventional SPECT pinhole systems, in our design each pinhole views only a portion of the object. However, the ensemble of projection data from all of the multipinhole detectors provide angular sampling that is sufficient to reconstruct tomographic data from the object. The performance of this multipinhole PSAPD imaging system was simulated using a ray tracing program that models the appropriate point spread functions and then was compared against the performance of a dual-headed pinhole SPECT system. The detection efficiency of both systems was simulated and projection data of a hot rod phantom were generated and reconstructed to assess spatial resolution. Appropriate Poisson noise was added to the data to simulate an acquisition time of 15 min and an activity of 18.5 MBq distributed in the phantom. Both sets of data were reconstructed with an ML-EM reconstruction algorithm. In addition, the imaging performance of both systems was evaluated with a uniformity phantom and a realistic digital mouse phantom. Simulations show that our proposed system produces a spatial resolution of 0.8 mm and an average detection efficiency of 630 cps/MBq. In contrast, simulations of the dual-headed pinhole SPECT system produce a spatial resolution of 1.1 mm and an average detection efficiency of 53 cps/MBq. These results suggest that our novel design will achieve high spatial resolution and will improve the detection efficiency by more than an order of magnitude compared to a dual-headed pinhole SPECT system. We expect that this system can perform SPECT with submillimeter spatial resolution, high throughput, and low radiation dose suitable for in vivo imaging of small animals.     

Three-dimensional SPECT reconstruction of combined cone beam and parallel beam data.

Single photon emission computed tomography (SPECT) using cone beam (CB) collimation exhibits increased sensitivity compared with acquisition geometries using parallel (P) hole collimation. However, CB collimation has a smaller field-of-view which may result in truncated projections and image artifacts. A primary objective of this work is to investigate maximum likelihood-expectation maximization (ML-EM) methods to reconstruct simultaneously acquired parallel and cone beam (P&CB) SPECT data. Simultaneous P&CB acquisition can be performed with commercially available triple camera systems by using two cone-beam collimators and a single parallel-hole collimator. The loss in overall sensitivity (relative to the use of three CB collimators) is about 15 to 20%. We have developed three methods to combine P&CB data using modified ML-EM algorithms. The first method consists of using both data sets to reconstruct a single intermediate image after each iteration using the ML-EM algorithm. The other two iterative algorithms combine intermediate parallel beam (PB) and CB source estimates to enhance image quality. For these methods, a PB estimate and a CB estimate are obtained for the first iteration. The second method consists of summing the PB and CB estimates for each subsequent iteration to obtain new PB and CB estimates. The third method is similar to the second method, with the exception that the new PB estimate simply is set equal to the PB estimate after each iteration. The combined source estimate is used in each subsequent iteration step of the EM algorithm. These algorithms are evaluated using projection data simulated using a Monte Carlo SPECT model. The P&CB SPECT images demonstrate marked improvement as compared with the CB-only reconstruction, particularly when the projections are truncated.     

Gaussian prefiltering of 123I DAT SPECT images when using depth-independent resolution recovery

Previously we have investigated a depth-independent compensation for collimator detector response (CDR) included in the OSEM reconstruction, intended for SPECT images that have been corrected for scatter and septal penetration using convolution-based methods. In this work, the aim was to study how different filtering strategies affect contrast as a function of noise when using Gaussian smoothing filters in combination with the above-described CDR compensation. The evaluation was performed for (123)I dopamine transporter (DAT) SPECT images. Prefiltering with 2D Gaussian filter kernels, where the deterioration in resolution is included in the depth-independent CDR compensation, was compared to conventional postfiltering with 3D Gaussian filter kernels. Images reconstructed without filtering are also included in the comparison. It was found that there is little benefit in noise reduction when using CDR compensation. However, this variant of prefiltering gives consistently higher contrasts as a function of noise compared with the postfiltering alternative, and that could be of interest when using other types of filters with contrast improving properties.     

Experimental results from a prototype slit-slat collimator with mixed multiplexed and non-multiplexed data

We have previously shown with simulations that a gain in signal-to-noise ratio (SNR) can be obtained by using mixed multiplexed (MX) and non-MX data in a slit-slat SPECT system as compared to using non-MX data only. We have now developed a prototype slit-slat collimator for a conventional gamma camera in order to validate these simulation results. The prototype collimator consists of seven slits and multiple parallel slats. Image reconstruction is performed using a modified OSEM algorithm, which takes into account geometric sensitivity variations and attenuation, but not scatter or resolution effects. Here, we first describe the calibration of the system and then we present the experimental validation with phantom experiments. SPECT acquisitions using different geometric and anthropomorphic phantoms were performed with and without multiplexing. The results show that reconstruction of the MX projections with the non-MX-projections eliminates artefacts caused by multiplexing. SNR gains obtained using the mixed MX and non-MX configurations were in the range of 26% to 51% for different phantoms. The results were in agreement with our previously published simulation work, proving that combining MX and non-MX data can result in artefact-free reconstructed images with improved SNR.     

System characteristics of SPECT with a slat collimated strip detector

In classical SPECT with parallel hole collimation, the sensitivity is constant over the field of view (FOV). This is no longer the case if a rotating slat collimator with planar photon collection is used: there will be a significant variation of the sensitivity within the FOV. Since not compensating for this inhomogeneous sensitivity distribution would result in non-quantitative images, an accurate knowledge of the sensitivity is mandatory to account for it during reconstruction. On the other hand, the spatial resolution versus distance dependency remains unaltered compared to parallel hole collimation. For deriving the sensitivity, different factors have to be taken into account: a first factor concerns the intrinsic detector properties and will be incorporated into the calculations as a detection efficiency term depending on the incident angle. The calculations are based on a second and more pronounced factor: the collimator and detector geometry. Several assumptions will be made for the calculation of the sensitivity formulae and it will be proven that these calculations deliver a valid prediction of the sensitivity at points far enough from the collimator. To derive a close field model which also accounts for points close to the collimator surface, a modified calculation method is used. After calculating the sensitivity in one plane it is easy to obtain the tomographic sensitivity. This is done by rotating the sensitivity maps for spin and camera rotation. The results derived from the calculations are then compared to simulation results and both show good agreement after including the aforementioned detection efficiency term. The validity of the calculations is also proven by measuring the sensitivity in the FOV of a prototype rotating slat gamma camera. An expression for the resolution of these planar collimation systems is obtained. It is shown that for equal collimator dimensions the same resolution-distance relationship is obtained as for parallel hole collimators. Although, a better spatial resolution can be obtained with our prototype camera due to the smaller pitch of the slats. This can be achieved without a major drop in system sensitivity due to the fact that the slats consist of less collimator material compared to a parallel hole collimator. The accuracy of the calculated resolution is proven by comparison with Monte Carlo simulation and measurement resolution values.