function resultat = calcul_I_moyen_elli(I,c,a,e,teta)
[nb_lignes,nb_colonnes] = size(I);

nb_pixels = 0;
somme_nvg = 0;
for j = max(1,floor(c(1)-a)):min(nb_colonnes,ceil(c(1)+a))
	for i = max(1,floor(c(2)-a)):min(nb_lignes,ceil(c(2)+a))
        f = a * e;
        f1 = f * [cos(teta) sin(teta)];
        f2 = f * [cos(teta) sin(teta)];
        d1 = norm([j i] - c - f1);
        d2 = norm([j i] - c - f2);
		if d1 + d2 < 2 * a
			nb_pixels = nb_pixels+1;
			somme_nvg = somme_nvg+I(i,j);
		end
	end
end
resultat = somme_nvg/nb_pixels;