Bei meiner Arbeit habe ich einige Wartungsarbeiten an interner Software durchgeführt, die behauptet, eine stereografische Projektion zu implementieren, deren Formeln vor langer Zeit von einem alten System kopiert wurden.
Soweit ich gesehen habe, scheint es das "Oblique and Equatorial Stereographic" von EPSG (Code 9809) zu implementieren. Aber es gibt einen Schritt, den ich nicht identifizieren konnte. Der in unserem Code verwendete "R" -Wert wird mit den folgenden Parametern berechnet:
(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
Vor diesem Hintergrund lautet die Formel für R:
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
Weiß jemand was hier los ist? Ich habe über stereografische Projektionen gesucht und gelesen, aber ich konnte so etwas nicht finden. Tatsächlich habe ich keine Projektionen gesehen, die diese "Null-Verzerrungskoordinate" verwenden.