yaogan.rar_site:www.pudn.com_对比度_对比度计算_计算 对比度_计算对比度

A=[2 6 10 14 2 6 10;6 10 14 2 6 10 14;10 14 2 6 10 14 2;14 2 6 10 14 2 6;2 6 10 14 2 6 10;6 10 14 2 6 10 14;10 14 2 6 10 14 2]; B=(A-2)/(14-2)*3; L=4; a=1;b=0; p=0;t=0; P=[];M=0; K=zeros(4,4);X=[]; Uy=zeros(4);Ux=zeros(4);Vx=zeros(4);Vy=zeros(4); for i=0:(L-1) for j=0:(L-1) for x=1:(7-a) for y=1:(7+b) if(B(x,y)==i)&&(B(x+a,y-b)==j) K(i+1,j+1)=K(i+1,j+1)+1; end end end end end for i=1:4 for j=1:4 p=K(i,j)+p; end end P=K()./p; for i=0:(L-1) for j=0:(L-1) M=(abs(i-j))^2*P(i+1,j+1)+M; end end for i=0:(L-1) for j=0:(L-1) Uy(i+1)=6*P(i+1,j+1) + Uy(i+1); t=P(i+1,j+1)+t; end Vy(i+1)=(i-Uy(i+1))^2*t + Vy(i+1); end for j=0:(L-1) for i=0:(L-1) Ux(j+1)=6*P(i+1,j+1)+Ux(j+1); t=P(i+1,j+1)+t; end Vx(j+1)=(j-Ux(j+1))^2*t+Vx(j+1); end for i=0:(L-1) for j=0:(L-1) X(i+1,j+1)=(i-Uy(i+1))*(j-Ux(j+1))*P(i+1,j+1)/( Vy(i+1)* Vx(j+1)); end end