function?[para]?=?findcircle（BWstepRstepAngleminRmaxRp）

[mn]?=?size（BW）;%?BW:二值圖像；
cntR?=?round（（maxR-minR）/stepR）+1;%?stepR:檢測的圓半徑步長，?minR:最小圓半徑，?maxR:最大圓半徑
cntAngle?=?round（2*pi/stepAngle）;%?stepAngle:角度步長，單位為弧度
hough_space?=?zeros（mncntR）;%?hough_space:參數空間，h（abr）表示圓心在（ab）半徑為r的圓上的點數
[rowscols]?=?find（BW）;
cntPoints?=?size（rows1）;
?
%?Hough變換將圖像空間（xy）對應到參數空間（abr）
%?a?=?x-r*cos（angle）?b?=?y-r*sin（angle）
for?i=1:cntPoints
????for?r=1:cntR
????????for?k=1:cntAngle
????????????a?=?round（rows（i）-（minR+（r-1）*stepR）*cos（k*stepAngle））;
????????????b?=?round（cols（i）-（minR+（r-1）*stepR）*sin（k*stepAngle））;
????????????if（a>0?&&?a<=m?&&?b>0?&&?b<=n）
????????????????hough_space（abr）?=?hough_space（abr）+1;
????????????end
????????end
????end
end
?
%?尋找滿足閾值的圓的參數
max_para?=?max（max（max（hough_space）））;
index?=?find（hough_space>=max_para*p）;?%?p:以p*hough_space的最大值為閾值，p取0，1之間的數
length?=?size（index1）;
hough_circle=zeros（mn）;
for?i=1:cntPoints
????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<（minR+（par3-1）*stepR）^2+5?&&?（rows（i）-par1）^2+（cols（i）-par2）^2>（minR+（par3-1）*stepR）^2-5）
????????????hough_circle（rows（i）cols（i））?=?1;%?hough_circl:二值圖像，檢測到的圓
????????end
????end
end
?
for?k=1:length
????par3?=?floor（index（k）/（m*n））+1;?????
????par2?=?floor（（index（k）-（par3-1）*（m*n））/m）+1;????%?圓心y坐標
????par1?=?index（k）-（par3-1）*（m*n）-（par2-1）*m;??????%?圓心x坐標
????par3?=?minR+（par3-1）*stepR;????????????????????%?圓的半徑
???%?fprintf（1‘Center?%d?%d?radius?%d\n‘par1par2par3）;
????para（k:）?=?[par1par2par3];??%?para:檢測到的圓的圓心、半徑
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????37906??2016-04-26?14:30??coins.png
?????文件????????1879??2016-04-26?15:15??findcircle.m
?????文件????????1338??2016-04-26?15:30??main.m
?????文件??????397455??2016-04-26?14:30??pillsetc.png
?????文件??????284401??2016-04-26?14:30??tape.png