Optimising rigid motion compensation for small animal brain PET imaging

Super-resolution in brain positron emission tomography using a real-time motion capture system

Super-resolution (SR) is a methodology that seeks to improve image resolution by exploiting the increased spatial sampling information obtained from multiple acquisitions of the same target with accurately known sub-resolution shifts. This work aims to develop and evaluate an SR estimation framework for brain positron emission tomography (PET), taking advantage of a high-resolution infra-red tracking camera to measure shifts precisely and continuously. Moving phantoms and non-human primate (NHP) experiments were performed on a GE Discovery MI PET/CT scanner (GE Healthcare) using an NDI Polaris Vega (Northern Digital Inc), an external optical motion tracking device. To enable SR, a robust temporal and spatial calibration of the two devices was developed as well as a list-mode Ordered Subset Expectation Maximization PET reconstruction algorithm, incorporating the high-resolution tracking data from the Polaris Vega to correct motion for measured line of responses on an event-by-event basis. For both phantoms and NHP studies, the SR reconstruction method yielded PET images with visibly increased spatial resolution compared to standard static acquisitions, allowing improved visualization of small structures. Quantitative analysis in terms of SSIM, CNR and line profiles were conducted and validated our observations. The results demonstrate that SR can be achieved in brain PET by measuring target motion in real-time using a high-resolution infrared tracking camera.

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

The linear parametric neurotransmitter positron emission tomography (lp-ntPET) kinetic model can be used to detect transient changes (activation) in endogenous neurotransmitter levels. Preclinical PET scans in awake animals can be performed to investigate neurotransmitter transient changes. Here we use the spatiotemporal kernel reconstruction (Kernel) for noise reduction in dynamic PET, and lp-ntPET kinetic modeling. Kernel is adapted for motion correction reconstruction, applied in awake rat PET scans. We performed 2D rat brain phantom simulation using the ntPET model at 3 different noise levels. Data was reconstructed with independent frame reconstruction (IFR), IFR with HYPR denoising, and Kernel, and lp-ntPET kinetic parameters (k 2a : efflux rate, γ: activation magnitude, t d : activation onset time, and t p : activation peak time) were calculated. Additionally, significant activation magnitude (γ) difference with respect to a region with no activation (rest) was calculated. Finally, [11C]raclopride experiments were performed in anesthetized and awake rats, injecting cold raclopride at 20 min after scan start to simulate endogenous neurotransmitter release. For simulated data at the regional level, IFR coefficient of variation (COV) of k 2a , γ, t d and t p was reduced with HYPR denoising, but Kernel showed the lowest COV (2 fold reduction compared with IFR). At the pixel level the same trend is observed for k 2a , γ, t d and t p COV, but reduction is larger with Kernel compared with IFR (10-14 fold). Bias in γ with respect with noise-free values was additionally reduced using Kernel (difference of 292, 72.4, and -6.92% for IFR, IFR+KYPR, and Kernel, respectively). Significant difference in activation between the rest and active region could be detected at a simulated activation of 160% for IFR and IFR+HYPR, and of 120% for Kernel. In rat experiments, lp-ntPET parameters have better confidence intervals using Kernel. In the γ, and t d parametric maps, the striatum structure can be identified with Kernel but not with IFR. Striatum voxel-wise γ, t d and t p values have lower variability using Kernel compared with IFR and IFR+HYPR. The spatiotemporal kernel reconstruction adapted for motion correction reconstruction allows to improve lp-ntPET kinetic modeling noise in awake rat studies, as well as detection of subtle neurotransmitter activations.

An iterative image-based inter-frame motion compensation method for dynamic brain PET imaging

As a non-invasive imaging tool, positron emission tomography (PET) plays an important role in brain science and disease research. Dynamic acquisition is one way of brain PET imaging. Its wide application in clinical research has often been hindered by practical challenges, such as patient involuntary movement, which could degrade both image quality and the accuracy of the quantification. This is even more obvious in scans of patients with neurodegeneration or mental disorders. Conventional motion compensation methods were either based on images or raw measured data, were shown to be able to reduce the effect of motion on the image quality. As for a dynamic PET scan, motion compensation can be challenging as tracer kinetics and relatively high noise can be present in dynamic frames. In this work, we propose an image-based inter-frame motion compensation approach specifically designed for dynamic brain PET imaging. Our method has an iterative implementation that only requires reconstructed images, based on which the inter-frame subject movement can be estimated and compensated. The method utilized tracer-specific kinetic modelling and can deal with simple and complex movement patterns. The synthesized phantom study showed that the proposed method can compensate for the simulated motion in scans with¹⁸F-FDG,¹⁸F-Fallypride and¹⁸F-AV45. Fifteen dynamic¹⁸F-FDG patient scans with motion artifacts were also processed. The quality of the recovered image was superior to the one of the non-corrected images and the corrected images with other image-based methods. The proposed method enables retrospective image quality control for dynamic brain PET imaging, hence facilitating the applications of dynamic PET in clinics and research.

Translation of Preclinical PET Imaging Findings: Challenges and Motion Correction to Overcome the Confounding Effect of Anesthetics

Preclinical brain positron emission tomography (PET) in animals is performed using anesthesia to avoid movement during the PET scan. In contrast, brain PET scans in humans are typically performed in the awake subject. Anesthesia is therefore one of the principal limitations in the translation of preclinical brain PET to the clinic. This review summarizes the available literature supporting the confounding effect of anesthesia on several PET tracers for neuroscience in preclinical small animal scans. In a second part, we present the state-of-the-art methodologies to circumvent this limitation to increase the translational significance of preclinical research, with an emphasis on motion correction methods. Several motion tracking systems compatible with preclinical scanners have been developed, each one with its advantages and limitations. These systems and the novel experimental setups they can bring to preclinical brain PET research are reviewed here. While technical advances have been made in this field, and practical implementations have been demonstrated, the technique should become more readily available to research centers to allow for a wider adoption of the motion correction technique for brain research.

Motion estimation and correction in SPECT, PET and CT

Patient motion impacts single photon emission computed tomography (SPECT), positron emission tomography (PET) and X-ray computed tomography (CT) by giving rise to projection data inconsistencies that can manifest as reconstruction artifacts, thereby degrading image quality and compromising accurate image interpretation and quantification. Methods to estimate and correct for patient motion in SPECT, PET and CT have attracted considerable research effort over several decades. The aims of this effort have been two-fold: to estimate relevant motion fields characterizing the various forms of voluntary and involuntary motion; and to apply these motion fields within a modified reconstruction framework to obtain motion-corrected images. The aims of this review are to outline the motion problem in medical imaging and to critically review published methods for estimating and correcting for the relevant motion fields in clinical and preclinical SPECT, PET and CT. Despite many similarities in how motion is handled between these modalities, utility and applications vary based on differences in temporal and spatial resolution. Technical feasibility has been demonstrated in each modality for both rigid and non-rigid motion, but clinical feasibility remains an important target. There is considerable scope for further developments in motion estimation and correction, and particularly in data-driven methods that will aid clinical utility. State-of-the-art machine learning methods may have a unique role to play in this context.

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

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.

Body motion detection and correction in cardiac PET: Phantom and human studies

PURPOSE

Patient body motion during a cardiac positron emission tomography (PET) scan can severely degrade image quality. We propose and evaluate a novel method to detect, estimate, and correct body motion in cardiac PET.

METHODS

Our method consists of three key components: motion detection, motion estimation, and motion-compensated image reconstruction. For motion detection, we first divide PET list-mode data into 1-s bins and compute the center of mass (COM) of the coincidences' distribution in each bin. We then compute the covariance matrix within a 25-s sliding window over the COM signals inside the window. The sum of the eigenvalues of the covariance matrix is used to separate the list-mode data into "static" (i.e., body motion free) and "moving" (i.e. contaminated by body motion) frames. Each moving frame is further divided into a number of evenly spaced sub-frames (referred to as "sub-moving" frames), in which motion is assumed to be negligible. For motion estimation, we first reconstruct the data in each static and sub-moving frame using a rapid back-projection technique. We then select the longest static frame as the reference frame and estimate elastic motion transformations to the reference frame from all other static and sub-moving frames using nonrigid registration. For motion-compensated image reconstruction, we reconstruct all the list-mode data into a single image volume in the reference frame by incorporating the estimated motion transformations in the PET system matrix. We evaluated the performance of our approach in both phantom and human studies.

RESULTS

Visually, the motion-corrected (MC) PET images obtained using the proposed method have better quality and fewer motion artifacts than the images reconstructed without motion correction (NMC). Quantitative analysis indicates that MC yields higher myocardium to blood pool concentration ratios. MC also yields sharper myocardium than NMC.

CONCLUSIONS

The proposed body motion correction method improves image quality of cardiac PET.

Markerless motion tracking and correction for PET, MRI, and simultaneous PET/MRI

OBJECTIVE

We demonstrate and evaluate the first markerless motion tracker compatible with PET, MRI, and simultaneous PET/MRI systems for motion correction (MC) of brain imaging.

METHODS

PET and MRI compatibility is achieved by careful positioning of in-bore vision extenders and by placing all electronic components out-of-bore. The motion tracker is demonstrated in a clinical setup during a pediatric PET/MRI study including 94 pediatric patient scans. PET MC is presented for two of these scans using a customized version of the Multiple Acquisition Frame method. Prospective MC of MRI acquisition of two healthy subjects is demonstrated using a motion-aware MRI sequence. Real-time motion estimates are accompanied with a tracking validity parameter to improve tracking reliability.

RESULTS

For both modalities, MC shows that motion induced artifacts are noticeably reduced and that motion estimates are sufficiently accurate to capture motion ranging from small respiratory motion to large intentional motion. In the PET/MRI study, a time-activity curve analysis shows image improvements for a patient performing head movements corresponding to a tumor motion of ±5-10 mm with a 19% maximal difference in standardized uptake value before and after MC.

CONCLUSION

The first markerless motion tracker is successfully demonstrated for prospective MC in MRI and MC in PET with good tracking validity.

SIGNIFICANCE

As simultaneous PET/MRI systems have become available for clinical use, an increasing demand for accurate motion tracking and MC in PET/MRI scans has emerged. The presented markerless motion tracker facilitate this demand.

Estimation of and correction for finite motion sampling errors in small animal PET rigid motion correction

Motion tracking with finite time sampling causing an associated unknown residual motion between two motion measurements is one of the factors contributing to resolution loss in small animal PET motion correction. The aim of this work is (i) to provide a means to estimate the effect of the finite motion sampling on the spatial resolution of the motion correction reconstructions and (ii) to correct for this residual motion thereby minimizing resolution loss. We calculate a tailored spatially variant deconvolution kernel from the measured motion data which is then used to deconvolve the motion corrected image using a 3D Richardson-Lucy algorithm. A simulation experiment of numerical phantoms as well as a microDerenzo phantom experiment wherein the phantom was manually moved at different speeds was performed to assess the performance of our proposed method. In the motion corrected images of the microDerenzo phantom there was an average rod FWHM differences between the slow and fast motion cases of 9.7%. This difference was reduced to 5.8% after applying the residual motion deconvolution. In awake animal experiments, the proposed method can serve to mitigate the finite sampling factor degrading the spatial resolution as well as the resolution differences between fast-moving and slow-moving animals.

Graphical abstract

Fast and Accurate Rat Head Motion Tracking With Point Sources for Awake Brain PET

To avoid the confounding effects of anesthesia and immobilization stress in rat brain positron emission tomography (PET), motion tracking-based unrestrained awake rat brain imaging is being developed. In this paper, we propose a fast and accurate rat headmotion tracking method based on small PET point sources. PET point sources (3-4) attached to the rat's head are tracked in image space using 15-32-ms time frames. Our point source tracking (PST) method was validated using a manually moved microDerenzo phantom that was simultaneously tracked with an optical tracker (OT) for comparison. The PST method was further validated in three awake [¹⁸F]FDG rat brain scans. Compared with the OT, the PST-based correction at the same frame rate (31.2 Hz) reduced the reconstructed FWHM by 0.39-0.66 mm for the different tested rod sizes of the microDerenzo phantom. The FWHM could be further reduced by another 0.07-0.13 mm when increasing the PST frame rate (66.7 Hz). Regional brain [¹⁸F]FDG uptake in the motion corrected scan was strongly correlated ( ) with that of the anesthetized reference scan for all three cases ( ). The proposed PST method allowed excellent and reproducible motion correction in awake in vivo experiments. In addition, there is no need of specialized tracking equipment or additional calibrations to be performed, the point sources are practically imperceptible to the rat, and PST is ideally suitable for small bore scanners, where optical tracking might be challenging.

The effect of isoflurane on 18F-FDG uptake in the rat brain: a fully conscious dynamic PET study using motion compensation

BACKGROUND

In preclinical positron emission tomography (PET) studies an anaesthetic is used to ensure that the animal does not move during the scan. However, anaesthesia may have confounding effects on the drug or tracer kinetics under study, and the nature of these effects is usually not known.

METHOD

We have implemented a protocol for tracking the rigid motion of the head of a fully conscious rat during a PET scan and performing a motion compensated list-mode reconstruction of the data. Using this technique we have conducted eight rat studies to investigate the effect of isoflurane on the uptake of ¹⁸F-FDG in the brain, by comparing conscious and unconscious scans.

RESULTS

Our results indicate that isoflurane significantly decreases the whole brain uptake, as well as decreasing the relative regional FDG uptake in the cortex, diencephalon, and inferior colliculi, while increasing it in the vestibular nuclei. No statistically significant changes in FDG uptake were observed in the cerebellum and striata.

CONCLUSION

The applied event-based motion compensation technique allowed for the investigation of the effect of isoflurane on FDG uptake in the rat brain using fully awake and unrestrained rats, scanned dynamically from the moment of injection. A significant effect of the anaesthesia was observed in various regions of the brain.