function ellipsedetect()
clc;
figure(1);
%----读原图并显示
RGB=ShowGRB();
%----将原图转换为二值图并显示
bw=Show2bw(RGB);
%---- 应用形态学函数去噪,
bw=quzao(bw);
%---- 闭运算后的图像
bw=fengbi(bw);
%---填充封闭区域或者空洞
bw=fillim(bw);
%---------------Step 4: 寻找边界
[B,L] = bwboundaries(bw,'noholes'); %bwboundaries寻找二值图像的边界,“noholes”为外边界,返回值B为边界像素所在的坐标,L??????
figure(2)
imshow(label2rgb(L, @jet, [0.5 0.5 0.5])) ; % label2rge把图像标记为@jet彩色,背景默认为灰色
hold on
for K = 5:length(B) % length行列中的最大值
boundary = B{K}; % B{k}每一个图形边界对应的坐标值
plot(boundary(:,2), boundary(:,1), 'k', 'LineWidth', 2)%描画每一个图形的边界
end
%-----------------Step 5: 决定那个图形为圆形
area_tru=0;
area_cal=0;
sta