--- a
+++ b/functions/functions_FeatExtr/computeGaborResponses.m
@@ -0,0 +1,26 @@
+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
+