A comprehensive numerical analysis of background phase correction with V-SHARP.
Sophisticated harmonic artifact reduction for phase data (SHARP) is a method to remove background field contributions in MRI phase images, which is an essential processing step for quantitative susceptibility mapping (QSM). To perform SHARP, a spherical kernel radius and a regularization parameter need to be defined. In this study, we carried out an extensive analysis of the effect of these two parameters on the corrected phase images and on the reconstructed susceptibility maps. As a result of the dependence of the parameters on acquisition and processing characteristics, we propose a new SHARP scheme with generalized parameters. The new SHARP scheme uses a high-pass filtering approach to define the regularization parameter. We employed the variable-kernel SHARP (V-SHARP) approach, using different maximum radii (Rm) between 1 and 15 mm and varying regularization parameters (f) in a numerical brain model. The local root-mean-square error (RMSE) between the ground-truth, background-corrected field map and the results from SHARP decreased towards the center of the brain. RMSE of susceptibility maps calculated with a spatial domain algorithm was smallest for Rmbetween 6 and 10 mm and f between 0 and 0.01 mm-1, and for maps calculated with a Fourier domain algorithm for Rmbetween 10 and 15 mm and f between 0 and 0.0091 mm-1. We demonstrated and confirmed the new parameter scheme in vivo. The novel regularization scheme allows the use of the same regularization parameter irrespective of other imaging parameters, such as image resolution. Copyright © 2016 John Wiley & Sons, Ltd.
Institute of Diagnostic and Interventional Radiology, University Hospital Zurich and University of Zurich, Zurich, Switzerland.|Institute for Biomedical Engineering, University of Zurich and ETH Zürich, Zurich, Switzerland.
Pinar Senay
PS
Özbay
Medical Physics Group, Institute of Diagnostic and Interventional Radiology, Jena University Hospital - Friedrich Schiller University Jena, Jena, Germany.
Andreas
A
Deistung
Medical Physics Group, Institute of Diagnostic and Interventional Radiology, Jena University Hospital - Friedrich Schiller University Jena, Jena, Germany.|Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, The State University of New York, Buffalo, NY, USA.
Xiang
X
Feng
Institute of Diagnostic and Interventional Radiology, University Hospital Zurich and University of Zurich, Zurich, Switzerland.
Daniel
D
Nanz
Medical Physics Group, Institute of Diagnostic and Interventional Radiology, Jena University Hospital - Friedrich Schiller University Jena, Jena, Germany.|Michael Stifel Center for Data-driven and Simulation Science Jena, Friedrich Schiller University Jena, Jena, Germany.
Jürgen Rainer
JR
Reichenbach
Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, The State University of New York, Buffalo, NY, USA.|MRI Clinical and Translational Research Center, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, The State University of New York, Buffalo, NY, USA.
Ferdinand
F
Schweser
D000465:Algorithms; D016477:Artifacts; D001921:Brain / Q000033:anatomy & histology*; D006801:Humans; D007089:Image Enhancement / Q000379:methods*; D007090:Image Interpretation, Computer-Assisted / Q000379:methods*; D008279:Magnetic Resonance Imaging / Q000379:methods*; D009716:Numerical Analysis, Computer-Assisted; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity
background field removal; phase; quantitative susceptibility mapping (QSM); sophisticated harmonic artifact reduction for phase data (SHARP)