latocella = 3;
passo = 10;
χ[S3medie_] := ( 1/2 ) ArcSin[Abs[S3medie]];
a[S3medie_] := 8.0*Cos[χ[S3medie]] ;
b[S3medie_] := 8.0*Sin[χ[S3medie]] ;
ψ[S1medie_, S2medie_] :=
  If[S1medie == 0 && S2medie == 0, 0, (1/2) ArcTan [S1medie, S2medie]];
Ellipse[i_, j_, S1medie_, S2medie_, S3medie_] :=
 Rotate[Circle[{i, j}, {a[S3medie], b[S3medie]}], ψ[S1medie,
   S2medie]]
fun[i_, j_] :=
 finalmatrix[[j*latocella + (latocella + 1)/2,
    900 - i*latocella - (latocella + 1)/2]] =
  Ellipse[j*latocella, 899 - i* latocella, S1medie[[i, j]],
   S2medie[[i, j]], S3medie[[i, j]]]

finalmatrix = Table[Null, {l, 1, 900}, {k, 1, 900}];
Outer[fun, Range[1, 299, passo], Range[1, 299, passo]];
Graphics[{Thick, finalmatrix}, ImageSize -> {900, 900}]

ParametricPlot[
 RotationMatrix[ψ[S1medie, S2medie]].{i +
    8.0*Cos[χ[S3medie]]*Cos[x],
   j + 8.0*Sin[χ[S3medie]] *Sin[x]} , {x, 0, 2 π}]