Different partial volume correction methods lead to different conclusions: An 18F-FDG-PET study of aging.

A cross-sectional group study of the effects of aging on brain metabolism as measured with 18F-FDG-PET was performed using several different partial volume correction (PVC) methods: no correction (NoPVC), Meltzer (MZ), Müller-Gärtner (MG), and the symmetric geometric transfer matrix (SGTM) using 99 subjects aged 65-87years from the Harvard Aging Brain study. Sensitivity to parameter selection was tested for MZ and MG. The various methods and parameter settings resulted in an extremely wide range of conclusions as to the effects of age on metabolism, from almost no changes to virtually all of cortical regions showing a decrease with age. Simulations showed that NoPVC had significant bias that made the age effect on metabolism appear to be much larger and more significant than it is. MZ was found to be the same as NoPVC for liberal brain masks; for conservative brain masks, MZ showed few areas correlated with age. MG and SGTM were found to be similar; however, MG was sensitive to a thresholding parameter that can result in data loss. CSF uptake was surprisingly high at about 15% of that in gray matter. The exclusion of CSF from SGTM and MG models, which is almost universally done, caused a substantial loss in the power to detect age-related changes. This diversity of results reflects the literature on the metabolism of aging and suggests that extreme care should be taken when applying PVC or interpreting results that have been corrected for partial volume effects. Using the SGTM, significant age-related changes of about 7% per decade were found in frontal and cingulate cortices as well as primary visual and insular cortices

Quantitative I-131 SPECT Reconstruction using CT Side Information from Hybrid Imaging

A penalized-likelihood (PL) SPECT reconstruction method using a modified regularizer that accounts for anatomical boundary side information was implemented to achieve accurate estimates of both the total target activity and the activity distribution within targets. In both simulations and experimental I-131 phantom studies, reconstructions from 1) penalized likelihood employing CT-side information based regularization (PL-CT); 2) penalized likelihood with edge preserving regularization (no CT); 3) penalized likelihood with conventional spatially invariant quadratic regularization (no CT) were compared with 4) Ordered Subset Expectation Maximization (OSEM), which is the iterative algorithm conventionally used in clinics for quantitative SPECT. Evaluations included phantom studies with perfect and imperfect (misregistered) side information and studies with uniform and non-uniform activity distributions in the target. For targets with uniform activity, the PL-CT images and profiles were closest to the `truth', avoided the edge offshoots evident with OSEM and minimized the blurring across boundaries evident with regularization without CT information. Apart from visual comparison, reconstruction accuracy was evaluated using the bias and standard deviation (STD) of the total target activity estimate and the root mean square error (RMSE) of the activity distribution within the target. PL-CT reconstruction reduced both bias and RMSE compared with regularization without side information. When compared with unregularized OSEM, PL-CT reduced RMSE and STD while bias was comparable. For targets with non-uniform activity, these improvements with PL-CT were observed only when the change in activity was matched by a change in the anatomical image and the corresponding inner boundary was also used to control the regularization. In summary, the present work demonstrates the potential of using CT side information to obtain improved estimates of the activity distribution in targets wi- - thout sacrificing the accuracy of total target activity estimation.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85862/1/Fessler243.pd

Recovery of total I-131 activity within focal volumes using SPECT and 3D OSEM

We experimentally investigated the SPECT recovery of I-131 activity in multiple spheres located simultaneously at different locations within a cylindrical phantom that had an elliptical cross section. The sphere volumes ranged from 209 cc down to 4.2 cc. A Prism 3000 camera and two types of parallel-hexagonal-hole collimation were employed: high energy (HE) and ultra high energy (UHE). Using appropriately-different 3D models of the point source response function for the two types of collimation, approximately the same recovery of activity could be achieved with either collimation by 3D OSEM reconstruction. The recovery coefficient was greater with no background activity in the phantom by 0.10, on average, compared to that with background. In the HE collimation case, the activity recovery was considerably better for all volumes using 3D OSEM reconstruction than it had been in the past using 1D SAGE reconstruction. Recovery-coefficient-based correction in a simulated patient case involving spherical tumours moderately improved the activity estimates (average error reduced from 14% to 9% for UHE collimation, and from 15% to 11% for HE collimation). For a test case with HE collimation, increasing the projection-image sampling density while decreasing the image voxel size increased the recovery coefficient by 0.075 on average, and, if used in a full set of calibration measurements of recovery coefficient versus volume, might lead to further improvement in accuracy for the patient case.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/58098/2/pmb7_3_017.pd

Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging

A penalized-likelihood (PL) SPECT reconstruction method using a modified regularizer that accounts for anatomical boundary side information was implemented to achieve accurate estimates of both the total target activity and the activity distribution within targets. In both simulations and experimental I-131 phantom studies, reconstructions from (1) penalized likelihood employing CT-side information-based regularization (PL-CT), (2) penalized likelihood with edge preserving regularization (no CT) and (3) penalized likelihood with conventional spatially invariant quadratic regularization (no CT) were compared with (4) ordered subset expectation maximization (OSEM), which is the iterative algorithm conventionally used in clinics for quantitative SPECT. Evaluations included phantom studies with perfect and imperfect side information and studies with uniform and non-uniform activity distributions in the target. For targets with uniform activity, the PL-CT images and profiles were closest to the 'truth', avoided the edge offshoots evident with OSEM and minimized the blurring across boundaries evident with regularization without CT information. Apart from visual comparison, reconstruction accuracy was evaluated using the bias and standard deviation (STD) of the total target activity estimate and the root mean square error (RMSE) of the activity distribution within the target. PL-CT reconstruction reduced both bias and RMSE compared with regularization without side information. When compared with unregularized OSEM, PL-CT reduced RMSE and STD while bias was comparable. For targets with non-uniform activity, these improvements with PL-CT were observed only when the change in activity was matched by a change in the anatomical image and the corresponding inner boundary was also used to control the regularization. In summary, the present work demonstrates the potential of using CT side information to obtain improved estimates of the activity distribution in targets without sacrificing the accuracy of total target activity estimation. The method is best suited for data acquired on hybrid systems where SPECT-CT misregistration is minimized. To demonstrate clinical application, the PL reconstruction with CT-based regularization was applied to data from a patient who underwent SPECT/CT imaging for tumor dosimetry following I-131 radioimmunotherapy.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85409/1/pmb10_9_007.pd

Towards accurate partial volume correction in (99m}^Tc oncology SPECT: perturbation for case-specific resolution estimation

BACKGROUND: Currently, there is no consensus on the optimal partial volume correction (PVC) algorithm for oncology imaging. Several existing PVC methods require knowledge of the reconstructed resolution, usually as the point spread function (PSF)-often assumed to be spatially invariant. However, this is not the case for SPECT imaging. This work aimed to assess the accuracy of SPECT quantification when PVC is applied using a case-specific PSF. METHODS: Simulations of SPECT [Formula: see text]Tc imaging were performed for a range of activity distributions, including those replicating typical clinical oncology studies. Gaussian PSFs in reconstructed images were estimated using perturbation with a small point source. Estimates of the PSF were made in situations which could be encountered in a patient study, including; different positions in the field of view, different lesion shapes, sizes and contrasts, noise-free and noisy data. Ground truth images were convolved with the perturbation-estimated PSF, and with a PSF reflecting the resolution at the centre of the field of view. Both were compared with reconstructed images and the root-mean-square error calculated to assess the accuracy of the estimated PSF. PVC was applied using Single Target Correction, incorporating the perturbation-estimated PSF. Corrected regional mean values were assessed for quantitative accuracy. RESULTS: Perturbation-estimated PSF values demonstrated dependence on the position in the Field of View and the number of OSEM iterations. A lower root mean squared error was observed when convolution of the ground truth image was performed with the perturbation-estimated PSF, compared with convolution using a different PSF. Regional mean values following PVC using the perturbation-estimated PSF were more accurate than uncorrected data, or data corrected with PVC using an unsuitable PSF. This was the case for both simple and anthropomorphic phantoms. For the simple phantom, regional mean values were within 0.7% of the ground truth values. Accuracy improved after 5 or more OSEM iterations (10 subsets). For the anthropomorphic phantoms, post-correction regional mean values were within 1.6% of the ground truth values for noise-free uniform lesions. CONCLUSION: Perturbation using a simulated point source could potentially improve quantitative SPECT accuracy via the application of PVC, provided that sufficient reconstruction iterations are used

Fully automated 3D segmentation of dopamine transporter SPECT images using an estimation-based approach

Quantitative measures of uptake in caudate, putamen, and globus pallidus in dopamine transporter (DaT) brain SPECT have potential as biomarkers for the severity of Parkinson disease. Reliable quantification of uptake requires accurate segmentation of these regions. However, segmentation is challenging in DaT SPECT due to partial-volume effects, system noise, physiological variability, and the small size of these regions. To address these challenges, we propose an estimation-based approach to segmentation. This approach estimates the posterior mean of the fractional volume occupied by caudate, putamen, and globus pallidus within each voxel of a 3D SPECT image. The estimate is obtained by minimizing a cost function based on the binary cross-entropy loss between the true and estimated fractional volumes over a population of SPECT images, where the distribution of the true fractional volumes is obtained from magnetic resonance images from clinical populations. The proposed method accounts for both the sources of partial-volume effects in SPECT, namely the limited system resolution and tissue-fraction effects. The method was implemented using an encoder-decoder network and evaluated using realistic clinically guided SPECT simulation studies, where the ground-truth fractional volumes were known. The method significantly outperformed all other considered segmentation methods and yielded accurate segmentation with dice similarity coefficients of ~ 0.80 for all regions. The method was relatively insensitive to changes in voxel size. Further, the method was relatively robust up to +/- 10 degrees of patient head tilt along transaxial, sagittal, and coronal planes. Overall, the results demonstrate the efficacy of the proposed method to yield accurate fully automated segmentation of caudate, putamen, and globus pallidus in 3D DaT-SPECT images

Quantitative imaging for targeted radionuclide therapy dosimetry : technical review

Targeted radionuclide therapy (TRT) is a promising technique for cancer therapy. However, in order to deliver the required dose to the tumor, minimize potential toxicity in normal organs, as well as monitor therapeutic effects, it is important to assess the individualized internal dosimetry based on patient-specific data. Advanced imaging techniques, especially radionuclide imaging, can be used to determine the spatial distribution of administered tracers for calculating the organ-absorbed dose. While planar scintigraphy is still the mainstream imaging method, SPECT, PET and bremsstrahlung imaging have promising properties to improve accuracy in quantification. This article reviews the basic principles of TRT and discusses the latest development in radionuclide imaging techniques for different theranostic agents, with emphasis on their potential to improve personalized TRT dosimetry

Applications of different machine learning approaches in prediction of breast cancer diagnosis delay

Background: The increasing rate of breast cancer (BC) incidence and mortality in Iran has turned this disease into a challenge. A delay in diagnosis leads to more advanced stages of BC and a lower chance of survival, which makes this cancer even more fatal. Objectives: The present study was aimed at identifying the predicting factors for delayed BC diagnosis in women in Iran. Methods: In this study, four machine learning methods, including extreme gradient boosting (XGBoost), random forest (RF), neural networks (NNs), and logistic regression (LR), were applied to analyze the data of 630 women with confirmed BC. Also, different statistical methods, including chi-square, p-value, sensitivity, specificity, accuracy, and area under the receiver operating characteristic curve (AUC), were utilized in different steps of the survey. Results: Thirty percent of patients had a delayed BC diagnosis. Of all the patients with delayed diagnoses, 88.5% were married, 72.1% had an urban residency, and 84.8% had health insurance. The top three important factors in the RF model were urban residency (12.04), breast disease history (11.58), and other comorbidities (10.72). In the XGBoost, urban residency (17.54), having other comorbidities (17.14), and age at first childbirth (>30) (13.13) were the top factors; in the LR model, having other comorbidities (49.41), older age at first childbirth (82.57), and being nulliparous (44.19) were the top factors. Finally, in the NN, it was found that being married (50.05), having a marriage age above 30 (18.03), and having other breast disease history (15.83) were the main predicting factors for a delayed BC diagnosis. Conclusion: Machine learning techniques suggest that women with an urban residency who got married or had their first child at an age older than 30 and those without children are at a higher risk of diagnosis delay. It is necessary to educate them about BC risk factors, symptoms, and self-breast examination to shorten the delay in diagnosis. Copyright © 2023 Dehdar, Salimifard, Mohammadi, Marzban, Saadatmand, Fararouei and Dianati-Nasab

A cycle-consistent adversarial network for brain PET partial volume correction without prior anatomical information

Purpose: Partial volume effect (PVE) is a consequence of the limited spatial resolution of PET scanners. PVE can cause the intensity values of a particular voxel to be underestimated or overestimated due to the effect of surrounding tracer uptake. We propose a novel partial volume correction (PVC) technique to overcome the adverse effects of PVE on PET images. Methods: Two hundred and twelve clinical brain PET scans, including 50 18F-Fluorodeoxyglucose (18F-FDG), 50 18F-Flortaucipir, 36 18F-Flutemetamol, and 76 18F-FluoroDOPA, and their corresponding T1-weighted MR images were enrolled in this study. The Iterative Yang technique was used for PVC as a reference or surrogate of the ground truth for evaluation. A cycle-consistent adversarial network (CycleGAN) was trained to directly map non-PVC PET images to PVC PET images. Quantitative analysis using various metrics, including structural similarity index (SSIM), root mean squared error (RMSE), and peak signal-to-noise ratio (PSNR), was performed. Furthermore, voxel-wise and region-wise-based correlations of activity concentration between the predicted and reference images were evaluated through joint histogram and Bland and Altman analysis. In addition, radiomic analysis was performed by calculating 20 radiomic features within 83 brain regions. Finally, a voxel-wise two-sample t-test was used to compare the predicted PVC PET images with reference PVC images for each radiotracer. Results: The Bland and Altman analysis showed the largest and smallest variance for 18F-FDG (95% CI: − 0.29, + 0.33 SUV, mean = 0.02 SUV) and 18F-Flutemetamol (95% CI: − 0.26, + 0.24 SUV, mean = − 0.01 SUV), respectively. The PSNR was lowest (29.64 ± 1.13 dB) for 18F-FDG and highest (36.01 ± 3.26 dB) for 18F-Flutemetamol. The smallest and largest SSIM were achieved for 18F-FDG (0.93 ± 0.01) and 18F-Flutemetamol (0.97 ± 0.01), respectively. The average relative error for the kurtosis radiomic feature was 3.32%, 9.39%, 4.17%, and 4.55%, while it was 4.74%, 8.80%, 7.27%, and 6.81% for NGLDM_contrast feature for 18F-Flutemetamol, 18F-FluoroDOPA, 18F-FDG, and 18F-Flortaucipir, respectively. Conclusion: An end-to-end CycleGAN PVC method was developed and evaluated. Our model generates PVC images from the original non-PVC PET images without requiring additional anatomical information, such as MRI or CT. Our model eliminates the need for accurate registration or segmentation or PET scanner system response characterization. In addition, no assumptions regarding anatomical structure size, homogeneity, boundary, or background level are required. © 2023, The Author(s)

Assessment of the sources of error affecting the quantitative accuracy of SPECT imaging in small animals

Small animal SPECT imaging systems have multiple potential applications in biomedical research. Whereas SPECT data are commonly interpreted qualitatively in a clinical setting, the ability to accurately quantify measurements will increase the utility of the SPECT data for laboratory measurements involving small animals. In this work, we assess the effect of photon attenuation, scatter and partial volume errors on the quantitative accuracy of small animal SPECT measurements, first with Monte Carlo simulation and then confirmed with experimental measurements. The simulations modeled the imaging geometry of a commercially available small animal SPECT system. We simulated the imaging of a radioactive source within a cylinder of water, and reconstructed the projection data using iterative reconstruction algorithms. The size of the source and the size of the surrounding cylinder were varied to evaluate the effects of photon attenuation and scatter on quantitative accuracy. We found that photon attenuation can reduce the measured concentration of radioactivity in a volume of interest in the center of a rat-sized cylinder of water by up to 50percent when imaging with iodine-125, and up to 25percent when imaging with technetium-99m. When imaging with iodine-125, the scatter-to-primary ratio can reach up to approximately 30percent, and can cause overestimation of the radioactivity concentration when reconstructing data with attenuation correction. We varied the size of the source to evaluate partial volume errors, which we found to be a strong function of the size of the volume of interest and the spatial resolution. These errors can result in large (>50percent) changes in the measured amount of radioactivity. The simulation results were compared with and found to agree with experimental measurements. The inclusion of attenuation correction in the reconstruction algorithm improved quantitative accuracy. We also found that an improvement of the spatial resolution through the use of resolution recovery techniques (i.e. modeling the finite collimator spatial resolution in iterative reconstruction algorithms) can significantly reduce the partial volume errors