function?[hough_spacehough_circlepara]?=?hough_circle（BWstep_rstep_angler_minr_maxp）;
%[HOUGH_SPACEHOUGH_CIRCLEPARA]?=?HOUGH_CIRCLE（BWSTEP_RSTEP_ANGLER_MAXP）
%------------------------------算法概述-----------------------------
%?該算法通過a?=?x-r*cos（angle），b?=?y-r*sin（angle）將圓圖像中的邊緣點
%?映射到參數空間（abr）中，由于是數字圖像且采取極坐標，angle和r都取
%?一定的范圍和步長，這樣通過兩重循環（angle循環和r循環）即可將原圖像
%?空間的點映射到參數空間中，再在參數空間（即一個由許多小立方體組成的
%?大立方體）中尋找圓心，然后求出半徑坐標。
%-------------------------------------------------------------------

%------------------------------輸入參數-----------------------------
%?BW:二值圖像；
%?step_r:檢測的圓半徑步長
%?step_angle:角度步長，單位為弧度
%?r_min:最小圓半徑
%?r_max:最大圓半徑
%?p:以p*hough_space的最大值為閾值，p取0，1之間的數
%-------------------------------------------------------------------

%------------------------------輸出參數-----------------------------
%?hough_space:參數空間，h（abr）表示圓心在（ab）半徑為r的圓上的點數
%?hough_circl:二值圖像，檢測到的圓
%?para:檢測到的圓的圓心、半徑
%-------------------------------------------------------------------

%?From?InternetModified?by?mhjerry2011-12-11

[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