1. A method of nuclear magnetic resonance imaging of an object, the method comprising: acquiring phase data, φ(r), from a gradient echo magnetic resonance scan of the object for a first matrix of N x by N y by N z data points, where N x , N y , and N z are each a positive integer; generating φ(r) new for a larger matrix of aN x by bN y by cN z data points, where aN, bN, and cN are each a positive integer, using φ(r) and by assigning a constant to data points that are not common to both the first matrix and the larger matrix; Fourier transforming φ(r) new and generating k-space data (corresponding to phase), herein referred to as φ(k), therefrom; multiplying φ(k) by a susceptibility mapping filter function and inverse Fourier transforming the product in a manner generating a susceptibility map estimate for the larger matrix; and generating a digital image file of the object using the susceptibility map estimate.