clc clear all I=imread('a003.bmp'); [m,n,p]=size(I); I=im2bw(I); se=strel('square',3); Ia=imerode(I,se); Iout=I-Ia; figure(1); imshow(Iout); Iout1=Iout; N=1; for i=1:m for j=1:n if Iout(i,j)==1 p(N,1)=i; <