图像边缘检测，轮廓拟合

在边缘检测中总会提取出不连续点，或伪轮廓。在这种情况下需要拟合出目标的轮廓，这样可以找到轮廓的数学表达式为后续的特征选取打下基础。博主用coins图像为例，用椭圆方程进行拟合，做出如下实验。

1、原图二值化

2、边缘检测（sobel算子）

3、填补孔洞

4、标记连通域

5、找到每个连通域坐标

6、用每个连通域坐标拟合出椭圆方程

7、在二值图像中画出每个椭圆函数

%%图像边缘检测和拟合轮廓
clc
clear
close all
%% 读取图像
I = imread('coins.png');
I = im2bw(I);                                   %二值化
I = imfill(I,'holes');                          %填补孔洞
[M,N] = size(I);
figure(1),imshow(I);title('原图');hold on
%% 选取待拟合坐标
% conicP = ginput(15);            
bw1 = edge(I,'sobel');                           %边缘检测
figure(2),imshow(bw1);title('边缘检测');
[L,num] = bwlabel(bw1);                           %标签
for i = 1:num
    [row,col] = find(L == i);
    conicP = zeros(length(row),2);
    conicP(:,1) = col;
    conicP(:,2) = row;
    figure(1),plot(conicP(:,1)', conicP(:,2)', 'xr');     %drawing sample points
%% 自定义椭圆函数拟合
    a0 = [1 1 1 1 1 1];
    f = @(a,x)a(1)*x(:,1).^2+a(2)*x(:,2).^2+a(3)*x(:,1).*x(:,2)+a(4)*x(:,1)+a(5)*x(:,2)+a(6);%建立方程
    p = nlinfit(conicP , zeros(size(conicP, 1), 1), f,[1 2 3 4 5 6]);
    syms x y
    conic = p(1)*x^2+p(2)*y^2+p(3)*x*y+p(4)*x+p(5)*y+p(6);
%% 在原图上显示拟合结果
    c = ezplot(conic,[0,N],[0,M]);
    figure(1),set(c, 'Color', 'Blue','LineWidth',2);
   
end
 hold off

实验结果图像如下，拟合结果为蓝线，原坐标为红点