function imfil=homomorfico(c,r,n,gh,gl)
clc;
t=size(c);
tfm=round(t(1)/2);
tcm=round(t(2)/2);
for i=1:t(1)
       for j=1:t(2)
           h=sqrt((i-tfm)^2+(j-tcm)^2);
           imfil(i,j)=((gh-gl)*(1-exp(-n*(h^2/r^2))))+gl;
       end
end