Minimizing dead-periods in flow-encoded or -compensated pulse sequences while imaging in oblique planes

A gradient optimization toolbox for general purpose time-optimal MRI gradient waveform design

Purpose

To introduce and demonstrate a software library for time‐optimal gradient waveform optimization with a wide range of applications. The software enables direct on‐the‐fly gradient waveform design on the scanner hardware for multiple vendors.

Methods

The open‐source gradient optimization (GrOpt) toolbox was implemented in C with both Matlab and Python wrappers. The toolbox enables gradient waveforms to be generated based on a set of constraints that define the features and encodings for a given acquisition. The GrOpt optimization routine is based on the alternating direction method of multipliers (ADMM). Additional constraints enable error corrections to be added, or patient comfort and safety to be adressed. A range of applications and compute speed metrics are analyzed. Finally, the method is implemented and tested on scanners from different vendors.

Results

Time‐optimal gradient waveforms for different pulse sequences and the constraints that define them are shown. Additionally, the ability to add, arbitrary motion (gradient moment) compensation or limit peripheral nerve stimulation is demonstrated. There exists a trade‐off between computation time and gradient raster time, but it was observed that acceptable gradient waveforms could be generated in 1‐40 ms. Gradient waveforms generated and run on the different scanners were functionally equivalent, and the images were comparable.

Conclusions

GrOpt is an open source toolbox that enables on‐the‐fly optimization of gradient waveform design, subject to a set of defined constraints. GrOpt was presented for a range of imaging applications, analyzed in terms of computational complexity, and implemented to run on the scanner for a multi‐vendor demonstration.

Optimization methods for magnetic resonance imaging gradient waveform design

The development and implementation of novel MRI pulse sequences remains challenging and laborious. Gradient waveforms are typically designed using a combination of analytical and ad hoc methods to construct each gradient waveform axis independently. This strategy makes coding the pulse sequence complicated, in addition to being time inefficient. As a consequence, nearly all commercial MRI pulse sequences fail to maximize use of the available gradient hardware or efficiently mitigate physiological effects. This results in expensive MRI systems that underperform relative to their inherent hardware capabilities. To address this problem, a software solution is proposed that incorporates numerical optimization methods into MRI pulse sequence programming. Examples are shown for rotational variant vs. invariant waveform designs, acceleration nulled velocity encoding gradients, and mitigation of peripheral nerve stimulation for diffusion encoding. The application of optimization methods to MRI pulse sequence design incorporates gradient hardware limits and the prescribed MRI protocol parameters (e.g. field-of-view, resolution, gradient moments, and/or b-value) to simultaneously construct time-optimal gradient waveforms. In many cases, the resulting constrained gradient waveform design problem is convex and can be solved on-the-fly on the MRI scanner. The result is a set of multi-axis time-optimal gradient waveforms that satisfy the design constraints, thereby increasing SNR-efficiency. These optimization methods can also be used to mitigate imaging artifacts (e.g. eddy currents) or account for peripheral nerve stimulation. The result of the optimization method is to enable easier pulse sequence gradient waveform design and permit on-the-fly implementation for a range of MRI pulse sequences.

Efficient implementation of hardware-optimized gradient sequences for real-time imaging

This work improves the performance of interactive real-time imaging with balanced steady-state free precession. The method employs hardware-optimized gradient pulses, together with a novel phase-encoding strategy that simplifies the design and implementation of the optimized gradient waveforms. In particular, the waveforms for intermediate phase-encode steps are obtained by simple linear combination, rather than separate optimized waveform calculations. Gradient waveforms are redesigned in real time as the scan plane is manipulated, and the resulting sequence operates at the specified limits of the MRI gradient subsystem for each new scan-plane orientation. The implementation provides 14-25% improvement in the sequence pulse repetition time over the vendor-supplied interactive real-time imaging sequence for similar scan parameters on our MRI scanner.

Time-optimal multidimensional gradient waveform design for rapid imaging

Magnetic resonance imaging (MRI) is limited in many cases by long scan times and low spatial resolution. Recent advances in gradient systems hardware allow very rapid imaging sequences, such as steady-state free precession (SSFP), which has repetition times (TRs) of 2-5 ms. The design of these rapid sequences demands time-optimal preparatory gradient waveforms to provide maximum readout duty-cycle, and preserve spatial resolution and SNR while keeping TRs low. Time-optimal gradient waveforms can be synthesized analytically for certain simple cases. However, certain problems, such as time-optimal 2D and 3D gradient design with moment constraints, either may not have a solution or must be solved numerically. We show that time-optimal gradient design is a convex-optimization problem, for which very efficient solution methods exist. These methods can be applied to solve gradient design problems for oblique gradient design, spiral imaging, and flow-encoding using either a constant slew rate or the more exact voltage-limited gradient models. Ultimately, these methods provide a time-optimal solution to many 2D and 3D gradient design problems in a sufficiently short time for interactive imaging.

In vivo measurement of pulsewave velocity in small vessels using intravascular MR

A one-dimensional intravascular MR (IVMR) technique for the measurement of pulsewave velocity in a single cardiac cycle is presented. The technique was used to measure pulsewave velocity in vivo in the intact rabbit model, where its sensitivity to different hemodynamic states was demonstrated using a pharmacological intervention with phenylephrine and nitroprusside. IVMR measurements of pulsewave velocity were found to increase with mean arterial pressure, as expected. Further, IVMR-based pulsewave velocity estimates were in agreement with those measured by pressure catheters and direct distensibility measurement. Because of their rapidity and highly localized nature, these measurements of vessel elasticity may complement the high-resolution vascular imaging information gained in an IVMR examination. This could allow assessment of atherosclerotic plaques and facilitate immediate treatment decisions. Magn Reson Med 45:53-60, 2001.