Accuracy of motion correction methods for PET brain imaging

Motion correction of neurological 3D PET data by the multiple acquisition frame (MAF) method presents some challenges that do not apply to the same extent in 2D. One of the most significant is that 3D PET scans have a much higher scatter component, and the inclusion of valid correction for scatter is essential for accurate quantification. This paper describes a technique for motion correction of PET data acquired in 3D that uses information provided by an optical motion tracking system, and includes corrections for temporal variations in attenuation and scatter. The technique is demonstrated in 3D dynamic and list mode scans of a Hoffman 3D brain phantom performed on a Siemens/CTI ECAT EXACT HR+ scanner in which multiple arbitrary movements were applied to the phantom during data acquisition. The MAF method cannot correct for motion that occurs during the acquisition of an individual frame. The feasibility of avoiding this limitation by acquiring the data in list mode, and re-framing into a dynamic scan in which a new frame is commenced whenever significant motion occurs was investigated. Both approaches provided a major reduction in motion distortion when compared to a conventional reconstruction without motion correction. It appears that effective motion compensation is achievable in 3D neurological PET.

Advances in PET technology have led to improved resolution capability in modern PET scanners, andhence the impact of motion has become increasingly significant in order to produce high resolutionPET images. This thesis aims to advance the utility of motion correction in combined PET and MR byimproving the capabilities of motion correction, and by easing the effort required in order to applythem. In particular, PET-based data-driven motion tracking methods, which are able to track motionwithout requiring additional hardware, are investigated and applied to the head. First, existing techniques for motion correction, especially in the context of PET/MR, are reviewedand a general framework for PET motion correction is presented. The framework provides a top-levelperspective on PET motion correction, simplifying the process of comparing and contrasting variousmotion correction techniques. Next, the implementation of a patient-based orientation system withinthe open-source PET reconstruction packages the Software for Tomographic Image Reconstruction(STIR) and the Synergistic Image Reconstruction Software (SIRF) is described. Additionally, animplementation of the general PET motion correction framework platform is introduced. To investigate current technological limitations in combined PET and MR motion correction, theefficacy of a commercial PET/MR motion correction product, Brain COMPASS (herein just COMPASSfor brevity), was assessed. An experiment introduced a methodology for measuring motion blur invivo and uses this to measure reduction in motion blur after applying COMPASS. N=10 subjects worea head-mounted apparatus with fiducial markers affixed. COMPASS reduced blur in the directionof motion, measured by the principal component of the full-width at half-maximum (FWHM) ofthe source point-spread-function (PSF) distribution, only when the motion was sufficient to causethat blur to exceed 1.9 mm. However, blur across all directions almost always increased (mean=0.3mm) when COMPASS was used. Results suggested that any error in estimating the patient posewas low and was not the cause of the blur. Increases in overall blur were therefore suspected to becaused by interpolation when resampling transformed reconstructions. A second experiment measuredneocortical standardised uptake value ratio (SUVr) in an 18F-Florbetaben (18F-FBB) study of N=42young, healthy, asymptomatic volunteers with high polygenetic risk of developing Alzheimer’s Disease.No significant changes in neocortical SUVr were observed. This cohort exhibited low motion, andthe neocortex is a large region-of-interest, so is less sensitive to motion corruption than other cohorts.Hence, it was concluded that COMPASS did not aid this investigation, but also does not have adiscernible negative impact despite low motion. During the COMPASS investigation, MR acquisitions were the time-limiting factor of the session,meaning that the MR acquisitions introduced for PET motion correction increased the overall sessiontime. Because MR is not always available for motion tracking, alternative PET-alone motion correctiontechniques were investigated. Several data-driven tracking techniques are available; however, it is notclear which are most appropriate. A phantom investigation was conducted on a Biograph mCT Flowto determine which data-driven head motion tracking techniques were most appropriate. Of the fourmethodologies investigated principal component analysis (PCA) was found to be the most promisingtechnique for PET data-driven motion tracking of head motion. The PET-alone motion correction technique outlined above had two key limitations. Firstly, theyprovide frame-based motion correction with low temporal resolution. However, the type of motionseen in the COMPASS investigation were continuous “drooping” of the head, is not modelled in adiscrete, frame-based scheme. Secondly, the PET reconstructions without attenuation provide noisypose estimates compared to those from MR volume registration. Therefore, a scheme was proposedthat, unlike previous techniques, models head motion directly from motion surrogates. Here, a radialbasis function (RBF) regression relates the endogenous tracking signal to the exogenous rigid poseparameters, which are supplied by temporal-sparse MR acquisitions. A proof-of-concept is presentedin which low temporal resolution MR pose estimates are upsampled. The technique is compared tolinear and nearest-neighbour interpolation, and for the regression, both linear and RBF are investigated.The proposed technique was found to introduce some negligible blurring when no motion is present (<0.1 mm mean-of-maximum displacement (MMD)) but reduce blurring by 0.7 mm MMD when motionis present. This demonstrates the capability of the technique to reduce error in motion estimation byimproving the temporal resolution. In a clinical setting, motion correction is likely to become increasingly important in the future. Thework in this thesis will make motion correction easier to apply in a clinical setting, so that it may beapplied more regularly, providing wider benefit.

PurposePatient motion during PET acquisition may affect measured time-activity curves, thereby reducing accuracy of tracer kinetic analyses. The aim of the present study was to evaluate different off-line frame-by-frame methods to correct patient motion, which is of particular interest when no optical motion tracking system is available or when older data sets have to be reanalysed.MethodsFour different motion correction methods were evaluated. In the first method attenuation-corrected frames were realigned with the summed image of the first 3 min. The second method was identical, except that non-attenuation-corrected images were used. In the third and fourth methods non-attenuation-corrected images were realigned with standard and cupped transmission images, respectively. Two simulation studies were performed, based on [11C]flumazenil and (R)-[11C]PK11195 data sets, respectively. For both simulation studies different types (rotational, translational) and degrees of motion were added. Simulated PET scans were corrected for motion using all correction methods. The optimal method derived from these simulation studies was used to evaluate two (one with and one without visible movement) clinical data sets of [11C]flumazenil, (R)-[11C]PK11195 and [11C]PIB. For these clinical data sets, the volume of distribution (VT) was derived using Logan analysis and values were compared before and after motion correction.ResultsFor both [11C]flumazenil and (R)-[11C]PK11195 simulation studies, optimal results were obtained when realignment was based on non-attenuation-corrected images. For the clinical data sets motion disappeared visually after motion correction. Regional differences of up to 433% in VT before and after motion correction were found for scans with visible movement. On the other hand, when no visual motion was present in the original data set, overall differences in VT before and after motion correction were <1.5 ± 1.3%.ConclusionFrame-by-frame motion correction using non-attenuation-corrected images improves the accuracy of tracer kinetic analysis compared to non-motion-corrected data.Electronic supplementary materialThe online version of this article (doi:10.1007/s00259-009-1193-y) contains supplementary material, which is available to authorised users.

Subject motion in positron emission tomography (PET) studies leads to image blurring and artifacts; simultaneously acquired magnetic resonance imaging (MRI) data provides a means for motion correction (MC) in integrated PET/MRI scanners. To assess the effect of realistic head motion and MR-based MC on static [18 F]-fluorodeoxyglucose (FDG) PET images in dementia patients. Observational study. Thirty dementia subjects were recruited. 3T hybrid PET/MR scanner where EPI-based and T1 -weighted sequences were acquired simultaneously with the PET data. Head motion parameters estimated from high temporal resolution MR volumes were used for PET MC. The MR-based MC method was compared to PET frame-based MC methods in which motion parameters were estimated by coregistering 5-minute frames before and after accounting for the attenuation-emission mismatch. The relative changes in standardized uptake value ratios (SUVRs) between the PET volumes processed with the various MC methods, without MC, and the PET volumes with simulated motion were compared in relevant brain regions. The absolute value of the regional SUVR relative change was assessed with pairwise paired t-tests testing at the P = 0.05 level, comparing the values obtained through different MR-based MC processing methods as well as across different motion groups. The intraregion voxelwise variability of regional SUVRs obtained through different MR-based MC processing methods was also assessed with pairwise paired t-tests testing at the P = 0.05 level. MC had a greater impact on PET data quantification in subjects with larger amplitude motion (higher than 18% in the medial orbitofrontal cortex) and greater changes were generally observed for the MR-based MC method compared to the frame-based methods. Furthermore, a mean relative change of ∼4% was observed after MC even at the group level, suggesting the importance of routinely applying this correction. The intraregion voxelwise variability of regional SUVRs was also decreased using MR-based MC. All comparisons were significant at the P = 0.05 level. Incorporating temporally correlated MR data to account for intraframe motion has a positive impact on the FDG PET image quality and data quantification in dementia patients. 3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:1288-1296.

BackgroundRespiratory motion in PET/CT leads to well-known image degrading effects commonly compensated using elastic motion correction approaches. Gate-to-gate motion correction techniques are promising tools for improving clinical PET data but suffer from relatively long reconstruction times. In this study, the performance of a fast elastic motion compensation approach based on motion deblurring (DEB-MC) was evaluated on patient and phantom data and compared to an EM-based fully 3D gate-to-gate motion correction method (G2G-MC) which was considered the gold standard.MethodsTwenty-eight patients were included in this study with suspected or confirmed malignancies in the thorax or abdomen. All patients underwent whole-body [18F]FDG PET/CT examinations applying hardware-based respiratory gating. In addition, a dynamic anthropomorphic thorax phantom was studied with PET/CT simulating tumour motion under controlled but realistic conditions. PET signal recovery values were calculated from phantom scans by comparing lesion activities after motion correction to static ground truth data. Differences in standardized uptake values (SUV) and metabolic volume (MV) between both reconstruction methods as well as between motion-corrected (MC) and non motion-corrected (NOMC) results were statistically analyzed using a Wilcoxon signed-rank test.ResultsPhantom data analysis showed high lesion recovery values of 91% (2 cm motion) and 98% (1 cm) for G2G-MC and 83% (2 cm) and 90% (1 cm) for DEB-MC. The statistical analysis of patient data found significant differences between NOMC and MC reconstructions for SUV max, SUV mean, MV, and contrast-to-noise ratio (CNR) for both reconstruction algorithms. Furthermore, both methods showed similar increases of 11–12% in SUV max and SUV mean after MC. The statistical analysis of the MC/NOMC ratio found no significant differences between the methods.ConclusionBoth motion correction techniques deliver comparable improvements of SUV max, SUV mean, and CNR after MC on clinical and phantom data. The fast elastic motion compensation technique DEB-MC may thereby be a valuable alternative to state-of-the art motion correction techniques.

Patient motion during dynamic PET imaging can induce errors in myocardial blood flow (MBF) estimation. Motion correction for dynamic cardiac PET is challenging because the rapid tracer kinetics of 82Rb leads to substantial tracer distribution change across different dynamic frames over time, which can cause difficulties for image registration-based motion correction, particularly for early dynamic frames. In this paper, we developed an automatic deep learning-based motion correction (DeepMC) method for dynamic cardiac PET. In this study we focused on the detection and correction of inter-frame rigid translational motion caused by voluntary body movement and pattern change of respiratory motion. A bidirectional-3D LSTM network was developed to fully utilize both local and nonlocal temporal information in the 4D dynamic image data for motion detection. The network was trained and evaluated over motion-free patient scans with simulated motion so that the motion ground-truths are available, where one million samples based on 65 patient scans were used in training, and 600 samples based on 20 patient scans were used in evaluation. The proposed method was also evaluated using additional 10 patient datasets with real motion. We demonstrated that the proposed DeepMC obtained superior performance compared to conventional registration-based methods and other convolutional neural networks (CNN), in terms of motion estimation and MBF quantification accuracy. Once trained, DeepMC is much faster than the registration-based methods and can be easily integrated into the clinical workflow. In the future work, additional investigation is needed to evaluate this approach in a clinical context with realistic patient motion.

Evaluation of functional magnetic resonance imaging (fMRI) as a reliable clinical imaging tool requires accurate assessment and correction of head motion artifacts. As the correction of bulk head motion is vital, the loss of signal strength from the confounding effect of head motion on spin magnetization may be an additional factor in activation analysis error. This study focuses on the evaluation and correction of the spin saturation artifact that occurs when parts of adjacent slices are selected due to changing head positions in single-shot multislice acquisitions. As a consequence of head movement, the acquired slices constituting a fMRI volume are no longer parallel to each other and the spin magnetization in fMRI voxels becomes dependent on head motion history. Motion corrections applying the same rigid motion estimates to all the slices in a volume may not be a reasonable approximation in cases where the magnitude of head motion exceeds a subvoxel range. For realistic ranges of motion in fMRI, an accurate estimate of rigid motion parameters for each echo planar imaging (EPI) slice is essential to correctly register voxel intensities. Previously we have implemented the map-slice-to-volume (MSV) motion correction method that maps each slice in a time series onto a reference anatomical volume, which proved to be effective in improving activation detection. To correctly evaluate the motion dependence of spin magnetization, each voxel is tracked with movement history that is available from MSV motion estimates. Relatively low in resolution, EPI voxels are composed of varying mixtures of white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF) and variations in the tissue composition give rise to voxel intensities that are functions of tissue T1 properties. We have developed a weighted-average spin saturation (WASS) correction method that can handle full rigid motion and account for the melange of different brain tissue isochromats at each EPI voxel location. We evaluated the effect of spin saturation artifacts and the performance of the WASS correction using simulated fMRI time series synthesized with known true activation, motion, and the associated spin saturation artifact. Two different ranges of head rotations, [-5,5] and [-2,2] deg, were introduced and the effect of the spin saturation artifact was quantified to show 18% and 13% reduction in activation detection rate, respectively. Following the MSV motion and WASS correction, results indicate that WASS correction can improve activation detection by 17% relative to MSV only correction.

Motion correction of echo-planar imaging (EPI) data used in functional MRI (fMRI) is an essential pre-processing step performed prior to statistical analysis. At ultra-high resolution fMRI, current requirements regarding translational and rotational motion may no longer be acceptable. This prompts the need for a systematic investigation of the effects of motion correction procedures with in-vivo fMRI data. Here we systematically evaluated the effect of retrospective motion correction with freely available fMRI analysis software packages (FSL, AFNI and SPM) on activation maps using fMRI data acquired with prospective motion detection, to identify and quantify confounding effects of retrospective motion correction, and to evaluate its dependence on spatial resolution and motion correction algorithms.Brain activation maps were obtained for two different resolutions, an ultra-high, i.e. 0.653 mm3, and a more widely used 2.03 mm3 isotropic resolutions at 7T. The EPI data were acquired using simultaneous non-image-based optical moiré phase tracking (MPT) of physical motion. The results showed that image-based motion detection, performed by SPM8 software package, may be erroneous in high-field fMRI data with partial brain coverage and can introduce spurious motion leading to false-positive and false-negative activation. Further analyses demonstrated that limited acquisition field of view has the dominant influence on the effect.

Motion correction in PET has become more important as system resolution has improved. The purpose of this study was to evaluate the accuracy of three motion compensation methods, event-by-event motion compensation with list-mode reconstruction (MOLAR), frame-based motion correction, and post-reconstruction image registration. Motion compensated image reconstructions were carried out with simulated HRRT data, using a range of motion information based on human motion data. ROI analyses in high contrast regions were performed to evaluate the accuracy of all the motion compensation methods, with particular attention to within-frame motion.Our study showed that MOLAR with list-mode based motion correction using accurate motion data can reliably correct for all reasonable head motions. Over all motions, the average ROI count was within 0.1±4.2% and 0.7±0.9% of the reference, no-motion value for two different ROIs. The location of the ROI centroid was found to be within 0.7±0.3mm of that of the reference image for the raphe nucleus. Frame-based motion compensation and post-reconstruction image registration were able to correct for small (<5mm), but the ROI intensity begins to deteriorate for medium motions (5-10mm), especially for small brain structures such as the raphe nucleus. For large (>10mm) motions, the average centroid locations of the raphe nucleus ROI had an offset error of 1.5±1.8mm and 1.8±1.8mm for each of the frame-based methods. For each frame-based method, the decrease in the average ROI intensity was 16.9±4.3% and 20.2±9.9% respectively for the raphe nucleus, and was 5.5±2.2% and 7.4±0.2% for putamen. Based on these data, we conclude that event-by-event based motion correction works accurately for all reasonable motions, whereas frame-based motion correction is accurate only when the within-frame motion is less than 10mm.

Degradation of image quality caused by cardiac and respiratory motions hampers the diagnostic quality of cardiac PET. It has been shown that improved diagnostic accuracy of myocardial defect can be achieved by tagged MR (tMR) based PET motion correction using simultaneous PET-MR. However, one major hurdle for the adoption of tMR-based PET motion correction in the PET-MR routine is the long acquisition time needed for the collection of fully sampled tMR data. In this work, the authors propose an accelerated tMR acquisition strategy using parallel imaging and/or compressed sensing and assess the impact on the tMR-based motion corrected PET using phantom and patient data. Fully sampled tMR data were acquired simultaneously with PET list-mode data on two simultaneous PET-MR scanners for a cardiac phantom and a patient. Parallel imaging and compressed sensing were retrospectively performed by GRAPPA and kt-FOCUSS algorithms with various acceleration factors. Motion fields were estimated using nonrigid B-spline image registration from both the accelerated and fully sampled tMR images. The motion fields were incorporated into a motion corrected ordered subset expectation maximization reconstruction algorithm with motion-dependent attenuation correction. Although tMR acceleration introduced image artifacts into the tMR images for both phantom and patient data, motion corrected PET images yielded similar image quality as those obtained using the fully sampled tMR images for low to moderate acceleration factors (<4). Quantitative analysis of myocardial defect contrast over ten independent noise realizations showed similar results. It was further observed that although the image quality of the motion corrected PET images deteriorates for high acceleration factors, the images were still superior to the images reconstructed without motion correction. Accelerated tMR images obtained with more than 4 times acceleration can still provide relatively accurate motion fields and yield tMR-based motion corrected PET images with similar image quality as those reconstructed using fully sampled tMR data. The reduction of tMR acquisition time makes it more compatible with routine clinical cardiac PET-MR studies.

Motion detection and correction for carotid MRI using a markerless optical system

The impact of a method for MR-based respiratory motion correction of PET data on lesion visibility and quantification in patients with oncologic findings in the lung was evaluated. Twenty patients with one or more lesions in the lung were included. Hybrid imaging was performed on an integrated PET/MR system using 18F-FDG as radiotracer. The standard thoracic imaging protocol was extended by a free-breathing self-gated acquisition of MR data for motion modelling. PET data was acquired simultaneously in list-mode for 5-10 mins. One experienced radiologist and one experienced nuclear medicine specialist evaluated and compared the post-processed data in consensus regarding lesion visibility (scores 1–4, 4 being best), image noise levels (scores 1–3, 3 being lowest noise), SUVmean and SUVmax. Motion-corrected (MoCo) images were additionally compared with gated images. Non-motion-corrected free-breathing data served as standard of reference in this study. Motion correction generally improved lesion visibility (3.19 ± 0.63) and noise ratings (2.95 ± 0.22) compared to uncorrected (2.81 ± 0.66 and 2.95 ± 0.22, respectively) or gated PET data (2.47 ± 0.93 and 1.30 ± 0.47, respectively). Furthermore, SUVs (mean and max) were compared for all methods to estimate their respective impact on the quantification. Deviations of SUVmax were smallest between the uncorrected and the MoCo lesion data (average increase of 9.1% of MoCo SUVs), while SUVmean agreed best for gated and MoCo reconstructions (MoCo SUVs increased by 1.2%). The studied method for MR-based respiratory motion correction of PET data combines increased lesion sharpness and improved lesion activity quantification with high signal-to-noise ratio in a clinical setting. In particular, the detection of small lesions in moving organs such as the lung and liver may thus be facilitated. These advantages justify the extension of the PET/MR imaging protocol by 5–10 minutes for motion correction.

To evaluate the feasibility of motion correction in glutamate chemical exchange saturation transfer (GluCEST) imaging, using a rat model of epileptic seizure. Epileptic seizure was induced in six male Wistar rats by intraperitoneal injection of kainic acid (KA). CEST data were obtained using a 7.0T Bruker MRI scanner before and 3h after KA injection. Retrospective motion correction was performed in CEST images using a gradient-based motion correction (GradMC) algorithm. GluCEST signals in the hippocampal regions were quantitatively evaluated with and without motion correction. Calculated GluCEST signals differed significantly between the pre-KA injection group, regardless of motion-correction implementation, and the post-KA injection group with motion correction (3.662 ± 1.393 % / 3.726 ± 1.982 % for pre-KA injection group with/without motion correction vs. 6.996 ± 1.684 % for post-KA injection group with motion correction; all P < 0.05). Our results clearly show that GradMC can be used in CEST imaging for efficient correction of seizure-like motion. The GradMC can be further implemented in various CEST imaging techniques to increase the accuracy of analysis.

Due to the limited temporal resolution and cardiac motion, coronary computed tomography angiography (CCTA) exam is one of the most challenging CT protocols which may require operating radiologist to apply additional phase adjustment or motion correction for image reconstruction. To evaluate image quality between automatic and manual CCTA reconstruction in a 0.25 second rotation time, 16 cm coverage, single-beat, CT scanner with automated phase selection and AI-assisted motion correction. CCTA exams of 535 consecutive patients were included. All exams were first reconstructed with an automatically selected phase. If there was an unacceptable motion artifact, a manual reconstruction process was performed by radiologists. Additionally, automatic image series which consist of auto-phase selection and a follow-up motion correction were reconstructed. For these two manual and automatic image series, a four-point Likert scale rating system was used to evaluate image quality of the coronary artery segment by two experienced radiologists, according to the 18-segment model. Fifty-one patients (9.5%) did not have satisfactory image quality after auto-phase selection. In these patients, the heart rate during scanning was higher (78.3±18.4 bpm) than in the remaining 484 patients (68.9±13.1 bpm). Overall, 734 out of the 918 vessel segments were identified for quality evaluation among 51 patients. Automatic and manual image series were rated as having average Likert scores of 3.48±0.62 and 3.32±0.67 (P < 0.001), respectively. Using a 0.25 second rotation speed, 16 cm z-coverage, CT scanner installed with an AI-assisted motion correction algorithm, the automatic image reconstruction with scanner equipped auto-phase-selection and motion correction algorithm outperforms manually controlled image reconstruction by radiologists. This suggests that the traditional CCTA exam reconstruction workflow could be altered allowing less radiologist involvement and becoming more efficient.

High-density neural devices are now offering the possibility to record from neuronal populations in vivo at unprecedented scale. However, the mechanical drifts often observed in these recordings are currently a major issue for “spike sorting,” an essential analysis step to identify the activity of single neurons from extracellular signals. Although several strategies have been proposed to compensate for such drifts, the lack of proper benchmarks makes it hard to assess the quality and effectiveness of motion correction. In this paper, we present a benchmark study to precisely and quantitatively evaluate the performance of several state-of-the-art motion correction algorithms introduced in the literature. Using simulated recordings with induced drifts, we dissect the origins of the errors performed while applying a motion correction algorithm as a preprocessing step in the spike sorting pipeline. We show how important it is to properly estimate the positions of the neurons from extracellular traces in order to correctly estimate the probe motion, compare several interpolation procedures, and highlight what are the current limits for motion correction approaches.