function ClassMap = recognition_phase(I, MODEL)
   ClassMap = zeros(size(I));
   N = 1/(15*15) * ones(15);
   BB=zeros(size(I,1),size(I,2),9);
      for k=1:9
          B=conv2(I,laws_kernel(k),'same');
          B_i=conv2(B.^2,N,'same');
          BB(:,:,k)=B_i(:,:);
      end
      
      for i = 1:size(B_i,1)
          for j = 1:size(B_i,2)
              for n = 1:size(MODEL,1)
                  buff(:)=BB(i,j,:);
                  sum_abs_diff(n)=sum(abs(buff(1,:)-MODEL(n,:)));
                  [M,ind]=min(sum_abs_diff);
                  ClassMap(i,j)=ind;
              end
          end
     end
end