Magnetic resonance imaging (MRI) being one of the most promising neuro-imaging modalities in medical diagnosis suffers from spatial inaccuracy when it comes to high precision processes such as neuro-stereotactic surgery. Spatial distortions of several millimeters are possible. The distortion originates from the fact that highly controllable magnetic fields are used to encode the positions of the spins. When the fields are perturbed by unknown interferences such as the para- and/or dia-magnetism of the patient or inhomogeneity in the main field the spatial encoding becomes inaccurate. However the distortions are not without helpful characteristics. For instance using an SE pulse sequence one may observe that 1) the distortion vectors are collinear 2) they are along an axis which is determined by the gradients and 3) the integral of image intensity between homologous coaxial points is unaffected by the inhomogeneity. With such an observation we propose a method that uses a pair of images from the same object that are acquired using a set of gradients having a prescribed relationship. While these two images are both distorted a set of independent first order ordinary differential equations have been derived to relate the two images. Based upon this relation we show that a distortion-compensated image can be derived. We also show some preliminary phantom study results as well as computer simulation results. The objective of this study is to arrive at spatially accurate MRI