twohistogram2.m

求灰度空间相关性图及计算标准差及标准差的平均值

clear;
I=imread('system4.bmp');
I2=imresize(I,0.2,'bilinear');
[m,n]=size(I2);
for x=1:256;
    for y=1:256;
        W(x,y)=0;
        V(x,y)=0;
    end;
end;
for i=1:m;
    for j=1:n;
        s=double(I2(i,j));
        R((s+1),i,j)=i;
        C((s+1),i,j)=j;
    end;
end;
a=double(max(max(I2)));
for s=1:(a+1);
    for i=1:m;
        for j=1:n;
            if ((R(s,i,j)~=0)&&(C(s,i,j)~=0));
                p=double(R(s,i,j));
                q=double(C(s,i,j));
                
                if (((p>1)&&(p<m))&&((q>1)&&(q<n)))
                    u=double(I2(p,(q+1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==1)&&(q==1));
                    u=double(I2(p,(q+1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==1)&&(q>1)&&(q<n));
                    u=double(I2(p,(q+1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==1)&&(q==n));
                    u=double(I2(p,(q-1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p>1)&&(p<m)&&(q==n));
                    u=double(I2((p+1),q));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==m)&&(q==n));
                    u=double(I2((p-1),(q-1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==m)&&(q>1)&&(q<n));
                    u=double(I2(p,(q+1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p==m)&&(q==1));
                    u=double(I2(p,(q+1)));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                    if ((p>1)&&(p<m)&&(q==1));
                    u=double(I2((p+1),q));
                    W(s,(u+1))=W(s,(u+1))+1;
                end;
                end;
            end;
        end;
    end;
    for i=1:256;
            for j=1:256;
                V(i,j)=W((257-i),j);
            end;
        end;
imshow(V);
    imwrite(V,'pingjia4.bmp');  
    BS=biaozhuncha(W);
    B=biaoarea(BS);