图像亚像素配准的乐队

应用背景在一个和第二部分的图像（二维刚性平移）寄存器由用户指定的像素%。而一零填充FFT计算%（快速傅里叶变换），此代码使用选择性采样的%矩阵乘以DFT（离散FT）极大地减少了计算时间和内存不牺牲精度。用这个程序所有的图像点都是用来%计算绕其峰一个很小的邻域采用互相关。这%的算法称为单步DFT算法。%曼努埃尔sicairos %归，塞缪尔·T·瑟曼，和杰姆斯R. Fienup，关键技术 ；代码接收参考FFT和移位的图像，和一个%（整数）采样的因素。代码期待与DC在FFT（1,1）所以不使用% fftshift。%%  ； ； ；输出= dftregistration（fft2（F），fft2（G），usfac）；%%的注册图像的一个像素在1 / usfac。%%（1）输出的均方根误差（NRMSE）[ 1 ]之间的F%%%的输出（2）是全球相位差的图像（应该是）%如果图像是真实值和非负的话。%%的输出（3）和输出（4）分别为。%%  ； ； ；[格雷戈] =输出dftregistration（fft2（F），fft2（G），usfac）；%%格雷戈是一个可选的输出，它返回的傅立叶变换的注册版本的克，%的全球相位差[输出（2）]也被补偿。

efficient_subpixel_registration_mine.m