speckle.rar

close all; clc; clear; %Iu = imread('U0V0.bmp彩色TO灰度.jpg'); %Id = imread('U10V10.BMP彩色TO灰度.jpg'); Iu = imread('平移前的图2.png'); Id = imread('平移后的图2.png'); [m,n]=size(Iu); %取H为以Iu中心为中心的41×41矩阵 H=Iu(m/2-20:m/2+20,n/2-20:n/2+20); tt=0; for i =21:m-21 for j =21:n-21 J=Id(i-20:i+20,j-20:j+20); %求矩阵H和J的相似度系数 R= corr2(H,J); %如果相关系数大于0，则搜索窗口高h1=i,长l1=j if R>tt tt=R; h1=i; l1=j; end end end %像素级坐标X,Y %meshgrid生成和H一样大小的网格矩阵，矩阵X和Y的大小为41×41 [X,Y] = meshgrid(m/2-20:m/2+20,n/2-20:n/2+20); Pos=[X(:), Y(:)]; %生成41×41全零矩阵 E=zeros(41,41); [Count, Width] = size(Pos); D= zeros(Count,4); tic for k =1: Count C=E; max=0; A=Iu(Pos(k,1)-20:Pos(k,1)+20,Pos(k,2)-20:Pos(k,2)+20); for i = Pos(k,1)+h1-m/2-1 : Pos(k,1)+h1-m/2+1 for j = Pos(k,2)+l1-n/2-1 : Pos(k,2)+l1-n/2+1 B=Id(i-20:i+20,j-20:j+20); R= corr2(A,B); C(i,j) = R; if R>max max=R; h=i; l=j; end end end Id=double(Id); [a,xx,yy]=fit2(x,y,z,n,m,p,q); %[subx,suby]=subpixel(h,l,C); D(k,:)=[Pos(k,1),Pos(k,2),suby-Pos(k,2),subx-Pos(k,1)]; %disp([Pos(k,1),Pos(k,2),suby-Pos(k,2),subx-Pos(k,1)]); end toc