20 lines
535 B
Mathematica
20 lines
535 B
Mathematica
|
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;
|