Вопрос по "геодезическому" коду на С

Имеется следующий код:

Code:

/// Пребразование геодезических координат(fiUT,laUT,H_UT)(радианы) в местные прямоугольные XZ (метры) с центром в т.с. Системы (fi,la,H)(радианы)
 
void IzGEOwXZ(double fi, double la, double , double fiUT, double laUT, double H_UT, double *x, double *y)
{
double R1WrtUT,R1Wrt,sinfi,sinfiUT,cosfi,cosfiUT,
Dlt_la,sinDlt_la,cosDlt_la,K2;
 
R1Wrt=R1Wert(fi); R1WrtUT=R1Wert(fiUT);
sinfi=sin(fi); sinfiUT=sin(fiUT); cosfi=cos(fi); cosfiUT=cos(fiUT);
Dlt_la=laUT-la; sinDlt_la=sin(Dlt_la); cosDlt_la=cos(Dlt_la);
 
K2=cosfi*sinfiUT-sinfi*cosfiUT*cosDlt_la;
 
*x=(R1WrtUT+H_UT)*K2+e2()*(R1Wrt*sinfi-R1WrtUT*sinfiUT)*cosfi;
*y=(R1WrtUT+H_UT)*cosfiUT*sinDlt_la;
}//IzGEOwXZ

При этом R1Wert() возвращает радиус кривизны первого вертикала, e2() возвращает квадрат эксцентриситета меридианного эллипса. Это все на эллипсоиде Крассовского. Соответственно, в функцию передается широта, долгота точки принятой за начало отсчета (x=0, y=0) и широта, долгота, высота некой другой точки. На выходе получаются координаты в метрах относительно начала отсчета.

Вопрос дурацкий, наверное и нечеткий, но все же : объясните, что стоит за этими формулами, это проекция какая-то картографическая, где можно про это почитать поподробнее?

 

Напоминает экваториальную проекцию (вроде не ошибся с названием) с "хитрым" центром в некой "центральной" точке. Но это ИМХО, догадки так сказать.

 

Show Spoiler Hide Spoiler

Кровь из глаз, я чуть глаза не сломал, читая назавания функций и переменных. Разобрать в этой каше ничего невозможно, для начала надо деобвфусцировать код Может после этого станет понятнее.