function res=epaicissement(I, f90, f45)
    res=zeros(size(I));
    for i=0:3
        res=addition(epaicissementF(I, f90), res);
        res=addition(epaicissementF(I, f45), res);
        f90=rot90(f90);
        f45=rot90(f45);
    end
end