1. A method of nuclear magnetic resonance imaging an object, the method comprising: (a) collecting k-space data from a gradient echo magnetic resonance scan of the object; (b) generating initial phase and initial magnitude values using the k-space data; (c) generating a revised set of phase values using at least some of the initial magnitude values and at least some of the initial phase values; (d) generating local change in magnetic field values using at least some of the revised set of phase values; (e) generating modified k-space data values (corresponding to phase), using at least some of the local change in magnetic field values; (f) generating regularized modified k-space data values by applying a regularized inverse filter to the modified k-space data values; (g) generating first susceptibility values using the regularized modified k-space data values; (h) generating second susceptibility values by applying a threshold to the first susceptibility values and/or replacing the susceptibility values by a set of pre-determined values; (i) generating third k-space data values using the first susceptibility values; (j) generating fourth k-space data values using the second susceptibility values; (k) merging the third k-space data values and the fourth k-space data values in a manner generating a fifth k-space data values; (k) generating third susceptibility values using the fifth k-space data values; and (l) generating a digital image file of the object using the third susceptibility values.