En mi trabajo, he estado haciendo un mantenimiento en el software interno que afirma implementar una proyección estereográfica cuyas fórmulas se copiaron de un sistema antiguo hace mucho tiempo.
Por lo que he visto, parece implementar el "Stereographic oblicuo y ecuatorial" de EPSG (código 9809). Pero hay un paso que no pude identificar. El valor "R" utilizado en nuestro código se calcula con los siguientes parámetros:
(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
Dado esto, la fórmula para obtener R es:
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
¿Alguien sabe qué está pasando aquí? Estuve buscando y leyendo sobre proyecciones estereográficas, pero no pude encontrar algo así. De hecho, no he visto ninguna proyección que use esta "coordenada de distorsión nula".
fuente