A complex-difference phase-contrast technique for measurement of volume flow rates

Magnetic resonance (MR) phase-difference methods work well for measuring volumetric flow rates when the vessel diameter is large compared with the in-plane voxel dimensions. For small vessels (eg, coronary arteries), partial-volume effects introduce substantial errors in the measured volume flow rate. To correctly measure flow rates through a voxel, both the fraction of the voxel containing moving spins and the phase shift imparted to those spins must be known. The authors propose a flow measurement method that combines information obtained with both the complex-difference and phase-difference processing techniques and thereby provides the fractional volume occupied by the moving spins and the phase of those spins. The complex-difference flow map method proposed results in improved accuracy of MR phase-contrast flow measurements in the presence of partial-volume effects.     

Rapid dark-blood carotid vessel-wall imaging with random bipolar gradients in a radial SSFP acquisition

PURPOSE: To investigate and evaluate a new rapid dark-blood vessel-wall imaging method using random bipolar gradients with a radial steady-state free precession (SSFP) acquisition in carotid applications. MATERIALS AND METHODS: The carotid artery bifurcations of four asymptomatic volunteers (28-37 years old, mean age = 31 years) were included in this study. Dark-blood contrast was achieved through the use of random bipolar gradients applied prior to the signal acquisition of each radial projection in a balanced SSFP acquisition. The resulting phase variation for moving spins established significant destructive interference in the low-frequency region of k-space. This phase variation resulted in a net nulling of the signal from flowing spins, while the bipolar gradients had a minimal effect on the static spins. The net effect was that the regular SSFP signal amplitude (SA) in stationary tissues was preserved while dark-blood contrast was achieved for moving spins. In this implementation, application of the random bipolar gradient pulses along all three spatial directions nulled the signal from both in-plane and through-plane flow in phantom and in vivo studies. RESULTS: In vivo imaging trials confirmed that dark-blood contrast can be achieved with the radial random bipolar SSFP method, thereby substantially reversing the vessel-to-lumen contrast-to-noise ratio (CNR) of a conventional rectilinear SSFP "bright-blood" acquisition from bright blood to dark blood with only a modest increase in TR (approximately 4 msec) to accommodate the additional bipolar gradients. CONCLUSION: Overall, this sequence offers a simple and effective dark-blood contrast mechanism for high-SNR SSFP acquisitions in vessel wall imaging within a short acquisition time.     

Single scan PC‐MRI by alternating the velocity encoding gradient polarity between phase encoding steps

The purpose of this study was to develop a faster approach to phase contrast magnetic resonance imaging. This article proposes a phase contrast imaging scheme called single scan phase contrast in which the polarity of the velocity‐encoding gradient is alternated between phase encoding steps. In single scan phase contrast, ghost images due to moving spins form. The signal intensity of the ghost images is modulated by the sine of the motion‐induce phase shift. Prior to image acquisition, the region of interest containing moving spins is identified, and the field of view is configured so to avoid overlap between the object in the image and the ghost image(s) due to motion in the region of interest. The image values of the region of interest and the ghost image are used to quantify velocity. At best, single scan phase contrast reduces the total acquisition time by a factor of two when compared to phase contrast. In this study, single scan phase contrast is validated against phase contrast in phantom and in vivo. Magn Reson Med, 2011. © 2011 Wiley‐Liss, Inc.     

Combined analysis of spatial and velocity displacement artifacts in phase contrast measurements of complex flows

MR phase contrast (PC) velocity imaging is a promising tool for quantifying blood flow velocity in vivo. PC velocity imaging is, however, susceptible to artifacts that result from the displacement of spins during the finite duration pulse sequences. Such displacement artifacts can lead to errors in velocity measurements, especially in the presence of oblique and accelerating flows, which are common throughout the cardiovascular system. By tracking particles (representing spins) through a computed velocity field, and assuming that spatial and velocity encodings occur at discrete times during the pulse sequence, we simulate the separate and combined effects of oblique and acceleration artifacts on PC velocity images. We demonstrate, both by simulation and MR measurement, the errors associated with such artifacts in PC velocity measurements in a representative flow geometry. Using example particle trajectories, we provide a fluid dynamic basis for characteristic phase-velocity image distortions that can arise when imaging complex, physiologically relevant flows.     

Time‐resolved absolute velocity quantification with projections

Quantitative information on time‐resolved blood velocity along the femoral/popliteal artery can provide clinical information on peripheral arterial disease and complement MR angiography as not all stenoses are hemodynamically significant. The key disadvantages of the most widely used approach to time‐resolve pulsatile blood flow by cardiac‐gated velocity‐encoded gradient‐echo imaging are gating errors and long acquisition time. Here, we demonstrate a rapid nontriggered method that quantifies absolute velocity on the basis of phase difference between successive velocity‐encoded projections after selectively removing the background static tissue signal via a reference image. The tissue signal from the reference image's center k‐space line is isolated by masking out the vessels in the image domain. The performance of the technique, in terms of reproducibility and agreement with results obtained with conventional phase contrast‐MRI was evaluated at 3 T field strength with a variable‐flow rate phantom and in vivo of the triphasic velocity waveforms at several segments along the femoral and popliteal arteries. Additionally, time‐resolved flow velocity was quantified in five healthy subjects and compared against gated phase contrast‐MRI results. To illustrate clinical feasibility, the proposed method was shown to be able to identify hemodynamic abnormalities and impaired reactivity in a diseased femoral artery. For both phantom and in vivo studies, velocity measurements were within 1.5 cm/s, and the coefficient of variation was less than 5% in an in vivo reproducibility study. In five healthy subjects, the average differences in mean peak velocities and their temporal locations were within 1 cm/s and 10 ms compared to gated phase contrast‐MRI. In conclusion, the proposed method provides temporally resolved arterial velocity with a temporal resolution of 20 ms with minimal post processing. Magn Reson Med, 2010. © 2010 Wiley‐Liss, Inc.     

Nonsubtractive spiral phase contrast velocity imaging.

Phase contrast velocity imaging is a standard method for accurate in vivo flow measurement. One drawback, however, is that it lengthens the scan time (or reduces the achievable temporal resolution) because one has to acquire two or more images with different flow sensitivities and subtract their phases to produce the final velocity image. Without this step, non-flow-related phase variations will give rise to an erroneous, spatially varying background velocity. In this paper, we introduce a novel phase contrast velocity imaging technique that requires the acquisition of only a single image. The idea is to estimate the background phase variation from the flow-encoded image itself and then have it removed, leaving only the flow-related phase to generate a corrected flow image. This technique is sensitive to flow in one direction and requires 50% less scan time than conventional phase contrast velocity imaging. Phantom and in vivo results were obtained and compared with those of the conventional method, demonstrating the new method's effectiveness in measuring flow in various vessels of the body. Magn Reson Med 42:704-713, 1999.     

Preferential arterial imaging using gated thick-slice gadolinium-enhanced phase-contrast acquisition in peripheral MRA

PURPOSE: To investigate the feasibility of preferential arterial imaging using gadolinium-enhanced thick-slice phase-contrast imaging. METHODS: Six healthy volunteers were studied using a peripheral-gated segmented k-space CINE phase-contrast pulse sequence using four views per RR interval with flow encoding in the superior-inferior direction. Images at the level of the popiteal trifurcation were acquired postcontrast with different section thicknesses (4-8 cm) and VENC values (20-150 cm/sec), and phase-difference processing. RESULTS: The post-gadolinium contrast-enhanced thick-slice phase-contrast acquisitions demonstrated the ability to visualize the tibio-peroneal (trifurcation) arteries, especially in systole. With MR contrast agents, the signal from blood is raised significantly above that of stationary tissue from T(1) shortening such that the partial volume artifact is reduced in thick-slice acquisitions. Furthermore, by selecting the VENC value as a function of the cardiac cycle, the noise floor can be raised to selectively suppress flow values less than that of the noise threshold, allowing better accentuation of arterial structures at systole. CONCLUSIONS: Thick-slice phase-contrast acquisition with phase-difference processing has been observed to reduce partial volume artifacts when an MR contrast agent substantially increases signal in the vasculature over that of normal background tissue. Preferential arterial images can be obtained by either increasing the VENC value to selectively suppress signal from slow flow in the veins or by subtracting the diastolic phase image from the peak systolic phase image. J. Magn. Reson. Imaging 2001;13:714-721.     

Rapid mapping of flow velocity using a new PARSE method.

A new method for flow velocity mapping is presented here. Instead of the conventional approach of employing two images (velocity sensitive and control) to generate velocity information, in the new method one determines the velocity directly from a single-shot acquisition by solving an inverse problem. This technique is a variant of single-shot parameter assessment by retrieval from signal encoding (SS-PARSE). The results of simulation and phantom studies show strong agreement with the actual velocities. The prototype method can measure velocities in the range of -50 to 50 cm/s, which is roughly appropriate for future applications in dynamic blood flow measurement in carotid arteries.     

Blood velocity assessment using 3D bright-blood time-resolved magnetic resonance angiography.

Blood velocity is a functional parameter that is not easily assessed noninvasively, especially in small animals. A new noninvasive method that uses magnetic resonance angiography (MRA) to measure blood flows is proposed. This method is based on the time-of-flight (TOF) phenomenon. By initially suppressing the signal from the stationary spins in the area of interest, it is possible to sequentially visualize only the signal from the moving spins entering a given volume. With this method, 3D cine images of the blood flow can be generated by positive contrast, with unparalleled spatial (<200 microm) and temporal resolutions (<10 ms/image). As a result, it is possible to measure flow in sinuous paths. The present method was applied in vivo to measure the blood velocity in mouse carotid arteries. Because of its robustness and simplicity of implementation, this method has numerous potential applications for fundamental studies in small animal models.     

Multistep phase difference phase contrast imaging

A new technique for multistep phase-contrast image processing is presented. The N-step method consists of simply forming the linear average of the N — 1 adjacent phase-difference signals. It has similar noise reduction properties as other multistep techniques, but the simplicity of the noise variance of the N-step technique allows intuitive insight into phase-difference phase-contrast processing and noise reduction, which can aid in the design of efficient and improved phase-contrast imaging sequences. As well, the computational simplicity of the N-step phase-difference technique compared with any other known multistep technique is advantageous. Like other multistep techniques, it has far more efficient noise reduction properties than simple two-step, multiple average phase-contrast imaging, even when normalized for total scan time. A three-step phase-difference velocity image has 50% less variance than an image acquired with two steps and two scans averaged but is obtained in 25% less scan time. Given its advantages, it should now be the chosen technique for increasing velocity-to-noise and contrast-to-noise ratios in all phase-difference phase-contrast clinical applications.     

High temporal resolution phase contrast MRI with multiecho acquisitions.

Velocity imaging with phase contrast (PC) MRI is a noninvasive tool for quantitative blood flow measurement in vivo. A shortcoming of conventional PC imaging is the reduction in temporal resolution as compared to the corresponding magnitude imaging. For the measurement of velocity in a single direction, the temporal resolution is halved because one must acquire two differentially flow-encoded images for every PC image frame to subtract out non-velocity-related image phase information. In this study, a high temporal resolution PC technique which retains both the spatial resolution and breath-hold length of conventional magnitude imaging is presented. Improvement by a factor of 2 in the temporal resolution was achieved by acquiring the differentially flow-encoded images in separate breath-holds rather than interleaved within a single breath-hold. Additionally, a multiecho readout was incorporated into the PC experiment to acquire more views per unit time than is possible with the single gradient-echo technique. A total improvement in temporal resolution by approximately 5 times over conventional PC imaging was achieved. A complete set of images containing velocity data in all three directions was acquired in four breath-holds, with a temporal resolution of 11.2 ms and an in-plane spatial resolution of 2 mm x 2 mm.     

Real-time blood flow imaging using autocalibrated spiral sensitivity encoding.

A novel spiral phase contrast (PC) technique was developed for high temporal resolution imaging of blood flow without cardiac gating. An autocalibrated spiral sensitivity encoding (SENSE) method is introduced and used to reconstruct PC images. Numerical simulations and a flow phantom study were performed to validate the technique. To study the accuracy of the flow measurement in vivo, a high-resolution cardiac experiment was performed and a subset of undersampled SENSE reconstructed data were reconstructed. Good agreement between the velocity measurement from the fully-sampled and undersampled data was achieved. Real-time experiments were performed to measure blood velocity in the ascending aorta and aortic valve, and during a Valsalva maneuver. The results demonstrate the potential of this technique for real-time flow imaging.     

Phase contrast using multiecho steady-state free precession.

Conventional phase-contrast (PC) MRI is limited in the temporal resolution (typically 50 ms) that can be achieved, due to the need to implement bipolar velocity encoding gradients. PC using steady-state free precession (SSFP) has recently been developed to acquire PC data at higher rates without sacrificing contrast-to-noise ratio (CNR). This work presents two multiecho SSFP PC implementations that can be used to increase the time efficiency of PCSSFP. Both approaches (extrinsic and intrinsic) enable reference image lines to be acquired within the same TR as the flow-encoded lines, thus minimizing the scan time and permitting TR-equivalent temporal resolutions. Both approaches have been implemented and tested successfully on human volunteers at 1.5T and 3T. While the intrinsic approach is useful for encoding higher velocity flows in-plane, the extrinsic implementation can be used for studying a wider range of encoding velocities for flow in the imaging plane and through the imaging plane.     

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.     

Dixon techniques in spiral trajectories with off-resonance correction: a new approach for fat signal suppression without spatial-spectral RF pulses.

Spiral imaging has recently gained acceptance in MR applications requiring rapid data acquisition. One of the main disadvantages of spiral imaging, however, is blurring artifacts that result from off-resonance effects. Spatial-spectral (SPSP) pulses are commonly used to suppress those spins that are chemically shifted from water and lead to off-resonance artifacts. However, SPSP pulses may produce nonuniform fat signal suppression or unwanted water signal suppression when applied in the presence of B(0) field inhomogeneities. Dixon techniques have been developed as methods for water-fat signal decomposition in rectilinear sampling schemes since they can produce unequivocal water-fat signal decomposition even in the presence of B(0) inhomogeneities. This article demonstrates that three-point and two-point Dixon techniques can be extended to conventional spiral and variable-density spiral data acquisitions for unambiguous water-fat decomposition with off-resonance blurring correction. In the spiral three-point Dixon technique, water-fat signal decomposition and image deblurring are performed based on the frequency maps that are directly derived from the acquired images. In the spiral two-point Dixon technique, several predetermined frequencies are tested to create a frequency map. The newly proposed techniques can achieve more effective and more uniform fat signal suppression when compared to the conventional spiral acquisition method with SPSP pulses.     

Adapting MRI acoustic radiation force imaging for in vivo human brain focused ultrasound applications

A variety of magnetic resonance imaging acoustic radiation force imaging (MR‐ARFI) pulse sequences as the means for image guidance of focused ultrasound therapy have been recently developed and tested ex vivo and in animal models. To successfully translate MR‐ARFI guidance into human applications, ensuring that MR‐ARFI provides satisfactory image quality in the presence of patient motion and deposits safe amount of ultrasound energy during image acquisition is necessary. The first aim of this work was to study the effect of motion on in vivo displacement images of the brain obtained with 2D Fourier transform spin echo MR‐ARFI. Repeated bipolar displacement encoding configuration was shown less sensitive to organ motion. The optimal signal‐to‐noise ratio of displacement images was found for the duration of encoding gradients of 12 ms. The second aim was to further optimize the displacement signal‐to‐noise ratio for a particular tissue type by setting the time offset between the ultrasound emission and encoding based on the tissue response to acoustic radiation force. A method for measuring tissue response noninvasively was demonstrated. Finally, a new method for simultaneous monitoring of tissue heating during MR‐ARFI acquisition was presented to enable timely adjustment of the ultrasound energy aimed at ensuring the safety of the MR‐ARFI acquisition. Magn Reson Med, 2013. © 2012 Wiley Periodicals, Inc.     

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.     

Blood attenuation with SSFP-compatible saturation (BASS)

PURPOSE: To investigate a rapid flow-suppression method for improving the contrast-to-noise ratio (CNR) between the vessel wall and the lumen for cardiovascular imaging applications. MATERIALS AND METHODS: In this study a new dark-blood steady-state free precession (SSFP) sequence utilizing two excitation pulses per TR was developed. The first pulse is applied immediately adjacent to the slice of interest, while the second is a conventional slice-selective pulse designed to excite an SSFP signal for the static spins in the slice of interest. The slice-selective pulse is followed by fully refocused gradients along all three imaging axes over each TR. The signal amplitude (SA) from the moving spins excited by the "saturation" pulse is attenuated since they are not fully refocused at the TE. RESULTS: This work provides confirmation, by both simulation and experiments, that modest adaptations of the basic True-FISP structure can limit unwanted "bright blood" signal within the vessels while simultaneously preserving the contrast and speed advantages of this well-established rapid imaging method. CONCLUSION: Animal imaging trials confirm that dark-blood contrast is achieved with the BASS sequence, which substantially reverses the lumen-to-muscle CNR of a conventional True-FISP "bright blood" acquisition from 14.77 (bright blood) to -13.96 (dark blood) with a modest increase (24.2% of regular TR of SSFP for this implementation) in acquisition time to accommodate the additional slab-selective excitation pulse and gradient pulses.     

Ultra-fast velocity imaging in stenotically produced turbulent jets using RUFIS

A method for rapidly producing velocity images is presented. This sequence combines a modified bipolar gradient pulse to magnitude encode the velocity with the rotating ultra-fast imaging sequence (RUFIS) to image the encoded spins. Velocity encoding is done in 3 msec, and RUFIS acquires 32 projections in 8 msec. The method is applied to turbulent jets associated with a 75% stenosis in a 15-mm inner diameter glass pipe. Data is acquired upstream and downstream from the stenosis for Reynolds numbers from 560 to 3750. In addition, a robust method of reconstructing the unobserved short time region of a free induction decay is presented and incorporated into the image processing.     

Three‐dimensional biexponential weighted 23Na imaging of the human brain with higher SNR and shorter acquisition time

A new method is presented for acquiring 3D biexponential weighted sodium images of the in vivo human brain with up to three times higher signal‐to‐noise ratio compared with conventional six‐step phase‐cycling triple‐quantum‐filtered imaging. To excite and detect multiple‐quantum coherences, a three‐pulse preparation is used. During the pulse train, two images are obtained. The first image is acquired with ultrashort echo time (0.3 ms) during preparation between the first two pulses to yield a spin‐density‐weighted image. After the last pulse, a single‐quantum‐filtered image is acquired with an echo time of 11 ms that maximizes the resulting signal. The biexponential weighted image is calculated by subtracting the single‐quantum‐filtered image from the spin‐density‐weighted image. The resulting image mainly shows signal from sodium ions with biexponential quadrupolar relaxation behavior. In isotropic environments, the resulting image mainly contains triple‐quantum‐filtered signal. The four‐step phase cycling yields similar signal‐to‐noise ratio in shorter acquisition time compared with six‐step phase‐cycling biexponential weighted imaging. Magn Reson Med 70:754–765, 2013. © 2012 Wiley Periodicals, Inc.