3,423
社区成员
发帖
与我相关
我的任务
分享
clear;close all;
im = imread('https://img-bbs.csdn.net/upload/201404/11/1397185779_271987.jpg');
BW=im2bw(im,0.5);
s = regionprops(BW, 'ConvexHull');
imshow(BW);
% hold on;
% plot(s(1).ConvexHull(:,1) ,s(1).ConvexHull(:,2) , 'r*');
%
% hold on;
% plot(s(2).ConvexHull(:,1) ,s(2).ConvexHull(:,2) , 'g*');
%
% hold on;
% plot(s(3).ConvexHull(:,1) ,s(3).ConvexHull(:,2) , 'b*');
[m1 ,~]=size(s(1).ConvexHull);
[m2 ,~]=size(s(2).ConvexHull);
len=length(s);
for i=1:len-1
for j=i+1:len
[m1 ,~]=size(s(i).ConvexHull);
[m2 ,~]=size(s(j).ConvexHull);
for h=1:m1
for k=1:m2
if pdist([s(i).ConvexHull(h,:) ;s(j).ConvexHull(k,:)],'euclidean')<100
hold on;
plot([s(i).ConvexHull(h,1) ,s(j).ConvexHull(k,1)] ,[s(i).ConvexHull(h,2) ,s(j).ConvexHull(k,2)],'r');
end
end
end
end
end