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

В том что многие пункты находятся за пределами Российской Федерации и для них известны
опубликованные современные геоцентрические координаты, что позволяет вычислить 7 параметров
и невязки этого преобразования. Где вы видели такие данные для СК-32 кроме Пулково и вокруг,
например Саратов или Ульяновск, поэтому обсуждение плавно и перешло на Шарнгорста.
Единственный такой пункт для СК-32 это Лапин, где известны координаты ГСК-2011,
хотя тоже не совсем ITRF, но для тестов подойдет.

 

Если каталог Шарнгорста не имеет никакого отношения ни к СК-10 и тем более к СК-32, то о каких невязках идёт речь? Ну хорошо, известны координаты неких пунктов, пусть даже всех, из Шарнгорста, также известны их координаты ITRF. Вычислили некие параметры... Если каталог не имеет никакого отношения к СК-10 и СК-32, то толку то от наших расчётов? Для перехода из чего во что мы расчитали некие параметры? Если он (каталог) не имел никакой практической реализации, в отличии от каталога 1935 года и 1937 года НКВД для СК-32?

Или вы имеете ввиду как-то поиграться с математикой и формулами, из-за неимения полных каталогов СК-32, и за счёт каталога Шарнгорста, сделать некое подобие перехода для СК-32 в WGS-84? Это и есть ваши хитрости?

А нельзя ли поступить так, как вы сами рассказали

что в СК-32 сеть восточнее Ульяновска/Саратова была "прицеплена" к полигонам Европейской части великого могучего? Зная каталожные координаты 42 пунктов из куска каталога, и зная координаты на Бесселе пунктов нужной области (из тех же карт) и зная их же в WGS 84, не получим ли мы СК-32 для нужной области?

Для тестов подойдут любые из 43 известных пунктов
Quality criteria
3.1326 Maximum spatial residual [m]
1.5776 Average spatial residual [m]
1.0103 Root Mean Square residual (RMS) [m]

 

Давайте по порядку: 7 параметров можно вычислить для любой пары пунктов и невязка (для них) будет равна нулю.
А вот при применении этих 7 параметров к любым другим пунктам там возникнут ненулевые невязки.
Если вычислить средние по больнице 7 параметров для всех известных пунктов
то невязки будут везде !
Я оцифровал все 703 пункта восьми полигонов СК-32 (без Уральского полигона и ручки к нему,
если кто-то закажет их скан в РГБ, то можем добавить и их).

 

Из 703 пунктов 1 класса СК-32 как минимум 150 находятся в Украине,
по крайней мере именно столько я нашел у них на карте.
Из них 23 с высотой "до 1934 года", ее сравнение с урезанной до метра БС1977 года тоже весьма поучительно:

Код:

 80 6395504.00 5544553.00 115.00 Войниця          |50  1 19.359 31 32 28.070 |+115.45
 81 6410930.00 5539304.00 141.00 Пологи           |49 58 38.509 31 45 27.223 |+139.59
 82 6385163.00 5535114.00 242.00 Григорівка       |49 56  7.080 31 23 58.819 |+242.78
 83 6404053.00 5533747.00 124.00 Хоцьки           |49 55 34.829 31 39 47.112 |+124.60
 84 6393779.00 5510819.00 225.00 Пекарі           |49 43  6.651 31 31 34.677 |+225.35
 85 6429948.00 5514889.00 136.00 Сеньківці        |49 45 37.432 32  1 37.584 |+135.26
 87 6452071.00 5485808.00 137.00 Іркліїв          |49 30  3.916 32 20 15.795 |+137.97
 89 6398339.00 5455508.00 204.00 Орловець         |49 13 19.440 31 36 13.431 |+204.3
 90 6413284.00 5437805.00 189.00 Сердюківка       |49  3 54.714 31 48 45.627 |+189.4
 92 6414775.00 5411173.00 226.00 Защита           |48 49 33.440 31 50 19.112 |+226.6
 94 6417052.00 5378663.00 239.00 Велика Виска     |48 33 44.752 31 55 44.804 |+232.83
 96 6452121.00 5355940.00 181.00 Губівка          |48 19 59.623 32 21 13.066 |+181.04
 99 6457866.00 5329294.00 167.00 Богодарівка      |48  5 38.377 32 26  1.580 |+167.47
108 6454429.00 5217306.00  62.00 Бармашове        |47  5 11.217 32 23 57.294 |+62.56
109 6443810.00 5213595.00  63.00 Калинівка        |47  3  8.076 32 15 35.495 |+63.16
553 6553527.00 5767600.00 208.00 Олександрівський |52  2  1.617 33 46 46.594 |+208.06
554 6548404.00 5739514.00 209.00 Собичев          |51 46 54.552 33 42  3.675 |+208.69 
185 7305235.00 5580004.00 218.00 Журавлівка       |50 19  4.008 36 15 54.242 |+215.76
276 7574760.00 5480685.00 217.00 Зоринівка        |49 27  8.222 40  1 53.022 |+217.34
279 7556600.00 5449470.00 189.00 Бараниківка      |49 10 24.835 39 46 35.412 |+188.82
286 7556677.00 5301634.00 249.00 Несвітай         |47 50 38.653 39 45 26.901 |+248.76
624 7464208.00 5321504.00 282.00 Ольховчик        |48  1 27.454 38 31 12.983 |+281.87
625 7472900.00 5314059.00 193.00 Мануйлівка       |47 57 27.919 38 38 14.248 |+191.76

 

У меня такой вопрос по таблице :
можно ли ввести указанные поправки арифметически в углы рамок карт, изданных до с/к 1942г. на эл.Бесселя.
Если по Таблицам Каврайского- 1941г. ВТУ ГШ КА определить прямоугольные координаты углов рамок трапеций топографической карты в системе координат "Пулково-32" на эл.Бесселя, а затем ввести соответствующие поправки.
Технически будет правильно, если не претендовать на исключительную точность?
Цель - привязка карты в с/к 1942г.

 

Найдите соответствующие координаты в WGS84 и привязывайте с помощью чего-нибудь типа:
* CPlaneLite — программа пересчёта координат по общим точкам
* GeoFindKey
Также, зная параметры конформного преобразования, можно и сами растры трансформировать:
* GeoConformImage
* PhotoQuick (сообщение #2)

 

А откуда у вас такой вопрос? Какие виды работ вы выполняете? Это госконтракт?

 

Это какая такая исключительная точность? Можете в цифрах написать? И для каких трапеций - какой масштаб?

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

Госконтракт не будет выполняться в СК-32. Вот грант - может быть.

 

СК-42 здесь вообще не при чем, для вычисления 7 параметров СК-32 в WGS84 достаточно иметь
2 точки в этих системах, я и предлагал взять в качестве примера "Пулково Сигнал А" и "Лапин".
Вычисления можно начать с вектора/хорды "на геоиде".
"Лапин" кстати есть и в каталоге Шарнгорста (под названием "Лапина", измеренный еще
при государе-императоре Николае I, но скорее всего он не совпадает физически с "Лапин" из СК-32).

Добавка к таблице Г.Н. Ефимова от Шарнгорста (долготы отъ Пулкова, для совместимости с СК-32
надо прибавить 30°19'38.55")

Код:

637       Окулова                                  55 52 22.040       +5 45 57.743
680       Лапина                                   55 42 47.931       +6 06 59.572

 

Чтобы продемонстрировать насколько это хитроумное дело
искать 7 параметров для классической АГС
возьмем пример дуги Струве из каталога Шарнгорста.
Разности геоцентрических координат (3 параметра) для известных 16 точек

Код:

2960 265 Старо-Некрасовка                 453.003 191.651 466.055 
1565 188 Белин                            456.349 211.207 458.941 
1528 186 Лесковичи                        455.888 210.691 460.096 
1476 185 Осовница                         455.951 210.154 460.702 
1510 -   Чекуцк                           455.981 210.591 460.182 
 894 157 Немеж                            455.689 203.309 468.092 
 840 155 Мешканцы                         455.603 203.216 469.271 
 899 159 Березняки                        455.568 203.517 468.088 
 853 158 Хорунжишки                       455.514 202.705 469.405 
 495 137 Якобштадт                        454.143 201.392 477.422 
 179 119 Юрьев,центр башни обсерватории   457.733 203.147 481.832 
 102 108 Хальяль, (Галльял)               459.489 203.697 485.481 
  72 106 Мекипелюс                        458.705 203.746 489.732 
  71 104 Свартвира                        465.162 204.826 486.518 
  39  70 Килпимяки                        464.462 207.475 490.181 
   2  36 Торнео                           470.916 206.133 497.829 

Вот какие средние по больнице 7 параметров получаются у helmparms3d
(надо только радианы в угловые секунды пересчитать)

Код:

0.9999999994 0.0000144493 0.0000311764 
-0.0000144489 0.9999999998 -0.0000127294 
-0.0000311766 0.0000127289 0.9999999994 
258.2462892763 311.1197073322 531.6367574653 
1.0000042997

Осталось только вычислить получившиеся невязки.
Исходные данные кого-нибудь интересуют ?

 

Так я же добавил helmdiff3d в пакет! Ты какую версию то пользуешь и откуда?

См. также https://github.com/dr-ni/helmert3d/issues/9

PS: Ежели нужен GUI, так вообще таблицами пользуйтесь: https://geodesist.ru/threads/osi-zdanija-ot-suschestvujuschego-gro.88950/page-3#post-1028887

 

Да, интересно

 

Координаты ETRS89 и высоты взяты из 6 выпуска "Вестника Санкт-Петербургского общества геодезии и картографии"
http://www.spbogik.ru/vestnik/80-6-2007.html
(страница 179 и текст, за исключением пункта Мякипялюс, его пришлось вычислять самому
триангуляцией по известным исправленным расстояниям
до Svartvira и Halljall).
Данные из каталога Шарнгорста в таблице (первая колонка номер в каталоге, вторая номер в дуге Струве):

Код:

   2  36 + Торнео                           65 49 54.4825      -6 10 12.6598 
  39  70 + Килпимяки                        62 38 08.0427      -3 33 32.7008 
  71 104 + Свартвира                        60 16 39.9643      -3 43 24.4265 
  72 106 + Остров Гохланд Мекипелюс         60 04 31.124       -3 21 25.643  
  79       Пулково, сигнал А                59 46 19.667       -0 00 13.772  
  87       Кабози S, пирамида               59 43 25.012       -0 02 24.051  
 102 108 + Хальяль, (Галльял)               59 25 57.600       -4 03 31.704  
 179 119 + Юрьев,центр башни обсерватории   58 22 47.560       -3 36 24.709  
 495 137 + Якобштадт                        56 30 07.739       -4 28 13.970  
 840 155 + Мешканцы                         54 55 53.846       -5 00 37.610  
 853 158 + Хорунжишки                       54 51 45.289       -4 41 53.365  
 894 157 + Немеж                            54 39 03.559       -5 00 37.492  
 899 159 + Березняки                        54 38 06.292       -4 53 53.229  
1476 185 + Осовница                         52 17 23.983       -4 40 39.355  
1510       Чекуцк                           52 12 28.719       -4 46 16.926  
1528 186 + Лесковичи                        52 09 40.196       -4 45 20.669  
1565 188 + Белин                            52 02 40.677       -5 06 35.214  
2960 265 + Старо-Некрасовка                 45 19 55.633       -1 23 55.929

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

Версия самая старая оригинальная, но результаты такие же.
Если умножить внедиагональные элементы матрицы поворота на 206264.8062
и поменять знаки поворота, получится

Код:

-rz=0.0000144493*206264.8062=2.98038
+ry=0.0000311764*206264.8062=6.43059
-rx=-0.0000127294*206264.8062=-2.62563

то есть +towgs84=258.246,311.12,531.64,2.62563,6.43059,-2.98038,4.2997,
невязки этого окончательного преобразования меня и интересуют
(есть то, что будет например в qgis).
Это все только проверка вычислительной технологии,
надо еще независимо проверять высоты геоида, пересчеты в геоцентрику и так далее.

Для Мякипелюс опубликованы такие координаты WGS84

Код:

60d04'27"N     26d58'11"E (страница 179)
60d04'26.8"N   26d58'11.2"E (координаты для навигатора из статьи про дугу Струве)

но это слишком грубо, поэтому я взял длины линий (+редукцию с геоида) и поправки со страницы 180,
и получил 60d04'26.843"E 26d58'11.197"N.
Раз нет ринекса, приходится использовать то, что есть.

 

Имеешь в виду разночтения в значениях осевых радиусов эллипсоидов или что?

 

cr2, опишу свою попытку вычисления 7 параметров.

В 6 выпуске «Вестника…» в таблице на с. 179 нашёл 10 пунктов из Вашей таблицы. Вычислил геоцентрические координаты от эллипсоида GRS 80. Эллипсоидальные высоты взял нулевые.

Код:

   2  36 Торнео                          2389026.657 1071460.181 5795689.820
  39  70 Килпимяки                       2624169.422 1323665.073 5640836.905
  71 104 Свартвира                       2834291.980 1419468.630 5515315.253
  72 106 Остров Гохланд Мекипелюс        2842597.640 1446434.703 5504098.507
 102 108 Хальяль, (Галльял)              2915554.203 1438883.390 5468039.118
 179 119 Юрьев, центр башни обсерватории 2993803.151 1506982.629 5407482.294
 495 137 Якобштадт                       3174762.291 1538536.142 5294959.995
 894 157 Немеж                           3342773.120 1581238.002 5178471.046
1565 188 Белин                           3556006.455 1674563.249 5005348.187
2960 265 Старо-Некрасовка                3930569.099 2172224.468 4512916.517

По данным Вашей таблицы для этих же пунктов вычислил долготы от Гринвича, приняв долготу Пулковского меридиана равной 30°19′33.79″, и вычислил геоцентрические координаты от эллипсоида Бесселя.

Код:

   2  36 Торнео                          2389480.271 1071724.776 5796205.784
  39  70 Килпимяки                       2624610.947 1323936.944 5641343.553
  71 104 Свартвира                       2834731.272 1419742.313 5515815.266
  72 106 Остров Гохланд Мекипелюс        2843029.828 1446707.544 5504601.639
 102 108 Хальяль, (Галльял)              2915988.323 1439158.237 5468539.396
 179 119 Юрьев, центр башни обсерватории 2994235.063 1507259.373 5407980.412
 495 137 Якобштадт                       3175191.606 1538815.973 5295455.346
 894 157 Немеж                           3343205.519 1581524.697 5178959.725
1565 188 Белин                           3556439.411 1674863.699 5005828.745
2960 265 Старо-Некрасовка                3930990.904 2172517.512 4513404.556

Вычислил 7 параметров.

Код:

∆X  368.318  ±  64.496 m
∆Y  216.109  ± 110.027 m
∆Z  512.306  ±  16.706 m
ωx  -0.02567 ± 2.92703 ″
ωy  -2.98756 ± 1.83240 ″
ωz  -3.74857 ± 2.37496 ″
s    5.5342  ± 2.5068 ppm

Строка +towgs84:
+towgs84=368.318,216.109,512.306,0.0256675,2.98756,3.74857,5.53422

Вычислил невязки.

Код:

   2  36 Торнео                          -7.602   0.139  -6.053   9.718
  39  70 Килпимяки                       -1.038  -1.448  -0.969   2.028
  71 104 Свартвира                       -1.201   1.104   1.940   2.535
  72 106 Остров Гохланд Мекипелюс         5.296   2.248  -1.359   5.912
 102 108 Хальяль, (Галльял)               3.383   1.531   0.238   3.721
 179 119 Юрьев, центр башни обсерватории  3.909   1.440   0.938   4.271
 495 137 Якобштадт                        5.304   1.830   0.467   5.631
 894 157 Немеж                            0.688  -1.730   4.065   4.471
1565 188 Белин                           -2.893 -11.071   8.151  14.050
2960 265 Старо-Некрасовка                -5.846   5.958  -7.418  11.167
 

Невязки меньше, чем в решении по трём параметрам. Однако уже по расположению пунктов вдоль линии ясно, что система плохо обусловлена. Корреляционная матрица это подтверждает.

Код:

 1.000  -0.958  -0.284   0.953   0.988  -0.958  -0.119
-0.958   1.000   0.289  -0.995  -0.949   0.981  -0.035
-0.284   0.289   1.000  -0.247  -0.411   0.285  -0.800
 0.953  -0.995  -0.247   1.000   0.944  -0.960   0.000
 0.988  -0.949  -0.411   0.944   1.000  -0.936   0.000
-0.958   0.981   0.285  -0.960  -0.936   1.000   0.000
-0.119  -0.035  -0.800   0.000   0.000   0.000   1.000

 

Параметры эллипсоидов проверять не стал, но определённую плохость обусловленности (cond(G)=109,8) и сами параметры подтверждаю.

PS: Только у меня почему то знак вращений противоположный.

 

Судя по формулам в ячейках $DATA.Q3:S22, Вы используете метод "Position Vector transformation", как PROJ.4 +towgs84.
Я же использую метод "Coordinate Frame rotation", как ГОСТ или MapInfo.

 

Принято и учтено.
Применив "научно-исследовательский подход" (метод тыка) и "поколдовав" с весами, удалось слегка снизить плохость обусловленности (cond(G)=93,9), при этом "слегка" "пожертвовав" последними двумя пунктами.

 

Почему именно это значение, а не 30°19'38.55"E как в Пулковской системе 1910 года
(или 30°19'38.558"E как в СК-32)?
Я причешу свои таблицы и тогда можно будет сравнить результаты с учетом геоида.

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

Есть вариант исключить точки в Финляндии, так как они нежестко связаны с южной частью.
По открытым финским данным сохранились еще несколько точек, но дьявол скрывается в центрировках,
к тому же если их добавить, то вес финских пунктов станет еще выше.
Также можно добавить Пулково, чтобы удалиться от меридиональной дуги.

 

А чего сразу исключать то? Я противник всяческих "исключений" без весомых доводов. Другое дело вес резануть "сомнительному" пункту. Только в приведённом мной примере я вообще ни шиша о пунктах не знаю, а резанул веса чисто из взаимной сходимости. Это имеет определённую плохость.