This paper presents a new magnetic resonance imaging (MRI) phase correction method. The linear phase correction method using autocorrelation proposed by Ahn and Cho (AC method) is extended to handle nonlinear terms, which are often important for polynomial expansion of phase variation in MRI. The polynomial coefficients are statistically determined from a cascade series of n-pixel-shift rotational differential fields (RDFs). The n-pixel-shift RDF represents local vector rotations of a complex field relative to itself after being shifted by n pixels. We have found that increasing the shift enhances the signal significantly and extends the AC method to handle higher order nonlinear phase error terms. The n-pixel-shift RDF can also be applied to improve other methods such as the weighted least squares phase unwrapping method proposed by Liang. The feasibility of the method has been demonstrated with two-dimensional (2-D) in vivo inversion-recovery MRI data.