?function???[xyrcount]=hough（BW1x0y0r0rr）???
??%???BW1為輸入的需檢測的圓圖象的矩陣[x0y0]為估算的定位中心的位置r0為估算的定位的半徑rr是檢測范圍的外擴值???
??[mn]=size（BW1）;???
??%確定搜索范圍，搜索的矩形的四個定點為[e3e1][e3e2][e4e1][e4e2]???
??e1=floor（x0-r0）-rr;???
??if???e1<1???
??????????e1=1;???
??end???
??e2=ceil（x0+r0）+rr;???
??if???e2>n???
??????????e2=n;???
??end???
??e3=floor（y0-r0）-rr;???
??if???e3<1???
??????????e3=1;???
??end???
??e4=ceil（y0+r0）+rr;???
??if???e4>m???
??????????e4=m;???
??end???
??%Hough變換求瞳孔中心???
??flag=1;count=0;?????????????????????????%flag為循環控制變量，count記錄落在所得圓上的點???
??while（flag）???
??????????flag=0;???
??????????BW2=zeros（mn）;?????????????????%設0矩陣BW2，以BW2為變換空間???
??????????for???j=e1:e2???????????????????????????
??????????????????for???i=e3:e4???
??????????????????????????if???BW1（ij）==1???????????????????????????????????????????%在搜索范圍內找到一個值為1點???
??????????????????????????????????k=（j-x0）^2+（i-y0）^2-r0^2;???????????????
??????????????????????????????????if???k>-1600&k<1600?????????????????????????????%這個點在有效區內???
??????????????????????????????????????????real=j-x0;imag=y0-i;???
??????????????????????????????????????????theta=atan2（imagreal）;??