clear all;
clc;
im=double(imread('c:\bmp\10\10bluembw.bmp'));
w=mean(im);

%plot(w);
%return

colthreshold=mean(w);
wc=find(w>=colthreshold);
k=2:size(wc,2);
wcmc=wc(k)-wc(k-1);
while std(wcmc)>0.2
    mc=mean(wcmc);
    wcmc(find(wcmc<mc))=[];
end
mc=mean(wcmc)
sc=std(wcmc)

w=mean(im');
rowthreshold=mean(w);
wc=find(w>=rowthreshold);
k=2:size(wc,2);
wcmr=wc(k)-wc(k-1);
while std(wcmr)>0.2
    mr=mean(wcmr);
    wcmr(find(wcmr<mr))=[];
end
mr=mean(wcmr)
sr=std(wcmr)