clc
clear all
load p
p=im2bw(p);
% BW = logical(p);
% L = bwlabel(BW,4);
% [r,c] = find(L == 2);
BW = logical(p);
[r,c] = find(BW == 1);
[m n]=size(p)
q=zeros(m,n);
count=1
for j=1:n;
count=1
for i=1:m-1;
if (xor(BW(i,j),BW(i+1,j))==0)&(BW(i,j)~=1)
count=count+1;
elseif (xor(BW(i,j),BW(i+1,j))~=0)
if BW(i,j)==1
q(i,j)=count;
end
end
end
end