for i=1:3
level=graythresh(I{i});
I1{i}=im2bw(I{i},level);
end
dirname='F:\moni2012\2013B\B\附件4';
files=dir([dirname,'\*.bmp']);
Num_file=numel(files);
M=cell(1,Num_file);
for j = 1:Num_file
M{j} = imread([dirname '\' files(j).name]);
end
for j=1:209
level=graythresh(M{j});
M1{j}=im2bw(M{j},level);
end
for i=1:3
Z{i}=I1{i}(:,[10,20,30,40,50,60]);
end
for i=1:3
ZS{i}=sum(Z{i},2)
end
for j=1:209
CS{j}=sum(M{j},2)
end
for i=1:3
for j=1:209
J{i,j}=ZS{i}-CS{j}
JS{i,j}=sum(abs(J{i,j}))
end
end
[X,Y]=sort(JS{i,j},2)