上一篇文章介绍了一种多边形最小外接圆的画法。今天给大家介绍一种基于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
测试结果
如有错误,欢迎指正。