На своей работе я занимался некоторым обслуживанием внутреннего программного обеспечения, которое претендует на реализацию стереографической проекции, формулы которой были скопированы из старой системы давным-давно.
Из того, что я видел, кажется, реализуется «Наклонная и экваториальная стереография» из EPSG (код 9809). Но есть один шаг, который я не смог определить. Значение «R», используемое в нашем коде, рассчитывается по следующим параметрам:
(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
Учитывая это, формула для получения 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
Кто-нибудь знает, что здесь происходит? Я искал и читал о стереографических проекциях, но я не смог найти ничего подобного. На самом деле, я не видел никаких проекций, использующих эту «нулевую координату искажения».
источник