多边形最大内接圆的简单画法

        上一篇文章介绍了一种多边形最小外接圆的画法。今天给大家介绍一种基于MATLAB的多边形最大内接圆的简单画法。话不多说，直接上代码。

函数

function [ X,Y,R ] = InsideCircle( BW )
%此函数用于计算多边形最大内接圆
%输入BW为含有一个多边形区域的二值图像，输出为最大内接圆的坐标和圆心
D = bwdist(~BW);
regionmax = imregionalmax(D);%寻找区域最大值
[x,y]=find(regionmax==1);
X=x(1);Y=y(1);%注：因为D中可能有两个极大值点
R=floor(D(X,Y));%向小取整
end

测试程序

 clc
img=zeros(20,20);%生成一个20×20的图
for i=8:15
 
    for j=8:i
        img(i,j)=1;
    end
end
img(16,9)=1;
img(17,9)=1;
img(18,9)=1;img(18,10)=1;%设置好多边形
figure(1)
subplot(121)
imshow(img)

[x,y,r]=InsideCircle(img);%调用函数，计算内接圆位置和大小
subplot(122)
imshow(img)
hold on
viscircles([y,x],r,'color','b');
text(y,x,'o');
hold off   

测试结果

如有错误，欢迎指正。