function A=imlog(I,sigma)
% A=imlog(I,sigma)
% Esegue un filtraggio dell'immagine I
% con filtro LOG di sigma assegnata. Restituisce in A il risultato

w=ceil(3*sigma)*2+1;
size=(w-1)/2;
x=[-size:size];
y=[-size:size];
for i= 1 : w
    for j= 1 : w
        arg=(x(i))^2+(y(j))^2;
        G(i,j)=((arg-2*sigma^2)/(2*pi*sigma^6))*exp(-arg/(2*sigma^2));
    end
end
somma=sum(sum(G));
Y=G/somma;

A=filter2(Y,I);