Position dependence of recovery coefficients in 177Lu-SPECT/CT reconstructions - phantom simulations and measurements.

BACKGROUND: Although the importance of quantitative SPECT has increased tremendously due to newly developed therapeutic radiopharmaceuticals, there are still no accreditation programs to harmonize SPECT imaging. Work is currently underway to develop an accreditation for quantitative 177Lu SPECT/CT. The aim of this study is to verify whether the positioning of the spheres within the phantom has an influence on the recovery and thus needs to be considered in SPECT harmonization. In addition, the effects of these recovery coefficients on a potential partial volume correction as well as absorbed-dose estimates are investigated. METHODS: Using a low-dose CT of a SPECT/CT acquisition, a computerized version of the NEMA body phantom was created using a semi-automatic threshold-based method. Based on the mass-density map, the detector orbit, and the sphere centers, realistic SPECT acquisitions of all possible 720 sphere configurations of both the PET and the SPECT versions of the NEMA Body Phantom were generated using Monte Carlo simulations. SPECT reconstructions with different numbers of updates were performed without (CASToR) and with resolution modeling (STIR). Recovery coefficients were calculated for all permutations, reconstruction methods, and phantoms, and their dependence on the sphere positioning was investigated. Finally, the simulation-based findings were validated using SPECT/CT acquisitions of six different sphere configurations. RESULTS: Our analysis shows that sphere positioning has a significant impact on the recovery for both of the reconstruction methods and the phantom type. Although resolution modeling resulted in significantly higher recovery, the relative variation in recovery within the 720 permutations was even larger. When examining the extreme values of the recovery, reconstructions without resolution modeling were influenced primarily by the sphere position, while with resolution modeling the volume of the two adjacent spheres had a larger influence. The SPECT measurements confirmed these observations, and the recovery curves showed good overall agreement with the simulated data. CONCLUSION: Our study shows that sphere positioning has a significant impact on the recovery obtained in NEMA sphere phantom measurements and should therefore be considered in a future SPECT accreditation. Furthermore, the single-measurement method normally performed for PVC should be reconsidered to account for the position dependency.

A Deep-Learning-Based Partial-Volume Correction Method for Quantitative 177Lu SPECT/CT Imaging.

With the development of new radiopharmaceutical therapies, quantitative SPECT/CT has progressively emerged as a crucial tool for dosimetry. One major obstacle of SPECT is its poor resolution, which results in blurring of the activity distribution. Especially for small objects, this so-called partial-volume effect limits the accuracy of activity quantification. Numerous methods for partial-volume correction (PVC) have been proposed, but most methods have the disadvantage of assuming a spatially invariant resolution of the imaging system, which does not hold for SPECT. Furthermore, most methods require a segmentation based on anatomic information. Methods: We introduce DL-PVC, a methodology for PVC of 177Lu SPECT/CT imaging using deep learning (DL). Training was based on a dataset of 10,000 random activity distributions placed in extended cardiac-torso body phantoms. Realistic SPECT acquisitions were created using the SIMIND Monte Carlo simulation program. SPECT reconstructions without and with resolution modeling were performed using the CASToR and STIR reconstruction software, respectively. The pairs of ground-truth activity distributions and simulated SPECT images were used for training various U-Nets. Quantitative analysis of the performance of these U-Nets was based on metrics such as the structural similarity index measure or normalized root-mean-square error, but also on volume activity accuracy, a new metric that describes the fraction of voxels in which the determined activity concentration deviates from the true activity concentration by less than a certain margin. On the basis of this analysis, the optimal parameters for normalization, input size, and network architecture were identified. Results: Our simulation-based analysis revealed that DL-PVC (0.95/7.8%/35.8% for structural similarity index measure/normalized root-mean-square error/volume activity accuracy) outperforms SPECT without PVC (0.89/10.4%/12.1%) and after iterative Yang PVC (0.94/8.6%/15.1%). Additionally, we validated DL-PVC on 177Lu SPECT/CT measurements of 3-dimensionally printed phantoms of different geometries. Although DL-PVC showed activity recovery similar to that of the iterative Yang method, no segmentation was required. In addition, DL-PVC was able to correct other image artifacts such as Gibbs ringing, making it clearly superior at the voxel level. Conclusion: In this work, we demonstrate the added value of DL-PVC for quantitative 177Lu SPECT/CT. Our analysis validates the functionality of DL-PVC and paves the way for future deployment on clinical image data.

Effect of kilovoltage and quality reference mAs on CT-based attenuation correction in 177Lu SPECT/CT imaging: a phantom study.

INTRODUCTION: CT-based attenuation correction (CT-AC) plays a major role in accurate activity quantification by SPECT/CT imaging. However, the effect of kilovoltage peak (kVp) and quality-reference mAs (QRM) on the attenuation coefficient image (µ-map) and volume CT dose index (CTDIvol) have not yet been systematically evaluated. Therefore, the aim of this study was to fill this gap and investigate the influence of kVp and QRM on CT-AC in 177Lu SPECT/CT imaging. METHODS: Seventy low-dose CT acquisitions of an Electron Density Phantom (seventeen inserts of nine tissue-equivalent materials) were acquired using various kVp and QRM combinations on a Siemens Symbia Intevo Bold SPECT/CT system. Using manufacturer reconstruction software, 177Lu µ-maps were generated for each CT image, and three low-dose CT related aspects were examined. First, the µ-map-based attenuation values (µmeasured) were compared with theoretical values (µtheoretical). Second, changes in 177Lu activity expected due to changes in the µ-map were calculated using a modified Chang method. Third, the noise in the µ-map was assessed by measuring the coefficient of variation in a volume of interest in the homogeneous section of the Electron Density Phantom. Lastly, two phantoms were designed to simulate attenuation in four tissue-equivalent materials for two different source geometries (1-mL and 10-mL syringes). 177Lu SPECT/CT imaging was performed using three different reconstruction algorithms (xSPECT Quant, Flash3D, STIR), and the SPECT-based activities were compared against the nominal activities in the sources. RESULTS: The largest relative errors between µmeasured and µtheoretical were observed in the lung inhale insert (range: 18%-36%), while it remained below 6% for all other inserts. The resulting changes in 177Lu activity quantification were -3.5% in the lung inhale insert and less than -2.3% in all other inserts. Coefficient of variation and CTDIvol ranged from 0.3% and 3.6 mGy (130 kVp, 35 mAs) to 0.4% and 0.9 mGy (80 kVp, 20 mAs), respectively. The SPECT-based activity quantification using xSPECT Quant reconstructions outperformed all other reconstruction algorithms. CONCLUSION: This study shows that kVp and QRM values in low-dose CT imaging have a minimum effect on quantitative 177Lu SPECT/CT imaging, while the selection of low values of kVp and QRM reduce the CTDIvol.

Quantitative 177Lu SPECT/CT imaging for personalized dosimetry using a ring-shaped CZT-based camera.

BACKGROUND: Dosimetry after radiopharmaceutical therapy with 177Lu (177Lu-RPT) relies on quantitative SPECT/CT imaging, for which suitable reconstruction protocols are required. In this study, we characterized for the first time the quantitative performance of a ring-shaped CZT-based camera using two different reconstruction algorithms: an ordered subset expectation maximization (OSEM) and a block sequential regularized expectation maximization (BSREM) combined with noise reduction regularization. This study lays the foundations for the definition of a reconstruction protocol enabling accurate dosimetry for patients treated with 177Lu-RPT. METHODS: A series of 177Lu-filled phantoms were acquired on a StarGuide™ (GE HealthCare), with energy and scatter windows centred at 208 (± 6%) keV and 185 (± 5%) keV, respectively. Images were reconstructed with the manufacturer implementations of OSEM (GE-OSEM) and BSREM (Q.Clear) algorithms, and various combinations of iterations and subsets. Additionally, the manufacturer-recommended Q.Clear-based reconstruction protocol was evaluated. Quantification accuracy, measured as the difference between the SPECT-based and the radionuclide calibrator-based activity, and noise were evaluated in a large cylinder. Recovery coefficients (RCs) and spatial resolution were assessed in a NEMA IEC phantom with sphere inserts. The reconstruction protocols considered suitable for clinical applications were tested on a cohort of patients treated with [177Lu]Lu-PSMA-I&T. RESULTS: The accuracy of the activity from the cylinder, although affected by septal penetration, was < 10% for all reconstructions. Both algorithms featured improved spatial resolution and higher RCs with increasing updates at the cost of noise build-up, but Q.Clear outperformed GE-OSEM in reducing noise accumulation. When the reconstruction parameters were carefully selected, similar values for noise (~0.15), spatial resolution (~1 cm) and RCs were found, irrespective of the reconstruction algorithm. Analogue results were found in patients. CONCLUSIONS: Accurate activity quantification is possible when imaging 177Lu with StarGuide™. However, the impact of septal penetration requires further investigations. GE-OSEM is a valid alternative to the recommended Q.Clear reconstruction algorithm, featuring comparable performances assessed on phantoms and patients.

PSMA-PET improves deep learning-based automated CT kidney segmentation.

For dosimetry of radiopharmaceutical therapies, it is essential to determine the volume of relevant structures exposed to therapeutic radiation. For many radiopharmaceuticals, the kidneys represent an important organ-at-risk. To reduce the time required for kidney segmentation, which is often still performed manually, numerous approaches have been presented in recent years to apply deep learning-based methods for CT-based automated segmentation. While the automatic segmentation methods presented so far have been based solely on CT information, the aim of this work is to examine the added value of incorporating PSMA-PET data in the automatic kidney segmentation. METHODS: A total of 108 PET/CT examinations (53 [68Ga]Ga-PSMA-I&T and 55 [18F]F-PSMA-1007 examinations) were grouped to create a reference data set of manual segmentations of the kidney. These segmentations were performed by a human examiner. For each subject, two segmentations were carried out: one CT-based (detailed) segmentation and one PET-based (coarser) segmentation. Five different u-net based approaches were applied to the data set to perform an automated segmentation of the kidney: CT images only, PET images only (coarse segmentation), a combination of CT and PET images, a combination of CT images and a PET-based coarse mask, and a CT image, which had been pre-segmented using a PET-based coarse mask. A quantitative assessment of these approaches was performed based on a test data set of 20 patients, including Dice score, volume deviation and average Hausdorff distance between automated and manual segmentations. Additionally, a visual evaluation of automated segmentations for 100 additional (i.e., exclusively automatically segmented) patients was performed by a nuclear physician. RESULTS: Out of all approaches, the best results were achieved by using CT images which had been pre-segmented using a PET-based coarse mask as input. In addition, this method performed significantly better than the segmentation based solely on CT, which was supported by the visual examination of the additional segmentations. In 80% of the cases, the segmentations created by exploiting the PET-based pre-segmentation were preferred by the nuclear physician. CONCLUSION: This study shows that deep-learning based kidney segmentation can be significantly improved through the addition of a PET-based pre-segmentation. The presented method was shown to be especially beneficial for kidneys with cysts or kidneys that are closely adjacent to other organs such as the spleen, liver or pancreas. In the future, this could lead to a considerable reduction in the time required for dosimetry calculations as well as an improvement in the results.

Analysis of a deep learning-based method for generation of SPECT projections based on a large Monte Carlo simulated dataset.

BACKGROUND: In recent years, a lot of effort has been put in the enhancement of medical imaging using artificial intelligence. However, limited patient data in combination with the unavailability of a ground truth often pose a challenge to a systematic validation of such methodologies. The goal of this work was to investigate a recently proposed method for an artificial intelligence-based generation of synthetic SPECT projections, for acceleration of the image acquisition process based on a large dataset of realistic SPECT simulations. METHODS: A database of 10,000 SPECT projection datasets of heterogeneous activity distributions of randomly placed random shapes was simulated for a clinical SPECT/CT system using the SIMIND Monte Carlo program. Synthetic projections at fixed angular increments from a set of input projections at evenly distributed angles were generated by different u-shaped convolutional neural networks (u-nets). These u-nets differed in noise realization used for the training data, number of input projections, projection angle increment, and number of training/validation datasets. Synthetic projections were generated for 500 test projection datasets for each u-net, and a quantitative analysis was performed using statistical hypothesis tests based on structural similarity index measure and normalized root-mean-squared error. Additional simulations with varying detector orbits were performed on a subset of the dataset to study the effect of the detector orbit on the performance of the methodology. For verification of the results, the u-nets were applied to Jaszczak and NEMA physical phantom data obtained on a clinical SPECT/CT system. RESULTS: No statistically significant differences were observed between u-nets trained with different noise realizations. In contrast, a statistically significant deterioration was found for training with a small subset (400 datasets) of the 10,000 simulated projection datasets in comparison with using a large subset (9500 datasets) for training. A good agreement between synthetic (i.e., u-net generated) and simulated projections before adding noise demonstrates a denoising effect. Finally, the physical phantom measurements show that our findings also apply for projections measured on a clinical SPECT/CT system. CONCLUSION: Our study shows the large potential of u-nets for accelerating SPECT/CT imaging. In addition, our analysis numerically reveals a denoising effect when generating synthetic projections with a u-net. Clinically interesting, the methodology has proven robust against camera orbit deviations in a clinically realistic range. Lastly, we found that a small number of training samples (e.g., ~ 400 datasets) may not be sufficient for reliable generalization of the u-net.

Reconstruction of Pulse Wave and Respiration From Wrist Accelerometer During Sleep.

OBJECTIVE: Nocturnal recordings of heart rate and respiratory rate usually require several separate sensors or electrodes attached to different body parts - a disadvantage for at-home screening tests and for large cohort studies. In this paper, we demonstrate that a state-of-the-art accelerometer placed at subjects' wrists can be used to derive reliable signal reconstructions of heartbeat (pulse wave intervals) and respiration during sleep. METHODS: Based on 226 full-night recordings, we evaluate the performance of our signal reconstruction methodology with respect to polysomnography. We use a phase synchronization analysis metrics that considers individual heartbeats or breaths. RESULTS: The quantitative comparison reveals that pulse-wave signal reconstructions are generally better than respiratory signal reconstructions. The best quality is achieved during deep sleep, followed by light sleep N2 and REM sleep. In addition, a suggested internal evaluation of multiple derived reconstructions can be used to identify time periods with highly reliable signals, particularly for pulse waves. Furthermore, we find that pulse-wave reconstructions are hardly affected by apnea and hypopnea events. CONCLUSION: During sleep, pulse wave and respiration signals can simultaneously be reconstructed from the same accelerometer recording at the wrist without the need for additional sensors. Reliability can be increased by internal evaluation if the reconstructed signals are not needed for the whole sleep duration. SIGNIFICANCE: The presented methodology can help to determine sleep characteristics and improve diagnostics and treatment of sleep disorders in the subjects' normal sleep environment.

Reconstruction of the respiratory signal through ECG and wrist accelerometer data.

Respiratory rate and changes in respiratory activity provide important markers of health and fitness. Assessing the breathing signal without direct respiratory sensors can be very helpful in large cohort studies and for screening purposes. In this paper, we demonstrate that long-term nocturnal acceleration measurements from the wrist yield significantly better respiration proxies than four standard approaches of ECG (electrocardiogram) derived respiration. We validate our approach by comparison with flow-derived respiration as standard reference signal, studying the full-night data of 223 subjects in a clinical sleep laboratory. Specifically, we find that phase synchronization indices between respiration proxies and the flow signal are large for five suggested acceleration-derived proxies with [Formula: see text] for males and [Formula: see text] for females (means ± standard deviations), while ECG-derived proxies yield only [Formula: see text] for males and [Formula: see text] for females. Similarly, respiratory rates can be determined more precisely by wrist-worn acceleration devices compared with a derivation from the ECG. As limitation we must mention that acceleration-derived respiration proxies are only available during episodes of non-physical activity (especially during sleep).