| Обсуждение |
| | пересчет координат из WGS84 в СК42 [ VAK Крюков ] Четверг, 31 августа 2006, 11:00
Решил опубликовать на уважаемом мной форуме один результат,
который может представлять интерес не только для меня...
В 2001 вышел очень нужный ГОСТ 51794
"Аппаратура радионавигационная глобальной навигационной спутниковой системы и глобальной системы позиционирования.
СИСТЕМЫ КООРДИНАТ.
Методы преобразований координат определяемый точек".
Таким образом, на мой взгляд:
1. Знания стали овладевать массами.
2. Государство устанавливало очень важное единообразие в области пересчета координат.
Пафосно говоря, фактически заложен геодезический камень в фундамент инфраструктуры пространственных данных страны.
В Интернете всплеск интереса к Госту наблюдался в 2003-2004гг., затем он сошел на нет, хотя некоторые вопросы по-моему остались открытыми.
Жизнь сложилась так, что по работе мне пришлось снимать точки на КПК Acer n35
и в несколько приемов рутинообразно конвертировать точки из формата в формат, из геодезических координат WGS84 в плоские координаты CR42.
Решил написать программку для автоматизации этого процесса.
Пришлось более детально изучить ГОСТ 51794.
Для перехода от "датума" WGS84 к эллипсоиду Красовского по материалам ГОСТ-а я перемножил две матрицы поворота
и матрицу поворота из ПЗ90 в референцную СК42 на вектор линейного смещения эллипсоида
(естественно, в однородных координатах будет только одно перемножение матриц),
затем привел знаки параметров как в формуле (20) и
получил следующие параметры трансформации из WGS84 в СК42:
Wx=0; Wy=1.8e-6; Wz=4.12e-6
dX=-23.9; dY=141.3; dZ=80.9
Затем я написал программку на c#, реализующей ГОСТ-ие формулы пересчета и
сравнил результаты с пересчетом в MapInfo: расхождение по углам после 7-ого знака,
в плоских координатах до 15 метров.
Найденные параметры поворота нигде не публиковались, поэтому решил попросить проверки у сообщества.
Может, в мои раскладки вкралась ошибка, ведь по образованию я не геодезист, просто немного помню институтский курс математики.
ВАК |
|
| | [ VAK ] Вторник, 19 сентября 2006, 10:55
Кто-нибудь может проверить приведенные параметры на каталожных координатах пунктов ГГС?
Не забудьте, что m = (– 0,12±0,06)•10e–6. |
|
ОтветитьЗнаком «*» отмечены обязательные для заполнения поля. |
| |
