for h=1:11
imageName=strcat(num2str(h),'.bmp');
H=imread(imageName);
[u,v]=size(H);
if rem(u,3)==0
m=u;
else m=u-rem(u,3);
end
if rem(v,3)==0
n=u;
else n=v-rem(v,3);
end %如果图片尺寸不是3*3的倍数，自动裁剪为3*3的倍数%
G=zeros(m,n);
for a=1:u;
for b=1:v
G(a,b)=H(a,b);
end
end
N=zeros(3,3,m*n/9);
I=zeros(1,9,m*n/9);
z=1;
for i=1:3:m
for j=1:3:n
N(:,:,z)=G(i:i+2,j:j+2);
z=z+1;
end
end %分块%
A=zeros(1,9,m*n/9);
for i=1:(m*n/9)
I(:,:,i)=reshape(N(:,:,i),1,9);
A(:,:,i)=sort(I(:,:,i));
end %升序排序%
for j=1:(m*n/9)
for i=1:9 z=1;
if I(1,i,j)~=A(1,5,j)
z=z+1;
else break
end
end
h1(1,j)=i;
end
count=unique(h1);
count1=hist(h1,unique(h1));
count2=zeros(1,9);
for j=1:length(unique(h1))
count2(1,count(1,j))=count1(1,j);
end
H1=count2/(m*n/9); %hDBM%
x=0;
for p=1:(m*n/9)
y=0;
x=x+1;
for i=1:9
if I(1,5,p)==I(1,i,p)
y=y+1;
h2(1,x)=y;
end
end
end
count3=unique(h2);
count4=hist(h2,unique(h2));
count5=zeros(1,9);
for i=1:length(unique(h2))
count5(1,i)=count4(1,i);
end
H2=count5/(m*n/9); %hOBC%
for i=1:(m*n/9)
h3(1,i)=length(unique(I(:,:,i)));
end
count6=unique(h3);
count7=hist(h3,unique(h3));
count8=zeros(1,9);
for i=1:length(unique(h3))
count8(1,count6(1,i))=count7(1,i);
end
H3=count8/(m*n/9); %hQGL%
for j=1:(m*n/9)
for i=1:9 z=1;
if A(1,i,j)~=I(1,5,j)
z=z+1;
else break
end
end
h4(1,j)=i;
end
count9=unique(h4);
count10=hist(h4,unique(h4));
count11=zeros(1,9);
for j=1:length(unique(h4))
count11(1,count9(1,j))=count10(1,j);
end
H4=count11/(m*n/9); %hDBC%
for j=1:(m*n/9)
for i=1:9
if I(1,5,j)~=A(1,i,j)
i=i;
else break
end
end
h5(1,j)=i;
end
count12=unique(h5);
count13=hist(h5,unique(h5));
count14=zeros(1,9);
for i=1:length(count13)
count14(1,count12(1,i))=count13(1,i);
end
H5=count14/(m*n/9); %hFBC%
f(1,h)=(H1(1,5)*H2(1,2)*H3(1,6)/(H4(1,1)*H3(1,9)))*(H4(1,3)+H4(1,7)-H4(1,2)-H4(1,8))*(H5(1,3)/(H5(1,2)*H5(1,9)))/(H4(1,2)+H4(1,8)-H4(1,1)-H4(1,9));%f%
end
评论0