Fast spiral SPECT with stationary γ-cameras and focusing pinholes

Small-animal SPECT systems with stationary detectors and focusing multiple pinholes can achieve excellent resolution-sensitivity trade-offs. These systems are able to perform fast total-body scans by shifting the animal bed through the collimator using an automated xyz stage. However, so far, a large number of highly overlapping central fields of view have been used, at the cost of overhead time needed for animal repositioning and long image reconstruction times due to high numbers of projection views.

METHODS

To improve temporal resolution and reduce image reconstruction time for such scans, we have developed and tested spiral trajectories (STs) of the animal bed requiring fewer steps. In addition, we tested multiplane trajectories (MPTs) of the animal bed, which is the standard acquisition method of the U-SPECT-II system that is used in this study. Neither MPTs nor STs require rotation of the animal. Computer simulations and physical phantom experiments were performed for a wide range of numbers of bed positions. Furthermore, we tested STs in vivo for fast dynamic mouse scans.

RESULTS

We found that STs require less than half the number of bed positions of MPTs to achieve sufficient sampling. The reduced number of bed positions made it possible to perform a dynamic total-body bone scan and a dynamic hepatobiliary scan with time resolutions of 60 s and 15 s, respectively.

CONCLUSION

STs open up new possibilities for high throughput and fast dynamic radio-molecular imaging.

Quantitative multi-pinhole small-animal SPECT: uniform versus non-uniform Chang attenuation correction

Attenuation of photon flux on trajectories between the source and pinhole apertures affects the quantitative accuracy of reconstructed single-photon emission computed tomography (SPECT) images. We propose a Chang-based non-uniform attenuation correction (NUA-CT) for small-animal SPECT/CT with focusing pinhole collimation, and compare the quantitative accuracy with uniform Chang correction based on (i) body outlines extracted from x-ray CT (UA-CT) and (ii) on hand drawn body contours on the images obtained with three integrated optical cameras (UA-BC). Measurements in phantoms and rats containing known activities of isotopes were conducted for evaluation. In (125)I, (201)Tl, (99m)Tc and (111)In phantom experiments, average relative errors comparing to the gold standards measured in a dose calibrator were reduced to 5.5%, 6.8%, 4.9% and 2.8%, respectively, with NUA-CT. In animal studies, these errors were 2.1%, 3.3%, 2.0% and 2.0%, respectively. Differences in accuracy on average between results of NUA-CT, UA-CT and UA-BC were less than 2.3% in phantom studies and 3.1% in animal studies except for (125)I (3.6% and 5.1%, respectively). All methods tested provide reasonable attenuation correction and result in high quantitative accuracy. NUA-CT shows superior accuracy except for (125)I, where other factors may have more impact on the quantitative accuracy than the selected attenuation correction.

Maximum-likelihood scintillation detection for EM-CCD based gamma cameras

Gamma cameras based on charge-coupled devices (CCDs) coupled to continuous scintillation crystals can combine a good detection efficiency with high spatial resolutions with the aid of advanced scintillation detection algorithms. A previously developed analytical multi-scale algorithm (MSA) models the depth-dependent light distribution but does not take statistics into account. Here we present and validate a novel statistical maximum-likelihood algorithm (MLA) that combines a realistic light distribution model with an experimentally validated statistical model. The MLA was tested for an electron multiplying CCD optically coupled to CsI(Tl) scintillators of different thicknesses. For (99m)Tc imaging, the spatial resolution (for perpendicular and oblique incidence), energy resolution and signal-to-background counts ratio (SBR) obtained with the MLA were compared with those of the MSA. Compared to the MSA, the MLA improves the energy resolution by more than a factor of 1.6 and the SBR is enhanced by more than a factor of 1.3. For oblique incidence (approximately 45°), the depth-of-interaction corrected spatial resolution is improved by a factor of at least 1.1, while for perpendicular incidence the MLA resolution does not consistently differ significantly from the MSA result for all tested scintillator thicknesses. For the thickest scintillator (3 mm, interaction probability 66% at 141 keV) a spatial resolution (perpendicular incidence) of 147 µm full width at half maximum (FWHM) was obtained with an energy resolution of 35.2% FWHM. These results of the MLA were achieved without prior calibration of scintillations as is needed for many statistical scintillation detection algorithms. We conclude that the MLA significantly improves the gamma camera performance compared to the MSA.

The effects of object activity distribution on multiplexing multi-pinhole SPECT

We aim to study the effects of activity distribution for multiplexing multi-pinhole (MPH) SPECT. Three digital phantoms, including a hot rod, a cold rod and a cold sphere phantom, were used. Different degrees of multiplexing were obtained by (i) adjusting the MPH pattern for the same 4-pinhole collimator (scheme 1) and (ii) increasing the number of pinholes (scheme 2). Noise-free and noisy projections were generated using a 3D analytical MPH projector based on the same acquisition time. Projections were reconstructed using OS-EM without resolution recovery. Normalized mean-square-error (NMSE), noise, image profiles and signal-to-background ratios (SBR) were assessed. For the hot rod phantom, the NMSE-noise trade-offs slightly improves for multiplexing designs in scheme 2. Substantial artifacts were observed and the NMSE-noise trade-offs slightly worsened for multiplexing designs for the cold phantoms. Resolutions slightly degraded for higher degrees of multiplexing (∼39-65%) for the cold rod phantom. For the cold sphere phantom, image profiles showed non-multiplexing designs better emulated the phantom, while ∼20% multiplexing performs similarly as compared to non-multiplexing in SBR. Our results indicate that multiplexing can help for sparse objects but leads to a significant image degradation in non-sparse distributions. Since many tracers are not highly specific, and the gain of detection efficiency by allowing multiplexing is fairly offset by image degradations, multiplexing needs to be kept to a minimum for optimum MPH collimator designs.

Targeted multi-pinhole SPECT

Purpose

Small-animal single photon emission computed tomography (SPECT) with focused multi-pinhole collimation geometries allows scanning modes in which large amounts of photons can be collected from specific volumes of interest. Here we present new tools that improve targeted imaging of specific organs and tumours, and validate the effects of improved targeting of the pinhole focus.

Methods

A SPECT system with 75 pinholes and stationary detectors was used (U-SPECT-II). An XYZ stage automatically translates the animal bed with a specific sequence in order to scan a selected volume of interest. Prior to stepping the animal through the collimator, integrated webcams acquire images of the animal. Using sliders, the user designates the desired volume to be scanned (e.g. a xenograft or specific organ) on these optical images. Optionally projections of an atlas are overlaid semiautomatically to locate specific organs. In order to assess the effects of more targeted imaging, scans of a resolution phantom and a mouse myocardial phantom, as well as in vivo mouse cardiac and tumour scans, were acquired with increased levels of targeting. Differences were evaluated in terms of count yield, hot rod visibility and contrast-to-noise ratio.

Results

By restricting focused SPECT scans to a 1.13-ml resolution phantom, count yield was increased by a factor 3.6, and visibility of small structures was significantly enhanced. At equal noise levels, the small-lesion contrast measured in the myocardial phantom was increased by 42%. Noise in in vivo images of a tumour and the mouse heart was significantly reduced.

Conclusion

Targeted pinhole SPECT improves images and can be used to shorten scan times. Scan planning with optical cameras provides an effective tool to exploit this principle without the necessity for additional X-ray CT imaging.

An enhanced high-resolution EMCCD-based gamma camera using SiPM side detection

Electron-multiplying charge-coupled devices (EMCCDs) coupled to scintillation crystals can be used for high-resolution imaging of gamma rays in scintillation counting mode. However, the detection of false events as a result of EMCCD noise deteriorates the spatial and energy resolution of these gamma cameras and creates a detrimental background in the reconstructed image. In order to improve the performance of an EMCCD-based gamma camera with a monolithic scintillation crystal, arrays of silicon photon-multipliers (SiPMs) can be mounted on the sides of the crystal to detect escaping scintillation photons, which are otherwise neglected. This will provide a priori knowledge about the correct number and energies of gamma interactions that are to be detected in each CCD frame. This information can be used as an additional detection criterion, e.g. for the rejection of otherwise falsely detected events. The method was tested using a gamma camera based on a back-illuminated EMCCD, coupled to a 3 mm thick continuous CsI:Tl crystal. Twelve SiPMs have been mounted on the sides of the CsI:Tl crystal. When the information of the SiPMs is used to select scintillation events in the EMCCD image, the background level for (99m)Tc is reduced by a factor of 2. Furthermore, the SiPMs enable detection of (125)I scintillations. A hybrid SiPM-/EMCCD-based gamma camera thus offers great potential for applications such as in vivo imaging of gamma emitters.

Absolute quantitative total-body small-animal SPECT with focusing pinholes

PURPOSE

In pinhole SPECT, attenuation of the photon flux on trajectories between source and pinholes affects quantitative accuracy of reconstructed images. Previously we introduced iterative methods that compensate for image degrading effects of detector and pinhole blurring, pinhole sensitivity and scatter for multi-pinhole SPECT. The aim of this paper is (1) to investigate the accuracy of the Chang algorithm in rodents and (2) to present a practical Chang-based method using body outline contours obtained with optical cameras.

METHODS

Here we develop and experimentally validate a practical method for attenuation correction based on a Chang first-order method. This approach has the advantage that it is employed after, and therefore independently from, iterative reconstruction. Therefore, no new system matrix has to be calculated for each specific animal. Experiments with phantoms and animals were performed with a high-resolution focusing multi-pinhole SPECT system (U-SPECT-II, MILabs, The Netherlands). This SPECT system provides three additional optical camera images of the animal for each SPECT scan from which the animal contour can be estimated.

RESULTS

Phantom experiments demonstrated that an average quantification error of -18.7% was reduced to -1.7% when both window-based scatter correction and Chang correction based on the body outline from optical images were applied. Without scatter and attenuation correction, quantification errors in a sacrificed rat containing sources with known activity ranged from -23.6 to -9.3%. These errors were reduced to values between -6.3 and +4.3% (with an average magnitude of 2.1%) after applying scatter and Chang attenuation correction.

CONCLUSION

We conclude that the modified Chang correction based on body contour combined with window-based scatter correction is a practical method for obtaining small-animal SPECT images with high quantitative accuracy.

Theoretical analysis of full-ring multi-pinhole brain SPECT

Presently used clinical brain SPECT suffers from limited spatio-temporal resolution. Here we investigate the feasibility of high-resolution and high-sensitivity full-ring multi-pinhole brain SPECT (MP-SPECT). Using an analytical model we optimized pinhole-detector geometries of MP-SPECT for different detector intrinsic resolutions R(i). System resolution and sensitivity of optimized MP-SPECT were compared to conventional clinical SPECT. The comparison of the system resolution of different systems was done at matched sensitivity, which was achieved by tuning pinhole diameters. Similarly, sensitivities were compared at matched system resolution. For MP-SPECT that uses detectors with intrinsic resolutions of 4 mm > R(i) 0.5 mm a sensitivity can be achieved that is 6.0 times higher than the sensitivity of conventional dual-head SPECT systems with parallel-hole collimators (DualPar), while system resolution can be improved by a factor of 2.4. To achieve these improvements a large detector-to-collimator distance is needed. In contrast, for detectors with intrinsic resolutions <0.2 mm, it is beneficial to place the detectors close to the pinholes, resulting in a high number of de-magnified projections. For a detector intrinsic resolution of 0.05 mm, a 14.5-fold improvement in sensitivity and a 3.8-fold improvement in system resolution compared to DualPar is predicted. Furthermore, we found that for optimized MP-SPECT the sensitivity scales proportionally to system resolution squared, with the proportionality constant depending on R(i). From our sensitivity-system resolution trade-off equations we deduced that MP-SPECT with an ideal detector (R(i) --> 0) can have a system resolution that is 2.0 times better than optimized MP-SPECT with a conventional detector (R(i) approximately 3 mm). The high performance of optimized MP-SPECT may open up completely new molecular imaging applications.

Theoretical analysis of full-ring multi-pinhole brain SPECT

Presently used clinical brain SPECT suffers from limited spatio-temporal resolution. Here we investigate the feasibility of high-resolution and high-sensitivity full-ring multi-pinhole brain SPECT (MP-SPECT). Using an analytical model we optimized pinhole-detector geometries of MP-SPECT for different detector intrinsic resolutions R(i). System resolution and sensitivity of optimized MP-SPECT were compared to conventional clinical SPECT. The comparison of the system resolution of different systems was done at matched sensitivity, which was achieved by tuning pinhole diameters. Similarly, sensitivities were compared at matched system resolution. For MP-SPECT that uses detectors with intrinsic resolutions of 4 mm > R(i) 0.5 mm a sensitivity can be achieved that is 6.0 times higher than the sensitivity of conventional dual-head SPECT systems with parallel-hole collimators (DualPar), while system resolution can be improved by a factor of 2.4. To achieve these improvements a large detector-to-collimator distance is needed. In contrast, for detectors with intrinsic resolutions <0.2 mm, it is beneficial to place the detectors close to the pinholes, resulting in a high number of de-magnified projections. For a detector intrinsic resolution of 0.05 mm, a 14.5-fold improvement in sensitivity and a 3.8-fold improvement in system resolution compared to DualPar is predicted. Furthermore, we found that for optimized MP-SPECT the sensitivity scales proportionally to system resolution squared, with the proportionality constant depending on R(i). From our sensitivity-system resolution trade-off equations we deduced that MP-SPECT with an ideal detector (R(i) --> 0) can have a system resolution that is 2.0 times better than optimized MP-SPECT with a conventional detector (R(i) approximately 3 mm). The high performance of optimized MP-SPECT may open up completely new molecular imaging applications.

Small field-of-view dedicated cardiac SPECT systems: impact of projection truncation

PURPOSE

Small field-of-view (FOV) dedicated cardiac SPECT systems suffer from truncated projection data. This results in (1) neglect of liver activity that otherwise could be used to estimate (and subsequently correct) the amount of scatter in the myocardium by model-based scatter correction, and (2) distorted attenuation maps. In this study, we investigated to what extent truncation impacts attenuation correction and model-based scatter correction in the cases of (99m)Tc, (201)Tl, and simultaneous (99m)Tc/(201)Tl studies. In addition, we evaluated a simple correction method to mitigate the effects of truncation.

METHODS

Digital thorax phantoms of different sizes were used to simulate the full FOV SPECT projections for (99m)Tc, (201)Tl, and simultaneous (99m)Tc/(201)Tl studies. Small FOV projections were obtained by artificially truncating the full FOV projections. Deviations from ideal heart positioning were simulated by axially shifting projections resulting in more severe liver truncation. Effects of truncation on SPECT images were tested for ordered subset (OS) expectation maximization reconstruction with (1) attenuation correction and detector response modelling (OS-AD), and (2) with additional Monte-Carlo-based scatter correction (OS-ADS). To correct truncation-induced artefacts, we axially extended truncated projections on both sides by duplicating pixel values on the projection edge.

RESULTS

For both (99m)Tc and (201)Tl, differences in the reconstructed myocardium between full FOV and small FOV projections were negligible. In the nine myocardial segments, the maximum deviations of the average pixel values were 1.3% for OS-AD and 3.5% for OS-ADS. For the simultaneous (99m)Tc/(201)Tl studies, reconstructed (201)Tl SPECT images from full FOV and small FOV projections showed clearly different image profiles due to truncation. The maximum deviation in defected segments was found to be 49% in the worst-case scenario. However, artificially extending projections reduced deviations in defected segments to a few percent.

CONCLUSION

Our results indicate that, for single isotope studies, using small FOV systems has little impact on attenuation correction and model-based scatter correction. For simultaneous (99m)Tc/(201)Tl studies, artificial projection extension almost fully eliminates the adverse effects of projection truncation.

A pinhole gamma camera with optical depth-of-interaction elimination

The performance of pinhole single photon emission computed tomography (SPECT) depends on the spatial resolution of the gamma-ray detectors used. Pinhole cameras suffer from strong resolution loss due to the varying depth-of-interaction (DOI) of gamma quanta that enter the detector material at an angle. We eliminate DOI effects in a scintillation gamma camera via a dedicated optic fiber bundle that acts as a focusing collimator for light generated in a scintillation crystal. A curved crystal is connected to a concavely shaped fiber-optic bundle such that the fibers connect perpendicularly to the crystal's convex surface and point straight at the pinhole opening. Limiting the fiber numerical apertures can be used to suppress resolution losses due to light spread. Here we demonstrate experimentally that this prototype position-sensitive gamma sensor successfully eliminates DOI effects, and has an intrinsic resolution of better than 280 microm full width at half maximum with an interaction probability of 67% for 140 keV photons. Therefore, the detector has great potential for increasing the resolution of pinhole SPECT.

A novel, SiPM-array-based, monolithic scintillator detector for PET

Silicon photomultipliers (SiPMs) are of great interest to positron emission tomography (PET), as they enable new detector geometries, for e.g., depth-of-interaction (DOI) determination, are MR compatible, and offer faster response and higher gain than other solid-state photosensors such as avalanche photodiodes. Here we present a novel detector design with DOI correction, in which a position-sensitive SiPM array is used to read out a monolithic scintillator. Initial characterization of a prototype detector consisting of a 4 x 4 SiPM array coupled to either the front or back surface of a 13.2 mm x 13.2 mm x 10 mm LYSO:Ce(3+) crystal shows that front-side readout results in significantly better performance than conventional back-side readout. Spatial resolutions <1.6 mm full-width-at-half-maximum (FWHM) were measured at the detector centre in response to an approximately 0.54 mm FWHM diameter test beam. Hardly any resolution losses were observed at angles of incidence up to 45 degrees , demonstrating excellent DOI correction. About 14% FWHM energy resolution was obtained. The timing resolution, measured in coincidence with a BaF(2) detector, equals 960 ps FWHM.

A micro-machined retro-reflector for improving light yield in ultra-high-resolution gamma cameras

High-resolution imaging of x-ray and gamma-ray distributions can be achieved with cameras that use charge coupled devices (CCDs) for detecting scintillation light flashes. The energy and interaction position of individual gamma photons can be determined by rapid processing of CCD images of individual flashes. Here we investigate the improvement of such a gamma camera when a micro-machined retro-reflector is used to increase the light output of a continuous scintillation crystal. At 122 keV we found that retro-reflectors improve the intrinsic energy resolution (full width at half maximum (FWHM)) by 32% (from 50% to 34%) and the signal-to-noise (SNR) ratio by 18%. The spatial resolution (FWHM) was improved by about 4%, allowing us to obtain a resolution of 159 microm. The full width at tenth maximum (FWTM) improvement was 13%. Therefore, this enhancement is a next step towards realizing compact high-resolution devices for imaging gamma emitters.

U-SPECT-II: An Ultra-High-Resolution Device for Molecular Small-Animal Imaging

We present a new rodent SPECT system (U-SPECT-II) that enables molecular imaging of murine organs down to resolutions of less than half a millimeter and high-resolution total-body imaging.

METHODS

The U-SPECT-II is based on a triangular stationary detector set-up, an XYZ stage that moves the animal during scanning, and interchangeable cylindric collimators (each containing 75 pinhole apertures) for both mouse and rat imaging. A novel graphical user interface incorporating preselection of the field of view with the aid of optical images of the animal focuses the pinholes to the area of interest, thereby maximizing sensitivity for the task at hand. Images are obtained from list-mode data using statistical reconstruction that takes system blurring into account to increase resolution.

RESULTS

For (99m)Tc, resolutions determined with capillary phantoms were smaller than 0.35 and 0.45 mm using the mouse collimator with 0.35- and 0.6-mm pinholes, respectively, and less than 0.8 mm using the rat collimator with 1.0-mm pinholes. Peak geometric sensitivity is 0.07% and 0.18% for the mouse collimator with 0.35- and 0.6-mm pinholes, respectively, and 0.09% for the rat collimator. Resolution with (111)In, compared with that with (99m)Tc, was barely degraded, and resolution with (125)I was degraded by about 10%, with some additional distortion. In vivo, kidney, tumor, and bone images illustrated that U-SPECT-II could be used for novel applications in the study of dynamic biologic systems and radiopharmaceuticals at the suborgan level.

CONCLUSION

Images and movies obtained with U-SPECT-II provide high-resolution radiomolecule visualization in rodents. Discrimination of molecule concentrations between adjacent volumes of about 0.04 microL in mice and 0.5 microL in rats with U-SPECT-II is readily possible.

Monolithic scintillator PET detectors with intrinsic depth-of-interaction correction

We developed positron emission tomography (PET) detectors based on monolithic scintillation crystals and position-sensitive light sensors. Intrinsic depth-of-interaction (DOI) correction is achieved by deriving the entry points of annihilation photons on the front surface of the crystal from the light sensor signals. Here we characterize the next generation of these detectors, consisting of a 20 mm thick rectangular or trapezoidal LYSO:Ce crystal read out on the front and the back (double-sided readout, DSR) by Hamamatsu S8550SPL avalanche photodiode (APD) arrays optimized for DSR. The full width at half maximum (FWHM) of the detector point-spread function (PSF) obtained with a rectangular crystal at normal incidence equals approximately 1.05 mm at the detector centre, after correction for the approximately 0.9 mm diameter test beam of annihilation photons. Resolution losses of several tenths of a mm occur near the crystal edges. Furthermore, trapezoidal crystals perform almost equally well as rectangular ones, while improving system sensitivity. Due to the highly accurate DOI correction of all detectors, the spatial resolution remains essentially constant for angles of incidence of up to at least 30 degrees . Energy resolutions of approximately 11% FWHM are measured, with a fraction of events of up to 75% in the full-energy peak. The coincidence timing resolution is estimated to be 2.8 ns FWHM. The good spatial, energy and timing resolutions, together with the excellent DOI correction and high detection efficiency of our detectors, are expected to facilitate high and uniform PET system resolution.

Multi-scale algorithm for improved scintillation detection in a CCD-based gamma camera

Gamma cameras based on charge-coupled devices (CCDs) and micro-columnar CsI scintillators can reach high spatial resolutions. However, the gamma interaction probability of these scintillators is low (typically <30% at 141 keV) due to the limited thickness of presently available micro-columnar scintillators. Continuous scintillators can improve the interaction probability but suffer from increased light spread compared to columnar scintillators. In addition, for both types of scintillators, gamma photons incident at an oblique angle reduce the spatial resolution due to the variable depth of interaction (DOI). To improve the spatial resolution and spectral characteristics of these detectors, we have developed a fast analytic scintillation detection algorithm that makes use of a depth-dependent light spread model and as a result is able to estimate the DOI in the scintillator. This algorithm, performing multi-scale frame analysis, was tested for an electron multiplying CCD (EM-CCD) optically coupled to CsI(Tl) scintillators of different thicknesses. For the thickest scintillator (2.6 mm) a spatial resolution of 148 microm full width half maximum (FWHM) was obtained with an energy resolution of 46% FWHM for perpendicularly incident gamma photons (interaction probability 61% at 141 keV). The multi-scale algorithm improves the spatial resolution up to 11%, the energy resolution up to 36% and the signal-to-background counts ratio up to 46% compared to a previously implemented algorithm that did not model the depth-dependent light spread. In addition, the multi-scale algorithm can accurately estimate DOI. As a result, degradation of the spatial resolution due to the variable DOI for gamma photons incident at a 45 degrees angle was improved from 2.0 10(3) to 448 microm FWHM. We conclude that the multi-scale algorithm significantly improves CCD-based gamma cameras as can be applied in future SPECT systems.

Intra-patient reproducibility of myocardial SPECT imaging with 201Tl

BACKGROUND

To define the physical and clinical reproducibility of (201)Tl myocardial perfusion SPECT (MPS), this study assesses the variation between two repeated rest (201)Tl MPS with repositioning only, with a two-hour time interval and with phantom measurements as a reference.

METHODS

Three anthropomorphic thorax phantoms were filled with (201)Tl. For each phantom five repeated (201)Tl MPS were obtained. In addition, in 20 patients repeated (201)Tl rest-MPS and in 26 patients early and delayed (201)Tl rest-MPS were performed. Quantitative analysis was done using MunichHeart. Statistical methods were used to calculate variability. Visual analysis was performed by 2 independent observers.

RESULTS

The average variation between repeated phantom MPS was 0.5% (95% confidence interval (CI): -0.4% to 1.4%). For patient scans this was -5.0% (95% CI: -2.5% to -7.5%) and between early and delayed (201)Tl MPS -15.5% (95% CI: -11.7% to -19.3%). Visual assessment revealed no clinical significant differences between rest (201)Tl and repeated or delayed (201)Tl MPS.

CONCLUSIONS

Repositioning in phantom (201)Tl MPS does not cause significant variation. Repeated (201)Tl MPS in patients shows 5.0% decrease of (201)Tl in 30 minutes, which increases to 15% during a two-hour time interval without quantitative or visual regional differences. This decrease indicates a time-related washout of (201)Tl, but does not change clinical diagnosis.

System calibration and statistical image reconstruction for ultra-high resolution stationary pinhole SPECT

For multipinhole single-photon emission computed tomography (SPECT), iterative reconstruction algorithms are preferred over analytical methods, because of the often complex multipinhole geometries and the ability of iterative algorithms to compensate for effects like spatially variant sensitivity and resolution. Ideally, such compensation methods are based on accurate knowledge of the position-dependent point spread functions (PSFs) specifying the response of the detectors to a point source at every position in the instrument. This paper describes a method for model-based generation of complete PSF lookup tables from a limited number of point-source measurements for stationary SPECT systems and its application to a submillimeter resolution stationary small-animal SPECT system containing 75 pinholes (U-SPECT-I). The method is based on the generalization over the entire object to be reconstructed, of a small number of properties of point-source responses which are obtained at a limited number of measurement positions. The full shape of measured point-source responses can almost be preserved in newly created PSF tables. We show that these PSFs can be used to obtain high-resolution SPECT reconstructions: the reconstructed resolutions judged by rod visibility in a micro-Derenzo phantom are 0.45 mm with 0.6-mm pinholes and below 0.35 mm with 0.3-mm pinholes. In addition, we show that different approximations, such as truncating the PSF kernel, with significant reduction of reconstruction time, can still lead to acceptable reconstructions.

Movies of dopamine transporter occupancy with ultra-high resolution focusing pinhole SPECT

A pivotal question in neuropharmacology is how the function of neurotransmitter systems relates to psychiatric diseases. In experimental neuropharmacology, we have dreamt about a looking glass that would allow us to see neurotransmitter systems in action, and about animals that would faithfully serve us as models for human psychiatric disease. Analysis of animal models has been limited by the availability of methods to study in vivo neurotransmitter dynamics. Now, a single photon emission computed tomography system called U-SPECT can localize dopamine transporters in sub-compartments of the mouse brain during a range of points in time. Applied to the midbrain dopamine system of different models of disease, this will aid the understanding of dynamic processes of this neurotransmitter that underlie brain functions and human brain pathology.

Optimizing multi-pinhole SPECT geometries using an analytical model

State-of-the-art multi-pinhole SPECT devices allow for sub-mm resolution imaging of radio-molecule distributions in small laboratory animals. The optimization of multi-pinhole and detector geometries using simulations based on ray-tracing or Monte Carlo algorithms is time-consuming, particularly because many system parameters need to be varied. As an efficient alternative we develop a continuous analytical model of a pinhole SPECT system with a stationary detector set-up, which we apply to focused imaging of a mouse. The model assumes that the multi-pinhole collimator and the detector both have the shape of a spherical layer, and uses analytical expressions for effective pinhole diameters, sensitivity and spatial resolution. For fixed fields-of-view, a pinhole-diameter adapting feedback loop allows for the comparison of the system resolution of different systems at equal system sensitivity, and vice versa. The model predicts that (i) for optimal resolution or sensitivity the collimator layer with pinholes should be placed as closely as possible around the animal given a fixed detector layer, (ii) with high-resolution detectors a resolution improvement up to 31% can be achieved compared to optimized systems, (iii) high-resolution detectors can be placed close to the collimator without significant resolution losses, (iv) interestingly, systems with a physical pinhole diameter of 0 mm can have an excellent resolution when high-resolution detectors are used.

Evaluation of 3D Monte Carlo-Based Scatter Correction for 201Tl Cardiac Perfusion SPECT

(201)Tl-Chloride ((201)Tl) is a myocardial perfusion SPECT agent with excellent biochemical properties commonly used for assessing tissue viability. However, cardiac (201)Tl SPECT images are severely degraded by photons scattered in the thorax. Accurate correction for this scatter is complicated by the nonuniform density and varied sizes of thoraxes, by the additional attenuation and scatter caused by female patients' breasts, and by the energy spectrum of (201)Tl. Monte Carlo simulation is a general and accurate method well suited to modeling this scatter.

METHODS

Statistical reconstruction that includes Monte Carlo modeling of scatter was compared with statistical reconstruction algorithms not corrected for scatter. In the ADS method, corrections for attenuation, detector response, and scatter (Monte Carlo-based) were implemented simultaneously via the dual-matrix ordered-subset expectation maximization algorithm with a Monte Carlo simulator as part of the forward projector. The ADS method was compared with the A method (ordered-subset expectation maximization with attenuation correction) and with the AD method (a method like the A method but with detector response modeling added). A dual-head SPECT system equipped with two (153)Gd scanning line sources was used for simultaneously acquiring transmission and emission data. Four clinically realistic phantom configurations (a large thorax and a small thorax, each with and without breasts) with a cardiac insert containing 2 cold defects were used to evaluate the proposed reconstruction algorithms. We compared the performance of the different algorithms in terms of noise properties, contrast-to-noise ratios, the contrast separability of perfusion defects, uniformity, and robustness to anatomic variations.

RESULTS

The ADS method provided images with clearly better visual defect contrast than did the other methods. The contrasts achieved with the ADS method were 10%-24% higher than those achieved with the AD method and 11%-37% higher than those achieved with the A method. For a typical contrast level, the ADS method exhibited noise levels around 27% lower than the AD method and 34% lower than the A method. Compared with the other 2 algorithms, the ADS reconstructions were less sensitive to anatomic variations and had better image uniformity in the homogeneously perfused myocardium. Finally, we found that the improvements that can be achieved with Monte Carlo-based scatter correction are stronger for (201)Tl than for (99m)Tc imaging.

CONCLUSION

Our results indicate that Monte Carlo-based scatter correction is suitable for (201)Tl cardiac imaging and that such correction simultaneously improves several image-quality metrics.

Front-illuminated versus back-illuminated photon-counting CCD-based gamma camera: important consequences for spatial resolution and energy resolution

Charge-coupled devices (CCDs) coupled to scintillation crystals can be used for high-resolution imaging with x-rays and gamma rays. When the CCD images can be read out fast enough, the energy and interaction position of individual gamma quanta can be estimated by a real-time image analysis of the scintillation light flashes ('photon-counting mode'). The electron-multiplying CCD (EMCCD) is well suited for fast read out, since even at high frame rates it has extremely low read-out noise. Back-illuminated (BI) EMCCDs have much higher quantum efficiency than front-illuminated (FI) EMCCDs. Here we compare the spatial and energy resolution of gamma cameras based on FI and BI EMCCDs. The CCDs are coupled to a 1000 microm thick columnar CsI(Tl) crystal for the purpose of Tc-99m and I-125 imaging. Intrinsic spatial resolutions of 44 microm for I-125 and 49 microm for Tc-99m were obtained when using a BI EMCCD, which is an improvement by a factor of about 1.2-2 over the FI EMCCD. Furthermore, in the energy spectrum of the BI EMCCD, the I-125 signal could be clearly separated from the background noise, which was not the case for the FI EMCCD. The energy resolution of a BI EMCCD for Tc-99m was estimated to be approximately 36 keV, full width at half maximum, at 141 keV. The excellent results for the BI EMCCD encouraged us to investigate the cooling requirements for our setup. We have found that for the BI EMCCD, the spatial and energy resolution, as well as image noise, remained stable over a range of temperatures from -50 degrees C to -15 degrees C. This is a significant advantage over the FI EMCCD, which suffered from loss of spatial and especially energy resolution at temperatures as low as -40 degrees C. We conclude that the use of BI EMCCDs may significantly improve the imaging capabilities and the cost efficiency of CCD-based high-resolution gamma cameras.

Statistical reconstruction for x-ray CT systems with non-continuous detectors

We analyse the performance of statistical reconstruction (SR) methods when applied to non-continuous x-ray detectors. Robustness to projection gaps is required in x-ray CT systems with multiple detector modules or with defective detector pixels. In such situations, the advantage of statistical reconstruction is that it is able to ignore missing or faulty pixels and that it makes optimal use of the remaining line integrals. This potentially obviates the need to fill the sinogram discontinuities by interpolation or any other approximative pre-processing techniques. In this paper, we apply SR to cone beam projections of (i) a hypothetical modular detector micro-CT scanner and of (ii) a system with randomly located defective detector elements. For the modular-detector system, SR produces reconstruction volumes free of noticeable gap-induced artefacts as long as the location of detector gaps and selection of the scanning range provide complete object sampling in the central imaging plane. When applied to randomly located faulty detector elements, SR produces images free of substantial ring artefacts even for cases where defective pixels cover as much as 3% of the detector area.

Evaluation of 3D Monte Carlo-based scatter correction for 99mTc cardiac perfusion SPECT

Cardiac SPECT images are degraded by photons scattered in the thorax. Accurate correction for scatter is complicated by the nonuniform density and varied sizes of thoraxes and by the additional attenuation and scatter caused by female patients' breasts. Monte Carlo simulation is a general and accurate method for detailed modeling of scatter. Hence, for 99mTc we compared statistical reconstruction based on Monte Carlo modeling of scatter with statistical reconstruction that incorporates the more commonly used triple-energy-window scatter correction. Both of these scatter correction methods were used in conjunction with attenuation correction and resolution recovery.

METHODS

Simultaneous attenuation, detector response, and Monte Carlo-based scatter correction were implemented via the dual-matrix ordered-subset expectation maximization algorithm with a Monte Carlo simulator as part of the forward projector (ADS-MC). ADS-MC was compared to (i) ordered-subset expectation maximization with attenuation correction and with detector response modeling (AD); (ii) like (i) but with scatter correction added using the triple-energy-window method (ADS-TEW). A dual-detector SPECT system equipped with 2 153Gd scanning line sources was used for acquiring 99mTc emission data and attenuation maps. Four clinically realistic phantom configurations (a large thorax and a small thorax, each with and without breasts) with a cardiac insert containing 2 cold defects were used to evaluate the proposed reconstruction algorithms. In these phantom configurations, we compared the performance of the algorithms in terms of noise properties, contrast-to-noise ratios, contrast separability of cold defects, and robustness to anatomic variation.

RESULTS

Noise was found to be approximately 14% lower in the ADS-MC images than in the ADS-TEW and AD images. Typically, the contrast obtained with the ADS-MC algorithm was 10%-20% higher than that obtained with the other 2 methods. Furthermore, compared with the other 2 algorithms, the ADS-MC method was less sensitive to anatomic variations.

CONCLUSION

Our results indicate that the imaging performance of 99mTc SPECT can be improved more by Monte Carlo-based scatter correction than by window-based scatter correction.

Hybrid scatter correction applied to quantitative holmium-166 SPECT

Ho-166 is a combined beta-gamma emitter of which the betas can be used therapeutically. From the 81 keV gammas of Ho-166, SPECT images can be obtained, which give opportunities to guide Ho-166 therapy. Accurate reconstruction of Ho-166 images is currently hampered by photopeak-scatter in the patient, down-scatter in the detector, collimator and patient caused by the 1.4 MeV photons and by bremsstrahlung. We developed and validated a method for quantitative SPECT of Ho-166 that involves correction for both types of scatter plus non-uniform attenuation correction using attenuation maps. Photopeak-scatter (S) is compensated for by a rapid 3D Monte Carlo (MC) method that is incorporated in ordered subset (OS) reconstruction of the emission data, together with simultaneous correction for attenuation (A) and detector response (D); this method is referred to as OS-ADS. Additionally, for correction of down-scatter, we use a 14 keV wide energy window centred at 118 keV (OS-ADSS). Due to a limited number of available energy windows, the same 118 keV energy window was used for down-scatter correction of the simultaneously acquired Gd-153 transmission data. Validations were performed using physical phantom experiments carried out on a dual-head SPECT system; Gd-153 transmission line sources were used for acquiring attenuation maps. For quantitative comparison of OS-ADS and OS-ADSS, bottles filled with Ho-166 were placed in both a cylindrical phantom and an anthropomorphic thorax phantom. Both OS-ADS and OS-ADSS were compared with an ordered subset reconstruction without any scatter correction (OS-AD). Underestimations of about 20% in the attenuation map were reduced to a few per cent after down-scatter correction. The average deviation from the true activity contained in the bottles was +72% with OS-AD. Using OS-ADS, this average overestimation was reduced to +28% and with OS-ADSS the deviation was further reduced to 16%. With OS-AD and OS-ADS, these numbers were more sensitive to the choice of volumes of interest than with OS-ADSS. For the reconstructed activity distributions, erroneous background activity found with OS-AD was reduced by a factor of approximately 2 by applying OS-ADS and reduced by a factor of approximately 4 by applying OS-ADSS. The combined attenuation, photopeak-scatter and down-scatter correction framework proposed here greatly enhanced the quantitative accuracy of Ho-166 imaging, which is of the uppermost importance for image-guided therapies. It is expected that the method, with adapted window settings, also can be applied to other isotopes with high energy peaks that contaminate the photopeak data, such as I-131 or In-111.

Efficient Monte Carlo based scatter artifact reduction in cone-beam micro-CT

Cupping and streak artifacts caused by the detection of scattered photons may severely degrade the quantitative accuracy of cone-beam X-ray computed tomography (CT) images. In order to overcome this problem, we propose and validate the following iterative scatter artifact reduction scheme. First, an initial image is reconstructed from the scatter-contaminated projections. Next, the scatter component of the projections is estimated from the initial reconstruction by a Monte Carlo (MC) simulation. The estimate obtained is then utilized during the reconstruction of a scatter-corrected image. The last two steps are repeated until an adequate correction is obtained. The estimation of the noise-free scatter projections in this scheme is accelerated in the following way: first, a rapid (i.e., based on a low number of simulated photon tracks) MC simulation is executed. The noisy result of this simulation is de-noised by a three-dimensional fitting of Gaussian basis functions. We demonstrate that, compared to plain MC, this method shortens the required simulation time by three to four orders of magnitude. Using simulated projections of a small animal phantom, we show that one cycle of the scatter correction scheme is sufficient to produce reconstructed images that barely differ from the reconstructions of scatter-free projections. The reconstructions of data acquired with a charge-coupled device based micro-CT scanner demonstrate a nearly complete removal of the scatter-induced cupping artifact. Quantitative errors in a water phantom are reduced from around 12% for reconstructions without the scatter correction to 1% after the proposed scatter correction has been applied. In conclusion, a general, accurate, and efficient scatter correction algorithm is developed that requires no mechanical modifications of the scanning equipment and results in only a moderate increase in the total reconstruction time.

Comparison of methods for suppressing edge and aliasing artefacts in iterative x-ray CT reconstruction

X-ray CT images obtained with iterative reconstruction (IR) can be hampered by the so-called edge and aliasing artefacts, which appear as interference patterns and severe overshoots in the areas of sharp intensity transitions. Previously, we have demonstrated that these artefacts are caused by discretization errors during the projection simulation step in IR. Although these errors are inherent to IR, they can be adequately suppressed by reconstruction on an image grid that is finer than that typically used for analytical methods such as filtered back-projection. Two other methods that may prevent edge artefacts are: (i) smoothing the projections prior to reconstruction or (ii) using an image representation different from voxels; spherically symmetric Kaiser-Bessel functions are a frequently employed example of such a representation. In this paper, we compare reconstruction on a fine grid with the two above-mentioned alternative strategies for edge artefact reduction. We show that the use of a fine grid results in a more adequate suppression of artefacts than the smoothing of projections or using the Kaiser-Bessel image representation.

Evaluation of accelerated iterative x-ray CT image reconstruction using floating point graphics hardware

Statistical reconstruction methods offer possibilities to improve image quality as compared with analytical methods, but current reconstruction times prohibit routine application in clinical and micro-CT. In particular, for cone-beam x-ray CT, the use of graphics hardware has been proposed to accelerate the forward and back-projection operations, in order to reduce reconstruction times. In the past, wide application of this texture hardware mapping approach was hampered owing to limited intrinsic accuracy. Recently, however, floating point precision has become available in the latest generation commodity graphics cards. In this paper, we utilize this feature to construct a graphics hardware accelerated version of the ordered subset convex reconstruction algorithm. The aims of this paper are (i) to study the impact of using graphics hardware acceleration for statistical reconstruction on the reconstructed image accuracy and (ii) to measure the speed increase one can obtain by using graphics hardware acceleration. We compare the unaccelerated algorithm with the graphics hardware accelerated version, and for the latter we consider two different interpolation techniques. A simulation study of a micro-CT scanner with a mathematical phantom shows that at almost preserved reconstructed image accuracy, speed-ups of a factor 40 to 222 can be achieved, compared with the unaccelerated algorithm, and depending on the phantom and detector sizes. Reconstruction from physical phantom data reconfirms the usability of the accelerated algorithm for practical cases.

U-SPECT-I: a novel system for submillimeter-resolution tomography with radiolabeled molecules in mice

A major advance in biomedical science and diagnosis was accomplished with the development of in vivo techniques to image radiolabeled molecules, but limited spatial resolution has slowed down applications to small experimental animals. Here, we present a SPECT system (U-SPECT-I) dedicated to radionuclide imaging of murine organs at a submillimeter resolution.

METHODS

The high performance of U-SPECT-I is based on a static triangular detector setup, with a cylindric imaging cavity in the center and 75 gold micropinhole apertures in the cavity wall. The pinholes are focused on a small volume of interest such as the mouse heart or spine to maximize the detection yield of gamma-photons. Three-dimensional molecular distributions are iteratively estimated using the detector data and a statistical reconstruction algorithm that takes into account system blurring and data noise to increase resolution and reduce image noise.

RESULTS

With 0.6-mm-diameter pinholes, the maximum fraction of detected photons emitted by a point source (peak sensitivity) is 0.22% for a 15%-wide energy window and remains higher than 0.12% in the central 12 mm of the central plane. In a resolution phantom, radioactively filled capillaries as small as 0.5 mm and separated by 0.5 mm can be distinguished clearly in reconstructions. Projection data needed for the reconstruction of cross sections of molecular distributions in mouse organs can readily be obtained without the need for any mechanical movements. Images of a mouse spine show 99mTc-hydroxymethylene diphosphonate uptake down to the level of tiny parts of vertebral processes. These are separated clearly from the vertebral and intervertebral foramina. Using another tracer, one can monitor myocardial perfusion in the left and right ventricular walls, even in structures as small as the papillary muscles.

CONCLUSION

U-SPECT-I allows discrimination between molecular concentrations in adjacent volumes of as small as about 0.1 muL, which is significantly smaller than can be imaged by any existing SPECT or PET system. Our initial in vivo images of the mouse heart and spine show that U-SPECT-I can be used for novel applications in the study of dynamic biologic systems with a clear projection to clinical applications. The combination of high resolution and detection efficiency of U-SPECT-I opens up new possibilities for the suborgan-level study of radiotracers in mouse models.

Photon-counting versus an integrating CCD-based gamma camera: important consequences for spatial resolution

Charge-coupled devices (CCDs) coupled to scintillation crystals can be used for high resolution imaging with x-rays and gamma-rays. When the CCD images can be read out fast enough, the energy and interaction position of individual gamma quanta can be estimated by real-time image analysis of scintillation light flashes ('photon counting mode'). We tested a set-up in which an electron-multiplying CCD was coupled to a 1 mm thick columnar CsI crystal by means of a fibre-optic taper. We found that, compared to light integration, photon counting improves the intrinsic spatial resolution by a factor of about 3 to 6. Applying our set-up to Tc-99m and I-125 imaging, we were able to obtain intrinsic resolutions below 60 microm (full width at half maximum). Counting losses due to overlapping of light flashes are negligible for event rates typical for biomedical radio-nuclide imaging and do strongly depend on energy window settings. Energy resolution was estimated to be approximately 35 keV FWHM for a 1:1 taper. We conclude that CCD-based gamma cameras have great potential for applications such as in vivo imaging of gamma emitters.

Parallel statistical image reconstruction for cone-beam x-ray CT on a shared memory computation platform

Statistical reconstruction methods offer possibilities of improving image quality as compared to analytical methods, but current reconstruction times prohibit routine clinical applications. To reduce reconstruction times we have parallelized a statistical reconstruction algorithm for cone-beam x-ray CT, the ordered subset convex algorithm (OSC), and evaluated it on a shared memory computer. Two different parallelization strategies were developed: one that employs parallelism by computing the work for all projections within a subset in parallel, and one that divides the total volume into parts and processes the work for each sub-volume in parallel. Both methods are used to reconstruct a three-dimensional mathematical phantom on two different grid densities. The reconstructed images are binary identical to the result of the serial (non-parallelized) algorithm. The speed-up factor equals approximately 30 when using 32 to 40 processors, and scales almost linearly with the number of cpus for both methods. The huge reduction in computation time allows us to apply statistical reconstruction to clinically relevant studies for the first time.

Evaluation of the ordered subset convex algorithm for cone-beam CT

Statistical methods for image reconstruction such as maximum likelihood expectation maximization (ML-EM) are more robust and flexible than analytical inversion methods and allow for accurate modelling of the photon transport and noise. Statistical reconstruction is prohibitively slow when applied to clinical x-ray cone-beam CT due to the large data sets and the high number of iterations required for reconstructing high resolution images. One way to reduce the reconstruction time is to use ordered subsets of projections during the iterations, which has been successfully applied to fan-beam x-ray CT. In this paper, we quantitatively analyse the use of ordered subsets in concert with the convex algorithm for cone-beam x-ray CT reconstruction, for the case of circular acquisition orbits. We focus on the reconstructed image accuracy of a 3D head phantom. Acceleration factors larger than 300 were obtained with errors smaller than 1%, with the preservation of signal-to-noise ratio. Pushing the acceleration factor towards 600 by using an increasing number of subsets increases the reconstruction error up to 5% and significantly increases noise. The results indicate that the use of ordered subsets can be extremely useful for cone-beam x-ray CT.

Design and simulation of a high-resolution stationary SPECT system for small animals

Exciting new SPECT systems can be created by combining pinhole imaging with compact high-resolution gamma cameras. These new systems are able to solve the problem of the limited sensitivity-resolution trade-off that hampers contemporary small animal SPECT. The design presented here (U-SPECT-III) uses a set of detectors placed in a polygonal configuration and a cylindrical collimator that contains 135 pinholes arranged in nine rings. Each ring contains 15 gold pinhole apertures that focus on the centre of the cylinder. A non-overlapping projection is acquired via each pinhole. Consequently, when a mouse brain is placed in the central field-of-view, each voxel in the cerebrum can be observed via 130 to 135 different pinholes simultaneously. A method for high-resolution scintillation detection is described that eliminates the depth-of-interaction problem encountered with pinhole cameras, and is expected to provide intrinsic detector resolutions better than 150 microm. By means of simulations U-SPECT-III is compared to a simulated dual pinhole SPECT (DP-SPECT) system with a pixelated array consisting of 2.0 x 2.0 mm NaI crystals. Analytic calculations indicate that the proposed U-SPECT-III system yields an almost four times higher linear and about sixty times higher volumetric system resolution than DP-SPECT, when the systems are compared at matching system sensitivity. In addition, it should be possible to achieve a 15 up to 30 times higher sensitivity with U-SPECT-III when the systems are compared at equal resolution. Simulated images of a digital mouse-brain phantom show much more detail with U-SPECT-III than with DP-SPECT. In a resolution phantom, 0.3 mm diameter cold rods are clearly visible with U-SPECT-III, whereas with DP-SPECT the smallest visible rods are about 0.6-0.8 mm. Furthermore, with U-SPECT-III, the image deformations outside the central plane of reconstruction that hamper conventional pinhole SPECT are strongly suppressed. Simulation results indicate that future pinhole SPECT systems are likely to bring about significant improvements in radio-molecular imaging of small animals.

Design and simulation of a high-resolution stationary SPECT system for small animals

Exciting new SPECT systems can be created by combining pinhole imaging with compact high-resolution gamma cameras. These new systems are able to solve the problem of the limited sensitivity-resolution trade-off that hampers contemporary small animal SPECT. The design presented here (U-SPECT-III) uses a set of detectors placed in a polygonal configuration and a cylindrical collimator that contains 135 pinholes arranged in nine rings. Each ring contains 15 gold pinhole apertures that focus on the centre of the cylinder. A non-overlapping projection is acquired via each pinhole. Consequently, when a mouse brain is placed in the central field-of-view, each voxel in the cerebrum can be observed via 130 to 135 different pinholes simultaneously. A method for high-resolution scintillation detection is described that eliminates the depth-of-interaction problem encountered with pinhole cameras, and is expected to provide intrinsic detector resolutions better than 150 microm. By means of simulations U-SPECT-III is compared to a simulated dual pinhole SPECT (DP-SPECT) system with a pixelated array consisting of 2.0 x 2.0 mm NaI crystals. Analytic calculations indicate that the proposed U-SPECT-III system yields an almost four times higher linear and about sixty times higher volumetric system resolution than DP-SPECT, when the systems are compared at matching system sensitivity. In addition, it should be possible to achieve a 15 up to 30 times higher sensitivity with U-SPECT-III when the systems are compared at equal resolution. Simulated images of a digital mouse-brain phantom show much more detail with U-SPECT-III than with DP-SPECT. In a resolution phantom, 0.3 mm diameter cold rods are clearly visible with U-SPECT-III, whereas with DP-SPECT the smallest visible rods are about 0.6-0.8 mm. Furthermore, with U-SPECT-III, the image deformations outside the central plane of reconstruction that hamper conventional pinhole SPECT are strongly suppressed. Simulation results indicate that future pinhole SPECT systems are likely to bring about significant improvements in radio-molecular imaging of small animals.

Characterization and suppression of edge and aliasing artefacts in iterative x-ray CT reconstruction

For the purpose of obtaining x-ray tomographic images, statistical reconstruction (SR) provides a general framework with possible advantages over analytical algorithms such as filtered backprojection (FBP) in terms of flexibility, resolution, contrast and image noise. However, SR images may be seriously affected by some artefacts that are not present in FBP images. These artefacts appear as aliasing patterns and as severe overshoots in the areas of sharp intensity transitions ('edge artefacts'). We characterize this inherent property of iterative reconstructions and hypothesize how discretization errors during reconstruction contribute to the formation of the artefacts. An adequate solution to the problem is to perform the reconstructions on an image grid that is finer than that typically employed for FBP reconstruction, followed by a downsampling of the resulting image to a granularity normally used for display. Furthermore, it is shown that such a procedure is much more effective than post-filtering of the reconstructions. Resulting SR images have superior noise-resolution trade-off compared to FBP, which may facilitate dose reduction during CT examinations.

Experimental validation of a rapid Monte Carlo based micro-CT simulator

We describe a newly developed, accelerated Monte Carlo simulator of a small animal micro-CT scanner. Transmission measurements using aluminium slabs are employed to estimate the spectrum of the x-ray source. The simulator incorporating this spectrum is validated with micro-CT scans of physical water phantoms of various diameters, some containing stainless steel and Teflon rods. Good agreement is found between simulated and real data: normalized error of simulated projections, as compared to the real ones, is typically smaller than 0.05. Also the reconstructions obtained from simulated and real data are found to be similar. Thereafter, effects of scatter are studied using a voxelized software phantom representing a rat body. It is shown that the scatter fraction can reach tens of per cents in specific areas of the body and therefore scatter can significantly affect quantitative accuracy in small animal CT imaging.