grayexitmonentr.m

尿沉渣图像识别

function T=GrayExitMonentR(IN)
[R,C]=size(IN);
gray=256;
p1=zeros(gray);
for M=1:R
    for N=1:(C-1)
        p1(IN(M,N)+1,IN(M,N+1)+1)= p1(IN(M,N)+1,IN(M,N+1)+1)+1;
        p1(IN(M,N+1)+1,IN(M,N)+1)= p1(IN(M,N+1)+1,IN(M,N)+1)+1;
    end
end

%NOMALIZATION
p1=double(p1)/numel(IN);

%COMPUTE SECONDMONENT
p=p1.^2;
secondmonent=sum(p(:));
T=secondmonent;
%COMPUTE CORRELATION
% colour=0:R-1;
% m1=sum(p1)*colour(:);
% m2=m1;
% pp=sum(p1);
% v1=((colour-m1).^2)*pp(:);
% % v2=v1;
% [m,n]=meshgrid(colour,colour);
% pp=m.*n.*p1;
% relation=(sum(pp(:))-m1*m2)/(v1+eps)*(v2+eps)
% %COMPUTE EMTROPY
% pp=(p1.*log2(p1+eps));
% emtropy=-sum(pp(:))
%COMPUTE CONTRAST
% M=zeros(R);   %INPUT IMAGE MUST BE SQUARE MONENT
% for k=1:R
%     M(k,k)=1;
% end
% % contrast=sum(sum(p1.*M));
% M=zeros(R);
% for k=1:R-1
%     M(k,k+1)=1;
% end
% for k=1:R-1
%     contrast=contrast+sum(sum(p1.*M^k))*k^2*2;
% end
% contrast