TY - JOUR
T1 - Magnetic Resonance RF Pulse Design by Optimal Control With Physical Constraints
AU - Rund, Armin
AU - Aigner, Christoph Stefan
AU - Kunisch, Karl
AU - Stollberger, Rudolf
PY - 2018/2
Y1 - 2018/2
N2 - Optimal control approaches have proved useful in designing RF pulses for large tip-angle applications. A typical challenge for optimal control design is the inclusion of constraints resulting from physiological or technical limitations that assure the realizability of the optimized pulses. In this paper, we show how to treat such inequality constraints, in particular, amplitude constraints on the B1 field, the slice-selective gradient, and its slew rate, as well as constraints on the slice profile accuracy. For the latter, a pointwise profile error and additional phase constraints are prescribed. Here, a penalization method is introduced that corresponds to a higher order tracking instead of the common quadratic tracking. The order is driven to infinity in the course of the optimization. We jointly optimize for the RF and slice-selective gradient waveform. The amplitude constraints on these control variables are treated efficiently by semismooth Newton or quasi-Newton methods. The method is flexible, adapting to many optimization goals. As an application, we reduce the power of refocusing pulses, which is important for spin echo-based applications with a short echo spacing. Here, the optimization method is tested in numerical experiments for reducing the pulse power of simultaneous multislice refocusing pulses. The results are validated by phantom and in-vivo experiments.
AB - Optimal control approaches have proved useful in designing RF pulses for large tip-angle applications. A typical challenge for optimal control design is the inclusion of constraints resulting from physiological or technical limitations that assure the realizability of the optimized pulses. In this paper, we show how to treat such inequality constraints, in particular, amplitude constraints on the B1 field, the slice-selective gradient, and its slew rate, as well as constraints on the slice profile accuracy. For the latter, a pointwise profile error and additional phase constraints are prescribed. Here, a penalization method is introduced that corresponds to a higher order tracking instead of the common quadratic tracking. The order is driven to infinity in the course of the optimization. We jointly optimize for the RF and slice-selective gradient waveform. The amplitude constraints on these control variables are treated efficiently by semismooth Newton or quasi-Newton methods. The method is flexible, adapting to many optimization goals. As an application, we reduce the power of refocusing pulses, which is important for spin echo-based applications with a short echo spacing. Here, the optimization method is tested in numerical experiments for reducing the pulse power of simultaneous multislice refocusing pulses. The results are validated by phantom and in-vivo experiments.
KW - Algorithms
KW - Brain/diagnostic imaging
KW - Humans
KW - Image Processing, Computer-Assisted/methods
KW - Magnetic Resonance Imaging/methods
KW - Phantoms, Imaging
KW - Radio Waves
KW - MRI
KW - Optimal Control
KW - RF Pulse
KW - Simultaneous multiple slice
U2 - 10.1109/TMI.2017.2758391
DO - 10.1109/TMI.2017.2758391
M3 - Article
C2 - 28981407
SN - 1558-254X
VL - 37
SP - 461
EP - 472
JO - IEEE Transactions on Medical Imaging
JF - IEEE Transactions on Medical Imaging
IS - 2
ER -