True myocardial motion tracking

Myocardial tagging is a powerful tool for the assessment of in-plane cardiac motion. However, for previous myocardial tagging techniques, the imaged slice is fixed with respect to the magnet coordinate system. Thus, images acquired at different heart phases do not always represent the same slice of the myocardium. A new myocardial tagging technique is presented, which takes the through-plane motion into consideration. It involves tagging of the desired myocardial slice and applying a subtraction imaging technique to image just that part of the myocardium. The examination time can be reduced considerably by the acquisition of two one-dimensionally tagged images. To increase the signal-to-noise ratio especially at later heart phases, variable imaging RF excitation flip angles are applied. To reduce motion artifacts a repetitive breathhold scheme was applied. in vivo results demonstrate that the tags can be accurately tracked within the entire heart period with a temporal resolution of 35 ms, even at a top basal level of the heart and right ventricle.     

Free-breathing radial acquisitions of the heart.

There is considerable interest in performing free-breathing acquisitions of the heart in order to obtain high-quality images without the need for multiple, long breathholds. In this article a 3D motion-correction method is described that is based on image registration of in-plane data and through-plane slice tracking. A number of fast radial undersampled images are acquired, each of which is free of motion artifacts. Initially, in-plane translational and rotational motion between each image was corrected before combining the data to give a fully sampled image. At the next stage, correction of in-plane deformation, in addition to translations and rotations, was performed in the image domain. Through-plane translational motion was compensated using a navigator echo to move the acquisition plane. Using this method, information on the motion of the heart was captured at the same time as acquiring the image data. No motion model, assumptions about the motion, or training data are required. The method is demonstrated on phantom data and cardiac images acquired on free-breathing volunteers.     

Simultaneous MRI tagging and through-plane velocity quantification: a three-dimensional myocardial motion tracking algorithm

A tracking algorithm was developed for calculation of three-dimensional point-specific myocardial motion. The algorithm was designed for images acquired with simultaneous magnetic resonance imaging (MRI) grid tagging and through-plane velocity quantification. The tagging grid provided the in-plane motion while the velocity quantification measured the through-plane motion. In four healthy volunteers, the in vivo performance was evaluated by comparing the systolic through-plane displacement with the displacement of tagging-grid intersections in long-axis images. The correlation coefficient was 0.93 (P < 0.001, N = 183). A t-test for paired samples revealed a small underestimation of the through-plane displacement by 0.04 +/- 0.09 cm (mean +/- SD, P < 0.001) on an average displacement of 0.77 +/- 0.23 cm toward the apex. The authors conclude that three-dimensional point-specific motion tracking based on simultaneous tagging and velocity quantification is competitive with other methods such as tagging in mutually orthogonal image planes or quantification of three orthogonal velocity components.     

Heart motion adapted cine phase-contrast flow measurements through the aortic valve.

A method for magnetic resonance cine velocity mapping through heart valves with adaptation of both slice offset and angulation according to the motion of the valvular plane of the heart is presented. By means of a subtractive labeling technique, basal myocardial markers are obtained and automatically extracted for quantification of heart motion at the valvular level. The captured excursion of the basal plane is used to calculate the slice offset and angulation of each required time frame for cine velocity mapping. Through-plane velocity offsets are corrected by subtracting velocities introduced by basal plane motion from the measured velocities. For evaluation of the method, flow measurements downstream from the aortic valve were performed both with and without slice adaptation in 11 healthy volunteers and in four patients with aortic regurgitation. Maximum through-plane motion at the aortic root level as calculated from the labeled markers averaged 8.9 mm in the volunteers and 6.5 mm in the patients. The left coronary root was visible in 2-4 (mean: 2.2) time frames during early diastole when imaging with a spatially fixed slice. Time frames obtained with slice adaptation did not contain the coronary roots. Motion correction increased the apparent regurgitant volume by 5.7 +/- 0.4 ml for patients with clinical aortic regurgitation, for an increase of approximately 50%. The proposed method provides flow measurements with correction for through-plane motion perpendicular to the aortic root between the valvular annulus and the coronary ostia throughout the cardiac cycle. Magn Reson Med 42:970-978, 1999.     

3D myocardial tissue tracking with slice followed cine DENSE MRI

PURPOSE: To track three-dimensional (3D) myocardial tissue motion using slice followed cine displacement encoded imaging with stimulated echoes (DENSE). MATERIALS AND METHODS: Slice following (SF) has previously been developed for 2D myocardial tagging to compensate for the effect of through-plane motion on 2D tissue tracking. By incorporating SF into a cine DENSE sequence, and applying displacement encoding in three orthogonal directions, we demonstrate the ability to track discrete elements of a slice of myocardium in 3D as the heart moves through the cardiac cycle. The SF cine DENSE tracking algorithm was validated on a moving phantom, and the effects of through-plane motion on 2D cardiac strain were investigated in six healthy subjects. RESULTS: A through-plane tracking accuracy of 0.46 +/- 0.32 mm was measured for a typical range of myocardial motion using a rotating phantom. In vivo 3D measurements of cardiac motion were consistent with prior myocardial tagging results. Through-plane rotation in a mid-ventricularshort-axis view was shown to decrease the magnitude of the 2D end-systolic circumferential strain by 3.91 +/- 0.43% and increase the corresponding radial strain by 6.01 +/- 1.07%. CONCLUSION: Slice followed cine DENSE provides an accurate method for 3D tissue tracking.     

Time-resolved undersampled projection reconstruction magnetic resonance imaging of the peripheral vessels using multi-echo acquisition.

The hybrid projection reconstruction (PR) imaging provides high temporal resolution through an undersampled PR acquisition for the in-plane dimensions and Cartesian slice encoding for the through-plane dimension. The undersampling of projection data introduces streak artifact, which may severely compromise image quality. This study reports on a combination of multi-echo acquisition with time-resolved undersampled PR imaging and its application to peripheral magnetic resonance angiography. Multi-echo acquisition improved imaging speed effectively, thereby reducing the undersampling streak artifact and improving the temporal resolution. The gradient distortion was reduced through gradient calibration and accurate k-space trajectory measurement.     

Artifact reduction in undersampled projection reconstruction MRI of the peripheral vessels using selective excitation.

The projection reconstruction (PR)-HyperTRICKS (time resolved imaging of contrast kinetics) acquisition integrates the benefits of through-plane Cartesian slice encoding and in-plane undersampled PR. It provides high spatial resolution both in-plane (about 1 mm(2)) and through-plane (1-2 mm), as well as relatively high temporal resolution (about 0.25 frames per second). However, undersampling artifacts that originate from anatomy superior or inferior to a coronal imaging FOV may severely degrade the image quality. In coronal MRA acquisitions, the slice coverage is limited in order to achieve high temporal resolution. In this report we describe an artifact reduction method that uses selective excitation in PR-HyperTRICKS. This technique significantly reduces undersampling streak artifacts while it increases the slice coverage.     

Improved turbo spin-echo imaging of the heart with motion-tracking

PURPOSE: To improve dark-blood and short tau inversion recovery (STIR) prepared turbo spin-echo (TSE) imaging of the heart, particularly in the basal short-axis plane where cardiac misregistration between the preparation and imaging phases is high. MATERIALS AND METHODS: In the first approach (tracked), the basal short-axis plane was labeled and tracked over the cardiac cycle. The slice-selective 180 degrees dark-blood and STIR preparation pulses were then independently positioned on the appropriately timed labeled images. In the second approach (offset), the preparation pulses were output in the same orientation as the imaging plane, but with a user-defined slice offset that was derived from the labeled data. Both approaches were compared with the standard untracked dark-blood STIR TSE sequence (7-mm slice thickness) in 10 healthy volunteers. RESULTS: For typical preparation slice thicknesses, tracked and offset TSE images were superior to the untracked images (both P < 0.01). For the more mobile right ventricle (RV), the image quality of the tracked images was superior to that of the offset images (P < 0.05). CONCLUSION: Tracking the through-plane motion of the heart between preparation and imaging phases improves the quality of thin-slice basal short-axis TSE images, particularly for the more mobile RV.     

Effect of through‐plane motion on left ventricular rotation: A study using slice‐following harmonic phase imaging

Noninvasive quantification of regional left ventricular rotation may improve understanding of cardiac function. Current methods used to quantify rotation typically acquire data on a set of prescribed short‐axis slices, neglecting effects due to through‐plane myocardial motion. We combine principles of slice‐following tagged imaging with harmonic phase analysis methods to account for through‐plane motion in regional rotation measurements. We compare rotation and torsion measurements obtained using our method to those obtained from imaging datasets acquired without slice‐following. Our results in normal volunteers demonstrate differences in the general trends of average and regional rotation‐time plots in midbasal slices and the rotation versus circumferential strain loops. We observe substantial errors in measured peak average rotation of the order of 58% for basal slices (due to change in the pattern of the curve), ?6.6% for midventricular slices, and ?8.5% for apical slices; and an average error in base‐to‐apex torsion of 19% when through‐plane motion is not considered. This study concludes that due to an inherent base‐to‐apex gradient in rotation that exists in the left ventricular, accounting for through‐plane motion is critical to the accuracy of left ventricular rotation quantification. Magn Reson Med, 2013. © 2011 Wiley Periodicals, Inc.     

Real-time adaptive motion correction in functional MRI

Functional magnetic resonance imaging (fMRI) of the brain is often degraded by bulk head motion. Algorithms that address this by retrospective re-registration of images in an fMRI time series are all fundamentally limited by any motion that occurs through-plane. Here, a technique is described that can account for such motion by prospective. correction in real time. A navigator echo is used before every image acquisition to detect superior/inferior displacements of the head. The displacement information is then used to adjust the plane of excitation of the ensuing single-shot echo-planar fMRI axial image. These correction updates can be completed in 100 ms with motion sensitivity at least as small as 0.5 mm. The efficacy of this method is documented in phantom and human studies.     

Quantification of in-plane motion of the coronary arteries during the cardiac cycle: Implications for acquisition window duration for MR flow quantification

Motion of the coronary arteries during the heart cycle can result in image blurring and inaccurate flow quantification by MR. This condition applies particularly for longer acquisition windows that are typical of breath-hold coronary flow measurements. To determine the sensitivity of the technique to in-plane motion of different coronary arteries, the temporal variation in coronary position was measured in a plane perpendicular to the proximal portion of the vessel. The results indicated the presence of substantial displacement of the coronary arteries within the cardiac cycle, with a magnitude of motion approximately twice as large for the right as for the left coronary arteries. An estimation of the resulting vessel blurring was calculated, showing that the duration of the acquisition window for high spatial resolution coronary flow acquisitions should be less than 25 to 120 msec, depending on the specific coronary artery studied. In addition, these data specify optimal acquisition window placement for high resolution coronary angiography.     

Aortic and mitral regurgitation: quantification using moving slice velocity mapping

Comprehensive assessment of the severity of valvular insufficiency includes quantification of regurgitant volumes. Previous methods lack reliable slice positioning with respect to the valve and are prone to velocity offsets due to through-plane motion of the valvular plane of the heart. Recently, the moving slice velocity mapping technique was proposed. In this study, the technique was applied for quantification of mitral and aortic regurgitation. Time-efficient navigator-based respiratory artifact suppression was achieved by implementing a prospective k-space reordering scheme in conjunction with slice position correction. Twelve patients with aortic insufficiency and three patients with mitral insufficiency were studied. Aortic regurgitant volumes were calculated from diastolic velocities mapped with a moving slice 5 mm distal to the aortic valve annulus. Mitral regurgitant flow was indirectly assessed by measuring mitral inflow at the level of the mitral annulus and net aortic outflow. Regurgitant fractions, derived from velocity data corrected for through-plane motion, were compared to data without correction for through-plane motion. In patients with mild and moderate aortic regurgitation, regurgitant fractions differed by 60% and 15%, on average, when comparing corrected and uncorrected data, respectively. Differences in severe aortic regurgitation were less (7%). Due to the large orifice area of the mitral valve, differences were still substantial in moderate-to-severe mitral regurgitation (19%). The moving slice velocity mapping technique was successfully applied in patients with aortic and mitral regurgitation. The importance of correction for valvular through-plane motion is demonstrated.     

Prospective motion correction using tracking coils

Intracavity imaging coils provide higher signal‐to‐noise than surface coils and have the potential to provide higher spatial resolution in shorter acquisition times. However, images from these coils suffer from physiologically induced motion artifacts, as both the anatomy and the coils move during image acquisition. We developed prospective motion‐correction techniques for intracavity imaging using an array of tracking coils. The system had <50 ms latency between tracking and imaging, so that the images from the intracavity coil were acquired in a frame of reference defined by the tracking array rather than by the system's gradient coils. Two‐dimensional gradient‐recalled and three‐dimensional electrocardiogram‐gated inversion‐recovery‐fast‐gradient‐echo sequences were tested with prospective motion correction using ex vivo hearts placed on a moving platform simulating both respiratory and cardiac motion. Human abdominal tests were subsequently conducted. The tracking array provided a positional accuracy of 0.7 ± 0.5 mm, 0.6 ± 0.4 mm, and 0.1 ± 0.1 mm along the X, Y, and Z directions at a rate of 20 frames‐per‐second. The ex vivo and human experiments showed significant image quality improvements for both in‐plane and through‐plane motion correction, which although not performed in intracavity imaging, demonstrates the feasibility of implementing such a motion‐correction system in a future design of combined tracking and intracavity coil. Magn Reson Med, 2013. © 2012 Wiley Periodicals, Inc.     

Fast method for correcting image misregistration due to organ motion in time-series MRI data.

Time-series MRI data often suffers from image misalignment due to patient movement and respiratory and other physiologic motion during the acquisition process. It is necessary that this misalignment be corrected prior to any automated quantitative analysis. In this article a fast and automated technique for removing in-plane misalignment from time-series MRI data is presented. The method is computationally efficient, robust, and fine-tuned for the clinical setting. The method was implemented and tested on data from 21 human subjects, including myocardial perfusion imaging, renal perfusion imaging, and blood-oxygen level-dependent cardiac T(2*) imaging. In these applications 10-fold or better reduction in image misalignment is reported. The improvement after registration on representative time-intensity curves is shown. Although the method currently corrects translation motion using image center of mass, the mathematical framework of our approach may be extended to correct rotation and other higher-order displacements.     

MR fluoroscopy using projection reconstruction multi-gradient-echo (prMGE) MRI.

A projection reconstruction multi-gradient-echo (prMGE) technique is presented. The introduced technique is an extension of a standard projection reconstruction steady-state gradient-echo technique allowing for the acquisition of several gradient echoes after each excitation of the spin system. Each echo train is used for acquiring data of a certain angular segment of k-space. By use of echo trains consisting of up to four echoes, the overall acquisition time for a 128(2) image can be reduced to 150 ms without sacrificing image quality. Results are presented for cardiac fluoroscopy, for the visualization of swallowing, and for the visualization of joint motion. For all investigated applications promising results have been obtained. Especially in parts of the body where motion on an even shorter time scale than the acquisition process or significant in-plane or through-plane flow are within the field of view, the introduced technique appears to be a promising technique for MR fluoroscopy. Magn Reson Med 42:324-334, 1999.     

Hybrid dante and phase-contrast imaging technique for measurement of three- dimensional myocardial wall motion

Characterization of myocardial stress and strain is necessary for a complete understanding of myocardial function. The precise quantification of regional myocardial strain is complicated by its time-varying pattern and regional variation resulting from the anisotropy of the myocardium and by complex torsional and shortening motions of the heart during the cardiac cycle. The authors have developed a technique for point-specific tracking of myocardial motion along all three axes in a constant selected section of myocardium by combining prospective section selection with in-plane DANTE (delays alternating with nutations for tailored excitation) tissue tagging and phase-contrast detection of motion perpendicular to the image plane. With this technique, it is possible to determine point-specific myocardial strain values in vivo.     

Accuracy of segmented MR velocity mapping to measure small vessel pulsatile flow in a phantom simulating cardiac motion

The purpose of this study was to investigate the accuracy of conventional, segmented, and echo-shared MR velocity mapping sequences to measure pulsatile flow in small moving vessels using a phantom with simulated cardiac motion. The phantom moved either cyclically in-plane, through-plane, in- and through-plane, or was stationary. The mean error in average flow was -2% +/- 3% (mean +/- SD) for all sequences under all conditions, with or without background correction, as long as the region of interest (ROI) size was equal to the vessel cross-sectional size. Overestimation of flow as a result of an oversized ROI was less than 20%, and independent of field of view (FOV) and matrix, as long as the offset in angle between the imaging plane and flow direction was less than 10 degrees. Segmented velocity mapping sequences are surprisingly accurate in measuring average flow and render flow profiles in small moving vessels despite the blurring in the images due to vessel motion. J. Magn. Reson. Imaging 2001;13:722-728.     

Assessment of flow in the right human coronary artery by magnetic resonance phase contrast velocity measurement: Effects of cardiac and respiratory motion

Flow in the human right coronary artery was determined using magnetic resonance phase contrast velocity quantification. Two methods were applied to reduce respiratory motion: imaging during breath holding, which is fast, and retrospective respiratory gating, which has a high temporal resolution (32 ms) in the cardiac cycle. Vessel cross-sectional area, through-plane velocity, and volume flow were determined in six healthy subjects. In-plane vessel displacement during the cardiac cycle, caused by cardiac contraction, was about 2–4 mm within a time frame of 32 ms in systole and early diastole. The motion resulted in blurring of images obtained during breath holding caused by the large acquisition time window (126 ms) within the cardiac cycle. Therefore, only with a high temporal resolution correct velocity images over the entire cardiac cycle could be obtained. The time- and cross-sectionally averaged velocity was 7 ± 2 cm/s, and the volume flow was 30 ± 10 ml/min.     

Correction of motion artifacts from cardiac cine magnetic resonance images

RATIONALE AND OBJECTIVES: An image registration method was developed to automatically correct motion artifacts, mostly from breathing, from cardiac cine magnetic resonance (MR) images. MATERIALS AND METHODS: The location of each slice in an image stack was optimized by maximizing a similarity measure of the slice with another image slice stack. The optimization was performed iteratively and both image stacks were corrected simultaneously. Two procedures to optimize the similarity were tested: standard gradient optimization and stochastic optimization in which one slice is chosen randomly from the image stacks and its location is optimized. In this work, cine short- and long-axis images were used. In addition to visual inspection results from real data, the performance of the algorithm was evaluated quantitatively by simulating the movements in four real MR data sets. The mean error and standard deviation were defined for 50 simulated movements as each slice was randomly displaced. The error rate, defined as the percentage of non-satisfactory registration results, was evaluated. The paired t-test was used to evaluate the statistical difference between the tested optimization methods. RESULTS: The algorithm developed was successfully applied to correct motion artifacts from real and simulated data. The results, where typical motion artifacts were simulated, indicated an error rate of about 3%. Subvoxel registration accuracy was also achieved. When different optimization methods were compared, the registration accuracy of the stochastic approach proved to be superior to the standard gradient technique (P < 10(-9)). CONCLUSIONS: The novel method was capable of robustly and accurately correcting motion artifacts from cardiac cine MR images.     

Model‐based reconstruction for cardiac cine MRI without ECG or breath holding

This paper describes an acquisition and reconstruction strategy for cardiac cine MRI that does not require the use of electrocardiogram or breath holding. The method has similarities with self‐gated techniques as information about cardiac and respiratory motion is derived from the imaging sequence itself; here, by acquiring the center k‐space line at the beginning of each segment of a balanced steady‐state free precession sequence. However, the reconstruction step is fundamentally different: a generalized reconstruction by inversion of coupled systems is used instead of conventional gating. By correcting for nonrigid cardiac and respiratory motion, generalized reconstruction by inversion of coupled systems (GRICS) uses all acquired data, whereas gating rejects data acquired in certain motion states. The method relies on the processing and analysis of the k‐space central line data: local information from a 32‐channel cardiac coil is used in order to automatically extract eigenmodes of both cardiac and respiratory motion. In the GRICS framework, these eigenmodes are used as driving signals of a motion model. The motion model is defined piecewise, so that each cardiac phase is reconstructed independently. Results from six healthy volunteers, with various slice orientations, show improved image quality compared to combined respiratory and cardiac gating. Magn Reson Med 63:1247–1257, 2010. © 2010 Wiley‐Liss, Inc.