function gaborResponses = computeGaborResponses(im, gabor, NumFilters)

%init
gaborResponses = zeros(size(im,1), size(im,2), NumFilters);

for oriIndex = 1 : NumFilters
    
    filter = real(gabor(oriIndex).filter);
    
    %normalize
    %%%%%%%%%%%%%%%%%%%%%
    filter = filter ./ max(filter(:));
    %%%%%%%%%%%%%%%%%%%%%
    
    gaborResponses(:, :, oriIndex) = imfilter(im, filter, 'replicate', 'same', 'conv');
    
    %gaborResponses(:, :, oriIndex) = abs(imfilter(im, gabor(oriIndex).filter, 'replicate', 'same', 'conv'));
    %gaborResponses(:, :, oriIndex) = gaborResponses(:, :, oriIndex) - mean2(gaborResponses(:, :, oriIndex));
    
%     figure,
%     imshow(gaborResponses(:, :, oriIndex), [])
%     pause

    
end %for oriIndex = 1 : param.divTheta