A novel method for harmonization of PET image spatial resolution without phantoms

• Published in: EJNMMI physics Vol. 12; no. 1; pp. 23 - 17
• Main Authors: Carbonell, Felix, Zijdenbos, Alex P., Hempel, Evan, Hajós, Mihály, Bedell, Barry J.
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 14.03.2025; Springer Nature B.V; SpringerOpen
• Subjects: Alzheimer’s disease; Applied and Technical Physics; Clinical trials; Computational Mathematics and Numerical Analysis; Correlation coefficients; Crystallography; Dependent variables; Engineering; Estimation; Estimators; FDG; Fourier transforms; Hoffman Phantom; Image acquisition; Image reconstruction; Imaging; Logarithms; Matching; Medical imaging; Medicine; Medicine & Public Health; Nuclear Medicine; Original Research; PET; Positron emission; Radiology; Scanners; Spatial resolution; Β-amyloid; Β-amyloid; Logarithmic intensity plots; Hoffman Phantom; Alzheimer’s disease; FWHM; FDG; Tau; Fourier transform; Spatial resolution; PET
• ISSN: 2197-7364; 2197-7364
• DOI: 10.1186/s40658-025-00740-9

Background Estimation of the spatial resolution in real images is extremely important in several fields, including crystallography, optics, microscopy, and tomography. In human PET imaging, estimating spatial resolution typically involves the acquisition of images from a physical phantom, typically a Hoffman phantom, which poses a logistical burden, especially in large multi-center studies. Indeed, phantom images may not always be readily available, and this method requires constant monitoring of scanner updates or replacements, scanning protocol changes, and image reconstruction guidelines to establish a equivalence with scans acquired from human subjects. Methods We propose a new computational approach that allows estimation of spatial resolution directly from human subject PET images. The proposed technique is based on the generalization of the logarithmic intensity plots in the 2D Fourier domain to the 3D case. The spatial resolution of the image is obtained through the estimated coefficients of a multiple linear regression problem having the logarithm of the squared norm of the Fourier transform as dependent variable and the squared 3D frequencies as multiple predictors. Results The proposed approach was applied to a cohort of subjects consisting of [18F]florbetapir amyloid PET images and matching phantoms from a Phase II clinical trial, and a second cohort including β-amyloid, FDG, and tau PET images from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) study. The resulting in-plane and axial resolution estimators varied between 3.5 mm and 8.5 mm for both PET and matching phantom images. They also yielded less than one voxel size across-subjects variability in groups of images sharing the same PET scanner model and reconstruction parameters. For human PET images, we also proved that the spatial resolution estimators showed: (1) a very high reproducibility, as measured by intraclass correlation coefficients (ICC > 0.985), (2) a strong cross-tracer linear correlations, and (3) a high within-subject longitudinal consistency, as measured by the maximum difference value between pairs of visits from the same subject. Conclusions Our novel approach does not only eliminate the need for surrogate phantom data, but also provides a general framework that can be applied to a wide range of tracers and other imaging modalities, such as SPECT. Clinical trial data Cognito Therapeutics’ OVERTURE clinical trial (NCT03556280, 2021-08-24), https://clinicaltrials.gov/study/NCT03556280 .