Premium
Monte Carlo SURE‐based parameter selection for parallel magnetic resonance imaging reconstruction
Author(s) -
Weller Daniel S.,
Ramani Sathish,
Nielsen JonFredrik,
Fessler Jeffrey A.
Publication year - 2014
Publication title -
magnetic resonance in medicine
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.696
H-Index - 225
eISSN - 1522-2594
pISSN - 0740-3194
DOI - 10.1002/mrm.24840
Subject(s) - undersampling , monte carlo method , regularization (linguistics) , iterative reconstruction , algorithm , mean squared error , computer science , mathematical optimization , mathematics , statistics , artificial intelligence
Purpose Regularizing parallel magnetic resonance imaging (MRI) reconstruction significantly improves image quality but requires tuning parameter selection. We propose a Monte Carlo method for automatic parameter selection based on Stein's unbiased risk estimate that minimizes the multichannel k‐space mean squared error (MSE). We automatically tune parameters for image reconstruction methods that preserve the undersampled acquired data, which cannot be accomplished using existing techniques. Theory We derive a weighted MSE criterion appropriate for data‐preserving regularized parallel imaging reconstruction and the corresponding weighted Stein's unbiased risk estimate. We describe a Monte Carlo approximation of the weighted Stein's unbiased risk estimate that uses two evaluations of the reconstruction method per candidate parameter value. Methods We reconstruct images using the denoising sparse images from GRAPPA using the nullspace method (DESIGN) and L 1 iterative self‐consistent parallel imaging (L 1 ‐SPIRiT). We validate Monte Carlo Stein's unbiased risk estimate against the weighted MSE. We select the regularization parameter using these methods for various noise levels and undersampling factors and compare the results to those using MSE‐optimal parameters. Results Our method selects nearly MSE‐optimal regularization parameters for both DESIGN and L 1 ‐SPIRiT over a range of noise levels and undersampling factors. Conclusion The proposed method automatically provides nearly MSE‐optimal choices of regularization parameters for data‐preserving nonlinear parallel MRI reconstruction methods. Magn Reson Med 71:1760–1770, 2014. © 2013 Wiley Periodicals, Inc.