High-precision MR velocity mapping by 3D-fourier phase encoding with a small number of encoding steps

The final result of Fourier velocity mapping is a set of images, each representing the spatial distribution of spins at a given velocity. To acquire data in a short time, the number of encoding gradient steps must be as small as possible, but this can mean sacrificing velocity resolution. We used interpolation methods to obtain high velocity resolution with a small number of encoding steps involving linear interpolation from 16 encoding steps or more and zero-filling interpolation from two to eight encoding steps. Velocity measured by interpolated Fourier-flow encoding agreed well with values obtained using a calibrated phantom. A simulation of noise on the images of the phantom showed that, for a given acquisition time, increasing number of encoding steps in the Fourier flow encoding gave better precision for velocity measurement than did averaging identical signals in phase-mapping methods.     

Acceleration mapping by fourier acceleration-encoding: in vitro study and initial results in the great thoracic vessels

Acceleration mapping can be conducted by replacing the bipolar gradient pulse of a velocity mapping sequence by a tripolar pulse. However, since the acceleration encoding pulse is longer, the image quality is altered by the requirement of a long echo time. Since Fourier encoding velocity imaging has been shown to be robust, this velocity mapping method was transformed into an acceleration mapping method. Four steps of the tripolar acceleration encoding gradient pulse were applied successively; acceleration was then obtained by Fourier transform after zero-filling. The accuracy of the method was assessed with a phantom giving a pulsatile flow. Acceleration maps of the ascending aorta and pulmonary artery were obtained in 10 healthy volunteers. The acceleration values measured were in the range of known physiologic values. The feasibility of Fourier encoding acceleration imaging was also demonstrated in four patients.     

Exercise-related changes in aortic flow measured with spiral echo-planar MR velocity mapping

Spiral echo-planar magnetic resonance (MR) velocity mapping was used to measure exercise-related changes in flow in the descending thoracic aorta in 10 healthy volunteers. Flow was measured at rest and Immediately after dynamic exercise, with a 0.5-T imager with a surface receiving coil and electrocardiographic triggering. Supine exercise was performed with a home-built pedaling apparatus. Spiral velocity mapping was performed in a transverse plane through the descending thoracic aorta with the subject at rest. The subject was then asked to perform maximum exercise, stop, and hold his breath during a four-heartbeat acquisition time. Eight cine frames with a temporal resolution of 50 msec were acquired through systole. Each image was acquired in 40 msec during spiral acquisition of k-space data, starting at the center, 6 msec after the excitation pulse. Reproduclbility of the technique was established by repeating the flow measurement in four consecutive heartbeats. At rest, the heart rate (mean ± standard deviation), mean aortic flow, peak aortic flow, and time to peak flow were 68 beats per minute ± 6, 41 mllliliters per beat ±8, 107 mL/sec ± 20, and 175 msec ± 25, respectively. After exercise, the heart rate and mean and peak aortic flow were significantly increased (P < 0.001), measuring 101 beats per minute ±12, 57 milliliters per beat ± 11, and 158 mL/sec ±29, respectively, while the time to peak flow (115 msec ±32) was significantly reduced (P < 0.001). The four sets of values obtained for the first four consecutive heartbeats measured at rest were similar, as were those obtained for the first four heartbeats after exercise.     

Multidimensional MR mapping of multiple components of velocity and acceleration by fourier phase encoding with a small number of encoding steps.

Previous studies have shown that the multi-step approach of velocity or acceleration encoding is highly efficient in terms of the signal-to-noise ratio per unit time. This work describes a multidimensional extension of this method for simultaneously measuring multiple components of velocity and acceleration with a few encoding steps. N flow dimensions were encoded with an ND-matrix, obtained by combining the various flow-encoding gradients. The small matrix obtained with as few as two encoding steps can be extended by zero-filling in all N dimensions and using ND-Fourier transformation to obtain the maximum of the resulting peak in the ND-matrix, which gives simultaneously all the components of velocity and/or acceleration. The processing time was shortened by using a method of phase computation that gives the same precision as Fourier transformation, but is much faster. A rotating disk was used to show that the velocity-to-noise ratio increases with the number of dimensions acquired, demonstrating the efficiency of multidimensional flow measurements. The feasibility of the method is illustrated by 3D maps of the myocardium velocity, and 2D measurement of velocity and acceleration in the ascending aorta-both obtained by multidimensional phase encoding in volunteers.     

Magnetic resonance phase velocity mapping through NiTi stents in a flow phantom model

PURPOSE: To assess constant and pulsatile flow velocity within the lumen of a peripheral NiTi stent using phase velocity mapping for comparison with independent assessments of flow velocity in a phantom model. MATERIALS AND METHODS: A 9 x 20-mm stent installed in flexible tubing was placed in a phantom filled with stationary fluid. Constant and pulsatile flow (produced by a pump programmed to produce a simulation of the carotid artery flow) was assessed using phase velocity mapping at 4.1 T (for constant flow) and at 1.5 T (for pulsatile flow). In all cases 256 x 256 gradient echo phase velocity maps were acquired. For the pulsatile flow condition, cine images with acquisition gated to the pump cycle were acquired with 40 msec temporal resolution across the simulated cardiac cycle. Computed flow volume rates were compared with fluid volume collection for the constant flow model, and with ultrasonic Doppler flow meter measurements for the pulsatile model. RESULTS: The data showed that volume flow rate assessments by phase velocity mapping agreed with independent measurements within 10% to 15%. CONCLUSION: Phase velocity mapping of the lumen of peripheral size NiTi stents is possible in an in vitro model.     

Real time blood flow imaging by spiral scan phase velocity mapping

The work describes the development of a novel sequence that uses rapid spiral k-space sampling, combined with phase velocity mapping, for real time flow velocity imaging. The performance of the technique is assessed on phantoms for both through-plane and in-plane flows. The flow measurements compared well with those measured using a bucket and stopwatch. One advantage of the technique is that flow related signal loss is minimal due to the early acquisition of the center of k-space data. Flow artifacts were observed for in-plane flow and these were understood with the aid of computer simulations. In vivo studies involved cine velocity mapping in normal volunteers; aortic blood flow waveforms acquired by spiral scanning in two cardiac cycles compared well with data from a conventional gradient-echo sequence. Potential applications of the method are demonstrated by studying the response of aortic flow to physical exercise and the real time monitoring of aortic flow during a valsalver maneuver.     

Mitral and aortic valvular flow: quantification with MR phase mapping.

When magnetic resonance phase mapping is used to quantitate valvular blood flow, the presence of higher-order-motion terms may cause a loss of phase information. To overcome this problem, a sequence with reduced encoding for higher-order motion was used, achieved by decreasing the duration of the flow-encoding gradient to 2.2 msec. Tested on a flow phantom simulating a severe valvular stenosis, the sequence was found to be robust for higher-order motion within the clinical velocity range. In eight healthy volunteers, mitral and aortic volume flow rates and peak velocities were quantified by means of phase mapping and compared with results of the indicator-dilution technique and Doppler echocardiography, respectively. Statistically significant correlations were found between phase mapping and the other two techniques. Similar studies in patients with valvular disease indicate that phase mapping is also valid for pathologic conditions. Phase mapping may be used as a noninvasive clinical tool for flow quantification in heart valve disease.     

Fast T(1) mapping with volume coverage.

Four different sequences which enable high-resolution, multislice T(1) relaxation-time mapping are presented. All these sequences are based on the Look-Locker method with differences arising from the use of either a saturation-recovery or inversion-recovery module prior to data acquisition with a full k-space or banded k-space acquisition scheme. The methods were implemented on a standard clinical scanner and the accuracy of the T(1) results was evaluated against spectroscopic measurements. The accuracy of the T(1) maps validated by phantom imaging measurements is around 1% for species which relax with T(1) times that mimic gray/white matter (T(1) < or = 1000 ms). Additionally, the inherent multislice, multipoint capability of the methods is demonstrated. Finally, in vivo results of the human brain obtained using the faster method are presented. The fastest data acquisition was achieved with a saturation-recovery, banded k-space method where k-space was divided into three segments; an overall acquisition time of around 5 min (for species with T(1) < or = 1 sec) was achieved for a T(1) map which can, in principle, provide whole-brain coverage with a matrix size of 256 x 256 at multiple time-points. Magn Reson Med 46:131-140, 2001.     

Strategies to improve contrast in turboFLASH imaging: reordered phase encoding and k-space segmentation.

TurboFLASH (fast low-angle shot) sequences enable the acquisition of an image in a fraction of a second. However, unique to T1-weighted ultrafast imaging, the magnetization variation during image acquisition can produce artifacts along the phase-encoding direction. In this study, the signal behavior and nature of these artifacts were analyzed with various acquisition schemes to improve image contrast. The magnetization variation during image acquisition and its filtering effect on the image were simulated for three different approaches to T1-weighted turboFLASH imaging: standard turboFLASH with (a) monotonically ascending phase-encoding steps, (b) reordered phase encoding, and (c) k-space segmentation. Each of the modified data acquisition schemes has advantages. However, for subsecond imaging, reordered phase encoding produced improved image contrast over that of standard turboFLASH, and segmented k-space imaging gave superior tissue contrast compared with that of both standard and reordered turboFLASH, with imaging time that permits breath-hold studies.     

Velocity encoding using both phase and magnitude

Imaging time constitutes a major limitation of phase-contrast (PC) angiography. It is possibly the main disadvantage of PC methods over the time-of-flight (TOF) methods that actually are used clinically. This relatively long imaging time comes from the fact that conventional PC methods require the acquisition of at least four images with different velocity sensitization to reconstruct a single angiogram (1, 2). However, more than one-half of the information gathered through the acquisition of these four images is either redundant or simply discarded. We propose a faster approach to making PC angiograms in which the quantity of data acquired is diminished by as much as a factor 2. This is made possible by encoding velocity information in both the phase and magnitude of the image. Due to the use of extra radiofrequency (RF) and gradient waveforms, decreases in data requirements do not translate in a direct manner into decreases in imaging time. Nevertheless, significant reductions in imaging time are achieved with the present approach.     

In vivo time-resolved quantitative motion mapping of the murine myocardium with phase contrast MRI.

Myocardial motion of healthy mice and mice with myocardial infarction was assessed in vivo by phase contrast (PC) cine MRI. The imaging module was a segmented fast low angle shot (FLASH) sequence with velocity compensation in all three gradient directions. To accomplish additional motion encoding, the spin phase was prepared using bipolar gradient pulses, which resulted in a linear dependence between the voxel velocity and spin phase. This method provided accurate quantification of the velocity magnitude and direction of the murine myocardium at a spatial resolution of 234 microm and a temporal resolution of about 10 ms. The acquisition was EKG-gated and the mice were anesthetized by inhalation of 1.5-4.0 vol.% isoflurane at 1.5 l/min oxygen flow. To validate the MRI measurements, an experiment with a calibrated rotating phantom was performed. Deviations between MR velocity measurements and optical assessment by a light detector were lower than 1.6%. During our study, myocardial motion velocities between 0.4 cm/s and 1.7 cm/s were determined for the healthy murine myocardium across the heart cycle. Areas with myocardial infarction were clearly segmented and showed a motion velocity which was significantly reduced. In conclusion, the method is an accurate technique for the assessment of murine myocardial motion in vivo.     

Flexible cardiac T1 mapping using a modified Look-Locker acquisition with saturation recovery

A modified Look-Locker acquisition using saturation recovery (MLLSR) for breath-held myocardial T(1) mapping is presented. Despite its reduced dynamic range, saturation recovery enables substantially higher imaging efficiency than conventional inversion recovery T(1) mapping because it does not require time for magnetization to relax to equilibrium. Therefore, MLLSR enables segmented readouts, shorter data acquisition windows, and shorter breath holds compared with inversion recovery. T(1) measurements in phantoms using MLLSR showed a high correlation with conventional single-point inversion recovery spin echo. In vivo T(1) measurements from normal and infarcted myocardium in 41 volunteers and patients were consistent with previously reported values. Twenty subjects were also scanned with MLLSR using an accelerated sampling scheme that required half the scan time (eight vs. 16 heartbeats) but yielded equivalent results. The flexibility afforded by the improved imaging efficiency of MLLSR allows the acquisition to be tailored to particular clinical needs and to individual patient's breath-holding abilities.     

Navigator gated high temporal resolution tissue phase mapping of myocardial motion.

Data acquisition for phase contrast velocity mapping of myocardial motion is typically based on multiple breath-held 2D measurements with limited acquisition duration and consequently relatively poor temporal resolution. In order to overcome the limitations of breath-hold acquisitions, an improved navigator-guided technique was implemented based on 2 navigator signals within each cardiac cycle in combination with paired acceptance and rejection criteria of successive navigator signals. Respiratory gated phase contrast measurements with 3-directional velocity encoding were performed in 12 healthy volunteers in basal, midventricular, and apical locations of the left ventricle during free breathing with a temporal resolution of 13.8 ms. Results were compared to standard breath-hold measurements with a temporal resolution of 69 ms. Data from the high temporal resolution study revealed details in left ventricular motion patterns that were previously not seen in phase contrast measurements and are only known from echocardiography. The proposed navigator gated technique for high temporal resolution velocity mapping is, therefore, highly promising for the detection of local and global motion abnormalities in patients with disturbed left ventricular performance, such as diastolic dysfunction.     

Rapid combined T1 and T2 mapping using gradient recalled acquisition in the steady state.

A novel, fully 3D, high-resolution T(1) and T(2) relaxation time mapping method is presented. The method is based on steady-state imaging with T(1) and T(2) information derived from either spoiling or fully refocusing the transverse magnetization following each excitation pulse. T(1) is extracted from a pair of spoiled gradient recalled echo (SPGR) images acquired at optimized flip angles. This T(1) information is combined with two refocused steady-state free precession (SSFP) images to determine T(2). T(1) and T(2) accuracy was evaluated against inversion recovery (IR) and spin-echo (SE) results, respectively. Error within the T(1) and T(2) maps, determined from both phantom and in vivo measurements, is approximately 7% for T(1) between 300 and 2000 ms and 7% for T(2) between 30 and 150 ms. The efficiency of the method, defined as the signal-to-noise ratio (SNR) of the final map per voxel volume per square root scan time, was evaluated against alternative mapping methods. With an efficiency of three times that of multipoint IR and three times that of multiecho SE, our combined approach represents the most efficient of those examined. Acquisition time for a whole brain T(1) map (25 x 25 x 10 cm) is less than 8 min with 1 mm(3) isotropic voxels. An additional 7 min is required for an identically sized T(2) map and postprocessing time is less than 1 min on a 1 GHz PIII PC. The method therefore permits real-time clinical acquisition and display of whole brain T(1) and T(2) maps for the first time.     

Assessment of regional left ventricular long-axis motion with MR velocity mapping in healthy subjects

The pattern of left ventricular long-axis motion during early diastole was assessed with magnetic resonance (MR) velocity mapping in 31 healthy volunteers. Regional long-axis velocity varied with time and position around the ventricle. During systole, the base descended toward the apex. The greatest magnitude of long-axis velocity occurred during early diastole. The lateral wall had the highest velocity (140 mm/sec ± 40 [mean ± standard deviation]); the anterior and inferior walls had lower velocities (96 mm/sec ± 27 and 92 mm/sec ± 34, respectively). The inferoseptal area consistently had the lowest velocities (87 mm/sec ± 40). Absolute values of peak early-diastolic velocity declined with age (r = ?.64, P <.001). Peak early-diastolic velocity was not dependent on heart rate (r =.014, P =.94). Regional variations in left ventricular wall motion were seen. MR velocity mapping is a useful technique for assessing regional left ventricular long-axis heart function.     

Interleaved pulsed MAMBA: a new parallel slice imaging method.

A method of acquiring slices in parallel is described which uses interleaved sets of pulsed B(0) field coils to generate discrete regions of uniform field within the main magnetic field known as interleaved MAMBA (multiple acquisition micro B(0) array). Simulations of a number of coil designs were performed using the Biot-Savart law. A six-step coil was built and interfaced to a 0.17 T Niche MRI system and the field steps measured using an imaging technique. Measured field steps were in good agreement with the values predicted by simulation. The coil design was then scaled up by a factor of three, interfaced to a 1.5 T whole-body MRI system, and scans of the hands and arms of volunteers were acquired from up to four field steps using standard spin and gradient echo sequences. Images were also acquired simultaneously from two field steps with no frequency encode aliasing and one excitation. The one-dimensional interleaved pulsed MAMBA step field technique shows great promise for enabling many slices to be acquired simultaneously along the axis of the coil for rapid volumetric studies without the need for multiple shot Hadamard encoding. Extension of interleaved coil design to two or three dimensions is feasible, which could provide full spatial coverage combined with ultra-rapid data acquisition.     

Time-resolved contrast-enhanced magnetic resonance angiography in pediatric patients using sensitivity encoding

PURPOSE: To evaluate the role of time-resolved contrast-enhanced magnetic resonance angiography (CE-MRA) using sensitivity encoding in imaging the thoraco-abdominal vessels in pediatric patients. MATERIALS AND METHODS: Thoraco-abdominal vessels of 22 pediatric patients (median age = 5 years) were evaluated with a 3D CE-MRA technique in combination with SENSE following a 0.2 mmol/kg injection of Gd-chelate. The acquisition parameters were as follows: TR/TE = 5/1.1 msec; flip angle = 40 degrees; in-plane phase encoding steps were reduced by a factor of 2 using sensitivity encoding (SENSE); 3D volume acquisition was repeated four to eight times consecutively during free breathing (four to eight dynamics) with a mean temporal resolution of 6.8 seconds/dynamic; and mean acquired voxel size = 1.4 x 1.7 x 3.1 mm (reconstructed as 1.4 x 1.4 x 1.55 mm). Arterial-to-venous signal intensity ratios (AVRs) were computed for each dynamic. RESULTS: All images were successfully reconstructed and were of diagnostic quality. The AVRs of prepeak, peak, and postpeak arterial volumes were 1.0 +/- 0.5, 6.1 +/- 3.3, and 1.3 +/-0.9, respectively, indicating good arterial-to-venous separation. The signal-to-noise ratio (SNR) of the peak arterial volume was 41 +/- 26. CONCLUSION: Our results suggest that it is feasible to apply SENSE to a conventional 3D CE-MRA technique in a time-resolved fashion for imaging the thoraco-abdominal vessels in pediatric patients during free breathing.     

Tracking of cyclic motion with phase-contrast cine MR velocity data

A method of computing trajectories of objects by using velocity data, particularly as acquired with phase-contrast magnetic resonance (MR) imaging, is presented. Starting from a specified location at one time point, the method recursively estimates the trajectory. The effects of measurement noise and eddy current-induced velocity offsets are analyzed. When the motion is periodic, trajectories can be computed by integrating in both the forward and backward temporal directions, and a linear combination of these trajectories minimizes the effect of velocity offsets and maximizes the precision of the combined trajectory. For representative acquisition parameters and signal-to- noise ratios, the limitations due to measurement noise are acceptable. In a phantom with reciprocal rotation, the measured and true trajectories agreed to within 3.3%. Sample trajectory estimates of human myocardial regions are encouraging.     

Phase-contrast pulsatility measurements: Preliminary results in normal and arteriosclerotic iliofemoral arteries work in progress

Magnetic resonance (MR) flow measurements were obtained in six healthy volunteers and 30 patients with arteriosclerotic disease with a 1.5-T imager and a pulse sequence for flow quantification based on flow-induced phase shifts. The iliac arteries were investigated in eight and the femoral arteries in 28 subjects. A trigger pulse, followed by the acquisition of 30 evenly distributed data sets, was applied every second heartbeat, thus eliminating any acquisition gap in a full heart cycle. For quantitative analysis, flow velocity was plotted as a function of time. Systolic acceleration, postsystolic deceleration, and pulsatility of flow were calculated and compared with stenosis grades determined from recent intraarterial digital subtraction angiograms. The flattening of the temporal flow patterns correlated with local severity of vascular occlusive disease.     

MR imaging of sodium in the human brain with a fast three-dimensional gradient-recalled-echo sequence at 4 T

RATIONALE AND OBJECTIVES: Sodium ions play a vital role in cellular homeostasis and electrochemical activity throughout the human body. However, the in vivo detection of sodium (23Na) with magnetic resonance (MR) techniques is hindered by the fast transverse relaxation, low tissue equivalent concentration, and small gyromagnetic ratio of sodium ions compared with protons (1H). The goals of this study were to acquire MR images of sodium in the whole human brain by using a fast three-dimensional gradient-recalled-echo sequence and to investigate the effect that restrictions on specific absorption ratio have on MR imaging of sodium at 4 T. MATERIALS AND METHODS: A three-dimensional gradient-recalled-echo sequence with short echo time was developed for MR imaging of sodium. Slab encoding was removed and a hard excitation pulse was used. Five healthy human volunteers were examined in a whole-body MR imager with the use of a custom transmit-and-receive birdcage coil. Fields of view were selected to cover the entire brain: 38 x 38 cm in the axial plane, with 24 sections of 5.8 mm each or 12 sections of 1.1 cm each. The in-plane acquisition matrix was 64 x 128, and voxel size was 0.2 cm(3). RESULTS: Sodium in white matter was depicted with an acceptable signal-to-noise ratio of 20-25. The echo time, and hence the signal-to-noise ratio, was limited by the MR imager's maximum allowable gradient strength. To keep the specific absorption ratio below 3 W/kg (the limit established by the Food and Drug Administration), it was necessary to prolong the repetition time to 30 msec. CONCLUSION: The MR imaging protocol used in this study provided acceptable visualization of sodium in the whole brain in a tolerable total acquisition time of 15 minutes.