相位相关matlab实现

Here is working code for an example that uses a single image and cuts out a 256x256 pixel portion of the image from one location 

即在同一副图像中，先切一个256*256的图片，然后位移，再切一个同样大小256*256的图像

and then cuts out the same size portion from a location that has been translated by 11 pixels in the x direction and 7 pixels in the y direction. 

This shows that the basic technique does work when you have the code formatting correct.

image1 = rgb2gray(imread('frame429.jpg')); % read in image and convert from RGB ro grayscale
rect1 = [1 1 255 255]; % define first cutout region [xmin ymin width height]
rect2 = [12 8 255 255]; % define second translated cutout region [xmin ymin width height]
sub_image1 = imcrop(image1,rect1); % 1st cut out portion of original image 
sub_image2 = imcrop(image1,rect2); % 2nd cut out portion of original image 
filter1=blackman(length(sub_image1),'symmetric');
filtimg1=(filter1*filter1');
Image1=double(sub_image1).*filtimg1;
Image2=double(sub_image2).*filtimg1
FFT1 = fft2(Image1);
FFT2 = conj(fft2(Image2));
FFTR = FFT1.*FFT2;
magFFTR = abs(FFTR);
FFTRN = (FFTR./magFFTR);
result = ifft2(double(FFTRN)); 
figure;
colormap('gray'); % choose colormap for plotting
imagesc(result);
figure;
colormap(jet);

mesh(result);

The code can easily be changed to use two different images of the same dimensions:

image1 = rgb2gray(imread('castle_color1.tif')); % read image 1 (convert to grayscale if RGB)
image2 = rgb2gray(imread('castle_color5.tif')); % read image 2 must have same dimensions as image 1
subplot(1,2,1);imshow(image1);
subplot(1,2,2);imshow(image2);
filter1=blackman(size(image1,1),'symmetric'); % setup window in one direction - I chose Blackman window - could use Hamming or others
filter2=blackman(size(image1,2),'symmetric'); % setup window in second direction
filtimg1=(filter1*filter2'); % make 2D windowing function that is sized to match image
Image1=double(image1).*filtimg1; % apply windowing function to image
Image2=double(image2).*filtimg1;
FFT1 = fft2(Image1); % 2d FFT
FFT2 = conj(fft2(Image2)); 
FFTR = FFT1.*FFT2; % G1*G2复共轭
magFFTR = abs(FFTR); %sqrt(real^2 + imag^2) 幅值
FFTRN = (FFTR./magFFTR);
result = ifft2(double(FFTRN)); 
figure;
colormap('gray'); % choose colormap for plotting
imagesc(result);
figure;
colormap(jet);
mesh(result);