To realize this goal, it is first shown that both the Laplace derivative of the z component of the static magnetic field in an object and the Laplacian of the corresponding phase distribution may be expected to be zero in regions of constant or linearly varying susceptibility and to be nonzero when there is an abrupt change in susceptibility, for instance, at a single point, a ridge, an interface, an edge or a boundary.
Next, a method is presented by which the Laplace derivative of a 3D phase map can be directly extracted from the complex data, without the need for phase unwrapping or subtraction of a reference image. The validity of this approach and of the theory behind it is subsequently demonstrated by simulations and phantom experiments with exactly known susceptibility distributions.
Finally, the potential of the Laplace derivative analysis is illustrated by simulations with a Shepp-Logan digital brain phantom and experiments with a gel phantom containing positive and negative focal susceptibility deviations.