По привычке для широты и долготы "proj4 string" Заодно тогда и Код: #define RAD2ARC 180/M_PI*60*60 /* number of arcseconds in 1 radian */ --- Сообщения объединены, 26 июн 2022, Оригинальное время сообщения: 26 июн 2022 --- proj4 всегда крутил в анти-гостовскую (анти-мапинфошную и также анти-RTCMовскую) сторону, также как и TOWGS84[] у ESRI WKT. Классическое определение Код: +towgs84=-96.063,-82.428,-121.751,4.801,0.345,-1.376,1.496 соответствует новомодной нотации Код: +proj=pipeline +step +proj=helmert +x=-96.063 +y=-82.428 +z=-121.751 +rx=4.801 +ry=0.345 +rz=-1.376 +s=1.496 +convention=position_vector
Немного дополню, у ESRI тоже есть выбор - или Coordinate Frame или Position Vector, тут просто от пользователя больше зависит. Так-то там и NTv2 и преобразования по Молоденскому...
@cr2 а не подсобишь? Какие радиусы эллипсоида ты использовал и как поступил с высотами? (не могу добиться твоего результата) PS: GIT: Geo-Linux-Calculations/helmparms3d : 1.0.6: Proj string output.
@cr2 . Разобрался я с этими финскими эллипсоидами: helmert3d-kkj-etrs89.xls: Proj string: +towgs84=-96.063,-82.429,-121.753,4.801,0.345,-1.376,1.496 Error: mX = 2.69 mY= 3.25 mZ= 2.05 Но решил посмотреть, что даст Wpress: helmert3d-kkj-etrs89-wpress.xls: Proj string: +towgs84=-94.507,-82.487,-120.096,4.836,0.347,-1.266,1.161 Error: mX = 2.39 mY = 2.67 mZ = 1.72 Всё "это" без остатка с помощью: Таблицы MS Excel для пересчёта координат по общим точкам.
Понедельник день тяжелый. Для EUREF-FIN брал +proj=latlong +ellps=GRS80 +to +proj=geocent, для классики +proj=latlong +ellps=intl +to +proj=geocent и высоту H_nat+N. --- Сообщения объединены, 27 июн 2022, Оригинальное время сообщения: 27 июн 2022 --- Файлы
Да уже разобрался (см. комент #225). Всё путём. ) PS: Горячку спорол. Не те высоты в файле для KKJ сохранил. Но в геоцентрике верные.
Имеются Таблицы прямоугольных координат Гаусса-Крюгера рамок трапеций, изд.ВТУ ГШ Красной армии на сфероиде Бесселя (Москва 1941 год) - для определения углов рамок трапеций + Таблицы с поправками перехода к с/к Пулк.-1942.
Местные системы координат старинных городов для создания планов, выполнялись в астрономических (геодезических) координатах на эллипсоиде Бесселя по астрономическому (геодезическому) азимуту, но не дирекционному углу , только после создания и принятия проекции, были получены местные прямоугольные координаты, при этом за начало координат принималась доминанта желательно в центре города. Проверка показала, что сеть местных координат старинных городов ориентирована по азимуту и расхождение с дирекционным углом равно сближению меридианов. --- Сообщения объединены, 3 авг 2022, Оригинальное время сообщения: 3 авг 2022 --- Поправки для перехода от Пулковской СК к системе СК 42 даны в таблицах для трапеций 1: 1 000 000 масштаба с точностью 2-3 м.
Vitali 43, ну вот у нас городу скоро 250, генплан с 180* года и СК города сориентирована по направлению основных улиц
Вспомнил молодость и забил данные СК-1932 по ИГД, базисам и (синтезированным) углам треугольников схемы сети в перфокарты для программы NGS ADJUST 6.4.2 Названия пунктов даны в КОИ-7, фортран с UTF-8 работать толком не умеет. На 149 странице распечатки АЦПУ вылетела ошибка "THE ADJUSTMENT IS SINGULAR", с 8 пунктами (из 703) что-то не в порядке: Код: 1 NATIONAL GEODETIC SURVEY PROGRAM ADJUST ADJUSTMENT PROGRAM PAGE 149 VERSION 6.4.2 0*********** COMMENCING ADJUSTMENT *********** THE AVERAGE BAND WIDTH FOR THE 2 DIM. ADJUSTMENT OF 703 STATIONS AND RANK 2109 IS 3.8%. D.P. WORDS NEEDED= 37783 0THE FOLLOWING 11 UNKNOWNS FALL BELOW TOLERANCE OF 1.00D-09 1A 373 %+bLAGODATNOE* (hWALYNSKIJ) EAST/WEST SHIFT GOOGE NUMBER IS 9.20D-10 2B 360 +*oHOTNI^IJ NORTH/SOUTH SHIFT GOOGE NUMBER IS 8.79D-10 3C 360 +*oHOTNI^IJ EAST/WEST SHIFT GOOGE NUMBER IS 4.15D-10 4D 109 +kALINOWKA EAST/WEST SHIFT GOOGE NUMBER IS 1.81D-10 5E 310 +*`w KONEC BAZISA NORTH/SOUTH SHIFT GOOGE NUMBER IS 1.73D-10 6F 310 +*`w KONEC BAZISA EAST/WEST SHIFT GOOGE NUMBER IS 1.00D-10 7G 113 +*sOKOLXNIKOWO NORTH/SOUTH SHIFT GOOGE NUMBER IS 2.87D-10 8H 465 dOLMATOWO EAST/WEST SHIFT GOOGE NUMBER IS 1.83D-10 9I 314 rQZANOWKA NORTH/SOUTH SHIFT GOOGE NUMBER IS 0.00D+00 10J 316 lEMVA NORTH/SOUTH SHIFT GOOGE NUMBER IS 0.00D+00 11K 316 lEMVA EAST/WEST SHIFT GOOGE NUMBER IS 1.66D-16 0THE ADJUSTMENT IS SINGULAR !! ******************************
Сам же снабдил версией 6.4.3: GNU Gama #23, а работаешь в 6.4.2? Подозрительно! GIT: adjust : 6.4.3-p
У меня в исходниках для 6.4.2 были заменены вычисления заокеанских долгот "западной долготы" на человеческие, поэтому ей и пользуюсь для отладки.
Исправил логическую ошибку в углах - и новогодний подарок Г. Н. Ефимову готов: метод развертывания, как его описывают в книгах заработал. Для метода проектирования надо добавить данные о высотах геоида из EGM2008 и высотах пунктов из какой-нибудь современной DEM. Остается технический вопрос о редукции длин в ADJUST и более принципиальный вопрос об азимутах Лапласа и среднеквадратичной минимизации УОЛ. По моему пониманию, в ADJUST азимутом Лапласа считается геодезический азимут вычисленный из астрономического с помощью априорно заданной модели УОЛ что не совсем соответствует идеям классической геодезии.
Буду искать, должны где-то лежать. Благодаря появлению "новых" отсканированных старых книг выявляются доселе совсем неизвестные исходники СК32, например гравиметрические данные для пунктов Саблино и Бугры в дополнение к их астрономическим данным.
Пока загвоздка в следующей проблеме курицы и яйца : для DEM (SRTM и иже с ними) нужно вычислить приближенные координаты WGS84, а чтобы их вычислить, надо перередуцировать длины базисов с эллипсоида Бесселя обратно на хорды, часто без высот. Каталог 1935 года не зря же назывался "предварительным". Астрономические данные я полностью ингнорирую, так как они (полноценные) есть только для Пулково и Саблинского базиса, а остальные очень фрагментарные (односторонние Ульяновск и Саратов из дополнительного каталога ГУГСК НКВД для Сибири, либо известны только астрономическая широта или только астрономическая долгота).
