Search

Purpose: We investigated the temporal variation of the induced magnetic field due to the transverse and the longitudinal gradient coils in tungsten collimators arranged in hexagonal and pentagonal geometries with and without gaps between the collimators., Methods: We modeled x-, y-, and z-gradient coils and different arrangements of single-photon emission computed tomography (SPECT) collimators using FEKO, a three-dimensional electromagnetic simulation tool. A time analysis approach was used to generate the pulsed magnetic field gradient. The approach was validated with measurements using a 7T MRI scanner., Results: Simulations showed an induced magnetic field representing 4.66% and 0.87% of the applied gradient field (gradient strength = 500 mT/m) for longitudinal and transverse gradient coils, respectively. These values can be reduced by 75% by adding gaps between the collimators for the pentagonal arrangement, bringing the maximum induced magnetic field to less than 2% of the applied gradient for all of the gradient coils., Conclusion: Characterization of the maximum induced magnetic field shows that by adding gaps between the collimators for an integrated SPECT/MRI system, eddy currents can be corrected by the MRI system to avoid artifact. The numerical model was validated and was proposed as a tool for studying the effect of a SPECT collimator within the MRI gradient coils., (© 2014 Wiley Periodicals, Inc.)

Parallel-hole SPECT collimators have traditionally been manufactured by stacking sheets of lead foil or by casting. These techniques significantly restrict our options in terms of collimator geometry. However, recent developments in metal additive manufacturing are making novel collimator designs possible, giving rise to new opportunities in SPECT imaging. In this paper we propose an innovative type of collimator for stationary SPECT, using parallel-holes whose collimation direction depends on their axial position. Its main advantage compared to current stationary SPECT systems (which are based on pinholes) is that, using only axial bed translations, we can achieve complete angular sampling of an increased portion of the transaxial area of the collimator bore. This allows the system to be much more compact than current stationary SPECT systems that image objects of the same size. We describe three possible designs, for full-body, brain and small-animal imaging, respectively, and test their feasibility using simulations. The system modeling method is validated against realistic Monte Carlo simulations, and then used in the evaluation of the systems' performances and reconstructions. The simulations show that the system is able to reconstruct objects occupying the predicted field of view ([Formula: see text] of the transaxial area of the bore) without sampling artifacts. In particular, we perform reconstructions from noisy projection data obtained for an activity and scanning time similar to standard protocols for the three applications, and the resulting images indicate the possibility of using the proposed systems in practice.

Purpose: Brain single photon emission computed tomography (SPECT) imaging is an important clinical tool, with unique tracers for studying neurological diseases. Nowadays, most commercial SPECT systems are combined with x-ray computed tomography (CT) in so-called SPECT/CT systems to obtain an anatomical background for the functional information. However, while CT images have a high spatial resolution, they have a low soft-tissue contrast, which is an important disadvantage for brain imaging. Magnetic resonance imaging (MRI), on the other hand, has a very high soft-tissue contrast and does not involve extra ionizing radiation. Therefore, the authors designed a brain SPECT insert that can operate inside a clinical MRI., Methods: The authors designed and simulated a compact stationary multipinhole SPECT insert based on digital silicon photomultiplier detector modules, which have shown to be MR-compatible and have an excellent intrinsic resolution (0.5 mm) when combined with a monolithic 2 mm thick LYSO crystal. First, the authors optimized the different parameters of the SPECT system to maximize sensitivity for a given target resolution of 7.2 mm in the center of the field-of-view, given the spatial constraints of the MR system. Second, the authors performed noiseless simulations of two multipinhole configurations to evaluate sampling and reconstructed resolution. Finally, the authors performed Monte Carlo simulations and compared the SPECT insert with a clinical system with ultrahigh-resolution (UHR) fan beam collimators, based on contrast-to-noise ratio and a visual comparison of a Hoffman phantom with a 9 mm cold lesion., Results: The optimization resulted in a stationary multipinhole system with a collimator radius of 150.2 mm and a detector radius of 172.67 mm, which corresponds to four rings of 34 diSPM detector modules. This allows the authors to include eight rings of 24 pinholes, which results in a system volume sensitivity of 395 cps/MBq. Noiseless simulations show sufficient axial sampling (in a Defrise phantom) and a reconstructed resolution of 5.0 mm (in a cold-rod phantom). The authors compared the 24-pinhole setup with a 34-pinhole system (with the same detector radius but a collimator radius of 156.63 mm) and found that 34 pinholes result in better uniformity but a worse reconstruction of the cold-rod phantom. The authors also compared the 24-pinhole system with a clinical triple-head UHR fan beam system based on contrast-to-noise ratio and found that the 24-pinhole setup performs better for the 6 mm hot and the 16 mm cold lesions and worse for the 8 and 10 mm hot lesions. Finally, the authors reconstructed noisy projection data of a Hoffman phantom with a 9 mm cold lesion and found that the lesion was slightly better visible on the multipinhole image compared to the fan beam image., Conclusions: The authors have optimized a stationary multipinhole SPECT insert for MRI and showed the feasibility of doing brain SPECT imaging inside a MRI with an image quality similar to the best clinical SPECT systems available.

The accurate determination of the local impulse response and the covariance in voxels from penalized maximum likelihood reconstructed images requires performing reconstructions from many noise realizations of the projection data. As this is usually a very time-consuming process, efficient analytical approximations based on the Fisher information matrix (FIM) have been extensively used in PET and SPECT to estimate these quantities. For 3D imaging, however, additional approximations need to be made to the FIM in order to speed up the calculations. The most common approach is to use the local shift-invariant (LSI) approximation of the FIM, but this assumes specific conditions which are not always necessarily valid. In this paper we take a single-pinhole SPECT system and compare the accuracy of the LSI approximation against two other methods that have been more recently put forward: the non-uniform object-space pixelation (NUOP) and the subsampled FIM. These methods do not assume such restrictive conditions while still increasing the speed of the calculations considerably. Our results indicate that in pinhole SPECT the NUOP and subsampled FIM approaches could be more reliable than the LSI approximation, especially when a high accuracy is required.

In single photon emission computed tomography, the choice of the collimator has a major impact on the sensitivity and resolution of the system. Traditional parallel-hole and fan-beam collimators used in clinical practice, for example, have a relatively poor sensitivity and subcentimeter spatial resolution, while in small-animal imaging, pinhole collimators are used to obtain submillimeter resolution and multiple pinholes are often combined to increase sensitivity. This paper reviews methods for production, sensitivity maximization, and task-based optimization of collimation for both clinical and preclinical imaging applications. New opportunities for improved collimation are now arising primarily because of (i) new collimator-production techniques and (ii) detectors with improved intrinsic spatial resolution that have recently become available. These new technologies are expected to impact the design of collimators in the future. The authors also discuss concepts like septal penetration, high-resolution applications, multiplexing, sampling completeness, and adaptive systems, and the authors conclude with an example of an optimization study for a parallel-hole, fan-beam, cone-beam, and multiple-pinhole collimator for different applications.

Multi-pinhole collimators are often used in pre-clinical SPECT systems because they have a better resolution-sensitivity tradeoff than parallel hole collimators when imaging small objects. Most multi-pinhole collimators are designed to allow no or only a limited amount of overlap between the different pinhole projections because the ambiguity introduced by multiplexing pinholes can result in artifacts. The origin of these artifacts is still not fully understood, but previous research has already shown that data incompleteness could be part of the explanation. Therefore, we developed a method to investigate data completeness in multiplexing multi-pinhole systems and showed that a certain activity distribution can be successfully reconstructed when the nonmultiplexed data is complete or when the overlap can be sufficiently de-multiplexed. We validated this method using computer simulated phantom data of different multiplexing systems. We also studied contrast-to-noise and nonprewhitening matched filter signal-to-noise ratio (NPW-SNR) to compare the image quality in a single pinhole system with multiplexing systems. We found that our method can indeed be used to evaluate data completeness in multiplexing systems and found no artifacts in the systems that had complete data. Sensitivity increased significantly with multiplexing but we found only small, nonsignificant differences in contrast-to-noise ratio. However, the NPW-SNR did slightly improve in the multiplexing setups. We conclude that more multiplexing does not necessarily result in more artifacts and that even a high amount of multiplexing can still result in artifact-free images if the nonmultiplexed data is complete or when the overlap can be sufficiently de-multiplexed.

MicroSPECT is one of the main functional imaging techniques used in the preclinical setting. Even though high-resolution images can be obtained with currently available systems, their sensitivity is often quite low due to the use of multi-pinhole collimation. This results in long acquisition times and hampers dynamic imaging. However, it has already been shown that this limited sensitivity can be overcome using high-resolution detectors. In this article, we therefore investigated the use of a digital photon counter (DPC) in combination with a 2 mm thick monolithic LYSO crystal for SPECT imaging. These light sensors contain arrays of avalanche photodiodes whose signals are directly digitised. The DPCs have the advantage that they are very compact, have a high intrinsic resolution, are MR compatible and allow disabling cells with a high dark count rate. In order to investigate the influence of the temperature dependent dark count rate on the detector performance, we compared it at 3 °C and 18 °C. At 3 °C, we observed an energy resolution of 28.8% and an intrinsic spatial resolution of 0.48 mm. Furthermore, the count rate at 10% loss is 60 kcps. Next, we looked at the event loss at 18 °C caused by the higher dark count rate and found a 5% loss compared to the 3 °C measurements. At this higher temperature the energy resolution becomes 29.2% and the intrinsic spatial resolution decreases to 0.52 mm. Due to the 5% count loss, the count rate at 10% loss increases to 63 kcps. A small degradation of the detector performance is thus observed at 18 °C.These results show the usefulness of this detector for SPECT imaging together with its excellent intrinsic spatial resolution. A drawback of the detector is its low, spatially varying energy resolution. Even though the detection efficiency and intrinsic spatial resolution are better at 3 °C, results are still acceptable at 18 °C.

Preclinical imaging with SPECT combined with CT or MRI is used more and more frequently and has proven to be very useful in translational research. In this article, an overview of current preclinical research applications and trends of SPECT combined with CT or MRI, mainly in tumour imaging and neuroscience imaging, is given and the advantages and disadvantages of the different approaches are described. Today SPECT and CT systems are often integrated into a single device (commonly called a SPECT/CT system), whereas at present combined SPECT and MRI is almost always carried out with separate systems and fiducial markers to combine the separately acquired images. While preclinical SPECT/CT is most widely applied in oncology research, SPECT combined with MRI (SPECT/MRI when integrated in one system) offers the potential for both neuroscience applications and oncological applications. Today CT and MRI are still mainly used to localize radiotracer binding and to improve SPECT quantification, although both CT and MRI have additional potential. Future technology developments may include fast sequential or simultaneous acquisition of (dynamic) multimodality data, spectroscopy, fMRI along with high-resolution anatomic MRI, advanced CT procedures, and combinations of more than two modalities such as combinations of SPECT, PET, MRI and CT all together. This will all strongly depend on new technologies. With further advances in biology and chemistry for imaging molecular targets and (patho)physiological processes in vivo, the introduction of new imaging procedures and promising new radiopharmaceuticals in clinical practice may be accelerated.

Single-photon emission computed tomography (SPECT) detectors with improved spatial resolution can be used to build multi-pinhole SPECT systems that have a higher sensitivity or a higher spatial resolution. In order to improve the spatial resolution we investigate the performance of a 2 mm thick continuous Lutetium Yttrium Orthosilicate (LYSO) scintillator and compare it to the performance of a 5 mm thick continuous NaI(Tl) scintillator. The advantages of LYSO are its high stopping power and its non-hygroscopicity. Drawbacks are the lower light output and the intrinsic radioactivity. The hypothesis of this study is that such a thin LYSO scintillator will have a small light spread and, as a consequence, will also have an improved spatial resolution when coupled to a Hamamatsu H8500 position sensitive photomultiplier tube. To optimize the spatial resolution and the useful detector area we used a mean nearest neighbor event-positioning method. Beam source measurements ((99m)Tc, 140 keV) were done to investigate the energy resolution and the spatial resolution of both detectors. The effect of the intrinsic radioactivity of the LYSO scintillator in the energy window was quantified. The mean energy resolution is 9.3% for the NaI(Tl) scintillator and 21.3% for the LYSO scintillator. The LYSO spectrum shows an X-ray escape peak which decreases the detection efficiency with 9.1%. The spatial resolution of the LYSO detector (0.93 mm full width at half maximum (FWHM)) is superior to the spatial resolution of the NaI(Tl) detector (1.37 mm FWHM). The intrinsic radioactivity in the energy window (42% window centered at 140 keV) is low (125.6 cps, 0.024 cps mm(-3)). LYSO is a promising scintillator for small-animal SPECT imaging, where spatial resolution is more important than energy resolution.

This work aimed to measure the in vivo quantification errors obtained when ray-based iterative reconstruction is used in micro-single-photon emission computed tomography (SPECT). This was investigated with an extensive phantom-based evaluation and two typical in vivo studies using 99mTc and 111In, measured on a commercially available cadmium zinc telluride (CZT)-based small-animal scanner. Iterative reconstruction was implemented on the GPU using ray tracing, including (1) scatter correction, (2) computed tomography-based attenuation correction, (3) resolution recovery, and (4) edge-preserving smoothing. It was validated using a National Electrical Manufacturers Association (NEMA) phantom. The in vivo quantification error was determined for two radiotracers: [99mTc]DMSA in naive mice (n  =  10 kidneys) and [111In]octreotide in mice (n  =  6) inoculated with a xenograft neuroendocrine tumor (NCI-H727). The measured energy resolution is 5.3% for 140.51 keV (99mTc), 4.8% for 171.30 keV, and 3.3% for 245.39 keV (111In). For 99mTc, an uncorrected quantification error of 28 ± 3% is reduced to 8 ± 3%. For 111In, the error reduces from 26 ± 14% to 6 ± 22%. The in vivo error obtained with 99mTc-dimercaptosuccinic acid ([99mTc]DMSA) is reduced from 16.2 ± 2.8% to -0.3 ± 2.1% and from 16.7 ± 10.1% to 2.2 ± 10.6% with [111In]octreotide. Absolute quantitative in vivo SPECT is possible without explicit system matrix measurements. An absolute in vivo quantification error smaller than 5% was achieved and exemplified for both [99mTc]DMSA and [111In]octreotide.

Purpose: Over the last ten years, there has been an extensive growth in the development of microSPECT imagers. Most of the systems are based on the combination of conventional, relatively large gamma cameras with poor intrinsic spatial resolution and multipinhole collimators working in large magnification mode. Spatial resolutions range from 0.58 to 0.76 mm while peak sensitivities vary from 0.06% to 0.4%. While pushing the limits of performance is of major importance, the authors believe that there is a need for smaller and less complex systems that bring along a reduced cost. While low footprint and low-cost systems can make microSPECT available to more researchers, the ease of operation and calibration and low maintenance cost are additional factors that can facilitate the use of microSPECT in molecular imaging. In this paper, the authors simulate the performance of a microSPECT imager that combines high space-bandwidth detectors and pinholes with truncated projection, resulting in a small and stationary system., Methods: A system optimization algorithm is used to determine the optimal SPECT systems, given our high resolutions detectors and a fixed field-of-view. These optimal system geometries are then used to simulate a Defrise disk phantom and a hot rod phantom. Finally, a MOBY mouse phantom, with realistic concentrations of Tc99m-tetrofosmin is simulated., Results: Results show that the authors can successfully reconstruct a Defrise disk phantom of 24 mm in diameter without any rotating system components or translation of the object. Reconstructed spatial resolution is approximately 800 μm while the peak sensitivity is 0.23%. Finally, the simulation of the MOBY mouse phantom shows that the authors can accurately reconstruct mouse images., Conclusions: These results show that pinholes with truncated projections can be used in small magnification or minification mode to obtain a compact and stationary microSPECT system. The authors showed that they can reach state-of-the-art system performance and can successfully reconstruct images with realistic noise levels in a preclinical context. Such a system can be useful for dynamic SPECT imaging.

Currently, clinical brain single photon emission computed tomography (SPECT) is mostly performed using rotating dual-head gamma cameras equipped with low-energy-high-resolution parallel-beam collimators (LEHR PAR). The resolution of these systems is rather poor (8-10 mm) and the rotation of the heavy gamma cameras can introduce misalignment errors. Therefore, we designed a static full-ring multi-lofthole brain SPECT insert for an existing ring of LaBr3 (5% Ce) detectors. The novelty of the design is found in the shutter mechanism that makes the system very flexible and eliminates the need for rotating parts. A stationary SPECT insert is not only more robust, it is also easier to integrate in a magnetic resonance imaging system (MRI) for simultaneous SPECT-MRI. The target spatial resolution of our design is 6 mm. In this study we used analytical calculations to optimize the collimator for an existing ring of LaBr3 (5% Ce) detectors. We fixed the target spatial resolution at 6 mm in the center of the field-of-view and maximized the volume sensitivity by changing the collimator radius, the aperture and the number of loftholes. Based on these optimal parameters we simulated phantom data and evaluated the image quality of our multi-lofthole system. We simulated a noiseless uniform and Defrise phantom to assess artifacts and sampling completeness and a noiseless hot-rod phantom to assess the reconstructed spatial resolution. We visually evaluated a simulated noisy Hoffman phantom with two lesions. Then, we evaluated the non-prewhitening matched filter signal-to-noise ratio (NPW-SNR) in two lesion detectability phantoms: one with hot lesions and one with cold lesions. Finally, a contrast-to-noise (CNR) study was performed on a phantom with both hot and cold lesions of different sizes (6-16 mm). All results were compared to a LEHR PAR system. The optimization resulted in a final collimator design with a volume sensitivity of 1.55 × 10(-4) cps Bq(-1), which is 2.5 times lower than the sensitivity of a dual-head system with LEHR PAR collimators. Spatial resolution, on the other hand, has clearly improved compared to LEHR PAR: with the multi-lofthole system we successfully reconstructed 4 mm hot rods. Although this improved resolution did not result in an unambiguous improvement in CNR or NPW-SNR, we believe that the flexibility of the shutter mechanism opens interesting perspectives toward time-multiplexing and integration with MRI.

Monolithic scintillation detectors for positron emission tomography and single-photon emission computed tomography (SPECT) imaging have many advantages over pixelated detectors. The use of monolithic crystals allows for reducing the scintillator cost per unit volume and increasing the sensitivity along with the energy and timing resolution of the detector. In addition, on thick detectors the depth-of-interaction can be determined without additional hardware. However, costly and complex calibration procedures have been proposed to achieve optimal detector performance for monolithic detectors. This hampers their use in commercial systems. There is thus, a need for simple calibration routines that can be performed on assembled systems. The main goal of this work is to develop a simplified calibration procedure based on acquired training data. In comparison with other methods that use training data acquired with beam sources attached to robotic stages, the proposed method uses a static un-collimated activity source with simple geometry acquiring in a reasonable time. Once the data are acquired, the calibration of the detector is accomplished in three steps: energy calibration based on the k-means clustering method, self-organization based on the self-organizing maps algorithm, and distortion correction based on the Monge-Kantorovich grid adaptation. The proposed calibration method was validated for 2D positioning using a SPECT detector. Similar results were obtained by comparison with an existing calibration method (maximum likelihood estimation). In conclusion, we proposed a novel calibration method for monolithic scintillation detectors that greatly simplifies their use with optimal performance in SPECT systems.

Purpose: We compared the performance of three commercial small-animal μSPECT scanners equipped with multipinhole general purpose (GP) and multipinhole high-resolution (HR) collimators designed for imaging mice., Methods: Spatial resolution, image uniformity, point source sensitivity and contrast recovery were determined for the U-SPECT-II (MILabs), the NanoSPECT-NSO (BioScan) and the X-SPECT (GE) scanners. The pinhole diameters of the HR collimator were 0.35 mm, 0.6 mm and 0.5 mm for these three systems respectively. A pinhole diameter of 1 mm was used for the GP collimator. To cover a broad field of imaging applications three isotopes were used with various photon energies: (99m)Tc (140 keV), (111)In (171 and 245 keV) and (125)I (27 keV). Spatial resolution and reconstructed image uniformity were evaluated in both HR and a GP mode with hot rod phantoms, line sources and a uniform phantom. Point source sensitivity and contrast recovery measures were additionally obtained in the GP mode with a novel contrast recovery phantom developed in-house containing hot and cold submillimetre capillaries on a warm background., Results: In hot rod phantom images, capillaries as small as 0.4 mm with the U-SPECT-II, 0.75 mm with the X-SPECT and 0.6 mm with the NanoSPECT-NSO could be resolved with the HR collimators for (99m)Tc. The NanoSPECT-NSO achieved this resolution in a smaller field-of-view (FOV) and line source measurements showed that this device had a lower axial than transaxial resolution. For all systems, the degradation in image resolution was only minor when acquiring the more challenging isotopes (111)In and (125)I. The point source sensitivity with (99m)Tc and GP collimators was 3,984 cps/MBq for the U-SPECT-II, 620 cps/MBq for the X-SPECT and 751 cps/MBq for the NanoSPECT-NSO. The effects of volume sensitivity over a larger object were evaluated by measuring the contrast recovery phantom in a realistic FOV and acquisition time. For 1.5-mm rods at a noise level of 8 %, the contrast recovery coefficient (CRC) was 42 %, 37 % and 34 % for the U-SPECT-II, X-SPECT and NanoSPECT-NSO, respectively. At maximal noise levels of 10 %, a CRCcold of 70 %, 52 % and 42 % were obtained for the U-SPECT-II, X-SPECT and NanoSPECT-NSO, respectively. When acquiring (99m)Tc with the GP collimators, the integral/differential uniformity values were 30 %/14 % for the U-SPECT-II, 50 %/30 % for the X-SPECT and 38 %/25 % for the NanoSPECT-NSO. When using the HR collimators, these uniformity values remained similar for U-SPECT-II and X-SPECT, but not for the NanoSPECT-NSO for which the uniformity deteriorated with larger volumes., Conclusion: We compared three μSPECT systems by acquiring and analysing mouse-sized phantoms including a contrast recovery phantom built in-house offering the ability to measure the hot contrast on a warm background in the submillimetre resolution range. We believe our evaluation addressed the differences in imaging potential for each system to realistically image tracer distributions in mouse-sized objects.

In single-photon emission computed tomography (SPECT), multi-pinhole collimation is often employed nowadays. Most multi-pinhole collimators avoid overlap (multiplexing) of the projections on the detector. This can be done by using additional shielding or by spacing the pinholes far enough apart. Using additional shielding has the drawback that it increases weight, design complexity and cost. Spacing the pinholes far enough apart results in sub-optimal detector usage, the valuable detector area is not entirely used. This is due to the circular projections of pinholes on the detector; these ellipses can not be tiled with high detector coverage. To overcome this we designed a new pinhole geometry, the lofthole, that has a rectangular projection on the detector. The lofthole has a circular aperture and a rectangular entrance/exit opening. Sensitivity formulae have been derived for pinholes and loftholes. These formulae take the penumbra effect into account; the proposed formulae do not take penetration into account. The derived formulae are valid for geometries where the field-of-view and the sensitivity of the aperture are solely limited by the exit window. A flood map measurement was performed to compare the rectangular projection of a lofthole with the circular projection of a pinhole. Finally, measurements were done to compare the amount of penetration of pinholes with the amount of penetration of a lofthole. A square lofthole collimator has less penetration than a knife-edge pinhole collimator that irradiates the same rectangular detector area with full coverage. A multi-lofthole collimator allows high detector coverage without using additional shielding. An additional advantage is the lower amount of penetration.

Parallel and converging hole collimators are most frequently used in nuclear medicine. Less common is the use of rotating slat collimators for single photon emission computed tomography (SPECT). The higher photon collection efficiency, inherent to the geometry of rotating slat collimators, results in much lower noise in the data. However, plane integrals contain spatial information in only one direction, whereas line integrals provide two-dimensional information. It is not a trivial question whether the initial gain in efficiency will compensate for the lower information content in the plane integrals. Therefore, a comparison of the performance of parallel hole and rotating slat collimation is needed. This study compares SPECT with rotating slat and parallel hole collimation in combination with MLEM reconstruction with accurate system modeling and correction for scatter and attenuation. A contrast-to-noise study revealed an improvement of a factor 2-3 for hot lesions and more than a factor of 4 for cold lesion. Furthermore, a clinically relevant case of heart lesion detection is simulated for rotating slat and parallel hole collimators. In this case, rotating slat collimators outperform the traditional parallel hole collimators. We conclude that rotating slat collimators are a valuable alternative for parallel hole collimators.

Today, new single photon emission computed tomography (SPECT) reconstruction techniques rely on accurate Monte Carlo (MC) simulations to optimize reconstructed images. However, existing MC scintillation camera models which usually include an accurate description of the collimator and crystal, lack correct implementation of the gamma camera's back compartments. In the case of dual isotope simultaneous acquisition (DISA), where backscattered photons from the highest energy isotope are detected in the imaging energy window of the second isotope, this approximation may induce simulation errors. Here, we investigate the influence of backscatter compartment modelling on the simulation accuracy of high-energy isotopes. Three models of a scintillation camera were simulated: a simple model (SM), composed only of a collimator and a NaI(Tl) crystal; an intermediate model (IM), adding a simplified description of the backscatter compartments to the previous model and a complete model (CM), accurately simulating the materials and geometries of the camera. The camera models were evaluated with point sources ((67)Ga, (99m)Tc, (111)In, (123)I, (131)I and (18)F) in air without a collimator, in air with a collimator and in water with a collimator. In the latter case, sensitivities and point-spread functions (PSFs) simulated in the photopeak window with the IM and CM are close to the measured values (error below 10.5%). In the backscatter energy window, however, the IM and CM overestimate the FWHM of the detected PSF by 52% and 23%, respectively, while the SM underestimates it by 34%. The backscatter peak fluence is also overestimated by 20% and 10% with the IM and CM, respectively, whereas it is underestimated by 60% with the SM. The results show that an accurate description of the backscatter compartments is required for SPECT simulations of high-energy isotopes (above 300 keV) when the backscatter energy window is of interest.

Purpose: Multiple investigators have recently reported the use of yttrium-90 (90Y) bremsstrahlung single photon emission computed tomography (SPECT) imaging for the dosimetry of targeted radionuclide therapies. Because Monte Carlo (MC) simulations are useful for studying SPECT imaging, this study investigates the MC simulation of 90Y bremsstrahlung photons in SPECT. To overcome the computationally expensive simulation of electrons, the authors propose a fast way to simulate the emission of 90Y bremsstrahlung photons based on prerecorded bremsstrahlung photon probability density functions (PDFs)., Methods: The accuracy of bremsstrahlung photon simulation is evaluated in two steps. First, the validity of the fast bremsstrahlung photon generator is checked. To that end, fast and analog simulations of photons emitted from a 90Y point source in a water phantom are compared. The same setup is then used to verify the accuracy of the bremsstrahlung photon simulations, comparing the results obtained with PDFs generated from both simulated and measured data to measurements. In both cases, the energy spectra and point spread functions of the photons detected in a scintillation camera are used., Results: Results show that the fast simulation method is responsible for a 5% overestimation of the low-energy fluence (below 75 keV) of the bremsstrahlung photons detected using a scintillation camera. The spatial distribution of the detected photons is, however, accurately reproduced with the fast method and a computational acceleration of approximately 17-fold is achieved. When measured PDFs are used in the simulations, the simulated energy spectrum of photons emitted from a point source of 90Y in a water phantom and detected in a scintillation camera closely approximates the measured spectrum. The PSF of the photons imaged in the 50-300 keV energy window is also accurately estimated with a 12.4% underestimation of the full width at half maximum and 4.5% underestimation of the full width at tenth maximum., Conclusions: Despite its limited accuracy, the fast bremsstrahlung photon generator is well suited for the simulation of bremsstrahlung photons emitted in large homogeneous organs, such as the liver, and detected in a scintillation camera. The computational acceleration makes it very useful for future investigations of 90Y bremsstrahlung SPECT imaging.

Purpose: Quantitative SPECT imaging is limited by many factors, including penetration of high energy photons through the collimator. Even small peaks of high energy yield extensive contamination in the photopeak energy window. Rotating slat collimators typically offer a 40- to 50-fold increase in geometric detection efficiency compared to parallel hole collimators, while the amount of high-energy contamination is expected to remain the same. The aim of this study is to show that SPECT for isotopes with high energy peaks benefits from the use of a rotating slat collimator by significantly reducing the effect of high energy contamination., Methods: Monte Carlo simulations were performed for I-123 and I-131. The simulation of a single point source in a scattering medium was used to investigate the history of each detected photon both for the rotating slat and the parallel hole collimator. Next, realistic simulations of an image quality phantom in combination with iterative image reconstruction enabled us to study the gain in image quality which can be expected for the rotating slat collimator., Results: The results show that the rotating slat collimator is able to reduce the contribution of high energy photons from about 60% to about 8% of all the detected photons in the photopeak energy window. This results in a recovered contrast gain of about 25%. For I-131, there was a reduction in contamination from 65% to 49% when using a rotating slat instead of a parallel hole collimator. This resulted in an average hot contrast improvement of 12% and average cold contrast improvement of 8% for tomographic images., Conclusions: Rotating slat collimators have an additional advantage regarding image quality for isotopes with high energy emissions, especially if the contamination arises from photons with energies higher than the main emission energy, due to the relatively lower number of measured collimator-penetrating photons with respect to geometrically collimated photons.

The GEANT4 application for tomographic emission (GATE) is one of the most detailed Monte Carlo simulation tools for SPECT and PET. It allows for realistic phantoms, complex decay schemes, and a large variety of detector geometries. However, only a fraction of the information in each particle history is available for postprocessing. In order to extend the analysis capabilities of GATE, a flexible framework was developed. This framework allows all detected events to be subdivided according to their type: In PET, true coincidences from others, and in SPECT, geometrically collimated photons from others. The framework of the authors can be applied to any isotope, phantom, and detector geometry available in GATE. It is designed to enhance the usability of GATE for the study of contamination and for the investigation of the properties of current and future prototype detectors. The authors apply the framework to a case study of Bexxar, first assuming labeling with 124I, then with 131I. It is shown that with 124I PET, results with an optimized window improve upon those with the standard window but achieve less than half of the ideal improvement. Nevertheless, 124I PET shows improved resolution compared to 131I SPECT with triple-energy-window scatter correction.

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.