In clinical brain SPECT, correction for photon attenuation in the patient is essential to obtain images which provide quantitative information on the regional activity concentration per unit volume ...(kBq.Formula: see text). This correction generally requires an attenuation map (Formula: see text map) denoting the attenuation coefficient at each voxel which is often derived from a CT or MRI scan. However, such an additional scan is not always available and the method may suffer from registration errors. Therefore, we propose a SPECT-only-based strategy for Formula: see text map estimation that we apply to a stationary multi-pinhole clinical SPECT system (G-SPECT-I) for
Tc-HMPAO brain perfusion imaging. The method is based on the use of a convolutional neural network (CNN) and was validated with Monte Carlo simulated scans. Data acquired in list mode was used to employ the energy information of both primary and scattered photons to obtain information about the tissue attenuation as much as possible. Multiple SPECT reconstructions were performed from different energy windows over a large energy range. Locally extracted 4D SPECT patches (three spatial plus one energy dimension) were used as input for the CNN which was trained to predict the attenuation coefficient of the corresponding central voxel of the patch. Results show that Attenuation Correction using the Ground Truth Formula: see text maps (GT-AC) or using the CNN estimated Formula: see text maps (CNN-AC) achieve comparable accuracy. This was confirmed by a visual assessment as well as a quantitative comparison; the mean deviation from the GT-AC when using the CNN-AC is within 1.8% for the standardized uptake values in all brain regions. Therefore, our results indicate that a CNN-based method can be an automatic and accurate tool for SPECT attenuation correction that is independent of attenuation data from other imaging modalities or human interpretations about head contours.
There is an increasing interest in iterative reconstruction (IR) as a key tool to improve quality and increase applicability of x-ray CT imaging. IR has the ability to significantly reduce patient ...dose; it provides the flexibility to reconstruct images from arbitrary x-ray system geometries and allows one to include detailed models of photon transport and detection physics to accurately correct for a wide variety of image degrading effects. This paper reviews discretization issues and modelling of finite spatial resolution, Compton scatter in the scanned object, data noise and the energy spectrum. The widespread implementation of IR with a highly accurate model-based correction, however, still requires significant effort. In addition, new hardware will provide new opportunities and challenges to improve CT with new modelling.
Background
Gamma camera imaging, including single photon emission computed tomography (SPECT), is crucial for research, diagnostics, and radionuclide therapy. Gamma cameras are predominantly based on ...arrays of photon multipliers tubes (PMTs) that read out NaI(Tl) scintillation crystals. In this way, standard gamma cameras can localize ɣ‐rays with energies typically ranging from 30 to 360 keV. In the last decade, there has been an increasing interest towards gamma imaging outside this conventional clinical energy range, for example, for theragnostic applications and preclinical multi‐isotope positron emission tomography (PET) and PET‐SPECT. However, standard gamma cameras are typically equipped with 9.5 mm thick NaI(Tl) crystals which can result in limited sensitivity for these higher energies.
Purpose
Here we investigate to what extent thicker scintillators can improve the photopeak sensitivity for higher energy isotopes while attempting to maintain spatial resolution.
Methods
Using Monte Carlo simulations, we analyzed multiple PMT‐based configurations of gamma detectors with monolithic NaI (Tl) crystals of 20 and 40 mm thickness. Optimized light guide thickness together with 2‐inch round, 3‐inch round, 60 × 60 mm2 square, and 76 × 76 mm2 square PMTs were tested. For each setup, we assessed photopeak sensitivity, energy resolution, spatial, and depth‐of‐interaction (DoI) resolution for conventional (140 keV) and high (511 keV) energy ɣ using a maximum‐likelihood algorithm. These metrics were compared to those of a “standard” 9.5 mm‐thick crystal detector with 3‐inch round PMTs.
Results
Estimated photopeak sensitivities for 511 keV were 27% and 53% for 20 and 40 mm thick scintillators, which is respectively, 2.2 and 4.4 times higher than for 9.5 mm thickness. In most cases, energy resolution benefits from using square PMTs instead of round ones, regardless of their size. Lateral and DoI spatial resolution are best for smaller PMTs (2‐inch round and 60 × 60 mm2 square) which outperform the more cost‐effective larger PMT setups (3‐inch round and 76 × 76 mm2 square), while PMT layout and shape have negligible (< 10%) effect on resolution. Best spatial resolution was obtained with 60 × 60 mm2 PMTs; for 140 keV, lateral resolution was 3.5 mm irrespective of scintillator thickness, improving to 2.8 and 2.9 mm for 511 keV with 20 and 40 mm thick crystals, respectively. Using the 3‐inch round PMTs, lateral resolutions of 4.5 and 3.9 mm for 140 keV and of 3.5 and 3.7 mm for 511 keV were obtained with 20 and 40 mm thick crystals respectively, indicating a moderate performance degradation compared to the 3.5 and 2.9 mm resolution obtained by the standard detector for 140 and 511 keV. Additionally, DoI resolution for 511 keV was 7.0 and 5.6 mm with 20 and 40 mm crystals using 60 × 60 mm2 square PMTs, while with 3‐inch round PMTs 12.1 and 5.9 mm were obtained.
Conclusion
Depending on PMT size and shape, the use of thicker scintillator crystals can substantially improve detector sensitivity at high gamma energies, while spatial resolution is slightly improved or mildly degraded compared to standard crystals.
Detection of scattered gamma quanta degrades image contrast and quantitative accuracy of single-photon emission computed tomography (SPECT) imaging. This paper reviews methods to characterize and ...model scatter in SPECT and correct for its image degrading effects, both for clinical and small animal SPECT. Traditionally scatter correction methods were limited in accuracy, noise properties and/or generality and were not very widely applied. For small animal SPECT, these approximate methods of correction are often sufficient since the fraction of detected scattered photons is small. This contrasts with patient imaging where better accuracy can lead to significant improvement of image quality. As a result, over the last two decades, several new and improved scatter correction methods have been developed, although often at the cost of increased complexity and computation time. In concert with (i) the increasing number of energy windows on modern SPECT systems and (ii) excellent attenuation maps provided in SPECT/CT, some of these methods give new opportunities to remove degrading effects of scatter in both standard and complex situations and therefore are a gateway to highly quantitative single- and multi-tracer molecular imaging with improved noise properties. Widespread implementation of such scatter correction methods, however, still requires significant effort.
Pinhole collimation is widely recognized for offering superior resolution-sensitivity trade-off in SPECT imaging of small subjects. The newly developed EXIRAD-3D autoradiography technique (MILabs ...B.V.) based on a highly focusing multi-pinhole collimator achieves micron-resolution SPECT for cryo-cooled tissue samples. For such high resolutions, the choice of pinhole material may have a significant impact on images. Therefore, this paper aims to compare the performance of EXIRAD-3D with lead, tungsten, gold, and depleted uranium pinhole collimators designed such that they achieve equal sensitivities. Performance in terms of resolution is characterized for several radioisotopes, namely
In (171 keV and 245 keV),
Tc (140 keV),
Tl (71 keV), and
I (27 keV). Using Monte Carlo simulation, point spread functions were generated and their profiles as well as their full-width-at-half-maximum and full-width-at-tenth-maximum were determined and evaluated for different materials and isotopes. Additionally, simulated reconstructions of a Derenzo resolution phantom, validated with experimental data, were judged by assessment of the resolvable rods as well as a contrast-to-noise ratio (CNR) analysis. Our results indicate that using materials with higher photon-stopping power yields images with better CNR for the studied isotopes with improvements ranging from 1.9% to 36.6%. Visual assessment on the reconstructed images suggests that for EXIRAD-3D, the tungsten collimator is generally a good choice for a wide range of SPECT isotopes. For relatively high-energy isotopes such as
In, using gold inserts can be beneficial.
Block-iterative image reconstruction methods, such as ordered subset expectation maximization (OSEM), are commonly used to accelerate image reconstruction. In OSEM, the speed-up factor over maximum ...likelihood expectation maximization (MLEM) is approximately equal to the number of subsets in which the projection data are divided. Traditionally, each subset consists of a couple of projection views, and the more subsets are used, the more the solution deviates from MLEM solutions. We found for multi-pinhole single photon emission computed tomography (SPECT) that even moderate acceleration factors in OSEM lead to inaccurate reconstructions. Therefore, we introduce pixel-based ordered subset expectation maximization (POSEM), which is based on an alternative subset choice. Pixels in each subset are spread out regularly over projections and are spatially separated as much as possible. We validated POSEM for data acquired with a focusing multi-pinhole SPECT system. Performance was compared with traditional OSEM and MLEM for a rat total body bone scan, a gated mouse myocardial perfusion scan and a Defrise phantom scan. We found that POSEM can be operated at acceleration factors that are often an order of magnitude higher than in traditional OSEM.
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.
Microscopic nuclear imaging down to spatial resolutions of a few hundred microns can already be achieved using low-energy gamma emitters (e.g.
I, ∼30 keV) and a basic single micro-pinhole gamma ...camera. This has been applied to
mouse thyroid imaging, for example. For clinically used radionuclides such as
Tc, this approach fails due to penetration of the higher-energy gamma photons through the pinhole edges. To overcome these resolution degradation effects, we propose a new imaging approach: scanning focus nuclear microscopy (SFNM). We assess SFNM using Monte Carlo simulations for clinically used isotopes. SFNM is based on the use of a 2D scanning stage with a focused multi-pinhole collimator containing 42 pinholes with narrow pinhole aperture opening angles to reduce photon penetration. All projections of different positions are used to iteratively reconstruct a three-dimensional image from which synthetic planar images are generated. SFNM imaging was tested using a digital Derenzo resolution phantom and a mouse ankle joint phantom containing
Tc (140 keV). The planar images were compared with those obtained using a single-pinhole collimator, either with matched pinhole diameter or with matched sensitivity. The simulation results showed an achievable
Tc image resolution of 0.04 mm and detailed
Tc bone images of a mouse ankle with SFNM. SFNM has strong advantages over single-pinhole imaging in terms of spatial resolution.
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.
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.
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.
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.