function?[hough_spacehough_circlepara]?=?hough_circle（BWstep_rstep_angler_minr_maxp）

%?%%%%%%%%%%%%%%%%%%%%%%%%%%
%?input
%?BW:二值圖像；
%?step_r:檢測的圓半徑步長
%?step_angle:角度步長，單位為弧度
%?r_min:最小圓半徑
%?r_max:最大圓半徑
%?p:閾值，0，1之間的數(shù)

%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?output
%?hough_space:參數(shù)空間，h（abr）表示圓心在（ab）半徑為r的圓上的點(diǎn)數(shù)
%?hough_circl:二值圖像，檢測到的圓
%?para:檢測到的圓的圓心、半徑

[mn]?=?size（BW）;
size_r?=?round（（r_max-r_min）/step_r）+1;
size_angle?=?round（2*pi/step_angle）;

hough_space?=?zeros（mnsize_r）;

[rowscols]?=?find（BW）;
ecount?=?size（rows）;

%?Hough變換
%?將圖像空間（xy）對(duì)應(yīng)到參數(shù)空間（abr）
%?a?=?x-r*cos（angle）
%?b?=?y-r*sin（angle）
for?i=1:ecount
????for?r=1:size_r
????????for?k=1:size_angle
????????????a?=?round（rows（i）-（r_min+（r-1）*step_r）*cos（k*step_angle））;
????????????b?=?round（cols（i）-（r_min+（r-1）*step_r）*sin（k*step_angle））;
????????????if（a>0&a<=m&b>0&b<=n）
????????????????hough_space（abr）?=?hough_space（abr）+1;
????????????end
????????end
????end
end

%?搜索超過閾值的聚集點(diǎn)
max_para?=?max（max（max（hough_space）））;
index?=?find（hough_space>=max_para*p）;
length?=?size（index）;
hough_circle?=?false（mn）;
for?i=1:ecount
????for?k=1:length
????????par3?=?floor（index（k）/（m*n））+1;
????????par2?=?floor（（index（k）-（par3-1）*（m*n））/m）+1;
????????par1?=?index（k）-（par3-1）*（m*n）-（par2-1）*m;
????????if（（rows（i）-par1）^2+（cols（i）-par2）^2<（r_min+（par3-1）*step_r）^2+5&...
????????????????（rows（i）-par1）^2+（cols（i）-par2）^2>（r_min+（par3-1）*step_r）^2-5）
????????????hough_circle（rows（i）cols（i））?=?true;
????????end
????end
end

%?打印檢測結(jié)果
for?k=1:length
????par3?=?floor（index（k）/（m*n））+1;
????par2?=?floor（（index（k）-（par3-1）*（m*n））/m）+1;
????par1?=?index（k）-（par3-1）*（m*n）-（par2-1）*m;
????par3?=?r_min+（par3-1）*step_r;
????fprintf（1‘Center?%d?%d?radius?%d\n‘par1par2par3）;
????para（:k）?=?[par1par2par3];
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????406??2005-12-17?12:39??main.m

?????文件??????36352??2005-12-17?13:37??Hough變換.doc

?????文件??????46926??2005-12-16?11:34??2.bmp

?????文件???????1973??2005-12-17?13:33??hough_circle.m

*????文件???????4862??2006-09-29?09:08??中國圖象圖形網(wǎng)下載說明.html

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

????????????????90519????????????????????5