Correction for collimator-detector response in SPECT using point spread function template

Image-guided patient-specific optimization of catheter placement for convection-enhanced nanoparticle delivery in recurrent glioblastoma

BACKGROUND

Proper catheter placement for convection-enhanced delivery (CED) is required to maximize tumor coverage and minimize exposure to healthy tissue. We developed an image-based model to patient-specifically optimize the catheter placement for rhenium-186 (186Re)-nanoliposomes (RNL) delivery to treat recurrent glioblastoma (rGBM).

METHODS

The model consists of the 1) fluid fields generated via catheter infusion, 2) dynamic transport of RNL, and 3) transforming RNL concentration to the SPECT signal. Patient-specific tissue geometries were assigned from pre-delivery MRIs. Model parameters were personalized with either 1) individual-based calibration with longitudinal SPECT images, or 2) population-based assignment via leave-one-out cross-validation. The concordance correlation coefficient (CCC) was used to quantify the agreement between the predicted and measured SPECT signals. The model was then used to simulate RNL distributions from a range of catheter placements, resulting in a ratio of the cumulative RNL dose outside versus inside the tumor, the "off-target ratio" (OTR). Optimal catheter placement) was identified by minimizing OTR.

RESULTS

Fifteen patients with rGBM from a Phase I/II clinical trial (NCT01906385) were recruited to the study. Our model, with either individual-calibrated or population-assigned parameters, achieved high accuracy (CCC > 0.80) for predicting RNL distributions up to 24 h after delivery. The optimal catheter placements identified using this model achieved a median (range) of 34.56 % (14.70 %-61.12 %) reduction on OTR at the 24 h post-delivery in comparison to the original placements.

CONCLUSIONS

Our image-guided model achieved high accuracy for predicting patient-specific RNL distributions and indicates value for optimizing catheter placement for CED of radiolabeled liposomes.

From SPECT/CT towards absolute quantification? - the case of unilateral condylar hyperplasia of the mandible

BACKGROUND

Unilateral condylar hyperplasia (UCH) of the mandible is a rare condition characterized by asymmetric growth of the mandibular condyles. Bone scintigraphy with SPECT(/CT) is commonly used to diagnose UCH and guide treatment. Still, varying results have been reported using the traditional threshold of 55%:45% in relative tracer uptake. While absolute quantification of uptake on SPECT/CT could improve results, optimal correction and reconstruction settings are currently unknown.

METHODS

Three anthropomorphic phantoms representing UCH were developed from patient CT volumes and produced using 3D printing technology. Fillable spherical inserts of different sizes (Ø: 8-15 mm) were placed in the condylar positions representing symmetrical and asymmetrical distributions. Recovery coefficients were determined for SPECT/CT using various reconstruction corrections, including attenuation and scatter correction (ACSC), resolution modeling (RM), and partial volume correction (PVC) using phantom measurements. Uptake ratios between condyles and condyle to clivus were evaluated. Finally, the impact of these correction techniques on absolute activity and diagnostic accuracy was assessed in a retrospective patient cohort for the diagnostic threshold of 55%:45%.

RESULTS

The activity was only partially recovered in all spherical inserts (range: 22.5-64.9%). However, RM improved relative recovery by 20.2-62.3% compared to ACSC. In the symmetric phantoms, the 95% confidence interval (CI) of condyle ratios included the diagnostic threshold (57.6%:42.4%) for UCH when using ACSC potentially leading to false positives, but not for ACSCRM datasets. Partial volume corrections coefficients from the NEMA IQ phantom was positionally dependent, with improvements seen performing PVC using coefficients derived from anthropomorphic phantoms. Retrospective application in a patient cohort showed only a weak linear correlation (R²: 0.25-0.67) and large limits of agreement (9.6-12.5%) between different reconstructions. Up to 44% of patients were reclassified using the 55%:45% threshold. Using clinical outcome data, ACSCRM had highest sensitivity (91%; 95% CI 59-100%) and specificity (66%; 95% CI 47-81%), significantly improving specificity (P = 0.038).

CONCLUSIONS

Anthropomorphic phantoms were shown to be essential in determining optimal settings for acquisition, reconstruction, and analysis. SPECT/CT reconstructions with attenuation and scatter correction and resolution modeling are recommended and could improve specificity when using the 55%:45% threshold to assess condylar growth.

Characterization of accurate 3D collimator-detector response function for single- and multi-lofthole collimated SPECT cameras

PURPOSE

Collimator-detector response function (CDRF) of a SPECT scanner refers to the image generated from a point source of activity. This research aims to characterize the CDRF of a breast-dedicated SPECT imager equipped with a lofthole collimator using GATE Monte Carlo simulation.

MATERIALS AND METHODS

To do so, a cylindrical multi-lofthole collimation system with lofthole apertures dedicated to breast imaging was modeled using the GATE Monte Carlo simulator. The dependency of the CDRF on the source-to-collimator distance of a single-lofthole as well as 8-lofthole collimations was assessed and then compared. In addition, the 3D-sensitivity map of the 8-lofthole collimation was derived. Finally, fair comparisons were conducted between the response of the 8-lofthole collimator and that of an 8-pinhole and also existing analytical derivations. In all cases, a data acquisition period of 5.0 min with an in-air 99mTc point source was considered.

RESULTS

For the single-lofthole collimator, 4.5 times increasing the magnification factor leads to a 16- and twofold improvement in the sensitivity and spatial resolution, respectively. In the single-lofthole collimator, the resolution and sensitivity are degraded as the source-to-aperture distance increases. For the cylindrical 8-lofthole collimator, the findings confirm that CDRF strongly depends on source-to-aperture distance and angle of photon incidence. For a 30 mm in-plane offset point, a 25% increase in sensitivity is observed compared to that of the center of the FOV. Increasing the angle from 0∘ to 34∘ results in a 50% reduction in sensitivity. Furthermore, the findings illustrate that spatial resolution follows a quadratic function as10 - 3 d 2 + 2 × 10 - 4 d + R 0 where d is an offset along the x-, y-, and z-axis, and R0 is the spatial resolution at the center of the FOV.

CONCLUSION

In conclusion, both spatial resolution and sensitivity of the lofthole collimation are considerably angle- and offset-dependent within the FOV of single- and multi-lofthole collimated SPECT imagers.

Resolution recovery on list mode MLEM reconstruction for dynamic cardiac SPECT system

The Dynamic Cardiac SPECT (DC-SPECT) system is being developed at the Massachusetts General Hospital, featuring a static cardio focus asymmetrical geometry enabling simultaneous high-resolution and high-sensitivity imaging. Among 14 design iterations of the DC-SPECT with varying number of detector heads, system sensitivity and resolution, the current version under development features 10 mm FWHM geometrical resolution (without resolution recovery) and 0.07% sensitivity at the center of the FOV, this is 1.5× resolution gain and 7× sensitivity gain compared to a conventional dual head gamma camera (0.01% sensitivity and 15-mm resolution). This work presents improvement in imaging resolution by implementing a spatially variant point spread function (SV-PSF) with list mode MLEM reconstruction. A resolution recovery method by PSF deconvolution is validated on list mode MLEM reconstruction for the DC-SPECT. A spatial invariant PSF is included as an additional test to show the influence of the PSF modelling accuracy on reconstructed image quality. We compare the MLEM reconstruction with and without PSF deconvolution; an analytic model is used for the calculation of system response, and the results are compared to the reconstruction with system modelling using Monte Carlo (MC) based methods. Results show that with PSF modelling applied, the quality of the reconstructed image is improved, and the DC-SPECT system can achieve a 4.5 mm central spatial resolution with average 795 counts/Mbq. Both the SV-PSF and the spatial-invariant PSF improve the image quality, and the reconstruction with SV-PSF generates line profiles closer to the ground truth. The results show substantial improvement over the GE Discovery 570c performance (7 mm spatial resolution with an average 460 counts/MBq, 5.8 mm resolution at the FOV center). The impact of PSF deconvolution is significant, improvement of the reconstructed image quality is evident in comparison to MC simulated system matrix with the same sampling size in the simulation.

In silico simulation of hepatic arteries: An open-source algorithm for efficient synthetic data generation

BACKGROUND

In silico testing of novel image reconstruction and quantitative algorithms designed for interventional imaging requires realistic high-resolution modeling of arterial trees with contrast dynamics. Furthermore, data synthesis for training of deep learning algorithms requires that an arterial tree generation algorithm be computationally efficient and sufficiently random.

PURPOSE

The purpose of this paper is to provide a method for anatomically and physiologically motivated, computationally efficient, random hepatic arterial tree generation.

METHODS

The vessel generation algorithm uses a constrained constructive optimization approach with a volume minimization-based cost function. The optimization is constrained by the Couinaud liver classification system to assure a main feeding artery to each Couinaud segment. An intersection check is included to guarantee non-intersecting vasculature and cubic polynomial fits are used to optimize bifurcation angles and to generate smoothly curved segments. Furthermore, an approach to simulate contrast dynamics and respiratory and cardiac motion is also presented.

RESULTS

The proposed algorithm can generate a synthetic hepatic arterial tree with 40 000 branches in 11 s. The high-resolution arterial trees have realistic morphological features such as branching angles (MAD with Murray's law= 1.2 ± 1 . 2 o $ = \;1.2 \pm {1.2^o}$ ), radii (median Murray deviation= 0.08 $ = \;0.08$ ), and smoothly curved, non-intersecting vessels. Furthermore, the algorithm assures a main feeding artery to each Couinaud segment and is random (variability = 0.98 ± 0.01).

CONCLUSIONS

This method facilitates the generation of large datasets of high-resolution, unique hepatic angiograms for the training of deep learning algorithms and initial testing of novel 3D reconstruction and quantitative algorithms designed for interventional imaging.

Activity recovery for differently shaped objects in quantitative SPECT

Objective.The aim was to theoretically and experimentally investigate recovery in SPECT images with objects of different shapes. Furthermore, the accuracy of volume estimation by thresholding was studied for those shapes.Approach.Nine spheres, nine oblate spheroids, and nine prolate spheroids phantom inserts were used, of which the six smaller spheres were part of the NEMA IEC body phantom and the rest of the inserts were 3D-printed. The inserts were filled with99mTc and177Lu. When filled with99mTc, SPECT images were acquired in a Siemens Symbia Intevo Bold gamma camera and when filled with177Lu in a General Electric NM/CT 870 DR gamma camera. The signal rate per activity (SRPA) was determined for all inserts and represented as a function of the volume-to-surface ratio and of the volume-equivalent radius using VOIs defined according to the sphere dimensions and VOIs defined using thresholding. Experimental values were compared with theoretical curves obtained analytically (spheres) or numerically (spheroids), starting from the convolution of a source distribution with a point-spread function. Validation of the activity estimation strategy was performed using four 3D-printed ellipsoids. Lastly, the threshold values necessary to determine the volume of each insert were obtained.Main results.Results showed that SRPA values for the oblate spheroids diverted from the other inserts, when SRPA were represented as a function of the volume-equivalent radius. However, SRPA values for all inserts followed a similar behaviour when represented as a function of the volume-to-surface ratio. Results for ellipsoids were in agreement with those results. For the three types of inserts the volume could be accurately estimated using a threshold method for volumes larger than 25 ml.Significance.Determination of SRPA independently of lesion or organ shape should decrease uncertainties in estimated activities and thereby, in the long term, be beneficial to patient care.

Validation of [Formula: see text]Tc and [Formula: see text]Lu quantification parameters for a Monte Carlo modelled gamma camera

PURPOSE

Monte Carlo (MC) simulation in Nuclear Medicine is a powerful tool for modeling many physical phenomena which are difficult to track or measure directly. MC simulation in SPECT/CT imaging is particularly suitable for optimizing the quantification of activity in a patient, and, consequently, the absorbed dose to each organ. To do so, validating MC results with real data acquired with gamma camera is mandatory. The aim of this study was the validation of the calibration factor (CF) and the recovery coefficient (RC) obtained with SIMIND Monte Carlo code for modeling a Siemens Symbia Intevo Excel SPECT-CT gamma camera to ensure optimal [Formula: see text]Tc and [Formula: see text]Lu SPECT quantification.

METHODS

Phantom experiments using [Formula: see text]Tc and [Formula: see text]Lu have been performed to measure spatial resolution and sensitivity, as well as to evaluate the CF and RC from acquired data. The geometries used for 2D planar imaging were (1) Petri dish and (2) capillary source while for 3D volumetric imaging were (3) a uniform filled cylinder phantom and (4) a Jaszczack phantom with spheres of different volumes. The experimental results have been compared with the results obtained from Monte Carlo simulations performed in the same geometries.

RESULTS

Comparison shows good accordance between simulated and experimental data. The measured planar spatial resolution was 8.3[Formula: see text] mm for [Formula: see text]Tc and 11.8±0.6 mm for [Formula: see text]Lu. The corresponding data obtained by SIMIND for [Formula: see text]Tc was 7.8±0.1 mm, while for [Formula: see text]Lu was 12.4±0.4 mm. The CF was 110.1±5.5 cps/MBq for Technetium and 18.3±1.0 cps/MBq for Lutetium. The corresponding CF obtained by SIMIND for [Formula: see text]Tc was 107.3±0.3 cps/MBq, while for [Formula: see text]Lu 20.4±0.7 cps/MBq. Moreover, a complete curve RCs vs Volume (ml) both for Technetium and Lutetium was determined to correct the PVE for all volumes of clinical interest. In none of the cases, a RC coefficient equal to 100 was found.

CONCLUSIONS

The validation of quantification parameters shows that SIMIND can be used for simulating both gamma camera planar and SPECT images of Siemens Symbia Intevo using [Formula: see text]Tc and [Formula: see text]Lu radionuclides for different medical purposes and treatments.

Biodistribution and dosimetry for combined [177Lu]Lu-PSMA-I&T/[225Ac]Ac-PSMA-I&T therapy using multi-isotope quantitative SPECT imaging

PURPOSE

Quantitative SPECT for patient-specific dosimetry is a valuable tool in the scope of radionuclide therapy, although its clinical application for ²²⁵Ac-based treatments may be limited due to low therapeutic activities. Therefore, the aim of this study was to demonstrate the feasibility of clinical quantitative low-count SPECT imaging during [¹⁷⁷Lu]Lu-PSMA-I&T/[²²⁵Ac]Ac-PSMA-I&T treatment.

METHODS

Eight prostate cancer patients (1000 MBq/8 MBq [¹⁷⁷Lu]Lu-PSMA-I&T/[²²⁵Ac]Ac-PSMA-I&T) received a single-bed quantitative ¹⁷⁷Lu/²²⁵Ac SPECT/CT acquisition (1 h) at 24 h post treatment (high-energy collimator, 16 projections p. head à 3.5 min, 128 × 128 pixel). The gamma peak at 440 keV (width: 10%) of the progeny ²¹³Bi was imaged along with the peak at 208 keV (width: 15%) of ¹⁷⁷Lu. Quantification included CT-based attenuation and window-based scatter correction plus resolution modelling. Gaussian post-filtering with a full-width-half-maximum of 30 mm and 40-45 mm was employed to match the signal-to-noise ratio of ²²⁵Ac and ¹⁷⁷Lu, respectively.

RESULTS

Kidney (r = 0.96, p < 0.01) and lesion (r = 0.94, p < 0.01) SUV for [¹⁷⁷Lu]Lu-PSMA-I&T and [²²⁵Ac]Ac-PSMA-I&T showed a strong and significant correlation. Kidney SUV were significantly higher (p < 0.01) for [²²⁵Ac]Ac-PSMA-I&T (2.5 ± 0.8 vs. 2.1 ± 0.9), while for [¹⁷⁷Lu]Lu-PSMA-I&T lesion SUV were significantly higher (p = 0.03; 1.8 ± 1.1 vs. 2.1 ± 1.5). For absorbed dose estimates, significant differences regarding the kidneys remained, while no significant differences for lesion dosimetry were found.

CONCLUSION

Quantitative low-count SPECT imaging of the peak at 440 keV during [²²⁵Ac]Ac-PSMA-I&T therapy is feasible. Multi-isotope imaging for [¹⁷⁷Lu]Lu-PSMA-I&T/[²²⁵Ac]Ac-PSMA-I&T therapy indicates accumulation of free ²¹³Bi in the kidneys.

Comparison of 90Y SIRT predicted and delivered absorbed doses using a PSF conversion method

PURPOSE

The aims of this study were to develop and apply a method to correct for the differences in partial volume effects of pre-therapy Technetium-99 m (^99mTc)-MAA SPECT and post-therapy Yttrium-90 (⁹⁰Y) bremsstrahlung SPECT imaging in selective internal radiation therapy, and to use this method to improve quantitative comparison of predicted and delivered ⁹⁰Y absorbed doses.

METHODS

The spatial resolution of ^99mTc SPECT data was converted to that of ⁹⁰Y SPECT data using a function calculated from ^99mTc and ⁹⁰Y point spread functions. This resolution conversion method (RCM) was first applied to ^99mTc and ⁹⁰Y SPECT phantom data to validate the method, and then to clinical data to assess the power of ^99mTc SPECT imaging to predict the therapeutic absorbed dose.

RESULTS

The maximum difference between absorbed doses to phantom spheres was 178%. This was reduced to 27% after the RCM was applied. The clinical data demonstrated differences within 38% for mean absorbed doses delivered to the normal liver, which were reduced to 20% after application of the RCM. Analysis of clinical data showed that therapeutic absorbed doses delivered to tumours greater than 100 cm³ were predicted to within 52%, although there were differences of up to 210% for smaller tumours, even after the RCM was applied.

CONCLUSIONS

The RCM was successfully verified using phantom data. Analysis of the clinical data established that the ^99mTc pre-therapy imaging was predictive of the ⁹⁰Y absorbed dose to the normal liver to within 20%, but had poor predictability for tumours smaller than 100 cm³.

Patient specific, imaging-informed modeling of rhenium-186 nanoliposome delivery via convection-enhanced delivery in glioblastoma multiforme

Convection-enhanced delivery of rhenium-186 (¹⁸⁶Re)-nanoliposomes is a promising approach to provide precise delivery of large localized doses of radiation for patients with recurrent glioblastoma multiforme. Current approaches for treatment planning utilizing convection-enhanced delivery are designed for small molecule drugs and not for larger particles such as¹⁸⁶Re-nanoliposomes. To enable the treatment planning for¹⁸⁶Re-nanoliposomes delivery, we have developed a computational fluid dynamics approach to predict the distribution of nanoliposomes for individual patients. In this work, we construct, calibrate, and validate a family of computational fluid dynamics models to predict the spatio-temporal distribution of¹⁸⁶Re-nanoliposomes within the brain, utilizing patient-specific pre-operative magnetic resonance imaging (MRI) to assign material properties for an advection-diffusion transport model. The model family is calibrated to single photon emission computed tomography (SPECT) images acquired during and after the infusion of¹⁸⁶Re-nanoliposomes for five patients enrolled in a Phase I/II trial (NCT Number NCT01906385), and is validated using a leave-one-out bootstrapping methodology for predicting the final distribution of the particles. After calibration, our models are capable of predicting the mid-delivery and final spatial distribution of¹⁸⁶Re-nanoliposomes with a Dice value of 0.69 ± 0.18 and a concordance correlation coefficient of 0.88 ± 0.12 (mean ± 95% confidence interval), using only the patient-specific, pre-operative MRI data, and calibrated model parameters from prior patients. These results demonstrate a proof-of-concept for a patient-specific modeling framework, which predicts the spatial distribution of nanoparticles. Further development of this approach could enable optimizing catheter placement for future studies employing convection-enhanced delivery.

Description of the methodology for dosimetric quantification in treatments with 177Lu-DOTATATE

Implementation of dosimetry calculations in the daily practice of Nuclear Medicine Departments is, at this time, a controversial issue, partly due to the lack of a standardized methodology that is accepted by all interested parties (patients, nuclear medicine physicians and medical physicists). However, since the publication of RD 601/2019 there is a legal obligation to implement it, despite the fact that it is a complex and high resource consumption procedure. The aim of this article is to review the theoretical bases of in vivo dosimetry in treatments with 177Lu-DOTATATE. The exposed methodology is the one proposed by the MIRD Committee (Medical Internal Radiation Dose) of the SNMMI (Society of Nuclear Medicine & Molecular Imaging). According to this method, the absorbed dose is obtained as the product of 2factors: the time-integrated activity of the radiopharmaceutical present in a source region and a geometrical factor S. This approach, which a priori seems simple, in practice requires several SPECT/CT acquisitions, several measurements of the whole body activity and taking several blood samples, as well as hours of image processing and computation. The systematic implementation of these calculations, in all the patients we treat, will allow us to obtain homogeneous data to correlate the absorbed doses in the lesions with the biological effect of the treatment. The final purpose of the dosimetry calculations is to be able to maximize the therapeutic effect in the lesions, controlling the radiotoxicity in the organs at risk.

Image quality characteristics of myocardial perfusion SPECT imaging using state-of-the-art commercial software algorithms: evaluation of 10 reconstruction methods

Myocardial perfusion imaging (MPI) is widely used as standard of care in patients with coronary artery disease. The availability of hybrid SPECT/CT imaging system and associated advanced reconstruction algorithms serve to improve diagnostic accuracy and enhances image quality of MPI. The aim of this work was to evaluate the relative performance of iterative reconstruction algorithms correcting for different combinations of image degrading factors versus filtered back projection on the quality of myocardial perfusion SPECT imaging. A standard cardiac phantom containing myocardial defects of different sizes and compositions was used to simulate myocardial perfusion SPECT/CT clinical studies. A clinically relevant activity was determined to avoid discordance with real data acquisition. Acquisition parameters including time per projection, angular rotation increment, and iterative reconstruction number of iterations and subsets were varied. The reconstruction was carried out applying different algorithms including 10 variants of analytical (e.g FBP) and iterative reconstructions with and without resolution recovery. Typical figures of merit were used to evaluate the image quality of MPI reconstructed with ten different reconstruction methods. OSEM-RR showed remarkable improvement of image quality of MPI in terms of SNR, CNR and defect contrast percentage compared to FBP algorithm. Full correction scheme IR-RR (i.e. IRACSCRR) provides clinically acceptable image quality of MPI compared to FBP.

Collimator and Energy Window Evaluation in Ga-67 Imaging by Monte Carlo Simulation

Objectives:

Gallium-67 (Ga-67) imaging is affected by collimator penetration and scatter components owing to the high-energy (HE) gamma-ray emissions. The characterization of penetration and scatter distribution is essential for the optimization of low-energy high-resolution (LEHR), medium energy (ME), and HE collimators and for the development of an effective correction technique. We compared the image quality that can be achieved by 3 collimators for different energy windows using the SIMIND Monte Carlo code.

Methods:

Simulation experiments were conducted for LEHR, ME, and HE collimators for Ga-67 point source placed at 12-cm distance from the detector surface using the Monte Carlo SIMIND simulation code. Their spectra point spread functions as well as the original, penetration, scattering, and X-rays curves were drawn and analyzed. The parameters full-width at half maximum and full-width at tenth maximum were also investigated.

Results:

The original, penetration, and scatter curves within 10% for LEHR were 34.46%, 33.52%, 17.29%, and 14.72%, respectively. Similarly, the original, penetration, scatter, and X-rays within 10% for ME and HE were 83.06%, 10.25%, 6.69%, and 0% and 81.44%, 11.51%, 7.05%, and 0%, respectively. The trade-off between spatial resolution and sensitivity was achieved by using the ME collimator at 185 photopeak of Ga-67.

Conclusion:

The Monte Carlo simulation outcomes can be applied for optimal collimator designing and for the development of new correction method in Ga-67 imaging.

Algorithms and Analyses for Joint Spectral Image Reconstruction in Y-90 Bremsstrahlung SPECT

Quantitative yttrium-90 (Y-90) SPECT imaging is challenging due to the nature of Y-90, an almost pure beta emitter that is associated with a continuous spectrum of bremsstrahlung photons that have a relatively low yield. This paper proposes joint spectral reconstruction (JSR), a novel bremsstrahlung SPECT reconstruction method that uses multiple narrow acquisition windows with accurate multi-band forward modeling to cover a wide range of the energy spectrum. Theoretical analyses using Fisher information and Monte-Carlo (MC) simulation with a digital phantom show that the proposed JSR model with multiple acquisition windows has better performance in terms of covariance (precision) than previous methods using multi-band forward modeling with a single acquisition window, or using a single-band forward modeling with a single acquisition window. We also propose an energy-window subset (ES) algorithm for JSR to achieve fast empirical convergence and maximum-likelihood based initialization for all reconstruction methods to improve quantification accuracy in early iterations. For both MC simulation with a digital phantom and experimental study with a physical multi-sphere phantom, our proposed JSR-ES, a fast algorithm for JSR with ES, yielded higher recovery coefficients (RCs) on hot spheres over all iterations and sphere sizes than all the other evaluated methods, due to fast empirical convergence. In experimental study, for the smallest hot sphere (diameter 1.6cm), at the 20th iteration the increase in RCs with JSR-ES was 66 and 31% compared with single wide and narrow band forward models, respectively. JSR-ES also yielded lower residual count error (RCE) on a cold sphere over all iterations than other methods for MC simulation with known scatter, but led to greater RCE compared with single narrow band forward model at higher iterations for experimental study when using estimated scatter.

Quantitative SPECT/CT-Technique and Clinical Applications

The continuous development of SPECT over the past 50 years has led to improved image quality and increased diagnostic confidence. The most influential developments include the realization of hybrid SPECT/CT devices, as well as the implementation of attenuation correction and iterative image reconstruction techniques. These developments have led to a preference for SPECT/CT devices over SPECT-only systems and to the widespread adoption of the former, strengthening the role of SPECT/CT as the workhorse of Nuclear Medicine imaging. New trends in the ongoing development of SPECT/CT are diverse. For example, whole-body SPECT/CT images, consisting of acquisitions from multiple consecutive bed positions in the manner of PET/CT, are increasingly performed. Additionally, in recent years, some interesting approaches in detector technology have found their way into commercial products. For example, some SPECT cameras dedicated to specific organs employ semiconductor detectors made of cadmium telluride or cadmium zinc telluride, which have been shown to increase the obtainable image quality by offering a higher sensitivity and energy resolution. However, the advent of quantitative SPECT/CT which, like PET, can quantify the amount of tracer in terms of Bq/mL or as a standardized uptake value could be regarded as most important development. It is a major innovation that will lead to increased diagnostic accuracy and confidence, especially in longitudinal studies and in the monitoring of treatment response. The current work comprises two main aspects. At first, physical and technical fundamentals of SPECT image formation are described and necessary prerequisites of quantitative SPECT/CT are reviewed. Additionally, the typically achievable quantitative accuracy based on reports from the literature is given. Second, an extensive list of studies reporting on clinical applications of quantitative SPECT/CT is provided and reviewed.

Fast quantitative reconstruction with focusing collimators for liver SPECT

BACKGROUND

Generation of a SPECT scan during procedure may aid in the optimization of treatments as liver radioembolization by offering image-guided dosimetry. This, however, requires both shortened acquisition times and fast quantitative reconstruction. Focusing collimators increase sensitivity and thus may speed up imaging. Monte Carlo-based iterative reconstruction has shown to provide quantitative results for parallel hole collimators but may be slow. The purpose of this work is to develop fast Monte Carlo-based reconstruction for focusing collimators and to evaluate the impact of reconstruction and collimator choice on quantitative accuracy of liver dosimetry by means of simulations.

RESULTS

The developed fast Monte Carlo simulator was found to accurately generate projections compared to a full Monte Carlo simulation, providing projections in several seconds instead of several days. Monte Carlo-based scatter correction was superior to other scatter correction methods in describing recovered activity and reached similar noise levels as dual-energy window scatter correction. Although truncation artifacts were present in the cone beam collimator (50 cm), the region inside the field of view (FOV) could be reconstructed without loss of accuracy. Provided the object to image is inside the FOV, the focusing collimator with 50 cm focal distance could retrieve the same noise levels as a parallel hole collimator in 68% of the total scanning time, the multifocal collimator in 73% of the time, and the 100-cm focal distance collimator in 84% of the time.

CONCLUSION

Focusing collimators combined with Monte Carlo-based reconstruction have the ability to enable quantitative imaging of the FOV in a significantly shorter timeframe. The proposed approach to the forward projector will additionally make it possible to reconstruct within minutes. These are crucial steps in moving toward real-time dosimetry during interventions.

Evaluation of quantitative 123I and 131I SPECT with Monte Carlo-based down-scatter compensation

OBJECTIVE

Quantitative I and I single-photon emission computed tomography (SPECT) is hampered by down-scatter from the high-energy peaks. This paper presents a down-scatter compensation method, where down-scatter generated in the patient and gamma camera collimator and detector is modelled using Monte Carlo simulation in the ordered subsets expectation maximization SPECT reconstruction algorithm.

MATERIALS AND METHODS

The new down-scatter compensation method was compared with conventional triple energy window (TEW) scatter compensation and Gaussian convolution-based forced detection Monte Carlo methods. The comparison was made with the NEMA-IEC phantom using six spherical inserts (diameters from 10 to 37 mm) and a lung compartment. The phantom was filled with I and I solutions to known sphere-to-background concentration ratios. Spherical volumes of interest with the same diameter as the inserts were drawn on the images, and recovery coefficients for the spheres were calculated in addition to lung-to-background ratio.

RESULTS

The new down-scatter compensation method provided higher recovery coefficients than the TEW scatter compensation or Gaussian convolution-based forced detection Monte Carlo algorithm for both isotopes. Background activity concentration could be accurately estimated with the new down-scatter compensation method and with the TEW scatter compensation, whereas activity concentration of the spheres was severely underestimated even with the new method.

CONCLUSION

Down-scatter compensation with Monte Carlo-simulation effectively reduces down-scatter effects in I and I SPECT imaging.

89 Zr-ImmunoPET companion diagnostics and their impact in clinical drug development

Therapeutic monoclonal antibodies have been used in cancer treatment for 30 years, with around 24 mAb and mAb:drug conjugates approved by the FDA to date. Despite their specificity, efficacy has remained limited, which, in part, derails nascent initiatives towards precision medicine. An image-guided approach to reinforce treatment decisions using immune positron emission tomography (immunoPET) companion diagnostic is warranted. This review provides a general overview of current translational research using Zr-89 immunoPET and opportunities for utilizing and harnessing this tool to its full potential. Patient case studies are cited to illustrate immunoPET probes as tools for profiling molecular signatures. Discussions on its utility in reinforcing clinical decisions as it relates to histopathological tumor assessment and standard diagnostic methods, and its potential as predictive biomarkers, are presented. We finally conclude with an overview of practical considerations to its utility in the clinic.

Penalized-Likelihood Reconstruction With High-Fidelity Measurement Models for High-Resolution Cone-Beam Imaging

We present a novel reconstruction algorithm based on a general cone-beam CT forward model, which is capable of incorporating the blur and noise correlations that are exhibited in flat-panel CBCT measurement data. Specifically, the proposed model may include scintillator blur, focal-spot blur, and noise correlations due to light spread in the scintillator. The proposed algorithm (GPL-BC) uses a Gaussian Penalized-Likelihood objective function, which incorporates models of blur and correlated noise. In a simulation study, GPL-BC was able to achieve lower bias as compared with deblurring followed by FDK as well as a model-based reconstruction method without integration of measurement blur. In the same study, GPL-BC was able to achieve better line-pair reconstructions (in terms of segmented-image accuracy) as compared with deblurring followed by FDK, a model-based method without blur, and a model-based method with blur but not noise correlations. A prototype extremities quantitative cone-beam CT test-bench was used to image a physical sample of human trabecular bone. These data were used to compare reconstructions using the proposed method and model-based methods without blur and/or correlation to a registered CT image of the same bone sample. The GPL-BC reconstructions resulted in more accurate trabecular bone segmentation. Multiple trabecular bone metrics, including trabecular thickness (Tb.Th.) were computed for each reconstruction approach as well as the CT volume. The GPL-BC reconstruction provided the most accurate Tb.Th. measurement, 0.255 mm, as compared with the CT derived value of 0.193 mm, followed by the GPL-B reconstruction, the GPL-I reconstruction, and then the FDK reconstruction (0.271 mm, 0.309 mm, and 0.335 mm, respectively).

Application of the 4-D XCAT Phantoms in Biomedical Imaging and Beyond

The four-dimensional (4-D) eXtended CArdiac-Torso (XCAT) series of phantoms was developed to provide accurate computerized models of the human anatomy and physiology. The XCAT series encompasses a vast population of phantoms of varying ages from newborn to adult, each including parameterized models for the cardiac and respiratory motions. With great flexibility in the XCAT's design, any number of body sizes, different anatomies, cardiac or respiratory motions or patterns, patient positions and orientations, and spatial resolutions can be simulated. As such, the XCAT phantoms are gaining a wide use in biomedical imaging research. There they can provide a virtual patient base from which to quantitatively evaluate and improve imaging instrumentation, data acquisition, techniques, and image reconstruction and processing methods which can lead to improved image quality and more accurate clinical diagnoses. The phantoms have also found great use in radiation dosimetry, radiation therapy, medical device design, and even the security and defense industry. This review paper highlights some specific areas in which the XCAT phantoms have found use within biomedical imaging and other fields. From these examples, we illustrate the increasingly important role that computerized phantoms and computer simulation are playing in the research community.

Improved quantitative 90 Y bremsstrahlung SPECT/CT reconstruction with Monte Carlo scatter modeling

PURPOSE

In ⁹⁰ Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative ⁹⁰ Y imaging by either SPECT or PET is challenging. In ⁹⁰ Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate ⁹⁰ Y bremsstrahlung SPECT/CT imaging.

METHODS

Since a fully Monte Carlo (MC) approach to ⁹⁰ Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration.

RESULTS

There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively.

CONCLUSIONS

Quantitative accuracy of ⁹⁰ Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.

The Use of Anatomical Information for Molecular Image Reconstruction Algorithms: Attenuation/Scatter Correction, Motion Compensation, and Noise Reduction

PET and SPECT are important tools for providing valuable molecular information about patients to clinicians. Advances in nuclear medicine hardware technologies and statistical image reconstruction algorithms enabled significantly improved image quality. Sequentially or simultaneously acquired anatomical images such as CT and MRI from hybrid scanners are also important ingredients for improving the image quality of PET or SPECT further. High-quality anatomical information has been used and investigated for attenuation and scatter corrections, motion compensation, and noise reduction via post-reconstruction filtering and regularization in inverse problems. In this article, we will review works using anatomical information for molecular image reconstruction algorithms for better image quality by describing mathematical models, discussing sources of anatomical information for different cases, and showing some examples.

Model-based iterative reconstruction for flat-panel cone-beam CT with focal spot blur, detector blur, and correlated noise

While model-based reconstruction methods have been successfully applied to flat-panel cone-beam CT (FP-CBCT) systems, typical implementations ignore both spatial correlations in the projection data as well as system blurs due to the detector and focal spot in the x-ray source. In this work, we develop a forward model for flat-panel-based systems that includes blur and noise correlation associated with finite focal spot size and an indirect detector (e.g. scintillator). This forward model is used to develop a staged reconstruction framework where projection data are deconvolved and log-transformed, followed by a generalized least-squares reconstruction that utilizes a non-diagonal statistical weighting to account for the correlation that arises from the acquisition and data processing chain. We investigate the performance of this novel reconstruction approach in both simulated data and in CBCT test-bench data. In comparison to traditional filtered backprojection and model-based methods that ignore noise correlation, the proposed approach yields a superior noise-resolution tradeoff. For example, for a system with 0.34 mm FWHM scintillator blur and 0.70 FWHM focal spot blur, using the correlated noise model instead of an uncorrelated noise model increased resolution by 42% (with variance matched at 6.9  ×  10(-8) mm(-2)). While this advantage holds across a wide range of systems with differing blur characteristics, the improvements are greatest for systems where source blur is larger than detector blur.

Investigation of Collimator Influential Parameter on SPECT Image Quality: a Monte Carlo Study

BACKGROUND

Obtaining high quality images in Single Photon Emission Tomography (SPECT) device  is the most important goal in nuclear medicine. Because if image quality is low, the possibility of making a mistake in diagnosing and treating the patient will rise. Studying effective factors in spatial resolution of imaging systems is thus deemed to be vital. One of the most important factors in SPECT imaging in nuclear medicine is the use of an appropriate collimator for a certain radiopharmaceutical feature in order to create the best image as it can be effective in the quantity of Full Width at Half Maximum (FWHM) which is the main parameter  in spatial resolution.

METHOD

In this research, the simulation of the detector and collimator of SPECT imaging device, Model HD3 made by Philips Co. and the investigation of important factors on the collimator were carried out using MCNP-4c code.

RESULTS

The results of the experimental measurments and simulation calculations revealed a relative difference of less than 5% leading to the confirmation of the accuracy of conducted simulation MCNP code calculation.

CONCLUSION

This is the first essential step in the design and modelling of new collimators used for creating high quality images in nuclear medicine.

Alternating direction method of multiplier for tomography with nonlocal regularizers

The ordered subset expectation maximization (OSEM) algorithm approximates the gradient of a likelihood function using a subset of projections instead of using all projections so that fast image reconstruction is possible for emission and transmission tomography such as SPECT, PET, and CT. However, OSEM does not significantly accelerate reconstruction with computationally expensive regularizers such as patch-based nonlocal (NL) regularizers, because the regularizer gradient is evaluated for every subset. We propose to use variable splitting to separate the likelihood term and the regularizer term for penalized emission tomographic image reconstruction problem and to optimize it using the alternating direction method of multiplier (ADMM). We also propose a fast algorithm to optimize the ADMM parameter based on convergence rate analysis. This new scheme enables more sub-iterations related to the likelihood term. We evaluated our ADMM for 3-D SPECT image reconstruction with a patch-based NL regularizer that uses the Fair potential function. Our proposed ADMM improved the speed of convergence substantially compared to other existing methods such as gradient descent, EM, and OSEM using De Pierro's approach, and the limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm.

Post-reconstruction non-local means filtering methods using CT side information for quantitative SPECT

Quantitative SPECT techniques are important for many applications including internal emitter therapy dosimetry where accurate estimation of total target activity and activity distribution within targets are both potentially important for dose–response evaluations. We investigated non-local means (NLM) post-reconstruction filtering for accurate I-131 SPECT estimation of both total target activity and the 3D activity distribution. We first investigated activity estimation versus number of ordered-subsets expectation–maximization (OSEM) iterations. We performed simulations using the XCAT phantom with tumors containing a uniform and a non-uniform activity distribution, and measured the recovery coefficient (RC) and the root mean squared error (RMSE) to quantify total target activity and activity distribution, respectively. We observed that using more OSEM iterations is essential for accurate estimation of RC, but may or may not improve RMSE. We then investigated various post-reconstruction filtering methods to suppress noise at high iteration while preserving image details so that both RC and RMSE can be improved. Recently, NLM filtering methods have shown promising results for noise reduction. Moreover, NLM methods using high-quality side information can improve image quality further. We investigated several NLM methods with and without CT side information for I-131 SPECT imaging and compared them to conventional Gaussian filtering and to unfiltered methods. We studied four different ways of incorporating CT information in the NLM methods: two known (NLM CT-B and NLM CT-M) and two newly considered (NLM CT-S and NLM CT-H). We also evaluated the robustness of NLM filtering using CT information to erroneous CT. NLM CT-S and NLM CT-H yielded comparable RC values to unfiltered images while substantially reducing RMSE. NLM CT-S achieved −2.7 to 2.6% increase of RC compared to no filtering and NLM CT-H yielded up to 6% decrease in RC while other methods yielded lower RCs than them: Gaussian filtering (up to 11.8% decrease in RC), NLM method without CT (up to 9.5% decrease in RC), and NLM CT-M and NLM CT-B (up to 19.4% decrease in RC). NLM CT-S and NLM CT-H achieved 8.2 to 33.9% and −0.9 to 36% decreased RMSE on tumors compared to no filtering respectively while other methods yielded less reduced or increased RMSE: Gaussian filtering (up to 7.9% increase in RMSE), NLM method without CT (up to 18.3% increase in RMSE), and NLM CT-M and NLM CT-B (up to 31.5% increase in RMSE). NLM CT-S and NLM CT-H also yielded images with tumor shapes that better-matched the true shapes than other methods. All NLM methods using CT information were robust to small misregistration between SPECT and CT, but NLM CT-S and NLM CT-H were more sensitive than NLM CT-M and NLM CT-B to missing CT information.