function?[ybeginxbegint2R]=SSDA（S?T）

%-------------------------------------------------------------------------%
%該函數(shù)用來查找原始圖像S中的子區(qū)塊T的位置，位置值取匹配區(qū)域的左上角坐標。
%函數(shù)使用互相關法實現(xiàn)。為了減小搜索區(qū)域，先將原圖像和子圖像抽取至較小范圍，獲取
%相應位置后，再在原圖像相應點附近詳細搜索。
%返回值xbegin、ybegin是相應的坐標值;t2是計算出來的匹配時間；R是累加次數(shù)矩陣。
%-------------------------------------------------------------------------%
%時間計算開始
t0=cputime;
hbar?=?waitbar（0‘Image?matching‘）;
tic;

%若原圖為彩色，轉換為灰度圖像（我用的是R2011b版本，isrgb函數(shù)不能使用）
if?size（S3）==3
????S=rgb2gray（S）;
end
%若模板為彩色，轉換為灰度圖像
if?size（T3）==3
????T=rgb2gray（T）;
end

%兩矩陣數(shù)值由整形轉換為雙精度型，便于后面的計算需要
S=double（S）;
T=double（T）;

%利用size函數(shù)得到兩矩陣的大小
[THeight?TWidth]?=?size（T）;
[imgH?imgW]?=?size（S）;

%求取模板圖片的灰度平均值TAverage
T1=0;
for?I=1:THeight
????for??J=1:TWidth
????????T1=T1+T（IJ）;
????end;
end;
TAverage=T1/（THeight*TWidth）;

%閾值Tk和累計次數(shù)矩陣R（Tk可調(diào)整大小，R初始設置為0）
R=zeros（imgH-THeight+1imgW-TWidth+1）;Tk=1;

%?從左上角開始運算
for?i=1:（imgH-THeight+1）
????for?j=1:（imgW-TWidth+1）
????????
%?求子圖片的平均值SijAverage
????????Sij=S（i:（i+THeight-1）j:（j+TWidth-1））;
????????SijAverage=0;
????????for?m=1:THeight
????????????for?n=1:TWidth
????????????????SijAverage=（SijAverage+Sij（mn））;
????????????end;
????????end;
????????SijAverage=SijAverage/（THeight*TWidth）;
????????
%?隨機選點，求誤差e和累積誤差Sum
????????Sum=0;
????????while?（（Sum????????????X=round（rand*（THeight-1）+1）;
????????????Y=round（rand*（TWidth-1）+1）;
????????????e=abs（Sij（XY）-SijAverage-（T（XY）-TAverage））;
????????????Sum=Sum+e;
????????????R（ij）=R（ij）+1;
????????end;
????end;
end;

%找出累計次數(shù)最大的點以及顯示目標點的坐標，并且保存
Max=max（max（R））?;
[rowcol]=find（R==Max）;
ybegin?=?col;
xbegin?=?row;

%計算結束，求取總的計算時間
t2=toc;
close（hbar）;
t1=cputime-t0

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????340746??2012-09-17?17:28??基于圖像灰度的SSDA匹配算法.pdf
?????文件??????294083??2012-09-18?09:14??SSDA圖像匹配算法的研究及實現(xiàn).pdf
?????文件????????3942??2012-06-04?21:00??Mark0.bmp
?????文件????????2094??2012-09-19?11:14??SSDA.m
?????文件?????????543??2012-09-19?11:14??test.m
?????文件??????954910??2006-12-24?17:32??MAIN0.bmp