Dans mon travail, j'ai effectué une maintenance sur un logiciel interne qui prétend implémenter une projection stéréographique dont les formules ont été copiées depuis un ancien système il y a longtemps.
D'après ce que j'ai vu, il semble implémenter la "stéréographie oblique et équatoriale" d'EPSG (code 9809). Mais il y a une étape que je n'ai pas pu identifier. La valeur "R" utilisée dans notre code est calculée avec les paramètres suivants:
(phi_n, lambda_n) - null distortion coordinates
(PHI_nc, LAMBDA_nc) - conformal null distortion coordinates
(PHI_o, LAMBDA_o) - conformal projection origin coordinates
e - eccenttricity
r_eq - earth radius at equator
Compte tenu de cela, la formule pour obtenir R est:
a = 1 + cos(PHI_nc) * cos(PHI_o) * cos(LAMBDA_nc - LAMBA_o)
b = sin(PHI_nc) * sin(PHI_o)
c = 2 * cos(PHI_nc) * [1 - (sin(phi_n) * e/2) ^ 2] ^ 1/2
R = r_eq * cos(phi_n) * (a + b) / c
Quelqu'un sait-il ce qui se passe ici? J'ai cherché et lu des projections stéréographiques, mais je n'ai rien trouvé de tel. En fait, je n'ai vu aucune projection utilisant cette "coordonnée de distorsion nulle".