Emilie Mussard1,2,3, Tom Hilbert1,2,3, Reto Meuli2, Jean-Philippe Thiran2,3, and Tobias Kober1,2,3
1Advanced Clinical Imaging Technology (HC CMEA SUI DI BM PI), Siemens Healthcare AG, Lausanne, Switzerland, 2Department of Radiology, University Hospital (CHUV), Lausanne, Switzerland, 3LTS5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
Synopsis
MP2RAGE is a T1 imaging
method providing greatly reduced B1 bias as well as less T2* and
PD-contributions. It requires, however, long acquisition time (standard
protocol with GRAPPAx3: ~8min) which hampers its clinical application. This
work proposes to use sparse iterative reconstruction techniques to shorten
MP2RAGE acquisition times. Resulting images are benchmarked using contrast
assessment, changes in obtained T1 values as well as evaluating the effect of
undersampling on an automated brain morphometry algorithm. Acceptable penalty
in image quality and morphometric outcome was achieved with up to 5-fold
acceleration, yielding a measurement time of 3.8min compared to fully sampled 20min.Introduction:
MP2RAGE
1 is a T1 imaging method that
greatly reduces the B1 bias field as well as T2* and PD contrast compared to
standard MPRAGE acquisitions and has the additional advantage to generate T1
maps from the obtained MP2RAGE contrast. To this end, two FLASH images are
sampled after inversion, resulting in a prolonged TR and thus long total
acquisition times (~8min using GRAPPA x3
2).
For clinical use, examinations of this duration are difficult to conduct. We
thus propose to apply sparse iterative reconstruction3 on MP2RAGE
images to reduce the required acquisition time. Results are benchmarked
calculating contrast figures for the different acceleration factors as well as
assessing the undersampling effects on an automated brain morphometry
algorithm.
Materials & Methods:
After obtaining written
consent, a fully
sampled MP2RAGE (TR 5s, TI1/TI2 0.7s/2.5s, flip angles 4
and 5 degrees, resolution 1mm isotropic, acq. matrix 256x240x176, TA=20mn) of one healthy volunteer was acquired at 3T (MAGNETOM
Skyra, Siemens Healthcare, Germany) using a 20-channel head/neck coil.
Artificial undersampling was performed using a variable-density Cartesian spiral
phyllotaxis pattern4 with different acceleration factors from 1 to
16. The reconstruction of the images from undersampled data was computed by
iteratively minimizing the following cost function enforcing both
consistency with acquired data and sparsity in the wavelet domain: $$\min_{X}\frac{1}{2}\sum_{i=1,2}\parallel PF\left\{S_{c}X_{i}\right\}-Y_{i}\parallel_2^2+\lambda\mid\Psi X_{i}\mid_{1}$$
with P being the
sampling mask, F the discrete Fourier transform, Sc complex sensitivities computed with ESPIRiT5, Y the undersampled
k-space, λ a regularization parameter and Ψ the wavelet-transform.
Both the MP2RAGE uniform contrast and the T1
map of the fully sampled dataset were reconstructed and served as ground truth
for comparing the obtained undersampling results. Furthermore, the Morphobox
prototype6 was applied on the fully sampled uniform image to obtain
six masks of structures of interest for further analysis, namely: thalamus,
caudate, putamen, hippocampus, global white matter and global grey matter.
The artificially
undersampled datasets were reconstructed using the procedure shown above,
obtaining uniform contrasts from which T1 maps were subsequently calculated
following [1]. After a first qualitative evaluation, contrast ratios (CR), contrast-to-noise ratios
(CNR), root-mean-square difference
(RMSD) as well as T1 map differences were quantitatively assessed in the six
brain structures defined above. CR and CNR were computed using the definition
given by Okubo’s comparison of MPRAGE and MP2RAGE7.
Results and Discussion:
Figure 1 shows reconstructed slices, relative
difference to the fully sampled conventional reconstruction and RMSD with
increasing acceleration factor R (R = 2.19, 5.23, 7.90). It can be seen that
with higher R, edges smooth out and that anatomical information is lost in some
structures, e.g. in the caudate or putamen. Similarly, the RMSD is increasing
with rising R. Changes in volume estimates of the structures of interest over
the different Rs are shown in Figure 2. To note, white-matter estimates
considerably drop with R>4. CR and CNR figures, however, remain stable as
can be seen in Table 1. CNR shows a slight increase for R>5 which is
probably due to the iterative denoising. T1 maps show that average T1 values
remain stable in all structures of interest; in some structures, however,
standard deviations increase to up to 10%.
Conclusion:
The application of sparse iterative
reconstruction on undersampled MP2RAGE acquisitions allows obtaining images
with only minor quality and CNR degradation for a not too high R. Increased
blurring might however impede visual reading with higher acceleration factors.
Our preliminary data suggests that acceleration factors up to 5.2,
corresponding to an acquisition time of 3.8 min, are feasible with acceptable
quality penalty. Notably, the iterative reconstruction proposed here has not
yet used the redundancy in the two inversion contrasts as proposed by Berkin et
al.
8; exploiting these, further improvements in image reconstruction
quality may be feasible.
Acknowledgements
No acknowledgement found.References
[1] Marques, José P., et al. "MP2RAGE, a self bias-field corrected
sequence for improved segmentation and T1-mapping at high field."
NeuroImage 49 (2010): 1271-1281.
[2]
Griswold, Mark A., et al. "Generalized Autocalibrating Partially Parallel
Acquisitions (GRAPPA)" Magnetic Resonance in Medicine 47.6 (2002):
1202-1210.
[3] Lustig,
Michael, et al. "Sparse MRI: The Application of Compressed Sensing for
Rapid MR Imaging." Magnetic Resonance in Medicine 58 (2007): 1182-1195.
[4] Forman,
Christoph, et al. "High-resolution 3D whole-heart coronary MRA: a study on
the combination of data acquisition in multiple breath-holds and 1D residual
respiratory motion compensation." Magn Reson Mater Phy 27 (2014):435-443.
[5] Uecker,
Martin, et al. "ESPIRiT—An Eigenvalue Approach to Autocalibrating Parallel
MRI: Where SENSE Meets GRAPPA." Magnetic Resonance in Medicine 71 (2014):
990-1001.
[6]
Schmitter, Daniel, et al. "An evaluation of volume-based morphometry for
prediction of mild cognitive impairment and Alzheimer's disease."
NeuroImage: Clinical 7 (2015): 7-17.
[7] Okubo,
Gosuke, et al. "MP2RAGE for Deep Gray Matter Measurement of the Brain: A
Comparative Study with MPRAGE." Journal of Magnetic Resonance Imaging
(2015).
[8] Bilgic,
Berkin, et al. "Multi-contrast reconstruction with Bayesian compressed
sensing." Magnetic Resonance in Medicine 66 (2011): 1601-1615.