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:參數(shù)空間，h（abr）表示圓心在（ab）半徑為r的圓上的點數(shù)
[rowscols]?=?find（BW）;
cntPoints?=?size（rows1）;
?
%?Hough變換將圖像空間（xy）對應到參數(shù)空間（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
?
%?尋找滿足閾值的圓的參數(shù)
max_para?=?max（max（max（hough_space）））;
index?=?find（hough_space>=max_para*p）;?%?p:以p*hough_space的最大值為閾值，p取0，1之間的數(shù)
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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????198919??2018-04-28?14:52??A題--吉林省第九屆大學生數(shù)學建模競賽.pdf
?????文件??????134895??2018-04-29?21:52??A題信息.docx
?????文件??????231102??2018-04-28?14:57??B題--吉林省第九屆大學生數(shù)學建模競賽.pdf
?????文件???????14538??2018-04-29?09:45??B題信息.docx
?????文件?????1126134??2018-05-01?13:47??DATA1.bmp
?????文件????????2370??2018-04-19?16:21??DATA1.mat
?????文件?????1126134??2018-05-01?13:49??DATA2.bmp
?????文件????????2580??2018-04-19?16:21??DATA2.mat
?????文件????????3925??2018-04-19?16:21??DATA3.mat
?????文件????????1879??2016-04-26?15:15??findcircle.m
?????文件?????????383??2018-04-29?14:04??findcircle_script.m
?????文件????????4517??2016-07-22?15:12??GeneralHough.py
?????目錄???????????0??2018-04-30?14:37??gh\
?????文件????????1204??2011-03-27?18:54??gh\1.jpg
?????文件???????10728??2011-03-26?19:10??gh\2.jpg
?????文件???????20331??2011-03-27?19:15??gh\3.jpg
?????文件????????1631??2011-03-27?19:16??gh\4.jpg
?????文件????????1051??2011-03-27?19:26??gh\5.jpg
?????文件????????1850??2011-03-27?19:28??gh\6.jpg
?????文件???????15243??2011-03-20?17:23??gh\7.jpg
?????文件????????2309??2011-03-20?17:24??gh\8.jpg
?????文件????????1106??2018-04-30?14:45??gh\ght.m
?????文件?????????837??2018-04-30?14:47??gh\main.m
?????文件????????1056??2011-03-27?21:14??gh\main1.m
?????文件???????39831??2011-03-27?16:27??gh\noise.jpg
?????文件????????6731??2011-03-27?16:29??gh\original.jpg
?????文件????????1059??2011-03-26?17:36??gh\rtable.m
?????文件????????2479??2011-03-26?17:48??gh\template1.jpg
?????文件????????2335??2011-03-27?15:13??gh\template2.jpg
?????文件????????2265??2011-03-27?15:17??gh\template3.jpg
?????目錄???????????0??2018-05-01?13:40??ght\
............此處省略229個文件信息