A fast and high precision imaging registration method is presented, which can be used to improve the speed and accuracy of spectral recovery in a hyperspectral imaging process. In this method, two local upsampling phase correlation methods are adopted. One is to estimate the rotation and scaling parameters in a log-polar space, and the other is to estimate the translation parameters in a Cartesian space. A gradient preprocessing process is introduced to make the method more robust. The principle of Fourier-Mellin transform is introduced and this transform can decouple the rotation and scaling parameters. Based on the traditional phase correlation method, the local upsampling phase correlation method is introduced. The necessity of the gradient preprocessing process is verified. The simulation results show that the proposed method can be used to achieve high precision image registration.