close?all;
clc;
clear;
Iu?=?imread（‘11162.mli.bmp‘）;
Id?=?imread（‘01272.mli.bmp‘）;

[mn]=size（Iu）;
?D=?zeros（40004）;
?k=1;
tic
for?m1=550:10:700;
????for?n1=900:10:1050;
????????H=Iu（m1-40:m1+40n1-40:n1+40）;
????????
????????tt=0;

for?i?=m1-41:m1+41
????
????for?j?=n1-41:n1+41
????????
????????J=Id（i-40:i+40j-40:j+40）;
????????
????????R=?corr2（HJ）;%求相關系數(shù)?
????????
????????if?R>tt
????????????tt=R;
????????????h1=i;
????????????l1=j;
????????end
????end
end??????????????????%到此步定出整像素的位移?

???E=zeros（8181）;
???????
??????C=E;
????????
??????max=0;
????????

?for?i?=?h1-1?:?h1+1
????
????for?j?=?l1-1?:?l1+1
????????
????????B=Id（i-40:i+40j-40:j+40）;
????????
????????R=?corr2（HB）;
????????
????????C（ij）?=?R;
????????
????????if?R>max

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1090??2014-05-13?21:27??曲面擬合亞像素位移\offset.m

?????文件???????1119??2010-04-25?20:43??曲面擬合亞像素位移\subpixel.m

?????目錄??????????0??2014-05-20?09:37??曲面擬合亞像素位移

-----------?---------??----------?-----??----

?????????????????2209????????????????????3