В том что многие пункты находятся за пределами Российской Федерации и для них известны опубликованные современные геоцентрические координаты, что позволяет вычислить 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 --- Есть вариант исключить точки в Финляндии, так как они нежестко связаны с южной частью. По открытым финским данным сохранились еще несколько точек, но дьявол скрывается в центрировках, к тому же если их добавить, то вес финских пунктов станет еще выше. Также можно добавить Пулково, чтобы удалиться от меридиональной дуги.
А чего сразу исключать то? Я противник всяческих "исключений" без весомых доводов. Другое дело вес резануть "сомнительному" пункту. Только в приведённом мной примере я вообще ни шиша о пунктах не знаю, а резанул веса чисто из взаимной сходимости. Это имеет определённую плохость.
