BW = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
B = bwboundaries(BW);
b = B{1};
imshow(BW,'InitialMagnification','fit');
hold on
[M,N] = size(BW);
x = [0.5 (N+0.5)];
for k = 1.5:(M-0.5)
y = [k k];
plot(x,y,'Color',[.7 .7 .7],'LineStyle','--');
end
y = [0.5 (M+0.5)];
for k = 1.5:(N-0.5)
x = [k k];
plot(x,y,'Color',[.7 .7 .7],'LineStyle','--');
end
plot(b(:,2),b(:,1),'r')
hold off
BW3 = imresize(Z,3,'nearest');
B3 = bwboundaries(BW3);
b3 = B3{1};
b3 = (b3 + 1)/3;
imshow(BW,'InitialMagnification','fit');
hold on
[M,N] = size(BW);
x = [0.5 (N+0.5)];
for k = 1.5:(M-0.5)
y = [k k];
plot(x,y,'Color',[.7 .7 .7],'LineStyle','--');
end
y = [0.5 (M+0.5)];
for k = 1.5:(N-0.5)
x = [k k];
plot(x,y,'Color',[.7 .7 .7],'LineStyle','--');
end
plot(b3(:,2),b3(:,1),'r')
hold off