Quantitative Susceptibility Mapping

The following MATLAB implementations and validation on testing datasets can be downloaded here. A link for download will be shown after registration. Since the manuscript is under review, the zip file is password protected by the Manuscript ID to limit access to reviewers. Anyone interested in a preview can also send a request to medi-users@googlegroups.com .

The code has been tested on 64-bit Windows 7 + MATLAB 2009a.

The testing dataset includes a numerical simulation, an experimental phantom, and a human brain.

For background field removal:



Laplacian Boundary Value (LBV)

Zhou et al. NMR in biomed 2014;27(3):312-9

Sophisticated Harmonic Artifact Reduction for Phase data (SHARP)

Schweser et al. Neuroimage 2011;54(4):2789-807

Regularization enabled SHARP (RESHARP)

Sun et al. MRM 2014;71(3):1151-7

Projection onto Dipole Fields (PDF)

Liu et al. NMR in biomed 2011;24(9):1129-36


For dipole inversion:



Truncated singular value decomposition (TSVD)

Wharton et al. MRM 2010;63(5):1292-304

Truncated k-space division (TKD)

Shmueli et al. MRM 2009;62(6):1510-1522

Iterative SWIM (iSWIM)

Tang et al. MRM 2013;69(5):1396-407

Haacke et al. JMRI 2010;32(3):663-76

Morphology Enabled Dipole Inversion (MEDI)

Liu et al. MRM 2011;66(3):777-83

Liu et al. Neuroimage 2012;59(3):2560-8

TV with Split-Bregmann (TVSB)

Bilgic et al. MRM epub

Results only


Liu et al. MRM 2009;61(1):196-204

Homogeneity enabled incremental dipole inversion (HEIDI)

Schweser et al. Neuroimage 2012;62(3):2083-100

Compressed Sensing Compensated (CSC)

Wu et al. MRM 2012;67(1):137-47