code: label = bwlabel(bw,4); stats = regionprops(label,'Centroid','Area'); >> for n = 1:length(stats) disp(sprintf('Blob number = %d, Area = %g, Centroid = (%g, %g)',... n,stats(n).Area,stats(n).Centroid)) end Blob number = 1, Area = 22, Centroid = (11.3182, 5.81818) Blob number = 2, Area = 80, Centroid = (31.825, 114.138) Blob number = 3, Area = 4, Centroid = (50, 144) Blob number = 4, Area = 3, Centroid = (58.6667, 54) Blob number = 5, Area = 166, Centroid = (108.042, 46.8675) Blob number = 6, Area = 186, Centroid = (73.2473, 94.0484) Blob number = 7, Area = 2, Centroid = (61, 114.5) Blob number = 8, Area = 2, Centroid = (63, 114.5) Blob number = 9, Area = 12, Centroid = (75.1667, 154.5) Blob number = 10, Area = 14, Centroid = (79.9286, 142.286) Blob number = 11, Area = 1, Centroid = (79, 139)