Проблема в том, что разные картографические проекции и разные эллипсоиды! (Не говоря о искажениях как СК42, так и WGS84). Естественно, координаты в СК42 сильно отличаются от координат в WGS84.
Пересчёт из одной системы в другую проводится так (например, СК42 в WGS84).
1. Координаты СК42, если они даны в формате x,y, приводятся к формату B,L (геодезические широта и долгота). Если необходима очень высокая точность пересчёта, учитывается высота (нормальная высота на европейской части РФ достаточно для наших целей близка к эллипсоидальной).
2. B,L, (H) СК42 пересчитываются в B,L WGS84, высота в WGS84 равна высоте в СК42 - нормальная высота.
3. B,L (WGS84) пересчитываем в x,y (WGS84), только масштабный коэффициент на центральном меридиане зоны для UTM равен 0,9996, а не 1, как в проекции Гаусса-Крюгера.
Пример (возьмём точку на Красной площади в Москве)
Дано (проекция Гаусса-Крюгера).
x42=6182000,0000
y42=7413600,0000
Найти:
x84,y84 в проекции UTM
Решение.
1. Найдём B42, L42 (геод. к-ты на элл. Красовского)
B42 = 55° 45' 09",890886
L42 = 37° 37' 26",293904
2. Найдём B84, L84 (эллипсоидальная высота принимается за 150 м)
B84 = 55° 45' 10",022431
L84 = 37° 37' 19",550035
Для пересчёта СК42 => WGS84 применены параметры 7-параметрического преобразования (Хелмерта)
DX=23,92м DY=-141,27м DZ=-80,90м WX=0" WY=-0,35" WZ=-0,82" M=-0,12*10^-6
Это условные параметры, взятые из открытых источников, на самом деле они могут значительно отличаться от приведённых. Примем также, что локальные искажения систем координат нулевые.
Если используются правильные параметры пересчёта, то точность конечного результата (x84, y84) не хуже 0,3 мм (каждый этап даёт ошибку не более 0,1 мм).
3. Координаты в UTM на эллипсоиде WGS84
x84 = 6179425,3349
y84 = _413518,5090 (номер зоны, если считать от 180° - 37)
Пересчёт B,L <=> x,y проводился с помощью 4-членных рядов Тейлора.
Voila!
P.S. Товарищ полковник, никаких секретов в таком пересчёте нет, все данные и алгоритмы взяты из открытых источников
