Motion-related Artifacts Research Articles

Coronary CT angiography: Diagnostic value and clinical challenges

Coronary computed tomography (CT) angiography has been increasingly used in the diagnosis of coronary artery disease due to improved spatial and temporal resolution with high diagnostic value being reported when compared to invasive coronary angiography. Diagnostic performance of coronary CT angiography has been significantly improved with the technological developments in multislice CT scanners from the early generation of 4-slice CT to the latest 320- slice CT scanners. Despite the promising diagnostic value, coronary CT angiography is still limited in some areas, such as inferior temporal resolution, motion-related artifacts and high false positive results due to severe calcification. The aim of this review is to present an overview of the technical developments of multislice CT and diagnostic value of coronary CT angiography in coronary artery disease based on different generations of multislice CT scanners. Prognostic value of coronary CT angiography in coronary artery disease is also discussed, while limitations and challenges of coronary CT angiography are highlighted.

Distinct neural signatures detected for ADHD subtypes after controlling for micro-movements in resting state functional connectivity MRI data

In recent years, there has been growing enthusiasm that functional magnetic resonance imaging (MRI) could achieve clinical utility for a broad range of neuropsychiatric disorders. However, several barriers remain. For example, the acquisition of large-scale datasets capable of clarifying the marked heterogeneity that exists in psychiatric illnesses will need to be realized. In addition, there continues to be a need for the development of image processing and analysis methods capable of separating signal from artifact. As a prototypical hyperkinetic disorder, and movement-related artifact being a significant confound in functional imaging studies, ADHD offers a unique challenge. As part of the ADHD-200 Global Competition and this special edition of Frontiers, the ADHD-200 Consortium demonstrates the utility of an aggregate dataset pooled across five institutions in addressing these challenges. The work aimed to (1) examine the impact of emerging techniques for controlling for “micro-movements,” and (2) provide novel insights into the neural correlates of ADHD subtypes. Using support vector machine (SVM)-based multivariate pattern analysis (MVPA) we show that functional connectivity patterns in individuals are capable of differentiating the two most prominent ADHD subtypes. The application of graph-theory revealed that the Combined (ADHD-C) and Inattentive (ADHD-I) subtypes demonstrated some overlapping (particularly sensorimotor systems), but unique patterns of atypical connectivity. For ADHD-C, atypical connectivity was prominent in midline default network components, as well as insular cortex; in contrast, the ADHD-I group exhibited atypical patterns within the dlPFC regions and cerebellum. Systematic motion-related artifact was noted, and highlighted the need for stringent motion correction. Findings reported were robust to the specific motion correction strategy employed. These data suggest that resting-state functional connectivity MRI (rs-fcMRI) data can be used to characterize individual patients with ADHD and to identify neural distinctions underlying the clinical heterogeneity of ADHD.

Novel insight into the detailed myocardial motion and deformation of the rodent heart using high-resolution phase contrast cardiovascular magnetic resonance.

BackgroundPhase contrast velocimetry cardiovascular magnetic resonance (PC-CMR) is a powerful and versatile tool allowing assessment of in vivo motion of the myocardium. However, PC-CMR is sensitive to motion related artifacts causing errors that are geometrically systematic, rendering regional analysis of myocardial function challenging. The objective of this study was to establish an optimized PC-CMR method able to provide novel insight in the complex regional motion and strain of the rodent myocardium, and provide a proof-of-concept in normal and diseased rat hearts with higher temporal and spatial resolution than previously reported.MethodsA PC-CMR protocol optimized for assessing the motion and deformation of the myocardium in rats with high spatiotemporal resolution was established, and ten animals with different degree of cardiac dysfunction underwent examination and served as proof-of-concept. Global and regional myocardial velocities and circumferential strain were calculated, and the results were compared to five control animals. Furthermore, the global strain measurements were validated against speckle-tracking echocardiography, and inter- and intrastudy variability of the protocol were evaluated.ResultsThe presented method allows assessment of regional myocardial function in rats with high level of detail; temporal resolution was 3.2 ms, and analysis was done using 32 circumferential segments. In the dysfunctional hearts, global and regional function were distinctly altered, including reduced global peak values, increased regional heterogeneity and increased index of dyssynchrony. Strain derived from the PC-CMR data was in excellent agreement with echocardiography (r = 0.95, p < 0.001; limits-of-agreement −0.02 ± 3.92%strain), and intra- and interstudy variability were low for both velocity and strain (limits-of-agreement, radial motion: 0.01 ± 0.32 cm/s and −0.06 ± 0.75 cm/s; circumferential strain: -0.16 ± 0.89%strain and −0.71 ± 1.67%strain, for intra- and interstudy, respectively).ConclusionWe demonstrate, for the first time, that PC-CMR enables high-resolution evaluation of in vivo circumferential strain in addition to myocardial motion of the rat heart. In combination with the superior geometric robustness of CMR, this ultimately provides a tool for longitudinal studies of regional function in rodents with high level of detail.

Free‐breathing MOLLI: Application to myocardial T1 mapping

Of the myocardial T(1) mapping techniques, the modified Look-Locker inversion recovery (MOLLI) sequence is accurate and highly reproducible. The MOLLI sequence requires patients to hold their breath for 17 heartbeats during the scanning process to minimize respiratory motion-related artifacts. However, some patients are unable to hold their breath because of illness or limited breath-hold capacity. This study, therefore, aimed to develop a robust myocardial T(1) mapping method based on the MOLLI sequence for patients unable to perform voluntary breath-holds. This study presents a free-breathing MOLLI (FB-MOLLI) sequence and an optimized reconstruction method to allow myocardial T(1) mapping in vivo without breath-hold. Nine healthy volunteers participated in this study after providing institutionally approved consent. The FB-MOLLI sequence acquires 20 images within 29 heartbeats. The reconstruction program employs a two-step automatic image registration technique and an image selection method inspired by the self-gating cardiac imaging method. Results indicate that the proposed reconstruction method increases the accuracy and reproducibility of free-breathing T(1) measurements significantly (p < 0.001). The FB-MOLLI method provides a robust tool for clinical application in free-breathing myocardial T(1) mapping, and could greatly facilitate acquisition procedures during routine examinations.

High-pitch Dual-source Computed Tomography Pulmonary Angiography in Freely Breathing Patients

To investigate pulmonary arterial (PA) enhancement, image noise, and artifacts related to breathing and heart motion in patients with suspected pulmonary embolism. Seventy-six consecutive patients underwent computed tomographic pulmonary angiography (CTPA) in dual-source high-pitch mode (pitch 3.0, 100 kV, 180 mAs, 50 mL contrast material) without breathing commands. PA enhancement, image noise, signal to noise ratio, overall image quality, incidence of total or partial interruption of the contrast column in the PAs, and heart motion-related and breathing-related artifacts of the diaphragm and pulmonary structures were recorded. Mean central and peripheral PA attenuation was 404 ± 104 and 453 ± 119 HU; mean image noise was 11 ± 2 HU; mean examination time was 0.67 ± 0.09 s; and mean dose-length product was 142 ± 31 mGy cm. There were no motion artifacts of the diaphragm or pulmonary vessels related to breathing or heart motion. There was no case of partial or total interruption of the contrast column in the PA tree. No examination was rated nondiagnostic. High-pitch dual-source CTPA in freely breathing patients effectively produces images that are free of artifacts related to breathing and cardiac motion. Hence, Valsalva-related artifacts can be eliminated using this technique.

Preoperative imaging evaluation of head and neck cancer: Comparison of 2D spin-echo and 3D THRIVE MRI techniques with resected tumours

To evaluate the usefulness of contrast-enhanced three-dimensional (3D) T1-weighted high-resolution isotropic volume examination (THRIVE) for the preoperative assessment of head and neck cancer, by comparison with spin-echo (SE) T1-weighted sequences and the pathology specimen. Thirty-seven consecutive patients who were diagnosed with oral cavity, oropharyngeal, and hypopharyngeal cancer and received surgical excision with preoperative magnetic resonance imaging (MRI) using both SE and 3D THRIVE sequences after gadolinium injection were studied. Tumour conspicuity, motion- and flow-related artefacts were evaluated, and the tumour size was measured based on both sequences, which were correlated with the surgical specimen. The population correlation coefficient (r(g)(2)) and Bland-Altman plots were used to assess measurement agreement. There was no statistical difference between SE and THRIVE in terms of tumour conspicuity and motion-related artefacts, however, flow-related artefacts significantly decreased with THRIVE. The measurement agreement of tumour size estimated on both SE and THRIVE was good (r(g)(2) = 0.87 for SE and 0.93 for THRIVE), and the difference was statistically significant. The 3D THRIVE sequence provided satisfactory image quality for the accurate measurement of tumour size with fewer artefacts and could be an acceptable alternative for SE T1-weighted images for preoperative staging in patients with head and neck cancer.

Motion Correction of Multi-b-value Diffusion-weighted Imaging in the Liver

Motion artifacts are a significant source of error in the acquisition and quantification of parameters from multi-b-value diffusion-weighted imaging (DWI). The objective of this article is to present a reliable method to reduce motion-related artifacts during free-breathing at higher b-values when signal levels are low. Twelve patients referred for magnetic resonance imaging of the liver underwent a clinical magnetic resonance imaging examination of the abdominal region that included DWI. Conventional single-shot spin-echo echo planar imaging acquisitions of the liver during free breathing were repeated in a "time-resolved" manner during a single acquisition to obtain data for multi-b-value analysis, alternating between low and high b-values. Image registration using a normalized mutual information similarity measure was used to correct for spatial misalignment of diffusion-weighted volumes caused by motion. Registration error was estimated indirectly by comparing the normalized root-mean-square error (NRMSE) values of data fitted to the biexponential intra-voxel incoherent motion model before and after motion correction. Regions of interest (ROIs) were selected in the liver close to the surface of the liver and close to internal structures such as large bile ducts and blood vessels. For the 12 patient datasets, the mean NRMSE value for the motion-corrected ROIs (0.38 ± 0.16) was significantly lower than the mean NRMSE values for the non-motion-corrected ROIs (0.41 ± 0.13) (P < .05). In cases where there was substantial respiratory motion during the acquisition, visual inspection verified that the algorithm markedly improved alignment of the liver contours between frames. The proposed method addresses motion-related artifacts to increase robustness in multi-b-value acquisitions.

Optimization of PROPELLER reconstruction for free‐breathing T1‐weighted cardiac imaging

Clinical cardiac MR imaging techniques generally require patients to hold their breath during the scanning process to minimize respiratory motion-related artifacts. However, some patients cannot hold their breath because of illness or limited breath-hold capacity. This study aims to optimize the PROPELLER reconstruction for free-breathing myocardial T1-weighted imaging. Eight healthy volunteers (8 men; mean age 26.4 years) participated in this study after providing institutionally approved consent. The PROPELLER encoding method can reconstruct a low-resolution image from every blade because of k-space center oversampling. This study investigated the feasibility of extracting a respiratory trace from the PROPELLER blades by implementing a fully automatic region of interest selection and introducing a best template index to account for the property of the human respiration cycle. Results demonstrated that the proposed algorithm significantly improves the contrast-to-noise ratio and the image sharpness (p < 0.05). The PROPELLER method is expected to provide a robust tool for clinical application in free-breathing myocardial T1-weighted imaging. It could greatly facilitate the acquisition procedures during such a routine examination.

Live animal myelin histomorphometry of the spinal cord with video-rate multimodal nonlinear microendoscopy

In vivo imaging of cellular dynamics can be dramatically enabling to understand the pathophysiology of nervous system diseases. To fully exploit the power of this approach, the main challenges have been to minimize invasiveness and maximize the number of concurrent optical signals that can be combined to probe the interplay between multiple cellular processes. Label-free coherent anti-Stokes Raman scattering (CARS) microscopy, for example, can be used to follow demyelination in neurodegenerative diseases or after trauma, but myelin imaging alone is not sufficient to understand the complex sequence of events that leads to the appearance of lesions in the white matter. A commercially available microendoscope is used here to achieve minimally invasive, video-rate multimodal nonlinear imaging of cellular processes in live mouse spinal cord. The system allows for simultaneous CARS imaging of myelin sheaths and two-photon excitation fluorescence microendoscopy of microglial cells and axons. Morphometric data extraction at high spatial resolution is also described, with a technique for reducing motion-related imaging artifacts. Despite its small diameter, the microendoscope enables high speed multimodal imaging over wide areas of tissue, yet at resolution sufficient to quantify subtle differences in myelin thickness and microglial motility.

Steps toward optimizing motion artifact removal in functional connectivity MRI; a reply to Carp

We recently described a motion-related artifact in resting state functional connectivity MRI (rs-fcMRI) data that arises from subject head motion (Power et al. (2012), see also Satterthwaite et al. (2012); Van Dijk et al. (2012)). Head motion produces well-known disruptions in BOLD signal (Friston et al., 1996), and these artifactual modulations of BOLD signal, which are similar at nearby voxels, create spurious patterns in correlations in rs-fcMRI. Specifically, head motion augments short-distance correlations and weakens long-distance correlations. Thus, for instance, a higher-motion dataset would typically display weakened correlations between (distantly spaced) default mode regions but enhanced correlations between (closely spaced) visual regions in comparison to a low-motion data-set. We described this artifact in cohorts of children, adolescents, and adults, and its severity (magnitude) was related to the prevalence of motion within a cohort. Motion-related functional connectivity artifact is thus a substantial confound in the examination of single rs-fcMRI datasets and in comparisons of multiple rs-fcMRI datasets.

Myocardial Perfusion Magnetic Resonance Imaging Using Sliding-Window Conjugate-Gradient Highly Constrained Back-Projection Reconstruction for Detection of Coronary Artery Disease

Myocardial perfusion magnetic resonance imaging (MRI) with sliding-window conjugate-gradient highly constrained back-projection reconstruction (SW-CG-HYPR) allows whole left ventricular coverage, improved temporal and spatial resolution and signal/noise ratio, and reduced cardiac motion-related image artifacts. The accuracy of this technique for detecting coronary artery disease (CAD) has not been determined in a large number of patients. We prospectively evaluated the diagnostic performance of myocardial perfusion MRI with SW-CG-HYPR in patients with suspected CAD. A total of 50 consecutive patients who were scheduled for coronary angiography with suspected CAD underwent myocardial perfusion MRI with SW-CG-HYPR at 3.0 T. The perfusion defects were interpreted qualitatively by 2 blinded observers and were correlated with x-ray angiographic stenoses ≥50%. The prevalence of CAD was 56%. In the per-patient analysis, the sensitivity, specificity, positive predictive value, negative predictive value, and accuracy of SW-CG-HYPR was 96% (95% confidence interval 82% to 100%), 82% (95% confidence interval 60% to 95%), 87% (95% confidence interval 70% to 96%), 95% (95% confidence interval 74% to100%), and 90% (95% confidence interval 82% to 98%), respectively. In the per-vessel analysis, the corresponding values were 98% (95% confidence interval 91% to 100%), 89% (95% confidence interval 80% to 94%), 86% (95% confidence interval 76% to 93%), 99% (95% confidence interval 93% to 100%), and 93% (95% confidence interval 89% to 97%), respectively. In conclusion, myocardial perfusion MRI using SW-CG-HYPR allows whole left ventricular coverage and high resolution and has high diagnostic accuracy in patients with suspected CAD.

Optimizing the order of operations for movement scrubbing: Comment on Power et al.

A recent study by Power and colleagues shows that BOLD artifacts induced by head movement can substantially alter patterns of resting-state functional connectivity and proposes a novel procedure for reducing these artifacts by deleting (or "scrubbing") movement-contaminated volumes. The authors acknowledge that this work is descriptive and not prescriptive, and note that future studies may refine the proposed scrubbing method. Nevertheless, it is worth pointing out that this method can be improved substantially by a single transposition in the order of operations. Temporal filtering is known to introduce ringing artifacts that emanate from sharp transitions in signal intensity. The method proposed in the target article applies temporal filtering before deleting contaminated volumes-in effect, spreading movement-related artifacts backwards and forwards in time, but deleting only the originally contaminated data. Using simulated data, we show that deleting and replacing contaminated volumes before temporal filtering removes a greater proportion of artifactual signal while retaining a greater proportion of the original data.

Hypoxia Imaging Using PET and SPECT: The Effects of Anesthetic and Carrier Gas on [64Cu]-ATSM, [99mTc]-HL91 and [18F]-FMISO Tumor Hypoxia Accumulation

BackgroundPreclinical imaging requires anaesthesia to reduce motion-related artefacts. For direct translational relevance, anaesthesia must not significantly alter experimental outcome. This study reports on the effects of both anaesthetic and carrier gas upon the uptake of [64Cu]-CuATSM, [99mTc]-HL91 and [18F]-FMISO in a preclinical model of tumor hypoxia.Methodology/Principal FindingsThe effect of carrier gas and anaesthetic was studied in 6 groups of CaNT-bearing CBA mice using [64Cu]-CuATSM, [99mTc]-HL91 or [18F]-FMISO. Mice were anaesthetised with isoflurane in air, isoflurane in pure oxygen, with ketamine/xylazine or hypnorm/hypnovel whilst breathing air, or in the awake state whilst breathing air or pure oxygen. PET or SPECT imaging was performed after which the mice were killed for organ/tumor tracer quantitation. Tumor hypoxia was confirmed. Arterial blood gas analysis was performed for the different anaesthetic regimes. The results demonstrate marked influences on tumor uptake of both carrier gas and anaesthetic, and show differences between [99mTc]-HL91, [18F]-FMISO and [64Cu]-CuATSM. [99mTc]-HL91 tumor uptake was only altered significantly by administration of 100% oxygen. The latter was not the case for [18F]-FMISO and [64Cu]-CuATSM. Tumor-to-muscle ratio (TMR) for both compounds was reduced significantly when either oxygen or anaesthetics (isoflurane in air, ketamine/xylazine or hypnorm/hypnovel) were introduced. For [18F]-FMISO no further decrease was measured when both isoflurane and oxygen were administered, [64Cu]-CuATSM did show an additional significant decrease in TMR. When using the same anaesthetic regimes, the extent of TMR reduction was less pronounced for [64Cu]-CuATSM than for [18F]-FMISO (40–60% versus 70% reduction as compared to awake animals breathing air).Conclusions/SignificanceThe use of anaesthesia can have profound effects on the experimental outcome. More importantly, all tested anaesthetics reduced tumor-hypoxia uptake. Anaesthesia cannot be avoided in preclinical studies but great care has to be taken in preclinical models of hypoxia as anaesthesia effects cannot be generalised across applications, nor disease states.

Self-gated PROPELLER-encoded cine cardiac imaging

This study uses the k-space center over-sampling property of PROPELLER encoding to detect cardiac and respiratory motion using raw k-space data. Using the motion information, cine cardiac imaging is self-gated. The data acquisition process requires neither electrocardiography triggering nor patient breath-holding. Using physiology motion information, this post-processing method rearranges the k-space data into groups corresponding to cardiac and respiratory phases. The PROPELLER blades of the same groups are combined to produce cine high-resolution images. This approach reduces the potential discrepancy of the k-space data and the motion-related artifacts in the reconstructed image because all blades in a group are acquired at consistent cardiac and respiratory phases. The study concludes that self-gated cine cardiac imaging is feasible using PROPELLER encoding. It is a potentially practical tool for cine cardiac imaging.

The NIRS Analysis Package: noise reduction and statistical inference.

Near infrared spectroscopy (NIRS) is a non-invasive optical imaging technique that can be used to measure cortical hemodynamic responses to specific stimuli or tasks. While analyses of NIRS data are normally adapted from established fMRI techniques, there are nevertheless substantial differences between the two modalities. Here, we investigate the impact of NIRS-specific noise; e.g., systemic (physiological), motion-related artifacts, and serial autocorrelations, upon the validity of statistical inference within the framework of the general linear model. We present a comprehensive framework for noise reduction and statistical inference, which is custom-tailored to the noise characteristics of NIRS. These methods have been implemented in a public domain Matlab toolbox, the NIRS Analysis Package (NAP). Finally, we validate NAP using both simulated and actual data, showing marked improvement in the detection power and reliability of NIRS.

Feasibility of fast MR‐thermometry during cardiac radiofrequency ablation

Online MR temperature monitoring during radiofrequency (RF) ablation of cardiac arrhythmias may improve the efficacy and safety of the treatment. MR thermometry at 1.5 T using the proton resonance frequency (PRF) method was assessed in 10 healthy volunteers under normal breathing conditions, using a multi-slice, ECG-gated, echo planar imaging (EPI) sequence in combination with slice tracking. Temperature images were post-processed to remove residual motion-related artifacts. Using an MR-compatible steerable catheter and electromagnetic noise filter, RF ablation was performed in the ventricles of two sheep in vivo. The standard deviation of the temperature evolution in time (TSD) was computed. Temperature mapping of the left ventricle was achieved at an update rate of approximately 1 Hz with a mean TSD of 3.6 ± 0.9 °C. TSD measurements at the septum showed a higher precision (2.8 ± 0.9 °C) than at the myocardial regions at the heart-lung and heart-liver interfaces (4.1 ± 0.9 °C). Temperature rose maximally by 9 °C and 16 °C during 5 W and 10 W RF applications, respectively, for 60 s each. Tissue temperature can be monitored at an update rate of approximately 1 Hz in five slices. Typical temperature changes observed during clinical RF application can be monitored with an acceptable level of precision.

Functional MRI assessment of orofacial articulators: neural correlates of lip, jaw, larynx, and tongue movements.

Compared with complex coordinated orofacial actions, few neuroimaging studies have attempted to determine the shared and distinct neural substrates of supralaryngeal and laryngeal articulatory movements when performed independently. To determine cortical and subcortical regions associated with supralaryngeal motor control, participants produced lip, tongue and jaw movements while undergoing functional magnetic resonance imaging (fMRI). For laryngeal motor activity, participants produced the steady-state/i/vowel. A sparse temporal sampling acquisition method was used to minimize movement-related artifacts. Three main findings were observed. First, the four tasks activated a set of largely overlapping, common brain areas: the sensorimotor and premotor cortices, the right inferior frontal gyrus, the supplementary motor area, the left parietal operculum and the adjacent inferior parietal lobule, the basal ganglia and the cerebellum. Second, differences between tasks were restricted to the bilateral auditory cortices and to the left ventrolateral sensorimotor cortex, with greater signal intensity for vowel vocalization. Finally, a dorso-ventral somatotopic organization of lip, jaw, vocalic/laryngeal, and tongue movements was observed within the primary motor and somatosensory cortices using individual region-of-interest (ROI) analyses. These results provide evidence for a core neural network involved in laryngeal and supralaryngeal motor control and further refine the sensorimotor somatotopic organization of orofacial articulators.

Signal quality and Bayesian signal processing in neurofeedback based on real-time fMRI

Real-time fMRI allows analysis and visualization of the brain activity online, i.e. within one repetition time. It can be used in neurofeedback applications where subjects attempt to control an activation level in a specified region of interest (ROI) of their brain. The signal derived from the ROI is contaminated with noise and artifacts, namely with physiological noise from breathing and heart beat, scanner drift, motion-related artifacts and measurement noise. We developed a Bayesian approach to reduce noise and to remove artifacts in real-time using a modified Kalman filter. The system performs several signal processing operations: subtraction of constant and low-frequency signal components, spike removal and signal smoothing. Quantitative feedback signal quality analysis was used to estimate the quality of the neurofeedback time series and performance of the applied signal processing on different ROIs. The signal-to-noise ratio (SNR) across the entire time series and the group event-related SNR (eSNR) were significantly higher for the processed time series in comparison to the raw data. Applied signal processing improved the t-statistic increasing the significance of blood oxygen level-dependent (BOLD) signal changes. Accordingly, the contrast-to-noise ratio (CNR) of the feedback time series was improved as well. In addition, the data revealed increase of localized self-control across feedback sessions. The new signal processing approach provided reliable neurofeedback, performed precise artifacts removal, reduced noise, and required minimal manual adjustments of parameters. Advanced and fast online signal processing algorithms considerably increased the quality as well as the information content of the control signal which in turn resulted in higher contingency in the neurofeedback loop.

Motion-related artefacts in EEG predict neuronally plausible patterns of activation in fMRI data

The simultaneous acquisition and subsequent analysis of EEG and fMRI data is challenging owing to increased noise levels in the EEG data. A common method to integrate data from these two modalities is to use aspects of the EEG data, such as the amplitudes of event-related potentials (ERP) or oscillatory EEG activity, to predict fluctuations in the fMRI data. However, this relies on the acquisition of high quality datasets to ensure that only the correlates of neuronal activity are being studied. In this study, we investigate the effects of head-motion-related artefacts in the EEG signal on the predicted T2*-weighted signal variation. We apply our analyses to two independent datasets: 1) four participants were asked to move their feet in the scanner to generate small head movements, and 2) four participants performed an episodic memory task. We created T2*-weighted signal predictors from indicators of abrupt head motion using derivatives of the realignment parameters, from visually detected artefacts in the EEG as well as from three EEG frequency bands (theta, alpha and beta). In both datasets, we found little correlation between the T2*-weighted signal and EEG predictors that were not convolved with the canonical haemodynamic response function (cHRF). However, all convolved EEG predictors strongly correlated with the T2*-weighted signal variation in various regions including the bilateral superior temporal cortex, supplementary motor area, medial parietal cortex and cerebellum. The finding that movement onset spikes in the EEG predict T2*-weighted signal intensity only when the time course of movements is convolved with the cHRF, suggests that the correlated signal might reflect a BOLD response to neural activity associated with head movement. Furthermore, the observation that broad-spectral EEG spikes tend to occur at the same time as abrupt head movements, together with the finding that abrupt movements and EEG spikes show similar correlations with the T2*-weighted signal, indicates that the EEG spikes are produced by abrupt movement and that continuous regressors of EEG oscillations contain motion-related noise even after stringent correction of the EEG data. If not properly removed, these artefacts complicate the use of EEG data as a predictor of T2*-weighted signal variation.

MR-Guided Thermotherapy of Abdominal Organs Using a Robust PCA-Based Motion Descriptor

Thermotherapies can now be guided in real-time using magnetic resonance imaging (MRI). This technique is rapidly gaining importance in interventional therapies for abdominal organs such as liver and kidney. An accurate online estimation and characterization of organ displacement is mandatory to prevent misregistration and correct for motion related thermometry artifacts. In addition, when the ablation is performed with an extracorporal heating device such as high intensity focused ultrasound (HIFU), the continuous estimation of the organ displacement is the basis for the dynamic adjustment of the focal point position to track the targeted pathological tissue. In this paper, we describe the use of an optimized principal component analysis (PCA)-based motion descriptor to characterize in real-time the complex organ deformation during the therapy. The PCA was used to detect, in a preparative learning step, spatio-temporal coherences in the motion of the targeted organ. During hyperthermia, incoherent motion patterns could be discarded, which enabled improvements in motion estimation robustness, the compensation of motion related errors in thermal maps, and the adjustment of the beam position. The suggested method was evaluated for a moving phantom, and tested in vivo in the kidney and the liver of 12 healthy volunteers under free breathing conditions. The ability to perform a MR-guided thermotherapy in vivo during HIFU intervention was finally demonstrated on a porcine kidney.