The disclosed method enables efficient registration and correction of random vertical and horizontal offsets of a plurality of images of a feature to enable combination thereof with maximum noise suppression. The method includes obtaining image cross correlations and estimating correlation peaks. If a correlation peak cannot be uniquely determined for a given correlation, all candidates are retained. An improved robust weighted regression is applied thereto to obtain estimates of the image shifts. Dimensional symmetry is exploited to remove redundancy from the weighted normal equation, and computational requirements are further reduced by analytically determining the coefficients matrix and the inhomogeneous matrix, thereby circumventing conventional computations. Solution of the resulting consistent, full rank, linear system requires only matrix inversion. For N images, the present method thereby increases computation speed and reduces storage by approximately N2.