The present invention provides a method of handling rapid phase aliasing in magnetic resonance images arising from local magnetic susceptibility differences. The methods of the present invention can be used to estimate the field effects within an object arising from the interfaces of regions having differences in magnetic susceptibilities, and to subtract out the resulting phase from the original or source phase data prior to any further phase processing. The methods of the present invention also include a process of accurately determining the susceptibility values of multiple voxel regions based on the geometry of such regions.