matlab做出自相关图,在matlab中实现图像的自相关和互相关

%图像的自相关?

clear

I1=imread('lenna.bmp','bmp');%输入图像1(参考图像)

I1=I1(:,:,1);

figure(1)%显示输入图像1

colormap(gray(255));

image(I1)

axis off

FI1=fft2(I1);

max1=max(FI1);

max2=max(max1);

scale=1.0/max2;

FI1=FI1.*scale;

I2=imread('lenna.bmp','bmp');%输入图像2(待识别图像)

I2=I2(:,:,1);

figure(2)%显示输入图像2

colormap(gray(255));

image(I2)

axis off

FI2=fft2(I2);

max1=max(FI2);

max2=max(max1);

scale=1.0/max2;

FI2=FI2.*scale;

FPR=FI1.*conj(FI2);%计算相关性

PR=ifft2(FPR);

PR=fftshift(PR);

max1=max(PR);

max2=max(max1);

scale=1.0/max2;

PR=PR.*scale;

figure(3)%空间域相关显示

colormap(gray(255));

image(abs(256*PR));

axis off

%图像的互相关?

clear

I1=imread('lenna.bmp','bmp');%输入图像1(参考图像)

I1=I1(:,:,1);

figure(1)%显示输入图像1

% colormap(gray(255));

imagesc(I1)

axis off

FI1=fft2(I1);

max1=max(FI1);

max2=max(max1);

scale=1.0/max2;

FI1=FI1.*scale;

I2=imread('lenna.bmp','bmp');%输入图像2(待识别图像)

A=zeros(256,256);%图像矩阵的大小

A=I2(1:256,1:256);%图像矩阵的大小

A=A(:,:,1);

figure(2)%显示输入图像2

% colormap(gray(255));

imagesc(A)

axis off

FI2=fft2(A);

max1=max(FI2);

max2=max(max1);

scale=1.0/max2;

FI2=FI2.*scale;

FPR=FI1.*conj(FI2);%计算相关性

PR=ifft2(FPR);

PR=fftshift(PR);

max1=max(PR);

max2=max(max1);

scale=1.0/max2;

PR=PR.*scale;

figure(3)%空间域相关显示

% colormap(gray(255));

imagesc(abs(256*PR));

axis off