[w,h]=size(a);
for i=1:w
for j=1:h
s(a((i-1)*h+j)+1)=0;
end
end
for i=1:w
for j=1:h
s(a((i-1)*h+j)+1)=s(a((i-1)*h+j)+1)+1;
end
end
for i=1:w
for j=1:h
p(a((i-1)*h+j)+1)=s(a((i-1)*h+j)+1)/h/w;
end
end
t=max(max(a))+1;
[w,h]=size(a);
for i=1:w
for j=1:h
s(a((i-1)*h+j)+1)=0;
end
end
for i=1:w
for j=1:h
s(a((i-1)*h+j)+1)=s(a((i-1)*h+j)+1)+1;
end
end
for i=1:w
for j=1:h
p(a((i-1)*h+j)+1)=s(a((i-1)*h+j)+1)/h/w;
end
end
t=max(max(a))+1;