Certain magnetic resonance balanced steady-state free precession (bSSFP) imaging applications require large and frequent changes in the gradient waveforms from one TR to the next. One such application of particular interest in our research group is bSSFP phase-contrast (PC) magnetic resonance imaging (MRI), which utilizes two or more images acquired with different gradient first moments at the echo time (TE). This can be achieved by switching between a positive and negative bipolar gradient pulse played before or after the data acquisition window, or by inverting the sign of the imaging gradients themselves. Unfortunately, different gradient wave-forms will in general produce slightly different residual eddy-current fields, which in turn give rise to small differences in the precession angle. As a result, switching between two waveforms on a TR-by-TR basis can produce an “oscillating” steady state for spins near on-resonance, resulting in severe image artifacts. For this reason, bSSFP PC-MRI has been limited to noninterleaved imaging, which limits the application of this technique to non real-time acquisitions, and makes the measurements more susceptible to patient motion.
It has been shown that eddy-current-induced bSSFP signal distortions for on-resonance spins can be mitigated by pairing the waveforms. See Bieri O., et al., “Analysis and compensation of eddy currents in balanced SSFP,” Magn. Reson. Med. 2005;54:129-137; and U.S. Pat. No. 7,046,004; the entire contents of both of which are incorporated herein by reference. Waveform pairing can ensure that any eddy-current-induced change in precession angle in one TR is simply reproduced in the following TR. This tends to bring spins that are near on-resonance back to the same starting point, regardless of the precise value of the small and unknown eddy-current-induced change in precession angle.
Waveform pairing has been applied to several bSSFP imaging applications including phase-contrast imaging, however, paired interleaved bSSFP introduces signal distortions approximately halfway between the center of the bSSFP passband and the signal nulls. In order to suppress within-band signal distortions in paired bSSFP, Markl et al. proposed a “double average” (dAVE) acquisition scheme in which every phase-encode is acquired twice and the complex signals are averaged. See, Markl M., et al., “Double average parallel steady-state free precession imaging: optimized eddy current and transient oscillation compensation,” Magn. Reson. Med., 2005;54:965-974. Such dAVE imaging has the advantage that it simultaneously suppresses signal distortions due to both unbalanced eddy-currents and off-resonance-induced signal oscillations during the transient toward steady state, at the cost of doubling the acquisition time (for a given phase-encode acquisition scheme), however, twice the number of acquisitions and thus double the time is required.
Embodiments of the present disclosure address shortcomings described previously by providing bSSFP sequences in which two gradient waveforms are interleaved in a “groupwise” fashion, i.e., each waveform is executed consecutively two or more times before switching to the other waveform. Such gradient waveform grouping provided by embodiments of the present disclosure can be utilized to mitigate steady-state eddy-current distortions in bSSFP sequences for various types of imaging applications, including both steady-state and transient imaging.
These, as well as other components, steps, features, objects, benefits, and advantages, will now become clear from a review of the following detailed description of illustrative embodiments, the accompanying drawings, and the claims.
The drawings disclose illustrative embodiments. They do not set forth all embodiments. Other embodiments may be used in addition or instead. Details that may be apparent or unnecessary may be omitted to save space or for more effective illustration. Conversely, some embodiments may be practiced without all of the details that are disclosed. When the same numeral appears in different drawings, it refers to the same or like components or steps.
Aspects of the disclosure may be more fully understood from the following description when read together with the accompanying drawings, which are to be regarded as illustrative in nature, and not as limiting. The drawings are not necessarily to scale, emphasis instead being placed on the principles of the disclosure. In the drawings:
While certain embodiments are depicted in the drawings, one skilled in the art will appreciate that the embodiments depicted are illustrative and that variations of, those shown, as well as other embodiments described herein, may be envisioned and practiced within the scope of the present disclosure.
Illustrative embodiments are now discussed. Other embodiments may be used in addition or instead. Details that may be apparent or unnecessary may be omitted to save space or for a more effective presentation. Conversely, some embodiments may be practiced without all of the details that are disclosed.
Embodiments of the present disclosure provide for magnetic resonance imaging utilizing bSSFP sequences in which two or more gradient waveforms are interleaved in a “groupwise” fashion, i.e., each waveform is executed consecutively two or more times before switching to the other waveform (where “N” counts the number of times each waveform is executed consecutively). As a result, embodiments of the present disclosure can mitigate steady-state signal distortions/artifacts in interleaved balanced steady-state free precession (bSSFP) caused by slightly unbalanced eddy-current fields.
Embodiments of the present disclosure were tested to analyze the effects of unbalanced precession on the steady-state bSSFP magnetization, with particular attention to bSSFP PC-MRI sequences that rely on accurate blood signal phase measurements. Using numerical Bloch simulations and phantom experiments, it was shown for exemplary embodiments that steady-state signal distortions in interleaved bSSFP can be quantitatively accounted for by means of a simple model that associates an effective resonance offset with each gradient waveform. Extending the idea of pairing by executing the waveforms in groups of two or more is seen as offering a systematic and general strategy for reducing steady-state distortions within the bSSFP pass-band. Each phase-encode can be acquired only once for each sign of the bipolar velocity-encoding gradient, and consequently the overall acquisition time is not increased.
The static field generator 110 generates a strong static magnetic field B0, which is used to line up nuclear spins in a target object (whose MRI image is being generated by the apparatus 100) along B0. The gradient field generator 120 generates a gradient field G(r), which is superimposed on the static field B0, so that nuclei in a selected plane can be excited by a proper choice of the frequency spectrum of the transverse RF excitation field. The RF excitation field generator 130 generates an RF excitation field B1. When B1 is applied to the object, typically as an RF excitation pulse transverse to B0, the nuclei become excited (due to the RF energy imparted by the RF excitation pulse), so that the nuclear spins rotate by a flip angle. Subsequently, the excited nuclei gradually return to alignment with the static field B0, giving up the excitation energy in the form of weak but detectable FID signals, which are processed by the controller 140 to produce images of the target object.
The controller 140 controls the operation of the MRI apparatus 100, including but not limited to the generation of the fields G(r), and B1, as well as the processing of the FID signals resulting from the de-excitation (precession and relaxation) of the nuclei in the object.
Numerical simulations were performed to investigate the effect of gradient waveform grouping on the steady-state magnetization profile over the full 1/TR SSFP band spacing. Simulations were performed in MatLab (MathWorks, Natick, Mass.) and assumed that SSFP imaging was performed using two different waveforms “A” and “B” that are switched every N TRs. For example, N=1 corresponds to directly interleaving the waveforms every TR (“ABABABAB . . . ” ordering), and N=2 corresponds to waveform pairing (“AABBAABB . . . ”). To simulate the effect of waveform-dependent changes in precession angle, the two waveforms were associated with two slightly different resonance offset frequencies fA and fB=ΔfAB. Simulations were performed using relaxation parameters corresponding to five different tissue types at 3T (12,13): gray matter (T1/T2 =1820/99 msec), fat (T1/T2=383/68 msec), muscle (T1/T2=1412/50 msec), liver (T1/T2=812/42 msec), and blood (T1/T2=1932/275 msec).
The steady-state magnetization mss for given values of fA and ΔfAB was obtained by letting the magnetization evolve over time according to the Bloch equations and requiring the magnetization to return to the same value every 2N TRs. To facilitate these calculations, spin evolution was expressed as a matrix multiplication, which allows the steady-state signal to be obtained by solving the linear system:
m
ss
=Am
ss
+b (EQ. 1)
where A (a 3×3 matrix) and b (a 3×1 column vector) express the accumulated influence of the RF excitations, T1 relaxation, T2 decay, and free precession that occur during the course of 2N TRs.
For exemplary embodiments, phantom imaging experiments were conducted. All experiments were performed on a GE Signa 3T Excite HD system (General Electric Medical Systems, Milwaukee, Wis.) using a “zoom” whole body gradient coil (peak gradient amplitude 40 mT/m, maximum slew rate 150 T/m/s). A spherical phantom (measured T1/T2=200/30 msec) was placed at the scanner isocenter and imaged with a bSSFP pulse sequence that interleaved the two different readout waveforms illustrated schematically in
The waveforms were switched every N TRs, with N=1 (direct interleaving), 2 (pairing), or 4 (an instance of grouping). In addition, a conventional noninterleaved reference image was acquired using waveform “A” only. Imaging parameters were: 1×1×3 mm voxel size; TR=8.0 msec; flip angle 40°; 180° RF phase cycling every TR; field-of-view (FOV)=16×20 cm (for noninterleaved bSSFP, and for N=1), 16×40 cm (for N=2), and 16×80 cm (for N=4). Gradient shims were manually adjusted such that the resonance offset varied linearly across the object along the phase-encode direction.
Signal profiles were measured along the phase-encode direction. To increase measurement signal-to-noise ratio (SNR), each point in the magnetization pro-file was obtained by averaging the signal from 40 adjacent pixels along the readout direction. It was empirically determined that reducing the number of averages did not affect the profile shape. Measured signal profiles were compared with simulation results by plotting the steady-state signal amplitude (i.e., the magnitude of the transverse magnetization) and phase-difference (between echoes 1 and N+1) over the 1/TR band spacing. The signal amplitude profiles were normalized by setting the peak signal in the noninterleaved case to 1 for both simulated and measured profiles. For the calculations pertaining to the embodiment, ΔfAB was not measured independently, but acted as a (constant) fitting parameter.
Exemplary embodiments of the present disclosure can be used for cardiac, e.g., CINE, imaging. For an in vivo imaging application, cardiac flow imaging was performed in a healthy volunteer using an interleaved bSSFP phase-contrast sequence with the following parameters: ECG-gated CINE; single-slice 2DFT; flip angle 50°; 180° RF phase cycling every TR; FOV 20×25 cm; voxel size 1.7×3.1×7 mm; TR=3.6 msec; six views per cardiac phase (temporal resolution 43.6 msec); linear phase-encode view ordering; flow-encoding along readout direction, with velocity-encoding value (VENC)=4.6 m/s. The data were acquired with an eight-channel cardiac receive-only coil array on the same scanner used for the phantom experiments. Two sequences were compared: paired interleaved bSSFP (N=2), and grouped interleaved bSSFP with N=6. It can be noted that in the grouped interleaved acquisition the value of N was dictated by the desired temporal resolution. In particular, the maximum achievable temporal resolution was seen as decreasing as N increases, since the acquisition window per cardiac phase is at least 2N TRs.
Simulation results are shown for ΔfAB=4.0 Hz, which corresponds to a precession imbalance of 11.5°.
With continued reference to
a) shows the maximum magnitude deviation over the passband of width 2/(3 TR), relative to the steady-state signal for a non-interleaved bSSFP sequence, as shown by the inset in
The images of
Although the analysis presented herein for some exemplary embodiments applies only to imaging in the steady state, additional simulation were performed that indicate that transient imaging can also benefit from gradient waveform grouping. Grouping the waveforms in groups of six resulted in reduced oscillations compared to N=2.
Embodiments of the present disclosure have employed a relatively simple model, namely, that the effect of different waveforms is to alter the precession angle by an amount 2πΔfABTR. Strictly speaking, since residual eddy-current fields have in fact been shown to be time-varying, it could be more accurate to include a time-varying ΔfAB(t) in a model. However, since eddy-fields have also been shown to be short-lived (<1 msec), the steady-state distortions may be expected to depend only on the net precession angle
over the entire TR, and thus to be independent of the precise form of ΔfAB(t). In other words, the “effect” of each waveform only lasts for at most one TR, and does not carry over into subsequent TRs. This further suggests that the analysis presented herein may apply to all bSSFP sequences that exhibit unbalanced precession, regardless of the physical cause of the imbalance.
In order to achieve reasonable image quality for embodiments including cardiac velocity-mapping experiments (even for N=6), it was necessary to choose a rather large VENC value (4.6 m/s). It was observed empirically that reducing the VENC, thereby making the two interleaved gradient waveforms increasingly dissimilar, leads to increased ghosting artifacts and signal loss. This is likely due to increased eddy-current-induced precession imbalance as the gradients become more dissimilar. A VENC value as high as 4-5 m/s may be suitable for imaging stenotic jets, but may be too high to allow meaningful myocardial tissue velocity mapping. In order to perform interleaved bSSFP velocity mapping with more typical VENC values of 1.5 m/s or less, it may be necessary to apply additional methods for eddy-current artifact suppression.
With continued reference to
Accordingly, embodiments of the present disclosure can provide gradient waveform grouping as a simple and general approach to mitigating steady-state signal distortions in interleaved balanced SSFP. When waveforms are executed in groups of N, signal profile variations tend to improve as N increases. Embodiments of the present disclosure can provide grouping of waveforms which may also benefit applications that acquire images during the transient approach to steady state.
The components, steps, features, objects, benefits and advantages that have been discussed are merely illustrative. None of them, nor the discussions relating to them, are intended to limit the scope of protection in any way. Numerous other embodiments are also contemplated. These include embodiments that have fewer, additional, and/or different components, steps, features, objects, benefits and advantages. These also include embodiments in which the components and/or steps are arranged and/or ordered differently.
In reading the present disclosure, one skilled in the art will appreciate that embodiments of the present disclosure can be implemented in hardware, software, firmware, or any combinations of such, and over one or more networks. Moreover, embodiments of the present disclosure can be included in or carried by various signals, e.g., as transmitted over a wireless RF or IR communications link or downloaded from the Internet.
Unless otherwise stated, all measurements, values, ratings, positions, magnitudes, sizes, and other specifications that are set forth in this specification, including in the claims that follow, are approximate, not exact. They are intended to have a reasonable range that is consistent with the functions to which they relate and with what is customary in the art to which they pertain.
All articles, patents, patent applications, and other publications which have been cited in this disclosure are hereby incorporated herein by reference.
The phrase “means for” when used in a claim is intended to and should be interpreted to embrace the corresponding structures and materials that have been described and their equivalents. Similarly, the phrase “step for” when used in a claim embraces the corresponding acts that have been described and their equivalents. The absence of these phrases means that the claim is not intended to and should not be interpreted to be limited to any of the corresponding structures, materials, or acts or to their equivalents.
Nothing that has been stated or illustrated is intended or should be interpreted to cause a dedication of any component, step, feature, object, benefit, advantage, or equivalent to the public, regardless of whether it is recited in the claims.
The scope of protection is limited solely by the claims that now follow. That scope is intended and should be interpreted to be as broad as is consistent with the ordinary meaning of the language that is used in the claims when interpreted in light of this specification and the' prosecution history that follows and to encompass all structural and functional equivalents.
This application is based upon and claims priority to U.S. Provisional Patent Application No. 61/106,782, entitled “Method for Eddy-Current Artifact Reduction in Balanced Steady-State Free Precession Magnetic Resonance Imaging (MRI),” filed 20 Oct. 2008, attorney docket number 028080-0403, client reference 09-079, the entire contents of which are incorporated herein by reference.
This invention was made with Government support under contract grant numbers R01-HL074332 and R21-HL079987 awarded by the National Institutes of Health. The Government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
61106782 | Oct 2008 | US |