function?[imaxjmaxxrou]=calculate（arbrmrnrabzmnleftright）
%計算
rou=zeros（mn）;
g1=mean2（left（a-（z-1）/2:a+（z-1）/2b-（z-1）/2:b+（z-1）/2））;
a=double（left（a-（z-1）/2:a+（z-1）/2b-（z-1）/2:b+（z-1）/2））-g1;
r1=mean2（a.*a）;
for?i=ar-（mr-1）/2:ar+（mr-1）/2
????for?j=br-（nr-1）/2:br+（nr-1）/2
????????g2=mean2（right（i-（z-1）/2:i+（z-1）/2j-（z-1）/2:j+（z-1）/2））;
????????b=double（right（i-（z-1）/2:i+（z-1）/2j-（z-1）/2:j+（z-1）/2））-g2;
????????r2=mean2（b.*b）;
????????r=mean2（a.*b）;
????????rou（ij）=r/sqrt（r1*r2）;
????end
end
x=max（max（rou））;
[imaxjmax]=find（rou==max（max（rou）））;

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

?????文件????????586??2012-01-12?10:30??影像匹配最終成果\calculate.m

?????文件???????1503??2012-01-12?10:37??影像匹配最終成果\practice3.m

?????文件????3367624??2012-01-10?08:33??影像匹配最終成果\原始影像右.tif

?????文件????3416416??2012-01-10?08:32??影像匹配最終成果\原始影像左.tif

?????文件?????????39??2012-01-12?09:54??影像匹配最終成果\右像特征點.txt

?????文件?????????40??2012-01-12?09:49??影像匹配最終成果\起始數據.txt

?????目錄??????????0??2012-01-12?16:31??影像匹配最終成果

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

??????????????6786208????????????????????7