An image alignment method includes computationally efficient methods of achieving high-accuracy local motion estimates by using phase correlation. The method also estimates motion reliability that allows a generic robust model fitting algorithm to produce more accurate results while operating much more efficiently. One of three methods are used to determine sub-pel motion estimation with improved accuracy. Each of the sub-pel motion estimation methods uses phase correlation, and are based on fitting computationally efficient 2-D quadratic surfaces to a phase correlation surface. A pre-filter is applied which shapes the phase correlation surface to enable appropriate fitting to the quadratic surface. Bias is also compensated for prior to applying a sub-pel motion estimation method. The method also estimates the reliability of the sub-pel motion estimates determined using phase correlation.