Motion Dependent and Spatially Variant Resolution Modeling for PET Rigid Motion Correction

Abstract

Recent advances in positron emission tomography (PET) have allowed to perform brain scans of freely moving animals by using rigid motion correction. One of the current challenges in these scans is that, due to the PET scanner spatially variant point spread function (SVPSF), motion corrected images have a motion dependent blurring since animals can move throughout the entire field of view (FOV). We developed a method to calculate the image-based resolution kernels of the motion dependent and spatially variant PSF (MD-SVPSF) to correct the loss of spatial resolution in motion corrected reconstructions. The resolution kernels are calculated for each voxel by sampling and averaging the SVPSF at all positions in the scanner FOV where the moving object was measured. In resolution phantom scans, the use of the MD-SVPSF resolution model improved the spatial resolution in motion corrected reconstructions and corrected the image deformation caused by the parallax effect consistently for all motion patterns, outperforming the use of a motion independent SVPSF or Gaussian kernels. Compared to motion correction in which the SVPSF is applied independently for every pose, our method performed similarly, but with more than two orders of magnitude faster computation time. Importantly, in scans of freely moving mice, brain regional quantification in motion-free and motion corrected images was better correlated when using the MD-SVPSF in comparison with motion independent SVPSF and a Gaussian kernel. The method developed here allows to obtain consistent spatial resolution and quantification in motion corrected images, independently of the motion pattern of the subject.

Development of a ligand for in vivo imaging of mutant huntingtin in Huntington’s disease

“…PET scans were followed by a 10-min 80-kV/500-μA CT scan on the same gantry for attenuation correction and coregistration purposes. Acquired PET data were reconstructed into 33 or 39 (for 60- or 90-min acquisition, respectively) frames of increasing length (12 × 10s, 3 × 20s, 3 × 30s, 3 × 60s, 3 × 150s, and 9 or 15 × 300s) using a list-mode iterative reconstruction with proprietary spatially variant resolution modeling in 8 iterations and 16 subsets of the three-dimensional (3D) ordered subset expectation maximization (OSEM 3D) algorithm ( 46 ). Normalization, dead time, and CT-based attenuation corrections were applied.…”

Validation, kinetic modeling, and test-retest reproducibility of [18F]SynVesT-1 for PET imaging of synaptic vesicle glycoprotein 2A in mice

“…Images were reconstructed using a list-mode iterative reconstruction with proprietary spatially variant resolution modeling in 8 iterations and 16 subsets of the 3D ordered subset expectation maximization (OSEM 3D) algorithm. 28 Normalization, dead time, and CT-based attenuation corrections were applied. PET image frames were reconstructed on a 128 × 128 × 159 grid with 0.776 × 0.776 × 0.796 mm 3 voxels.…”

Spatiotemporal Kernel Reconstruction for Linear Parametric Neurotransmitter PET Kinetic Modeling in Motion Correction Brain PET of Awake Rats

“…The attenuation map was calculated using the binary image of the activity body outline with an uniform attenuation factor for soft tissue (0.096 cm –1 ) ( Angelis et al, 2013 ). Motion dependent and spatially variant resolution modeling was implemented as well ( Miranda et al, 2020 ). Dynamic images were reconstructed with a framing of 12 frames × 10 s, 6 ×20 s, 2 ×60 s, and 27 ×120 s.…”