Параметры (датум) для Системы координат 1932 года (СК-32)

По привычке для широты и долготы

"proj4 string"
Заодно тогда и

Code:

#define RAD2ARC 180/M_PI*60*60 /* number of arcseconds in 1 radian */

--- Сообщения объединены, Jun 26, 2022, Оригинальное время сообщения: Jun 26, 2022 ---

proj4 всегда крутил в анти-гостовскую (анти-мапинфошную и также анти-RTCMовскую) сторону,
также как и TOWGS84[] у ESRI WKT.

Классическое определение

Code:

+towgs84=-96.063,-82.428,-121.751,4.801,0.345,-1.376,1.496

соответствует новомодной нотации

Code:

+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

 

Нда. "Не верь глазам своим". Cast @ErnieBoyd

 

Немного дополню, у 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.

--- Сообщения объединены, Jun 27, 2022, Оригинальное время сообщения: Jun 27, 2022 ---

Файлы

 

Да уже разобрался (см. комент #225). Всё путём. )

PS: Горячку спорол. Не те высоты в файле для KKJ сохранил. Но в геоцентрике верные.

 

Имеются Таблицы прямоугольных координат Гаусса-Крюгера рамок трапеций, изд.ВТУ ГШ Красной армии на сфероиде Бесселя (Москва 1941 год) - для определения углов рамок трапеций + Таблицы с поправками перехода к с/к Пулк.-1942.

 

Не поможет ничем. "Это" самое дно качества.

 

Местные системы координат старинных городов для создания планов, выполнялись в астрономических (геодезических) координатах на эллипсоиде Бесселя по астрономическому (геодезическому) азимуту, но не дирекционному углу , только после создания и принятия проекции, были получены местные прямоугольные координаты, при этом за начало координат принималась доминанта желательно в центре города. Проверка показала, что сеть местных координат старинных городов ориентирована по азимуту и расхождение с дирекционным углом равно сближению меридианов.

--- Сообщения объединены, Aug 3, 2022, Оригинальное время сообщения: Aug 3, 2022 ---

Поправки для перехода от Пулковской СК к системе СК 42 даны в таблицах для трапеций 1: 1 000 000 масштаба с точностью 2-3 м.

 

Вспомнил молодость и забил данные СК-1932 по ИГД, базисам и (синтезированным) углам треугольников схемы сети
в перфокарты для программы NGS ADJUST 6.4.2
Названия пунктов даны в КОИ-7, фортран с UTF-8 работать толком не умеет.
На 149 странице распечатки АЦПУ вылетела ошибка "THE ADJUSTMENT IS SINGULAR",
с 8 пунктами (из 703) что-то не в порядке:

Code:

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 года не зря же назывался "предварительным".
Астрономические данные я полностью ингнорирую, так как они (полноценные) есть только для Пулково
и Саблинского базиса, а остальные очень фрагментарные (односторонние Ульяновск и Саратов
из дополнительного каталога ГУГСК НКВД для Сибири,
либо известны только астрономическая широта или только астрономическая долгота).