www.pudn.com > ICA.rar > ContrastRatio.m, change:2005-07-11,size:427b


function [valueCR,valueENL] = ContrastRatio(I1) 
 
% EQIM --- Return the ContrastRatio of input vector X 
% Date: April,20,2004 
%***********************X is a vector  
%***********************Yis ContrastRatio 
 
 
for i=1:100 
    for j=1:128 
        temp(i,j)=I1(i+64,j+126); 
    end 
end 
vecorX=reshape(temp',1,128*100); 
 
 
temp1=std(vecorX); 
temp2=mean(vecorX); 
valueCR=temp1/temp2; 
valueENL=(1/valueCR)^2;