A comprehensive numerical analysis of background phase correction with V‐SHARP. (3rd June 2016)

• Record Type: Journal Article
• Authors: Özbay, Pinar Senay; Deistung, Andreas; Feng, Xiang; Nanz, Daniel; Reichenbach, Jürgen Rainer; Schweser, Ferdinand
• Other Names: Liu Chunlei guestEditor.; Bowtell Richard guestEditor.; Schenck John guestEditor.
• Is Part Of: NMR in biomedicine. Volume 30:Number 4(2017:Apr.)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2016-06-03
• Subjects: background field removal -- phase -- quantitative susceptibility mapping (QSM) -- sophisticated harmonic artifact reduction for phase data (SHARP); Nuclear magnetic resonance -- Periodicals; Magnetic Resonance Spectroscopy -- Periodicals; 574
• DOI: 10.1002/nbm.3550
• Languages: English
• ISSNs: 0952-3480

Abstract

Abstract : 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 ( R m ) 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 R m between 6 and 10 mm and f between 0 and 0.01 mm −1, and for maps calculated with a Fourier domain algorithm for R m between 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.