Image-based EPI ghost correction using an algorithm based on projection onto convex sets (POCS).

This work describes the use of a method, based on the projection onto convex sets (POCS) algorithm, for reduction of the N/2 ghost in echo-planar imaging (EPI). In this method, ghosts outside the parent image are set to zero and a model k-space is obtained from the Fourier transform (FT) of the resulting image. The zeroth- and first-order phase corrections for each line of the original k-space are estimated by comparison with the corresponding line in the model k-space. To overcome problems of phase wrapping, the first-order phase corrections for the lines of the original k-space are estimated by registration with the corresponding lines in the model k-space. It is shown that applying these corrections will result in a reduction of the ghost, and that iterating the process will result in a convergence towards an image in which the ghost is minimized. The method is tested on spin-echo EPI data. The results show that the method is robust and remarkably effective, reducing the N/2 ghost to a level nearly comparable to that achieved with reference scans.     

RINGLET motion correction for 3D MRI acquired with the elliptical centric view order.

A new rigid-body motion correction algorithm is described that is compatible with 3D image sets acquired with the elliptical centric (EC) view order. With this view order, an annular ring of k-space data is acquired in the ky-kz plane during any short time interval. Images for tracking motion can be reconstructed in the yz-plane from any ring of the acquisition data. In these tracking images, a point source (such as an external marker) shows a characteristic bull's-eye pattern that permits motion monitoring and correction. The true position of the point object is located at the center of the bull's-eye pattern. Cross correlation can be performed to automatically track the positions of markers reconstructed from adjacent rings of k-space. To increase the marker signal, the markers are encased in inductively coupled RF coils. Rigid-body motion in the yz-plane is calculated directly with the Euclidean group for rotation and translation, and corrected by rotating and applying phase shifts to any corrupted rings of data. In the current work we present a theoretical analysis of this method, as well as results of volunteer and controlled phantom experiments that demonstrate its initial feasibility. Although the EC view order has mainly been used for MR angiography (MRA), it can also be used for most 3D acquisitions.     

IMPACT: image-based physiological artifacts estimation and correction technique for functional MRI.

Functional MRI (fMRI) signal variation induced by respiratory and cardiac motion affects the activation signal and limits the accuracy of analysis. Current physiological motion correction methods require either synchronization with external monitoring of respiration and heartbeat, specialized pulse sequence design, or k-space data. The IMage-based Physiological Artifacts estimation and Correction Technique (IMPACT), which is free from these constraints, is described. When images are acquired fast enough to sample physiological motion without aliasing, respiratory and cardiac signals can be directly estimated from magnitude images. Physiological artifacts are removed by reordering images according to the estimated respiratory and cardiac phases and then subtracting the Fourier-fitted variation from magnitude images. Compared with the k-space-based method, this method can efficiently and effectively reduce the overall signal fluctuation in the brain and increase the activated area. With this new technique, physiological artifacts can be reduced using traditional fMRI pulse sequences, and existing data can be corrected and reanalyzed without additional experiments.     

k-Space weighted least-squares algorithm for accurate and fast motion extraction from magnetic resonance navigator echoes.

Navigator echoes provide an effective method to monitor physiological motion, allowing motion artifacts to be suppressed by modifying the data acquisition accordingly. The displacement can be measured directly from the navigator phase using an algorithm suggested by Ahn and Cho (IEEE Trans Med Imaging 1987;MI-6:32-36). Although computationally efficient, it is susceptible to noise, particularly contributions from points at the edges of k-space. A k-space weighted least-squares (kLS) algorithm is proposed which fits a straight line to the motion-induced phase shift. The linear fit is weighted strongly to high SNR points near the k-space center and only weakly to the poor SNR points at the edges of k-space. This algorithm was found to be as efficient as Ahn's algorithm but more robust against noise.     

Pulsatile motion effects on 3D magnetic resonance angiography: implications for evaluating carotid artery stenoses.

In-plane carotid artery motion during a 3D MR angiography (MRA) scan can significantly degrade the resulting image resolution. This study characterizes the effect of cardiac pulsatility on 3D contrast-enhanced (CE) MRA with elliptical centric acquisitions using a point-spread function (PSF) analysis. Internal carotid artery (ICA) motion was collected from volunteers and patients using both MR and ultrasound (US) scans. After measuring the carotid artery motion displacement, a simulation was performed which calculated the blurring effects for three different protocols: nongated and two different cardiac gating schemes. The motion sensitivity of each protocol was evaluated for different spatial resolutions. The selection of optimal imaging parameters for a given scan time was investigated. The final results showed that cardiac-gated acquisitions only over a limited region of k-space high spatial frequencies are more time-efficient than cardiac gating for the entire k-space, as it allows for higher resolutions to be achieved and for capturing the arterial phase with low spatial frequencies. Selecting the optimal gating parameters depends directly on the motion characteristics of each individual. Our initial clinical experience is presented, and the need for a real-time tool that characterizes motion behavior for each individual as a prescan protocol is discussed.     

TSE with average-specific phase encoding ordering for motion detection and artifact suppression

PURPOSE: To detect motion-corrupted measurements in multi-average turbo-spin-echo (TSE) acquisitions and reduce motion artifacts in reconstructed images. MATERIALS AND METHODS: An average-specific phase encoding (PE) ordering scheme was developed for multi-average TSE sequences in which each echo train is assigned a unique PE pattern for each pre-averaged image (PAI). A motion detection algorithm is developed based on this new PE ordering to identify which echo trains in which PAIs are motion-corrupted. The detected PE views are discarded and replaced by uncorrupted k-space data of the nearest PAI. Both phantom and human studies were performed to investigate the effectiveness of motion artifact reduction using the proposed method. RESULTS: Motion-corrupted echo trains were successfully detected in all phantom and human experiments. Significant motion artifact suppression has been achieved for most studies. The residual artifacts in the reconstructed images are mainly caused by residual inconsistencies that remain after the corrupted k-space data is corrected. CONCLUSION: The proposed method combines a novel data acquisition scheme, a robust motion detection algorithm, and a simple motion correction algorithm. It is effective in reducing motion artifacts for images corrupted by either bulk motion or local motion that occasionally happens during data acquisition.     

Direct monitoring of coronary artery motion with cardiac fat navigator echoes.

Navigator echoes (NAVs) provide an effective means of monitoring physiological motion in magnetic resonance imaging (MRI). Motion artifacts can be suppressed by adjusting the data acquisition accordingly. The standard pencil-beam NAV has been used to detect diaphragm motion; however, it does not monitor cardiac motion effectively. Here we report a navigator approach that directly measures coronary artery motion by exciting the surrounding epicardial fat and sampling the signal with a k-space trajectory sensitized to various motion parameters. The present preliminary human study demonstrates that superior-inferior (SI) respiratory motion of the coronary arteries detected by the cardiac fat NAV highly correlates with SI diaphragmatic motion detected by the pencil-beam NAV. In addition, the cardiac fat navigator gating is slightly more effective than the diaphragmatic navigator gating in suppressing motion artifacts in free-breathing 3D coronary MR angiography (MRA).     

Motion artifact reduction and vessel enhancement for free-breathing navigator-gated coronary MRA using 3D k-space reordering.

Breathing-induced bulk motion of the myocardium during data acquisition may cause severe image artifacts in coronary magnetic resonance angiography (MRA). Current motion compensation strategies include breath-holding or free-breathing MR navigator gating and tracking techniques. Navigator-based techniques have been further refined by the applications of sophisticated 2D k-space reordering techniques. A further improvement in image quality and a reduction of relative scanning duration may be expected from a 3D k-space reordering scheme. Therefore, a 3D k-space reordered acquisition scheme including a 3D navigator gated and corrected segmented k-space gradient echo imaging sequence for coronary MRA was implemented. This new zonal motion-adapted acquisition and reordering technique (ZMART) was developed on the basis of a numerical simulation of the Bloch equations. The technique was implemented on a commercial 1.5T MR system, and first phantom and in vivo experiments were performed. Consistent with the results of the theoretical findings, the results obtained in the phantom studies demonstrate a significant reduction of motion artifacts when compared to conventional (non-k-space reordered) gating techniques. Preliminary in vivo findings also compare favorably with the phantom experiments and theoretical considerations. Magn Reson Med 45:645-652, 2001.     

Spherical navigator echoes for full 3D rigid body motion measurement in MRI.

We developed a 3D spherical navigator (SNAV) echo technique that can measure rigid body motion in all six degrees of freedom simultaneously by sampling a spherical shell in k-space. 3D rotations of an imaged object simply rotate the data on this shell and can be detected by registration of k-space magnitude values. 3D translations add phase shifts to the data on the shell and can be detected with a weighted least-squares fit to the phase differences at corresponding points. MRI pulse sequences were developed to study k-space sampling strategies on such a shell. Data collected with a computer-controlled motion phantom with known rotational and translational motions were used to evaluate the technique. The accuracy and precision of the technique depend on the sampling density. Roughly 2000 sample points were necessary for accurate detection to within the error limits of the motion phantom when using a prototype time-intensive sampling method. This number of samples can be captured in an approximately 27-ms double excitation SNAV pulse sequence with a 3D helical spiral trajectory. Preliminary results with the helical SNAV are encouraging and indicate that accurate motion measurement suitable for retrospective or prospective correction should be feasible with SNAV echoes.     

Partial k-space reconstruction in single-shot diffusion-weighted echo-planar imaging.

Partial k-space sampling is frequently used in single-shot diffusion-weighted echo-planar imaging (DW-EPI) to reduce the TE and thereby improve the SNR. However, it increases the sensitivity of the technique to bulk rotational motion, which introduces a phase gradient across the tissue that shifts the echo in k-space. If the echo is displaced into the high spatial frequencies, conventional homodyne reconstruction fails, causing intensity oscillations across the image. Zero-padding, on the other hand, compromises the image resolution and may cause truncation artifacts. We present an adaptive version of the homodyne algorithm that detects the location of the echo in k-space and adjusts the center and width of the homodyne filters accordingly. The adaptive algorithm produces artifact-free images when the echo is shifted into the high positive k-space range, and reduces to the standard homodyne algorithm in the absence of bulk motion.     

Coronary artery imaging using three-dimensional breath-hold steady-state free precession with two-dimensional iterative partial fourier reconstruction

PURPOSE: To assess the feasibility of using a two-dimensional partial Fourier (PF) reconstruction scheme to reduce the acquisition time of magnetic resonance imaging (MRI) of coronary arteries. MATERIALS AND METHODS: Symmetric k-space data sets of coronary arteries were collected in seven volunteers using a three-dimensional breath-hold steady-state free precession (SSFP) sequence. Partial, asymmetric k-space data sets were generated by removing 25% of the data in the readout direction and 25% of the data in the phase encoding direction. The missing data were then estimated using a two-dimensional projection-onto-convex-sets (POCS) algorithm or filled with zeroes. Images were reconstructed from the full data set, the PF data set, and the zero-filled (ZF) data set, respectively. Coronary artery sharpness was evaluated quantitatively and qualitatively. RESULTS: Coronary artery sharpness in PF images was comparable to that in full k-space images and significantly better than that in ZF images. CONCLUSION: Two-dimensional POCS PF reconstruction is a potentially useful technique for reducing acquisition time or improving spatial resolution for breath-hold coronary MR angiography.     

Motion correction using the k-space phase difference of orthogonal acquisitions.

Rigid body translations of an object in MRI create image artifacts along the phase-encode (PE) direction in standard 2DFT imaging. If two images are acquired with swapped PE direction, it is possible to determine and correct for arbitrary in-plane translational interview motions in both images directly from phase differences in the k-space acquisitions by solving a large system of linear equations. For example, if one assumes two N x N 2D acquisitions with in-plane translational interview motion, 4N unknown motions may corrupt the two images, but the phase difference at each point in k-space yields a system of N(2) equations in these 4N unknowns. If the acquisitions have orthogonal PE directions, this highly overdetermined system of equations can be solved to provide the motion records, which in turn can be used to correct the motion artifacts in each image. The theory of this orthogonal k-space phase difference (ORKPHAD) technique is described, and results are presented for synthetic and in vivo motion-corrupted data sets. In all cases, the data showed clear improvement of translation-induced artifacts. These methods do not require special pulse sequences and are theoretically generalizable to partial Fourier imaging and 3D acquisitions.     

Floating navigator echo (FNAV) for in-plane 2D translational motion estimation.

A modification of the classical navigator echo (NAV) technique is presented whereby both 2D translational motion components are computed from a single navigator line. Instead of acquiring the NAV at the center of the k-space, a kx line is acquired off-center in the phase-encoding (ky) direction as a floating NAV (FNAV). It is shown that the translational motion in both the readout and phase-encoding directions can be computed from this line. The algorithm used is described in detail and verified experimentally. The new technique can be readily implemented to replace classic NAV in MRI sequences, with little to no additional cost or complexity. The new method can help suppress 2D translational motion and provide more accurate motion estimates for other motion-suppression techniques, such as the diminishing variance algorithm.     

Impact of bulk cardiac motion on right coronary MR angiography and vessel wall imaging

The purpose of this study was to investigate the impact of in-plane coronary artery motion on coronary magnetic resonance angiography (MRA) and coronary MR vessel wall imaging. Free-breathing, navigator-gated, 3D-segmented k-space turbo field echo ((TFE)/echo-planar imaging (EPI)) coronary MRA and 2D fast spin-echo coronary vessel wall imaging of the right coronary artery (RCA) were performed in 15 healthy adult subjects. Images were acquired at two different diastolic time periods in each subject: 1) during a subject-specific diastasis period (in-plane velocity <4 cm/second) identified from analysis of in-plane coronary artery motion, and 2) using a diastolic trigger delay based on a previously implemented heart-rate-dependent empirical formula. RCA vessel wall imaging was only feasible with subject-specific middiastolic acquisition, while the coronary wall could not be identified with the heart-rate-dependent formula. For coronary MRA, RCA border definition was improved by 13% (P < 0.001) with the use of subject-specific trigger delay (vs. heart-rate-dependent delay). Subject-specific middiastolic image acquisition improves 3D TFE/EPI coronary MRA, and is critical for RCA vessel wall imaging.     

Application and assessment of a robust elastic motion correction algorithm to dynamic MRI

The purpose of this study was to assess the performance of a new motion correction algorithm. Twenty-five dynamic MR mammography (MRM) data sets and 25 contrast-enhanced three-dimensional peripheral MR angiographic (MRA) data sets which were affected by patient motion of varying severeness were selected retrospectively from routine examinations. Anonymized data were registered by a new experimental elastic motion correction algorithm. The algorithm works by computing a similarity measure for the two volumes that takes into account expected signal changes due to the presence of a contrast agent while penalizing other signal changes caused by patient motion. A conjugate gradient method is used to find the best possible set of motion parameters that maximizes the similarity measures across the entire volume. Images before and after correction were visually evaluated and scored by experienced radiologists with respect to reduction of motion, improvement of image quality, disappearance of existing lesions or creation of artifactual lesions. It was found that the correction improves image quality (76% for MRM and 96% for MRA) and diagnosability (60% for MRM and 96% for MRA).     

POCS‐enhanced correction of motion artifacts in parallel MRI

A new method for correction of MRI motion artifacts induced by corrupted k‐space data, acquired by multiple receiver coils such as phased arrays, is presented. In our approach, a projections onto convex sets (POCS)‐based method for reconstruction of sensitivity encoded MRI data (POCSENSE) is employed to identify corrupted k‐space samples. After the erroneous data are discarded from the dataset, the artifact‐free images are restored from the remaining data using coil sensitivity profiles. The error detection and data restoration are based on informational redundancy of phased‐array data and may be applied to full and reduced datasets. An important advantage of the new POCS‐based method is that, in addition to multicoil data redundancy, it can use a priori known properties about the imaged object for improved MR image artifact correction. The use of such information was shown to improve significantly k‐space error detection and image artifact correction. The method was validated on data corrupted by simulated and real motion such as head motion and pulsatile flow. Magn Reson Med 63:1104–1110, 2010. © 2010 Wiley‐Liss, Inc.     

Respiratory self-gated four-dimensional coronary MR angiography: a feasibility study.

The four-dimensional (4D) coronary MR angiography (MRA) approach has been developed to eliminate the need for accurate determination of the acquisition window and trigger delay time. Diaphragm navigator (NAV) has been the conventional respiratory gating method for free-breathing coronary MRA. However, NAV echo acquisition interrupts the continuous radiofrequency pulse application required for 4D steady-state free precession coronary MRA. The objective of this work was to investigate the feasibility of a respiratory self-gating (RSG) technique for 4D coronary MRA and its effectiveness by comparing with retrospective NAV gating. Data were acquired continuously throughout the cardiac cycle and retrospectively remapped to cardiac phases based on the electrocardiogram signal simultaneously recorded. An RSG signal extracted from a direct measurement of the heart position was used for retrospective respiratory gating and motion correction. In seven healthy volunteers, 4D MRA images were reconstructed, allowing retrospective assessment of the cardiac motion of the coronary artery and selection of the images with the best vessel delineation. Statistical analysis shows that 4D RSG provides coronary artery delineation comparable to mid-diastole images acquired using NAV. Respiratory self-gating is an effective method for eliminating respiratory motion artifacts and allows 4D coronary MRA during free breathing.     

Free-breathing renal magnetic resonance angiography with steady-state free-precession and slab-selective spin inversion combined with radial k-space sampling and water-selective excitation.

The impact of radial k-space sampling and water-selective excitation on a novel navigator-gated cardiac-triggered slab-selective inversion prepared 3D steady-state free-precession (SSFP) renal MR angiography (MRA) sequence was investigated. Renal MRA was performed on a 1.5-T MR system using three inversion prepared SSFP approaches: Cartesian (TR/TE: 5.7/2.8 ms, FA: 85 degrees), radial (TR/TE: 5.5/2.7 ms, FA: 85 degrees) SSFP, and radial SSFP combined with water-selective excitation (TR/TE: 9.9/4.9 ms, FA: 85 degrees). Radial data acquisition lead to significantly reduced motion artifacts (P < 0.05). SNR and CNR were best using Cartesian SSFP (P < 0.05). Vessel sharpness and vessel length were comparable in all sequences. The addition of a water-selective excitation could not improve image quality. In conclusion, radial k-space sampling reduces motion artifacts significantly in slab-selective inversion prepared renal MRA, while SNR and CNR are decreased. The addition of water-selective excitation could not improve the lower CNR in radial scanning.     

High-resolution diffusion-weighted imaging with interleaved variable-density spiral acquisitions

PURPOSE: To develop a multishot magnetic resonance imaging (MRI) pulse sequence and reconstruction algorithm for diffusion-weighted imaging (DWI) in the brain with submillimeter in-plane resolution. MATERIALS AND METHODS: A self-navigated multishot acquisition technique based on variable-density spiral k-space trajectory design was implemented on clinical MRI scanners. The image reconstruction algorithm takes advantage of the oversampling of the center k-space and uses the densely sampled central portion of the k-space data for both imaging reconstruction and motion correction. The developed DWI technique was tested in an agar gel phantom and three healthy volunteers. RESULTS: Motions result in phase and k-space shifts in the DWI data acquired using multishot spiral acquisitions. With the two-dimensional self-navigator correction, diffusion-weighted images with a resolution of 0.9 x 0.9 x 3 mm3 were successfully obtained using different interleaves ranging from 8-32. The measured apparent diffusion coefficient (ADC) in the homogenous gel phantom was (1.66 +/- 0.09) x 10(-3) mm2/second, which was the same as measured with single-shot methods. The intersubject average ADC from the brain parenchyma of normal adults was (0.91 +/- 0.01) x 10(-3) mm2/second, which was in a good agreement with the reported literature values. CONCLUSION: The self-navigated multishot variable-density spiral acquisition provides a time-efficient approach to acquire high-resolution diffusion-weighted images on a clinical scanner. The reconstruction algorithm based on motion correction in the k-space data is robust, and measured ADC values are accurate and reproducible.