www.pudn.com > ICA.rar > logeqim.m, change:2005-03-10,size:505b


function Y = logeqim(X) 
 
xmax = max(max(X)); 
xmin = min(min(X)); 
[N,M]=size(X); 
tmpmean1=mean(X); 
tmpstd1=std(X); 
tmpvalue1=tmpstd1/tmpmean1; 
for (k =1:N) 
    for(j=1:M) 
        if (xmin<=0) 
           X(k,j)=X(k,j)-xmin+1; 
        elseif (xmin>=0) & (xmin<1) 
            X(k,j)=X(k,j)+1; 
        else 
            X(k,j)=X(k,j); 
       end 
    end 
end 
x1max = max(max(X)); 
x1min = min(min(X)); 
 
 
for (r=1:N) 
    Y(r,:)=((log(X(r,:))-log(x1min))*255/(log(x1max)-log(x1min))); 
end