Nazar написав:Вообще побороли, но по науке было бы лучше. Если что накопаешь пиши, будем пробовать.
...
Всего, основных, есть десятка ~два проекций, вычитывать где какая и чем отличается - даже не столько лень, сколько времени нет.
Судя по общим описаниям лучше всего использовать координаты Гаусса-Крюгера.
Эдакая компиляция:
=== Cut ===
Алгоритм вычисления координат Гаусса-Крюгера из геодезических координат (ведется внутри 6-ти градусной зоны):
X=S+ (L^2*r*sin(B))/2.+ (L^4*r*cos(B)^2*sin(B))*(5-tan(B)^2)/24; (1)
Y=L*r+ L^3*r*cos(B)^2*(1-tan(B)^2)/6 (2)
где:
B- параллель точки в радианах.
L- отклонение от центрального меридиана ( в радианах).
S=6367558.5*B-16036.5*sin(2*B)+16.8*sin(4*B)
r=Re*cos(B)/sqrt(1-E2*sin(B)*sin(Bp)) - радиус кривизны параллели.
Re=6378245.0 - экваториальный радиус эллипсоида Красовского.
Rp=6356863.0 - полярный радиус эллипсоида Красовского.
E2==0.00669343
Расч тные и справочные значения совпадают с точностью до 0.25 метра.
=== Cut ===
Есть некие сомнения насчет вычисления "r" - там где "sin(B)*sin(Bp)" - либо "sin(B)*sin(B)", либо "sin(B)*sin(Rp)" - без вникания в суть определить точно не получится - нужно или попробовать посчитать или разобраться что считается. Если не обломаешься посчитать - хорошо, т.к. вникать лень.