Continuous three-dimensional imaging for magnetic resonance elastography
Described here are systems and methods for a robust magnetic resonance elastography (“MRE”) imaging platform for rapid dynamic 3D MRE imaging. The imaging platform includes an MRE pulse sequence and advanced image reconstruction framework that work synergistically in order to greatly expand the domains where MRE can be deployed successfully.
1 . A method for generating a mechanical property map using a magnetic resonance imaging (MRI) system, the method comprising:
(a) acquiring data from a subject using the MRI system while a vibratory stimulus is provided to the subject in order to induce motion within a tissue of the subject, wherein:
the data are acquired using a pulse sequence that samples k-space in a plurality of planes using an echo planar readout,
each plane is radially rotated about a phase encoding axis in successive repetition time periods of the pulse sequence, and
the pulse sequence includes motion encoding gradients that are generated before each echo planar readout in order to encode the motion within the tissue into an image phase of the data;
wherein the data are acquired with uniform undersampling in a through-plane dimension;
(b) reconstructing a plurality of images from the data while accounting for the uniform undersampling in the through-plane dimension, wherein the reconstructed images are indicative of the image phase into which the motion within the tissue has been encoded; and
(c) generating a mechanical property map from the reconstructed images.
2 . The method of claim 1 , wherein the plurality of images are reconstructed using an iterative image reconstruction technique that accounts for the uniform undersampling in the through-plane dimension.
3 . The method of claim 2 , wherein the iterative image reconstruction technique comprises a regularized iterative image reconstruction technique.
4 . The method of claim 3 , wherein the regularized iterative image reconstruction technique implements a Tikhonov regularizer.
5 . The method of claim 4 , wherein the regularized iterative image reconstruction technique is a Tikhonov-regularized SENSE reconstruction.
6 . The method of claim 2 , wherein the iterative image reconstruction technique implements a signal model that includes a forward operator that accounts for the uniform undersampling in the through-plane dimension.
7 . The method of claim 6 , wherein the forward operator also accounts for each plane being radially rotated about the phase encoding axis in successive repetition time periods of the pulse sequence by using in-plane radial Fourier encoding.
8 . The method of claim 1 , wherein second data are acquired with the MRI system while the vibratory stimulus is provided to the subject, wherein the vibratory stimulus is triggered to be temporally shifted relative to the motion encoding gradients in the pulse sequence by a time delay such that the second data are acquired during a different phase of the motion within the tissue.
9 . The method of claim 1 , wherein each plane is radially rotated about the phase encoding axis in successive repetition time periods of the pulse sequence according to a uniform distribution of rotation angles over an angular range.
10 . The method of claim 9 , wherein the angular range is 180 degrees.
11 . The method of claim 1 , wherein each plane is radially rotated about the phase encoding axis in successive repetition time periods of the pulse sequence according to a golden angle increment.
12 . The method of claim 1 , wherein an echo time of the pulse sequence is reduced using fractional motion encoding.
13 . The method of claim 2 , wherein field effects are incorporated into a signal model used in the iterative image reconstruction such that the field effects are prospectively accounted for during reconstruction.
14 . The method of claim 1 , wherein the mechanical property map comprises a stiffness map that depicts tissue stiffness in the tissue of the subject.