Accurate Monte Carlo modelling of the back compartments of SPECT cameras

Monte Carlo simulation of SPECT characterization for 223 Ra post-injection scintigraphy

OBJECTIVES

223 Ra is a promising α-emitting radionuclide for prostate cancer metastasis palliative treatment. Post-injection scintigraphy is of major importance to verify the concentration of the radiopharmaceutical in the targeted sites. Given the low activity administered to patients, the choice of acquisition parameters, including the collimator type, the energy window's width and the photopeak energy to be used, is primordial for the image quality. The purpose of our work was to select the SPECT configuration suitable for 223 Ra post-injection scintigraphy.

METHODS

We conducted simulation studies with a Symbia T6 Siemens SPECT-CT, available in our department. 223 Ra photons energy spectra were assessed for low energy high resolution (LEHR), medium energy (ME) and high energy (HE) collimators. Then, depending on the energy window, we calculated the scatter fraction, the sensitivity and the spatial resolution.

RESULTS

Scatter fraction was low for all collimators; however, the contribution of photons that scattered more than twice under the low energy photopeaks was important in the case of LEHR. Sensitivity's best values were obtained in the case of the LEHR collimator; nevertheless, the spatial resolution was very low for this collimator. The latter was best for ME and HE collimators.

CONCLUSION

A combination between a good sensitivity, a high spatial resolution and a low scatter fraction has been determined in the case of the ME collimator, followed by HE collimator as an alternative. To increase the image acquisition statistics with ME collimator, we recommend to use simultaneous energy windows: 20% centered at 82 keV, 20% centered at 154 keV and 20% centered at 270 keV.

Quantitative validation of Monte Carlo SPECT simulation: application to a Mediso AnyScan GATE simulation

BACKGROUND

Monte Carlo (MC) simulations are used in nuclear medicine imaging as they provide unparalleled insight into processes that are not directly experimentally measurable, such as scatter and attenuation in an acquisition. Whilst MC is often used to provide a 'ground-truth', this is only the case if the simulation is fully validated against experimental data. This work presents a quantitative validation for a MC simulation of a single-photon emission computed tomography (SPECT) system.

METHODS

An MC simulation model of the Mediso AnyScan SCP SPECT system installed at the UK National Physical Laboratory was developed in the GATE (Geant4 Application for Tomographic Emission) toolkit. Components of the detector head and two collimator configurations were modelled according to technical specifications and physical measurements. Experimental detection efficiency measurements were collected for a range of energies, permitting an energy-dependent intrinsic camera efficiency correction function to be determined and applied to the simulation on an event-by-event basis. Experimental data were collected in a range of geometries with [Formula: see text]Tc for comparison to simulation. The procedure was then repeated with [Formula: see text]Lu to determine how the validation extended to another isotope and set of collimators.

RESULTS

The simulation's spatial resolution, sensitivity, energy spectra and the projection images were compared with experimental measurements. The simulation and experimental uncertainties were determined and propagated to all calculations, permitting the quantitative agreement between simulated and experimental SPECT acquisitions to be determined. Statistical agreement was seen in sinograms and projection images of both [Formula: see text]Tc and [Formula: see text]Lu data. Average simulated and experimental sensitivity ratios of ([Formula: see text]) were seen for emission and scatter windows of [Formula: see text]Tc, and ([Formula: see text]) and ([Formula: see text]) for the 113 and 208 keV emissions of [Formula: see text]Lu, respectively.

CONCLUSIONS

MC simulations will always be an approximation of a physical system and the level of agreement should be assessed. A validation method is presented to quantify the level of agreement between a simulation model and a physical SPECT system.

Optimized SPECT Imaging of 224Ra α-Particle Therapy by 212Pb Photon Emissions

In preparation for an α-particle therapy trial using 1-7 MBq of ²²⁴Ra, the feasibility of tomographic SPECT/CT imaging was of interest. The nuclide decays in 6 steps to stable ²⁰⁸Pb, with ²¹²Pb as the principle photon-emitting nuclide. ²¹²Bi and ²⁰⁸Tl emit high-energy photons up to 2,615 keV. A phantom study was conducted to determine the optimal acquisition and reconstruction protocol. Methods: The spheres of a body phantom were filled with a ²²⁴Ra-RaCl₂ solution, and the background compartment was filled with water. Images were acquired on a SPECT/CT system. In addition, 30-min scans were acquired for 80- and 240-keV emissions, using triple-energy windows, with both medium-energy and high-energy collimators. Images were acquired at 90-95 and 29-30 kBq/mL, plus an explorative 3-min acquisition at 20 kBq/mL (using only the optimal protocol). Reconstructions were performed with attenuation correction only, attenuation plus scatter correction, 3 levels of postfiltering, and 24 levels of iterative updates. Acquisitions and reconstructions were compared using the maximum value and signal-to-scatter peak ratio for each sphere. Monte Carlo simulations were performed to examine the contributions of key emissions. Results: Secondary photons of the 2,615-keV ²⁰⁸Tl emission produced in the collimators make up most of the acquired energy spectrum, as revealed by Monte Carlo simulations, with only a small fraction (3%-6%) of photons in each window providing useful information for imaging. Still, decent image quality is possible at 30 kBq/mL, and nuclide concentrations are imageable down to approximately 2-5 kBq/mL. The overall best results were obtained with the 240-keV window, medium-energy collimator, attenuation and scatter correction, 30 iterations and 2 subsets, and a 12-mm gaussian postprocessing filter. However, all combinations of the applied collimators and energy windows were capable of producing adequate results, even though some failed to reconstruct the 2 smallest spheres. Conclusion: SPECT/CT imaging of ²²⁴Ra in equilibrium with daughters is possible, with sufficient image quality to provide clinical utility for the current trial of intraperitoneally administrated activity. A systematic scheme for optimization was designed to select acquisition and reconstruction settings.

Mesh modeling of system geometry and anatomy phantoms for realistic GATE simulations and their inclusion in SPECT reconstruction

Objective.Monte-Carlo simulation studies have been essential for advancing various developments in single photon emission computed tomography (SPECT) imaging, such as system design and accurate image reconstruction. Among the simulation software available, Geant4 application for tomographic emission (GATE) is one of the most used simulation toolkits in nuclear medicine, which allows building systems and attenuation phantom geometries based on the combination of idealized volumes. However, these idealized volumes are inadequate for modeling free-form shape components of such geometries. Recent GATE versions alleviate these major limitations by allowing users to import triangulated surface meshes.Approach.In this study, we describe our mesh-based simulations of a next-generation multi-pinhole SPECT system dedicated to clinical brain imaging, called AdaptiSPECT-C. To simulate realistic imaging data, we incorporated in our simulation the XCAT phantom, which provides an advanced anatomical description of the human body. An additional challenge with the AdaptiSPECT-C geometry is that the default voxelized XCAT attenuation phantom was not usable in our simulation due to intersection of objects of dissimilar materials caused by overlap of the air containing regions of the XCAT beyond the surface of the phantom and the components of the imaging system.Main results.We validated our mesh-based modeling against the one constructed by idealized volumes for a simplified single vertex configuration of AdaptiSPECT-C through simulated projection data of123I-activity distributions. We resolved the overlap conflict by creating and incorporating a mesh-based attenuation phantom following a volume hierarchy. We then evaluated our reconstructions with attenuation and scatter correction for projections obtained from simulation consisting of mesh-based modeling of the system and the attenuation phantom for brain imaging. Our approach demonstrated similar performance as the reference scheme simulated in air for uniform and clinical-like123I-IMP brain perfusion source distributions.Significance.This work enables the simulation of complex SPECT acquisitions and reconstructions for emulating realistic imaging data close to those of actual patients.

Accurate Monte Carlo Modeling of Small-Animal Multi-Pinhole SPECT for Non-Standard Multi-Isotope Applications

Recent advances in preclinical SPECT instrumentation enable non-standard multi-isotope acquisitions at the edge of physical feasibility to improve efficiency of pharmaceutical research. Due to the variety of applications, optimization of imaging hardware, acquisition protocols and reconstruction algorithms is a central and recurring task. For this purpose, we developed a Monte Carlo simulation model of a preclinical state-of-the-art multi-pinhole SPECT system, the NanoSPECT/CT^PLUS, with emphasis on high accuracy for multi-isotope experiments operating near the system range limits. The GATE/ GEANT4 model included an accurate description of multi-pinhole collimators and all substructures of the detector back compartment. The readout electronics was modeled with a variety of signal processors partially extended to incorporate non-simplified measured response functions. The final model was able to predict energy spectra, planar images and tomographic reconstructions with high accuracy for both standard and non-standard multi-isotope experiments. Complex activity distributions could be reproduced for a wide range of noise levels and different modes of angular undersampling. Using the example of a dual-isotope triple-tracer experiment, the model has proven to be a powerful tool for protocol optimization and quantitative image correction at the performance range limits of multi-isotope multi-pinhole SPECT.

Generation of clinical 177Lu SPECT/CT images based on Monte Carlo simulation with GATE

PURPOSE

Patient-specific dosimetry in MRT relies on quantitative imaging, pharmacokinetic assessment and absorbed dose calculation. The DosiTest project was initiated to evaluate the uncertainties associated with each step of the clinical dosimetry workflow through a virtual multicentric clinical trial. This work presents the generation of simulated clinical SPECT datasets based on GATE Monte Carlo modelling with its corresponding experimental CT image, which can subsequently be processed by commercial image workstations.

METHODS

This study considers a therapy cycle of 6.85 GBq ¹⁷⁷Lu-labelled DOTATATE derived from an IAEA-Coordinated Research Project (E23005) on "Dosimetry in Radiopharmaceutical therapy for personalised patient treatment". Patient images were acquired on a GE Infinia-Hawkeye 4 gamma camera using a medium energy (ME) collimator. Simulated SPECT projections were generated based on experimental time points and validated against experimental SPECT projections using flattened profiles and gamma index. The simulated projections were then incorporated into the patient SPECT/CT DICOM envelopes for processing and their reconstruction within a commercial image workstation.

RESULTS

Gamma index passing rate (2% - 1 pixel criteria) between 95 and 98% and average gamma between 0.28 and 0.35 among different time points revealed high similarity between simulated and experimental images. Image reconstruction of the simulated projections was successful on HERMES and Xeleris workstations, a major step forward for the initiation of a multicentric virtual clinical dosimetry trial based on simulated SPECT/CT images.

CONCLUSIONS

Realistic ¹⁷⁷Lu patient SPECT projections were generated in GATE. These modelled datasets will be circulated to different clinical departments to perform dosimetry in order to assess the uncertainties in the entire dosimetric chain.

Evaluation of a variable-aperture full-ring SPECT system using large-area pixelated CZT modules: A simulation study for brain SPECT applications

PURPOSE

Single photon emission computed tomography (SPECT) scanners using cadmium zinc telluride (CZT) offer compact, lightweight, and improved imaging capability over conventional NaI(Tl)-based SPECT scanners. The main purpose in this study is to propose a full-ring SPECT system design with eight large-area CZT detectors that can be used for a broad spectrum of SPECT radiopharmaceuticals and demonstrate the performance of our system in comparison to the reference conventional NaI(Tl)-based two-head Anger cameras.

METHODS

A newly designed full-ring SPECT system is composed of eight large-area CZT cameras (128 mm × 179.2 mm effective area) that can be independently swiveled around their own axes of rotation independently and can have radial motion for varying aperture sizes that can be adapted to different sizes of imaging volume. Extended projection data were generated by conjoining projections of two adjacent detectors to overcome the limited field-of-view (FOV) by each CZT camera. Using Monte Carlo simulations, we evaluated this new system design with digital phantoms including a Derenzo hot rod phantom and a Zubal brain phantom. Comparison of performance metrics such as spatial resolution, sensitivity, contrast-to-noise ratio (CNR), and contrast-recovery ratio was made between our design and conventional SPECT scanners having different pixel sizes and radii of rotation (one clinically well-known type and two arbitrary types matched to our proposed CZT-SPECT geometries).

RESULTS

The proposed scanner could result in up to about three times faster in acquisition time over conventional scan time at same acquisition time per step. The spatial resolution improvement, or deterioration, of our proposed scanner compared to the clinical-type scanner was dependent upon the location of the point source. However, there were overall performance improvements over the three different setups of the conventional scanner particularly in volume sensitivity (approximately up to 1.7 times). Overall, we successfully reconstructed the phantom image for both ^99m Tc-based perfusion and ¹²³ I-based dopamine transporter (DaT) brain studies simulated for our new design. In particular, the striatal/background contrast-recovery ratio in 3-to-1 reference ratio was over 0.8 for the ¹²³ I-based DaT study.

CONCLUSIONS

We proposed a variable-aperture full-ring SPECT system using combined pixelated CZT and energy-optimized parallel-hole collimator modules and evaluated the performance of this scanner using relevant digital phantoms and MC simulations. Our studies demonstrated the potential of our new full-ring CZT-SPECT design, showing reduced acquisition time and improved sensitivity with acceptable CNR and spatial resolution.

Modelling of energy-dependent spectral resolution for SPECT Monte Carlo simulations using SIMIND

PURPOSE

Monte Carlo (MC) modelling techniques have been used extensively in Nuclear Medicine (NM). The theoretical energy resolution relationship ( ∝ 1 / E ), does not accurately predict the gamma camera detector response across all energies. This study aimed to validate the accuracy of an energy resolution model for the SIMIND MC simulation code emulating the Siemens Symbia T16 dual-head gamma camera.

METHODS

Measured intrinsic energy resolution data (full width half maximum (FWHM) values), for Ba-133, Lu-177, Am-241, Ga-67, Tc-99m, I-123, I-131 and F-18 sources in air, were used to create a fitted model of the energy response of the gamma camera. Both the fitted and theoretical models were used to simulate intrinsic and extrinsic energy spectra using three different scenarios (source in air; source in simple scatter phantom and a clinical voxel-based digital patient phantom).

RESULTS

The results showed the theoretical model underestimated the FWHM values at energies above 160.0 keV up to 23.5 keV. In contrast, the fitted model better predicted the measured FWHM values with differences less than 3.3 keV. The I-131 in-scatter energy spectrum simulated with the fitted model better matched the measured energy spectrum. Higher energy photopeaks, (I-123: 528.9 keV and I-131: 636.9 keV) simulated with the fitted model, more accurately resembled the measured photopeaks. The voxel-based digital patient phantom energy spectra, simulated with the fitted and theoretical models, showed the potential impact of an incorrect energy resolution model when simulating isotopes with multiple photopeaks.

CONCLUSION

Modelling of energy resolution with the proposed fitted model enables the SIMIND user to accurately simulate NM images. A great improvement was seen for high-energy photon emitting isotopes (e.g. I-131), as well as isotopes with multiple photopeaks (e.g. Lu-177, I-131 and Ga-67) in comparison to the theoretical model. This will result in accurate evaluation of radioactivity quantification, which is vital for dosimetric purposes.

Inclusion of quasi-vertex views in a brain-dedicated multi-pinhole SPECT system for improved imaging performance

With brain-dedicated multi-detector systems employing pinhole apertures the usage of detectors facing the top of the patient's head (i.e. quasi-vertex (QV) views) can provide the advantage of additional viewing from close to the brain for improved detector coverage. In this paper, we report the results of simulation and reconstruction studies to investigate the impact of the QV views on the imaging performance of AdaptiSPECT-C, a brain-dedicated stationary SPECT system under development. In this design, both primary and scatter photons from regions located inferior to the brain can contribute to SPECT projections acquired by the QV views, and thus degrade AdaptiSPECT-C imaging performance. In this work, we determined the proportion, origin, and nature (i.e. primary, scatter, and multiple-scatter) of counts emitted from structures within the head and throughout the body contributing to projections from the different AdaptiSPECT-C detector rings, as well as from a true vertex view detector. We simulated phantoms used to assess different aspects of image quality (i.e. uniform activity concentration sphere, and Derenzo), as well as anthropomorphic phantoms with different count levels emulating clinical ¹²³I activity distributions (i.e. DaTscan and perfusion). We determined that attenuation and scatter in the patient's body greatly diminish the probability of the photons emitted outside the volume of interest reaching to detectors and being recorded within the 15% photopeak energy window. In addition, we demonstrated that the inclusion of the residual of such counts in the system acquisition does not degrade visual interpretation or quantitative analysis. The addition of the QV detectors improves volumetric sensitivity, angular sampling, and spatial resolution leading to significant enhancement in image quality, especially in the striato-thalamic and superior regions of the brain. Besides, the use of QV detectors improves the recovery of clinically relevant metrics such as the striatal binding ratio and mean activity in selected cerebral structures. Our findings proving the usefulness of the QV ring for brain imaging with ¹²³I agents can be generalized to other commonly used SPECT imaging agents labelled with isotopes, such as ^99mTc and likely ¹¹¹In.

SPECT performance evaluation on image of Yttrium 90 - Bremsstrahlung using Monte Carlo simulation

Yttrium-90 (⁹⁰Y) is one of the most widely used radionuclides in Nuclear Medicine practice. However, characteristic energy of this beta emitter constitutes a difficulty for dose planning using SPECT imaging. This work aimed to study bremsstrahlung X-rays effects produced by ⁹⁰Y beta particles during SPECT image acquisition using Monte Carlo code MCNPX. Several simulations were carried out to evaluate different aspects that could affect SPECT image quality, such as: collimator type, source-collimator distance and composition of each interacting material. Two configurations of ⁹⁰Y sources were simulated: a point source in several spheres of different materials (soft tissue, water, articular cartilage, and bone) and dimensions with radius ranging from 1 to 20 mm; and a uniformly distributed source in a Lucite cylindrical phantom filled with water. It was evaluated the bremsstrahlung photon emission generated inside different materials; for this was considered the number photons that passing through every different sphere's surface for each radii and material. In case of cylindrical phantom filled with water, in order to obtain the energy deposited over NaI (Tl) crystal detector; there was considered Median Energy General Purpose (MEGP) and Low Energy High Resolution (LEHR) collimators. Moreover, using TMESH routine available in the MCNPX Monte Carlo code, energy distribution images according to the collimator type and the source-collimator distance were obtained. The simulation was validated by comparing with the spectral distribution of the ⁹⁰Y bremsstrahlung X-rays obtained experimentally from an acrylic cylindrical phantom. Results corroborated the importance of Monte Carlo simulation method to evaluate a performance of SPECT image acquisition with ⁹⁰Y. The best resolution was obtained with MEGP collimator independent of source-collimator distance.

Multi-Isotope Capabilities of a Small-Animal Multi-Pinhole SPECT System

The quantitative accuracy and image quality of multi-isotope SPECT is affected by various hardware-related perturbations. The present study evaluates the simultaneous acquisition of multiple isotopes using a multiplexed multi-pinhole SPECT system, assesses the extent of different error sources, and proposes experimental procedures for its objective characterization. Methods: Phantom measurements with single-, dual- and triple-isotope combinations of ^99mTc, ¹¹¹In, ¹²³I, ¹⁷⁷Lu, and ²⁰¹Tl were performed with the NanoSPECT/CT^PLUS to evaluate system energy resolution, count rate performance, sensitivity, collimator penetration, hardware versus object scatter, spectral crosstalk, spatial resolution, spatial registration accuracy, image uniformity, image noise, and image quality. Results: The intrinsic detector properties were suitable for the simultaneous acquisition of up to 3 isotopes with limitations for count rates exceeding 104 kcps and γ-energies lower than 75 keV. Spectral crosstalk between isotopes was more likely mediated by hardware than by source scatter and was strongly dependent on the isotope combination. Simultaneous multi-isotope acquisitions slightly degraded spatial resolution and image uniformity for spatially superimposed but not for spatially separated activity distributions while the background noise level was increased for all multi-isotope studies. For particular isotopes, collimator penetration and x-ray fluorescence contributed a significant portion of error. Conclusion: The NanoSPECT/CT^PLUS enables the simultaneous acquisition of 3 radioisotopes with high quantitative accuracy and only little loss of image quality when the activity ratio is adapted to isotope-specific count rate sensitivities and when the system calibration is performed with phantoms of appropriate size.

Validation of a Monte Carlo modelled gamma camera for Lutetium-177 imaging

This study validated a model of the Siemens Symbia T16 dual-head SPECT/CT gamma camera created using the Monte Carlo program SIMIND for ¹⁷⁷Lu. The validation was done by comparing experimental and simulated gamma camera performance criteria tests for the ¹⁷⁷Lu 208 keV photopeak with a medium-energy collimator. Results showed good agreement between the experimental and simulated values. These results illustrated that SIMIND could emulate the Symbia T16 successfully and therefore, can be used with confidence to model ¹⁷⁷Lu images.

Primary, scatter, and penetration characterizations of parallel-hole and pinhole collimators for I-123 SPECT

Multi-pinhole (MPH) collimators are known to provide better trade-off between sensitivity and resolution for preclinical, as well as for smaller regions in clinical SPECT imaging compared to conventional collimators. In addition to this geometric advantage, MPH plates typically offer better stopping power for penetration than the conventional collimators, which is especially relevant for I-123 imaging. The I-123 emits a series of high-energy (>300 keV, ~2.5% abundance) gamma photons in addition to the primary emission (159 keV, 83% abundance). Despite their low abundance, high-energy photons penetrate through a low-energy parallel-hole (LEHR) collimator much more readily than the 159 keV photons, resulting in large downscatter in the photopeak window. In this work, we investigate the primary, scatter, and penetration characteristics of a single pinhole collimator that is commonly used for I-123 thyroid imaging and our two MPH collimators designed for I-123 DaTscan imaging for Parkinson's Disease, in comparison to three different parallel-hole collimators through a series of experiments and Monte Carlo simulations. The simulations of a point source and a digital human phantom with DaTscan activity distribution showed that our MPH collimators provide superior count performance in terms of high primary counts, low penetration, and low scatter counts compared to the parallel-hole and single pinhole collimators. For example, total scatter, multiple scatter, and collimator penetration events for the LEHR were 2.5, 7.6 and 14 times more than that of MPH within the 15% photopeak window. The total scatter fraction for LEHR was 56% where the largest contribution came from the high-energy scatter from the back compartments (31%). For the same energy window, the total scatter for MPH was 21% with only 1% scatter from the back compartments. We therefore anticipate that using MPH collimators, higher quality reconstructions can be obtained in a substantially shorter acquisition time for I-123 DaTscan and thyroid imaging.

Development and validation of a full model of a four-headed neuroimaging single-photon emission computed tomography scanner

Objective

The Nucline X-Ring 4R is a four-headed gamma camera dedicated to neuroimaging. In this paper, we describe and validate a GATE (Geant4 Application for Tomographic Emission) model of the Nucline X-Ring 4R.

Materials and methods

Images produced during model simulations were compared with those acquired experimentally to confirm the model was an accurate representation of the scanner. The most commonly reported measurements used to validate a GATE model include energy resolution, spatial resolution and sensitivity. In addition to the commonly reported static imaging measures, single-photon emission computed tomography (SPECT) spatial resolution was investigated to confirm that the model produces similar SPECT images to the experimental output.

Results

The experimental full-width at half-maximum was calculated to be 12.3 keV, which corresponds to an energy resolution of 8.8%. The simulated full-width at half-maximum was measured to be 12 keV, giving an energy resolution of 8.6%. The average spatial resolutions were found to be well matched (5.69 mm – simulated and 5.64 mm – experimental). However, the sensitivity was overestimated using the GATE model (47.8 and 54.3 cps/MBq) compared with the values obtained experimentally (42.7 and 44.3 cps/MBq). Finally, the simulated SPECT spatial resolution images were found to produce qualitatively comparable results.

Conclusion

The model developed has been shown to produce similar results and images to those obtained experimentally. This model has the potential to simulate patient scans with the aim of improving patient care by optimizing scanner protocols.

Molecular imaging using the theranostic agent 197(m)Hg: phantom measurements and Monte Carlo simulations

Background

Radiomercury ^197mHg and ¹⁹⁷Hg, henceforth referred to as ^197(m)Hg, is a promising theranostic radionuclide endowed with properties that allow diagnostic and therapeutic applications. The aim of this work was to investigate the capabilities of ^197(m)Hg for nuclear medicine imaging. Therefore measurements were performed by using a Philips BrightView SPECT camera. Furthermore, Monte Carlo simulations using the GATE software were performed to theoretically explore the imaging contribution from the various gamma and X-ray emissions from ^197(m)Hg for a commercial clinical camera with low-energy high-resolution (LEHR) and high-energy general-purpose (HEGP) collimators. We estimated the spatial resolution by using a four-quadrant bar phantom, and we evaluated the planar and tomographic images from an abdominal phantom containing three cylindrical sources of ^197(m)Hg solution.

Results

A good accordance between measurements and simulations was found for planar and SPECT imaging. Simulations allowed the decomposition of the detected energy spectrum into photon origins. Measurements and simulations for the bar phantom revealed that for the LEHR collimator, the 6-mm pattern could be resolved, whereas for the HEGP collimator, the resolution is about 10 mm. Furthermore, we found that no significant image distortion results from high-energy photons when using the LEHR collimator.

Conclusions

We demonstrated the imaging capabilities of ^197(m)Hg which is essential both for diagnostic applications and to determine the in vivo biodistribution for dose calculations in therapeutic applications.

Fast GPU-based Monte Carlo code for SPECT/CT reconstructions generates improved 177Lu images

Background

Full Monte Carlo (MC)-based SPECT reconstructions have a strong potential for correcting for image degrading factors, but the reconstruction times are long. The objective of this study was to develop a highly parallel Monte Carlo code for fast, ordered subset expectation maximum (OSEM) reconstructions of SPECT/CT images. The MC code was written in the Compute Unified Device Architecture language for a computer with four graphics processing units (GPUs) (GeForce GTX Titan X, Nvidia, USA). This enabled simulations of parallel photon emissions from the voxels matrix (128³ or 256³). Each computed tomography (CT) number was converted to attenuation coefficients for photo absorption, coherent scattering, and incoherent scattering. For photon scattering, the deflection angle was determined by the differential scattering cross sections. An angular response function was developed and used to model the accepted angles for photon interaction with the crystal, and a detector scattering kernel was used for modeling the photon scattering in the detector. Predefined energy and spatial resolution kernels for the crystal were used. The MC code was implemented in the OSEM reconstruction of clinical and phantom ¹⁷⁷Lu SPECT/CT images. The Jaszczak image quality phantom was used to evaluate the performance of the MC reconstruction in comparison with attenuated corrected (AC) OSEM reconstructions and attenuated corrected OSEM reconstructions with resolution recovery corrections (RRC).

Result

The performance of the MC code was 3200 million photons/s. The required number of photons emitted per voxel to obtain a sufficiently low noise level in the simulated image was 200 for a 128³ voxel matrix. With this number of emitted photons/voxel, the MC-based OSEM reconstruction with ten subsets was performed within 20 s/iteration. The images converged after around six iterations. Therefore, the reconstruction time was around 3 min. The activity recovery for the spheres in the Jaszczak phantom was clearly improved with MC-based OSEM reconstruction, e.g., the activity recovery was 88% for the largest sphere, while it was 66% for AC-OSEM and 79% for RRC-OSEM.

Conclusion

The GPU-based MC code generated an MC-based SPECT/CT reconstruction within a few minutes, and reconstructed patient images of ¹⁷⁷Lu-DOTATATE treatments revealed clearly improved resolution and contrast.

131I activity quantification of gamma camera planar images

A procedure to estimate the activity in target tissues in patients during the therapeutic administration of ¹³¹I radiopharmaceutical treatment for thyroid conditions (hyperthyroidism and differentiated thyroid cancer) using a gamma camera (GC) with a high energy (HE) collimator, is proposed. Planar images are acquired for lesions of different sizes r, and at different distances d, in two HE GC systems. Defining a region of interest (ROI) on the image of size r, total counts n _g are measured. Sensitivity S (cps MBq^-1) in each acquisition is estimated as the product of the geometric G and the intrinsic efficiency η ₀. The mean fluence of 364 keV photons arriving at the ROI per disintegration G, is calculated with the MCNPX code, simulating the entire GC and the HE collimator. Intrinsic efficiency η ₀ is estimated from a calibration measurement of a plane reference source of ¹³¹I in air. Values of G and S for two GC systems-Philips Skylight and Siemens e-cam-are calculated. The total range of possible sensitivity values in thyroidal imaging in the e-cam and skylight GC measure from 7 cps MBq^-1 to 35 cps MBq^-1, and from 6 cps MBq^-1 to 29 cps MBq^-1, respectively. These sensitivity values have been verified with the SIMIND code, with good agreement between them. The results have been validated with experimental measurements in air, and in a medium with scatter and attenuation. The counts in the ROI can be produced by direct, scatter and penetration photons. The fluence value for direct photons is constant for any r and d values, but scatter and penetration photons show different values related to specific r and d values, resulting in the large sensitivity differences found. The sensitivity in thyroidal GC planar imaging is strongly dependent on uptake size, and distance from the GC. An individual value for the acquisition sensitivity of each lesion can significantly alleviate the level of uncertainty in the measurement of thyroid uptake activity for each patient.

Characteristics of Bremsstrahlung emissions of (177)Lu, (188)Re, and (90)Y for SPECT/CT quantification in radionuclide therapy

PURPOSE

Beta particles emitted by radioisotopes used in targeted radionuclide therapies (TRT) create Bremsstrahlung (BRS) which may affect SPECT quantification when imaging these isotopes. The purpose of the current study was to investigate the characteristics of Bremsstrahlung produced in tissue by three β-emitting radioisotopes used in TRT.

METHODS

Monte Carlo simulations of (177)Lu, (188)Re, and (90)Y sources placed in water filled cylinders were performed. BRS yields, mean energies and energy spectra for (a) all photons generated in the decays, (b) photons that were not absorbed and leave the cylinder, and (c) photons detected by the camera were analyzed. Next, the results of simulations were compared with those from experiments performed on a clinical SPECT camera using same acquisition conditions and phantom configurations as in simulations.

RESULTS

Simulations reproduced relatively well the shapes of the measured spectra, except for (90)Y which showed an overestimation in the low energy range. Detailed analysis of the results allowed us to suggest best collimators and imaging conditions for each of the investigated isotopes. Finally, our simulations confirmed that the BRS contribution to the energy spectra in quantitative imaging of (177)Lu and (188)Re could be ignored.

CONCLUSIONS

For (177)Lu and (188)Re, BRS contributes only marginally to the total spectra recorded by the camera. Our analysis shows that MELP and HE collimators are the best for imaging these two isotopes. For (90)Y, HE collimator should be used.

Yttrium-90-labeled microsphere tracking during liver selective internal radiotherapy by bremsstrahlung pinhole SPECT: feasibility study and evaluation in an abdominal phantom

BACKGROUND

The purpose of the study is to evaluate whether a pinhole collimator is better adapted to bremsstrahlung single photon emission computed tomography [SPECT] than parallel-hole collimators and in the affirmative, to evaluate whether pinhole bremsstrahlung SPECT, including a simple model of the scatter inside the patient, could provide a fast dosimetry assessment in liver selective internal radiotherapy [SIRT].

MATERIALS AND METHODS

Bremsstrahlung SPECT of an abdominal-shaped phantom including one cold and five hot spheres was performed using two long-bore parallel-hole collimators: a medium-energy general-purpose [MEGP] and a high-energy general-purpose [HEGP], and also using a medium-energy pinhole [MEPH] collimator. In addition, ten helical MEPH SPECTs (acquisition time 3.6 min) of a realistic liver-SIRT phantom were also acquired.

RESULTS

Without scatter correction for SPECT, MEPH SPECT provided a significantly better contrast recovery coefficient [CRC] than MEGP and HEGP SPECTs. The CRCs obtained with MEPH SPECT were still improved with the scatter correction and became comparable to those obtained with positron-emission tomography [PET] for the 36-, 30- (cold), 28-, and 24-mm-diameter spheres: CRC = 1.09, 0.59, 0.91, and 0.69, respectively, for SPECT and CRC = 1.07, 0.56, 0.84, and 0.63, respectively, for PET. However, MEPH SPECT gave the best CRC for the 19-mm-diameter sphere: CRC = 0.56 for SPECT and CRC = 0.01 for PET. The 3.6-min helical MEPH SPECT provided accurate and reproducible activity estimation for the liver-SIRT phantom: relative deviation = 10 ± 1%.

CONCLUSION

Bremsstrahlung SPECT using a pinhole collimator provided a better CRC than those obtained with parallel-hole collimators. The different designs and the better attenuating material used for the collimation (tungsten instead of lead) explain this result. Further, the addition of an analytical modeling of the scattering inside the phantom resulted in an almost fully recovered contrast. This fills the gap between the performance of90Y-PET and bremsstrahlung pinhole SPECT which is a more affordable technique and could even be used during the catheterization procedure in order to optimize the90Y activity to inject.

Fast Monte Carlo-simulator with full collimator and detector response modelling for SPECT

OBJECTIVE

Monte Carlo (MC)-simulations have proved to be a valuable tool in studying SPECT-reconstruction algorithms. Despite their popularity, the use of Monte Carlo-simulations is still often limited by their large computation demand. This is especially true in situations where full collimator and detector modelling with septal penetration, scatter and X-ray fluorescence needs to be included. This paper presents a rapid and simple MC-simulator, which can effectively reduce the computation times.

METHODS

The simulator was built on the convolution-based forced detection principle, which can markedly lower the number of simulated photons. Full collimator and detector response look-up tables are pre-simulated and then later used in the actual MC-simulations to model the system response. The developed simulator was validated by comparing it against (123)I point source measurements made with a clinical gamma camera system and against (99m)Tc software phantom simulations made with the SIMIND MC-package.

RESULTS

The results showed good agreement between the new simulator, measurements and the SIMIND-package. The new simulator provided near noise-free projection data in approximately 1.5 min per projection with (99m)Tc, which was less than one-tenth of SIMIND's time.

CONCLUSION

The developed MC-simulator can markedly decrease the simulation time without sacrificing image quality.