function?Center=CenterDetect（In_cn_rrange）
%range表示霍夫檢測圓的半徑范圍
%n_c表示水平方向上，圓的個數
%n_r表示豎直方向上，圓的個數

figure（1）
imshow（I）
[BWxiyi]?=?roipoly;
%構建直線方程，用于后面圓點排序（重要）
flag=1;?%表示斜率存在
if?xi（1）~=xi（4）
????lk=（yi（4）-yi（1））/（xi（4）-xi（1））;??%斜率
????lb=yi（1）-lk*xi（1）;
else
????LX=xi（1）;
????flag=0;
end
BW=double（BW）;
NewI=double（I）.*BW;
NewI=uint8（NewI）;
level?=?graythresh（NewI）;
A=im2bw（NewIlevel）;

%centers為圓心坐標，radii為半徑
[centersradii]?=?imfindcircles（~Arange）;

%對圓心進行排序
figure（2）
imshow（I）
hold?on
pp=centers;

%計算中心點到左端直線的x方向的（由原點的分布特性，可以不用求點到直線的距離，從而提高算法效率）
if?flag
????xe=（pp（:2）-lb）/lk;??%中間變量
????pp（:3）=pp（:1）-xe;
else
????pp（:3）=pp（:1）-LX;
end

Sort_R=sortrows（pp3）;
psort_row=Sort_R（:1:2）;

%再按照y的大小，對每一列的圓心進行排序
Center=[];
for?k=1:n_c
????exm=psort_row（（k-1）*n_r+1:k*n_r:）;
????exm=sortrows（exm2）;
????Center=[Center;exm];
end

for?k=1:length（Center（:1））
????plot（Center（k1）Center（k2）‘*‘）;
????%?????text（centers（i1）+5centers（i2）+5num2str（i））;
????text（Center（k1）+5Center（k2）+5num2str（k）‘FontSize‘10‘Color‘‘red‘‘FontWeight‘‘Bold‘）;
end
viscircles（centers?radii‘EdgeColor‘‘b‘）;
Center=Center‘;??%轉置
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????1921078??2014-11-23?19:09??三種圓心檢測精度對比\001.bmp

?????文件????1921078??2014-11-23?19:10??三種圓心檢測精度對比\007.bmp

?????文件???????1351??2017-07-05?19:21??三種圓心檢測精度對比\CenterDetect.m

?????文件???????1756??2018-11-28?20:42??三種圓心檢測精度對比\CenterDetectNew.m

?????文件???????1338??2017-07-06?22:43??三種圓心檢測精度對比\circle_fitting.m

?????文件???????4236??2019-02-26?10:25??三種圓心檢測精度對比\Compare_Three_CenterMethod.m

?????文件???????5303??2019-01-13?23:22??三種圓心檢測精度對比\兩種方法檢測的圓心比較.mat

?????目錄??????????0??2019-02-26?10:28??三種圓心檢測精度對比

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

??????????????3856140????????????????????8