MRI with TRELLIS: a novel approach to motion correction

A motion-correcting pulse sequence and reconstruction algorithm, termed TRELLIS, is presented. k-Space is filled using orthogonal overlapping strips and the directions for phase- and frequency-encoding are alternated such that the frequency-encode direction always runs lengthwise along each strip. The overlap between strips is used both for signal averaging and to produce a system of equations that, when solved, quantifies the rotational and translational motion of the object. Results obtained from simulations with computer-generated phantoms, a purpose-built moving phantom, and in human subjects show the method is effective. TRELLIS offers some advantages over existing techniques in that k-space is sampled uniformly and all acquired data are used for both motion detection and image reconstruction.     

Trajectory correction for free-breathing radial cine MRI

Radial acquisitions can suffer from trajectory errors leading to reduced image quality. Here we present a new method of trajectory correction that uses all spokes of a radial acquisition and compare it to an existing method that uses a two-spoke pre-scan calibration. For both methods, estimates of the necessary shifts were made using magnitude or phase data and the performances were compared. The additional effect of B0 correction was considered in all cases. Mouse cardiac scans were used for the comparisons and we also compared the quality of navigator signals obtained from the radial data with each technique.     

Retrospective intra-scan motion correction

This paper analyzes the effects of intra-scan motion and demonstrates the possibility of correcting them directly in k-space with a new automatic retrospective method. The method is presented for series of 2D acquisitions with Cartesian sampling. Using a reference k-space acquisition (corrected for translations) within the series, intra-scan motion parameters are accurately estimated for each trajectory in k-space of each data set in the series resulting in pseudo-random sample positions. The images are reconstructed with a Bayesian estimator that can handle sparse arbitrary sampling in k-space and reduces intra-scan rotation artefacts to the noise level. The method has been assessed by means of a Monte Carlo study on axial brain images for different signal-to-noise ratios. The accuracy of motion estimates is better than 0.1 degrees for rotation, and 0.1 and 0.05 pixel, respectively, for translation along the read and phase directions for signal-to-noise ratios higher than 6 of the signals on each trajectory. An example of reconstruction from experimental data corrupted by head motion is also given.     

Manifestation and post hoc correction of gradient cross-term artifacts in DTI

Cross-terms between imaging and diffusion gradients, unaccounted for during tensor calculations, can lead to erroneous estimation of diffusivity and fractional anisotropy (FA) in regions of isotropic and anisotropic diffusion. Cross-term of magnitude 136.8±1.6 s/mm(2), artificially introduced in the slice-encode direction, caused an increase in FA in isotropic phantom from 0.0546±0.0001 to 0.0996±0.0001, while the change in chimpanzee brain depended on the orientation of the white matter (WM). Mean diffusivity (MD) remained unchanged in isotropic phantom, but increased by ～20% in the WM due to cross-terms. A bias was observed in the principal eigenvectors in both phantom and chimpanzee brain, resulting in significant increase in midline crossing fibers along the bias than perpendicular to it in tractography in chimpanzee brain. Post hoc correction of these artifacts was achieved by estimating the cross-term factors using calibration scans on an isotropic phantom and modifying the b-matrix before tensor calculation. Upon correction, the FA and MD values closely resembled the values obtained from sequence without cross-terms, and the bias in principal eigenvectors was eliminated. Customized sequences involving large b-values, high-resolution imaging, or long diffusion or echo times should therefore be evaluated and any residual cross-terms corrected before implementation.     

Intra-voxel incoherent motion MRI in rodent model of diethylnitrosamine-induced liver fibrosis

Rationale and Objectives

To compare the apparent diffusion coefficient (ADC) and the perfusion fraction measured by intra-voxel incoherent motion (IVIM) Magnetic Resonance Imaging (MRI) with liver fibrosis degrees in a rodent model.

Materials and Methods

All experiments received approval from our institutional animal care and use committee. Liver fibrosis was induced in 13 rats by oral gavage with diethylnitrosamine; 4 untreated rats with normal livers were used as controls. Diffusion Weighted MRI was performed and 8 gradient factors (0, 50, 100, 150, 200, 300, 400 and 500 s/mm²) were acquired. The values of ADC, true diffusion coefficient D and perfusion fraction f were measured based on Li Bihan’s method. The percentage of liver fibrosis was assessed via quantitative analysis of Masson trichrome staining using an average of 30 fields per section. The MRI measurements were compared to the histological fibrotic grade to evaluate the correlation between them.

Results

ADC contained the contribution of diffusion and perfusion. The ADC and f values decreased significantly with the increasing fibrosis level (correlation coefficient: ADC: ρ = − 0.781, p < 0.001; f: ρ = − 0.720, p = 0.001); but D was poorly correlated with fibrosis level (ρ = − 0.502, p = 0.040).

Conclusion

The hepatic ADC and the perfusion fraction f were significantly correlated with the liver fibrosis level; however, D was not. This might suggest that hepatic perfusion is altered during the progression of hepatic fibrosis.     

Characteristics of geometric distortion correction with increasing field-of-view in open-configuration MRI

Open-configuration magnetic resonance imaging (MRI) systems are becoming increasingly desirable for musculoskeletal imaging and image-guided radiotherapy because of their non-claustrophobic configuration. However, geometric image distortion in large fields-of-view (FOV) due to field inhomogeneity and gradient nonlinearity hinders the practical applications of open-type MRI. We demonstrated the use of geometric distortion correction for increasing FOV in open MRI. Geometric distortion was modeled and corrected as a global polynomial function. The appropriate polynomial order was identified as the minimum difference between the coordinates of control points in the distorted MR image space and those predicted by polynomial modeling. The sixth order polynomial function was found to give the optimal value for geometric distortion correction. The area of maximum distortion was < 1 pixel with an FOV of 285 mm. The correction performance error was increased at most 1.2% and 2.9% for FOVs of 340 mm and ~ 400 mm compared with the FOV of 285 mm. In particular, unresolved distortion was generated by local deformation near the gradient coil center.     

Improving cerebral blood flow quantification for arterial spin labeled perfusion MRI by removing residual motion artifacts and global signal fluctuations

Denoising is critical to improving the quality and stability of cerebral blood flow (CBF) quantification in arterial spin labeled (ASL) perfusion magnetic resonance imaging (MRI) due to the intrinsic low signal-to-noise-ratio (SNR) of ASL data. Previous studies have been focused on reducing the spatial or temporal noise using standard filtering techniques, and less attention has been paid to two global nuisance effects, the residual motion artifacts and the global signal fluctuations. Since both nuisances affect the whole brain, removing them in advance should enhance the CBF quantification quality for ASL MRI. The purpose of this paper was to assess this potential benefit. Three methods were proposed to suppress each or both of the two global nuisances. Their performances for CBF quantification were validated using ASL data acquired from 13 subjects. Evaluation results showed that covarying out both global nuisances significantly improved temporal SNR and test-retest stability of CBF measurement. Although the concept of removing both nuisances is not technically novel per se, this paper clearly showed the benefits for ASL CBF quantification. Dissemination of the proposed methods in a free ASL data processing toolbox should be of interest to a broad range of ASL users.     

Elimination of motion and pulsation artifacts using BLADE sequences in knee MR imaging

The purpose of this study is to evaluate the ability of proton density (PD)-BLADE sequences in reducing or even eliminating motion and pulsatile flow artifacts in knee magnetic resonance imaging examinations. Eighty consecutive patients, who had been routinely scanned for knee examination, participated in the study. The following pairs of sequences with and without BLADE were compared: (a) PD turbo spin echo (TSE) sagittal (SAG) fat saturation (FS) in 35 patients, (b) PD TSE coronal (COR) FS in 19 patients, (c) T2 TSE axial in 13 patients and (d) PD TSE SAG in 13 patients. Both qualitative and quantitative analyses were performed based on the signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and relative contrast (ReCon) measures of normal anatomic structures. The qualitative analysis was performed by experienced radiologists. Also, the presence of image motion and pulsation artifacts was evaluated. Based on the results of the SNR, CRN and ReCon for the different sequences and anatomical structures, the BLADE sequences were significantly superior in 19 cases, whereas the corresponding conventional sequences were significantly superior in only 6 cases. BLADE sequences eliminated motion artifacts in all the cases. However, motion artifacts were shown in (a) six PD TSE SAG FS, (b) three PD TSE COR FS, (c) three PD TSE SAG and (d) two T2 TSE axial conventional sequences. In our results, it was found that, in PD FS sequences (sagittal and coronal), the differences between the BLADE and conventional sequences regarding the elimination of motion and pulsatile flow artifacts were statistically significant. In all the comparisons, the PD FS BLADE sequences (coronal and sagittal) were significantly superior to the corresponding conventional sequences regarding the classification of their image quality. In conclusion, this technique appears to be capable to potentially eliminate motion and pulsatile flow artifacts in MR images.     

Concurrent correction of geometric distortion and motion using the map-slice-to-volume method in echo-planar imaging

The accuracy of measuring voxel intensity changes between stimulus and rest images in fMRI echo-planar imaging (EPI) data is severely degraded in the presence of head motion. In addition, EPI is sensitive to susceptibility-induced geometric distortions. Head motion causes image shifts and associated field map changes that induce different geometric distortion at different time points. Conventionally, geometric distortion is "corrected" with a static field map independently of image registration. That approach ignores all field map changes induced by head motion. This work evaluates the improved motion correction capability of mapping slice to volume with concurrent iterative field corrected reconstruction using updated field maps derived from an initial static field map that has been spatially transformed and resampled. It accounts for motion-induced field map changes for translational and in-plane rotation motion. The results from simulated EPI time series data, in which motion, image intensity and activation ground truths are available, show improved accuracy in image registration, field corrected image reconstruction and activation detection.     

An investigation of motion correction algorithms for pediatric spinal cord DTI in healthy subjects and patients with spinal cord injury

Patient and physiological motion can cause artifacts in DTI of the spinal cord which can impact image quality and diffusion indices. The purpose of this investigation was to determine a reliable motion correction method for pediatric spinal cord DTI and show effects of motion correction on DTI parameters in healthy subjects and patients with spinal cord injury. Ten healthy subjects and ten subjects with spinal cord injury were scanned using a 3 T scanner. Images were acquired with an inner field-of-view DTI sequence covering cervical spine levels C1 to C7. Images were corrected for motion using two types of transformation (rigid and affine) and three cost functions. Corrected images and transformations were examined qualitatively and quantitatively using in-house developed code. Fractional anisotropy (FA) and mean diffusivity (MD) indices were calculated and tested for statistical significance pre- and post- motion correction. Images corrected using rigid methods showed improvements in image quality, while affine methods frequently showed residual distortions in corrected images. Blinded evaluation of pre and post correction images showed significant improvement in cord homogeneity and edge conspicuity in corrected images (p < 0.0001). The average FA changes were statistically significant (p < 0.0001) in the spinal cord injury group, while healthy subjects showed less FA change and were not significant. In both healthy subjects and subjects with spinal cord injury, quantitative and qualitative analysis showed the rigid scaled-least-squares registration technique to be the most reliable and effective in improving image quality.     

Quantitative assessment of motion correction for high angular resolution diffusion imaging

Several methods have been proposed for motion correction of high angular resolution diffusion imaging (HARDI) data. There have been few comparisons of these methods, partly due to a lack of quantitative metrics of performance. We compare two motion correction strategies using two figures of merit: displacement introduced by the motion correction and the 95% confidence interval of the cone of uncertainty of voxels with prolate tensors. What follows is a general approach for assessing motion correction of HARDI data that may have broad application for quality assurance and optimization of postprocessing protocols. Our analysis demonstrates two important issues related to motion correction of HARDI data: (1) although neither method we tested was dramatically superior in performance, both were dramatically better than performing no motion correction, and (2) iteration of motion correction can improve the final results. Based on the results demonstrated here, iterative motion correction is strongly recommended for HARDI acquisitions.     

Usefulness of slice encoding for metal artifact correction (SEMAC) for reducing metallic artifacts in 3-T MRI

Purpose

The purpose of the study was to assess the usefulness of slice encoding for metal artifact correction (SEMAC) in 3.0-T magnetic resonance (MR) in minimizing metallic artifacts in patients with spinal prostheses.

Materials and Methods

Institutional review board approval and informed consent were obtained for this study. Twenty-seven spine MR scans were performed with metal artifact reduction SEMAC between May 2011 and July 2012 in patients with metallic devices. The MR scans were performed on a 3-T MR system (Achieva; Philips Healthcare, Best, the Netherlands) including SEMAC-corrected T2-weighted axial/sagittal images and two-dimensional fast spin echo (FSE) axial/sagittal images. The SEMAC-corrected images were compared to conventional T2-weighted FSE images. Two musculoskeletal radiologists qualitatively analyzed the images in terms of visualization of the pedicle, vertebral body, dural sac, intervertebral disc, intervertebral neural foramina, screws and metallic artifacts. The paired images were rated using a 5-point scale. P values less than .05 were considered to indicate statistically significant differences.

Results

The SEMAC-corrected MR images significantly reduced the metal-related artifacts. The T2-weighted images with SEMAC sequences enabled significantly improved periprosthetic visualizations of the pedicle, vertebral body, dural sac and neural foramina, with the exception of the intervertebral disc (P < .05). In addition, there was significant improvement in prosthesis visualization (P < .05).

Conclusion

MR images with SEMAC can reduce metal-related artifacts, providing improved delineation of the prosthesis and periprosthetic region. However, for the evaluation of the intervertebral disc, the SEMAC-corrected MR images showed no significant benefits.     

Elimination of motion,pulsatile flow and cross-talk artifacts using blade sequences in lumbar spine MR imaging

The purpose of this study is to evaluate the ability of T2 turbo spin echo (TSE) axial and sagittal BLADE sequences in reducing or even eliminating motion, pulsatile flow and cross-talk artifacts in lumbar spine MRI examinations. Forty four patients, who had routinely undergone a lumbar spine examination, participated in the study. The following pairs of sequences with and without BLADE were compared: a) T2 TSE Sagittal (SAG) in thirty two cases, and b) T2 TSE Axial (AX) also in thirty two cases. Both quantitative and qualitative analyses were performed based on measurements in different normal anatomical structures and examination of seven characteristics, respectively. The qualitative analysis was performed by experienced radiologists. Also, the presence of image motion, pulsatile flow and cross-talk artifacts was evaluated. Based on the results of the qualitative analysis for the different sequences and anatomical structures, the BLADE sequences were found to be significantly superior to the conventional ones in all the cases. The BLADE sequences eliminated the motion artifacts in all the cases. In our results, it was found that in the examined sequences (sagittal and axial) the differences between the BLADE and conventional sequences regarding the elimination of motion, pulsatile flow and cross-talk artifacts were statistically significant. In all the comparisons, the T2 TSE BLADE sequences were significantly superior to the corresponding conventional sequences regarding the classification of their image quality. In conclusion, this technique appears to be capable of potentially eliminating motion, pulsatile flow and cross-talk artifacts in lumbar spine MR images and producing high quality images in collaborative and non-collaborative patients.     

Ring artifacts correction based on the projection-field in neutron CT

Ring artifacts will happen mostly when the detector has inconsistent response among the detector channels, and the characteristic produced rings centered in the iso-center in the reconstructed slices inevitably affect the recognition and analysis of the corresponding sample structures in neutron computed tomography(CT). In this work, a ring correction method based on the projection-field(RCP) is proposed, it is a pre-processing method and provides the corrected projection data directly, which is also conducive to efficient data storage and other algorithmic researches. Simulation and physical experiments are performed for verifying the effect of the method, and one of the correction methods based on the imagefield is used for comparison. The results demonstrate that the RCP can correct the ring artifacts well without reducing the image resolution or over-correction.     

Evaluation of MRI artifacts at 3 Tesla for 38 commonly used cosmetics

Purpose

To evaluate MRI artifacts at 3-Tesla for 38 commonly used cosmetics.

Materials and Methods

Thirty-eight cosmetics (16, nail polishes; 5, eyeliners; 3, mascaras; 10, eye shadows; 1, lip gloss; 1, body lotion; 1, body glitter, and 1, hair loss concealer) underwent evaluation for MRI artifacts at 3-Tesla. The cosmetics were applied a copper-sulfate-filled, phantom and initially assessed using a “screening” gradient echo (GRE) pulse sequence. Of the 38 different cosmetics, 14 (37%) exhibited artifacts. For these 14 cosmetics, additional characterization of artifacts was performed using a GRE pulse sequence. A qualitative scale was applied to characterize the artifact size.

Results

Artifacts were observed, as follows: 2, nail polishes; 5, eyeliners; 3, mascaras; 3, eye shadows; 1, hair loss concealer. Artifact size ranged from small (eye shadow) to very large (hair loss concealer) and tended to be associated with the presence of iron oxide or other metal-based ingredient.

Conclusions

Commonly used cosmetics caused artifacts that may create issues if the area of interest is the same as where the cosmetic was applied or if its presence was unknown, thus, potentially causing it to be construed as pathology. Therefore, these findings have important implications for patients referred for MRI examinations.     

Beat-to-beat respiratory motion correction with near 100% efficiency: a quantitative assessment using high-resolution coronary artery imaging

This study quantitatively assesses the effectiveness of retrospective beat-to-beat respiratory motion correction (B2B-RMC) at near 100% efficiency using high-resolution coronary artery imaging. Three-dimensional (3D) spiral images were obtained in a coronary respiratory motion phantom with B2B-RMC and navigator gating. In vivo, targeted 3D coronary imaging was performed in 10 healthy subjects using B2B-RMC spiral and navigator gated balanced steady-state free-precession (nav-bSSFP) techniques. Vessel diameter and sharpness in proximal and mid arteries were used as a measure of respiratory motion compensation effectiveness and compared between techniques. Phantom acquisitions with B2B-RMC were sharper than those acquired with navigator gating (B2B-RMC vs. navigator gating: 1.01±0.02 mm⁻¹ vs. 0.86±0.08 mm⁻¹, P<.05). In vivo B2B-RMC respiratory efficiency was significantly and substantially higher (99.7%±0.5%) than nav-bSSFP (44.0%±8.9%, P<.0001). Proximal and mid vessel sharpnesses were similar (B2B-RMC vs. nav-bSSFP, proximal: 1.00±0.14 mm⁻¹ vs. 1.08±0.11 mm⁻¹, mid: 1.01±0.11 mm⁻¹ vs. 1.05±0.12 mm⁻¹; both P=not significant [ns]). Mid vessel diameters were not significantly different (2.85±0.39 mm vs. 2.80±0.35 mm, P=ns), but proximal B2B-RMC diameters were slightly higher (2.85±0.38 mm vs. 2.70±0.34 mm, P<.05), possibly due to contrast differences. The respiratory efficiency of B2B-RMC is less variable and significantly higher than navigator gating. Phantom and in vivo vessel sharpness and diameter values suggest that respiratory motion compensation is equally effective.     

Correction of concomitant gradient artifacts in experimental microtesla MRI

Magnetic resonance imaging (MRI) suffers from artifacts caused by concomitant gradients when the product of the magnetic field gradient and the dimension of the sample becomes comparable to the static magnetic field. To investigate and correct for these artifacts at very low magnetic fields, we have acquired MR images of a 165-mm phantom in a 66-microT field using gradients up to 350 microT/m. We prepolarize the protons in a field of about 100 mT, apply a spin-echo pulse sequence, and detect the precessing spins using a superconducting gradiometer coupled to a superconducting quantum interference device (SQUID). Distortion and blurring are readily apparent at the edges of the images; by comparing the experimental images to computer simulations, we show that concomitant gradients cause these artifacts. We develop a non-perturbative, post-acquisition phase correction algorithm that eliminates the effects of concomitant gradients in both the simulated and the experimental images. This algorithm assumes that the switching time of the phase-encoding gradient is long compared to the spin precession period. In a second technique, we demonstrate that raising the precession field during phase encoding can also eliminate blurring caused by concomitant phase-encoding gradients; this technique enables one to correct concomitant gradient artifacts even when the detector has a restricted bandwidth that sets an upper limit on the precession frequency. In particular, the combination of phase correction and precession field cycling should allow one to add MRI capabilities to existing 300-channel SQUID systems used to detect neuronal currents in the brain because frequency encoding could be performed within the 1-2 kHz bandwidth of the readout system.     

Efficacy of motion artifact reduction in neonatal DW segmented EPI at 3 T using phase correction by numerical optimization and segment data swapping

Segmented echoplanar imaging (EPI) is a potentially valuable acquisition method for neonatal diffusion-weighted imaging (DWI) due to the lower acoustic noise levels as well as reduced blurring and distortion associated with it, as compared with single-shot EPI. Reduced acoustic noise may be important for the safety of neonates. However, little information regarding the efficacy of segmented EPI motion correction schemes is available for the neonatal population. We quantitatively assessed the efficacy of a postprocessing technique for motion artifact reduction involving phase correction by nonlinear optimization, alone and in combination with a novel method of utilizing a second data set (referred to as segment data swapping). These methods were applied to three-directional eight-segment echoplanar DW images obtained from 13 sedated neonates and to nine-directional DW images from 3 unsedated neonates. For comparison, the efficacy of the nonlinear optimization method was also evaluated in four adults. Motion correction efficacy was quantified using the motion artifact-to-signal ratio (ASR). The median, 70th percentile and 90th percentile ASR values obtained from neonatal three-directional DWI using nonlinear optimization alone were 2.8%, 4.6% and 9.6%, respectively. Efficacy improved (P<.005), particularly in dealing with the images most difficult to correct, when the phase correction by numerical optimization was combined with segment data swapping (median ASR=1.9%, 70th percentile ASR=2.7%, 90th percentile ASR=4.3%). Similar results were obtained for nine-directional diffusion tensor imaging. Nonlinear optimization alone applied to adult images showed significantly (P<.001) lower ASR values (median ASR=0.9%, 70th percentile ASR=2.1%, 90th percentile ASR=4.1%), demonstrating the greater challenge in DWI of neonates with segmented EPI. In conclusion, phase correction by nonlinear optimization provides effective motion correction for neonatal DW eight-segment EPI, especially when used in conjunction with segment data swapping.     

Improved reduction of motion artifacts in diffusion imaging using navigator echoes and velocity compensation

Navigator echoes provide a means with which to remove motion artifacts from diffusion-weighted images obtained using any multishot imaging technique. However, residual motion artifact is often present in the corrected images rendering the technique unreliable. It is shown that velocity-compensated diffusion sensitization when used in tandem with a navigator echo further reduces the degree of residual motion artifacts present in the corrected images and improves the reliability and clinical utility of the technique. This is demonstrated by applying a method for quantification of motion artifact to brain images of healthy volunteers scanned using both conventional (Stejskal-Tanner) and velocity-compensated gradient sensitization. Other factors affecting the efficacy of the navigator echo technique, such as brain pulsatile motion, gradient b factor, and navigator echo signal-to-noise ratio, are also discussed.     

A quality control method for detecting and suppressing uncorrected residual motion in fMRI studies

Motion correction is an important step in the functional magnetic resonance imaging (fMRI) analysis pipeline. While many studies simply exclude subjects who are estimated to have moved beyond an arbitrary threshold, there exists no objective method for determining an appropriate threshold. Furthermore, any criterion based only upon motion estimation ignores the potential for proper realignment. The method proposed here uses unsupervised learning (specifically k-means clustering) on features derived from the mean square derivative (MSD) of the signal before and after realignment to identify problem data. These classifications are refined through analysis of correlation between subject activation maps and the mean activation map, as well as the relationship between tasking and motion as measured through regression of the canonical hemodynamic response functions to fit both estimated motion parameters and MSD. The MSD is further used to identify specific scans containing residual motion, data which is suppressed by adding nuisance regressors to the general linear model; this statistical suppression is performed for identified problem subjects, but has potential for use over all subjects. For problem subjects, our results show increased hemodynamic activity more consistent with group results; that is, the addition of nuisance regressors resulted in a doubling of the correlation between the activation map for the problem subjects and the activation map for all subjects. The proposed method should be useful in helping fMRI researchers make more efficient use of their data by reducing the need to exclude entire subjects from studies and thus collect new data to replace excluded subjects.