function N_features=cp_gabor(x,FeatureFN,NonlinId,SmoothId,NormId,Decimation)
% CP_GABOR  Comparative texture measure study: Gabor features
%	N_features=cp_gabor(x,FeatureFN,NonlinId,SmoothId,NormId,Decimation)

N_features=0;
[M_x,N_x]=size(x);

for CtrFreq=2.^(-6:-2)*sqrt(2)
  for Angle=(0:3)*pi/4
    N_features=N_features+1;
    fprintf(1,'.');
    y=sepgabor(x,CtrFreq,Angle);
    v=cp_energy(y,NonlinId,SmoothId,NormId,[1 1]*CtrFreq); % Local energy func.
    catframe(v,FeatureFN,[],'Interm',Decimation); % Write feature file
  end
end

