function ClassMap = recognition_phase(I, MODEL)
      ClassMap=zeros(size(I,1),size(I,2));
      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');
          Bcs=conv2(B.^2,N,'same');
          BB(:,:,k)=Bcs(:,:);
      end
      
      for x=1:size(Bcs,1)
          for y=1:size(Bcs,2)
     
              for n=1:size(MODEL,1)
                  buff(:)=BB(x,y,:);
                  sum_abs_diff(n)=sum(abs(buff(1,:)-MODEL(n,:)));
                  [M,ind]=min(sum_abs_diff);
                  ClassMap(x,y)=ind;
         end
       end
     end
end