close?all;
clc;
clear;
%Iu?=?imread（‘U0V0.bmp彩色TO灰度.jpg‘）;
%Id?=?imread（‘U10V10.BMP彩色TO灰度.jpg‘）;
Iu?=?imread（‘平移前的圖2.png‘）;
Id?=?imread（‘平移后的圖2.png‘）;

[mn]=size（Iu）;
%取H為以Iu中心為中心的41×41矩陣
H=Iu（m/2-20:m/2+20n/2-20:n/2+20）;

tt=0;

for?i?=21:m-21
????
????for?j?=21:n-21
????????
????????J=Id（i-20:i+20j-20:j+20）;
???????%求矩陣H和J的相似度系數?
????????R=?corr2（HJ）;
???????%如果相關系數大于0，則搜索窗口高h1=i長l1=j?
????????if?R>tt
????????????tt=R;
????????????h1=i;
????????????l1=j;
????????end
????end
end

%像素級坐標XY
%meshgrid生成和H一樣大小的網格矩陣，矩陣X和Y的大小為41×41
[XY]?=?meshgrid（m/2-20:m/2+20n/2-20:n/2+20）;

?Pos=[X（:）?Y（:）];
?%生成41×41全零矩陣
?E=zeros（4141）;
?
?[Count?Width]?=?size（Pos）;
?
D=?zeros（Count4）;

tic

for?k?=1:?Count
????????
????????C=E;
????????
????????max=0;
????????
????A=Iu（Pos（k1）-20:Pos（k1）+20Pos（k2）-20:Pos（k2）+20）;
????
?for?i?=?Pos（k1）+h1-m/2-1?:?Pos（k1）+h1-m/2+1
????
????for?j?=?Pos（k2）+l1-n/2-1?:?Pos（k2）+l1-n/2+1
????????
????????B=Id（i-20:i+20j-20:j+20）;
????????
????????R=?corr2（AB）;
????????
????????C（ij）?=?R;
????????
????????if?R>max
???????????max=R;
????????????h=i;
????????????l=j;
????????end
????????
????end
?end

Id=double（Id）;

[axxyy]=fit2（xyznmpq）
%[subxsuby]=subpixel（hlC）;

D（k:）=[Pos（k1）Pos（k2）suby-Pos（k2）subx-Pos（k1）];

%disp（[Pos（k1）Pos（k2）suby-Pos（k2）subx-Pos（k1）]）;

end


toc

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-11-09?14:24??speckle\
?????文件????????1548??2017-11-07?16:06??speckle\SANBAN.m
?????文件????????1460??2017-10-26?11:31??speckle\subpixel.m
?????文件???????15745??2010-04-29?09:38??speckle\x0y0.bmp彩色TO灰度.jpg
?????文件???????16530??2010-04-29?09:45??speckle\x0y1.BMP彩色TO灰度.jpg
?????文件???????16752??2010-04-29?09:38??speckle\x11y11.BMP彩色TO灰度.jpg
?????文件???????16751??2010-04-29?09:38??speckle\x17y5.BMP彩色TO灰度.jpg
?????文件???????16381??2010-04-29?09:45??speckle\x1y0.BMP彩色TO灰度.jpg
?????文件???????16603??2010-04-29?09:38??speckle\x3y19.BMP彩色TO灰度.jpg
?????文件???????16551??2017-10-31?14:34??speckle\平移前的模擬散斑圖2.png
?????文件???????27829??2017-10-31?14:34??speckle\平移后的模擬散斑圖2.png