function features
% Get the feature image
global IMAGE DATA CLASSIFY_PARAMS

%[fmsk,fmsk_horizon] = cu_farmsk(CLASSIFY_PARAMS.minDst,CLASSIFY_PARAMS.maxDst);
%CLASSIFY_PARAMS.fmsk_horizon = fmsk_horizon;
%CLASSIFY_PARAMS.fmsk = fmsk;


fim = double(colrHist(IMAGE.rgb,5,0,255,5,5));
%fim = double(IMAGE.rgb);


[image_height, image_width, num_features_per_pixel] = size(fim);
DATA.X = double(reshape(fim, image_height * image_width, num_features_per_pixel));
DATA.image_height = image_height;
DATA.image_width = image_width;

if (CLASSIFY_PARAMS.save_images == 1) && (CLASSIFY_PARAMS.img_cnt < CLASSIFY_PARAMS.max_imag)
   curr_imag = IMAGE.rgb;
   str = ['save img' num2str(CLASSIFY_PARAMS.img_cnt) '  curr_imag'];
   CLASSIFY_PARAMS.img_cnt = CLASSIFY_PARAMS.img_cnt + 1;
   eval(str);
end


return;