Igray = rgb2gray(I) Ibw = im2bw(Igray,graythresh(Igray)) Iedge = edge(uint8(Ibw)) se = strel( square ,3) Iedge2 = imdilate(Iedge, se) Ifill= imfill(Iedge2, holes ) [Ilabel num] = bwlabel(Ifill) Iprops = r