Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Есть. Называются "определение параметров локальной системы координат". На сайте компании Ракурс есть программа, которая именно это и делает. Причем бесплатная.Skywalker писал(а):Имеется не очень точно-привязанное графическое приложение в системе координат Пулково-42, 9 зона, с него оцифрованы точки, для них вычислены прямоугольные координаты в СК Пулково-42, 9 зона. Также имеется приложение с координатами для этих точек в неизвестной СК. Есть ли какие-нибудь методы для определения неизвестной СК?
Не поможет. Ибо исходные координаты такой кривизны, а территория настолько мала, что ничего не поможет.Boris писал(а):Есть. Называются "определение параметров локальной системы координат". На сайте компании Ракурс есть программа, которая именно это и делает. Причем бесплатная.
Если поменять местами местные X и Y, невязки аффинного преобразования не кажутся страшными:stout писал(а):...не очень точное для этих целей простое аффинное преобразование даёт всего лишь метровые невязки. А тут невязки в десятки и сотни метров.
Код: Выделить всё
2026 -2.718 0.259
2027 -6.934 5.174
2028 5.965 5.284
2029 2.392 -9.281
2031 5.121 10.540
2032 -2.956 -5.981
2033 -6.385 -0.278
2226 5.516 -5.717
Код: Выделить всё
1.1684850297031111, -0.050508779240812984, 9326030.6515536346, 0.0049954036126287557, 0.96381045319532099, 6576472.9219575021Эта программа непригодна для решения подобных задач: она не выдаст результата, если координаты хотя бы одной точки не согласуются с прочими с геодезической точностью.Boris писал(а):На сайте компании Ракурс есть программа, которая именно это и делает. Причем бесплатная.
Код: Выделить всё
1.1684850297031111, -0.050508779240812984, 9326030.6515536346, 0.0049954036126287557, 0.96381045319532099, 6576472.9219575021Эти параметры можно использовать в MapInfo, в ArcGIS они, по-моему, бесполезны. Нужно в ArcMap выполнить векторную трансформацию. О методах преобразований, включая аффинное, подробнее здесь.Skywalker писал(а):Могу ли я на основе этих параметров ... создать новую систему координат X в ArcGis?
Ну, во-первых, менять надо не местные X и Y, а то, что ТС называет "Пулково-42". Потому как есть его указание "9 зона"ErnieBoyd писал(а):Если поменять местами местные X и Y, невязки аффинного преобразования не кажутся страшными:
Код: Выделить всё
2026 9507976.00 7144902.00 181161.00 588834.00 1.0
2027 9508373.00 7145482.00 181523.30 589439.00 1.0
2028 9508754.00 7146081.00 181887.19 590058.72 1.0
2029 9509155.00 7146688.00 182253.80 590671.50 1.0
2031 9507048.00 7147542.00 180492.53 591587.26 1.0
2032 9507433.00 7147926.00 180831.51 591966.78 1.0
2033 9506603.00 7147081.00 180080.78 591099.86 1.0
2226 9506067.00 7145271.00 179550.95 589219.00 1.0Код: Выделить всё
2026 9507976.000 7144902.000 181161.000 588834.000 1 -42.349 -123.207
2027 9508373.000 7145482.000 181523.300 589439.000 1 -76.930 -52.981
2028 9508754.000 7146081.000 181887.190 590058.720 1 -95.460 13.395
2029 9509155.000 7146688.000 182253.800 590671.500 1 -130.496 65.981
2031 9507048.000 7147542.000 180492.530 591587.260 1 72.238 96.884
2032 9507433.000 7147926.000 180831.510 591966.780 1 32.794 126.769
2033 9506603.000 7147081.000 180080.780 591099.860 1 98.507 29.851
2226 9506067.000 7145271.000 179550.950 589219.000 1 141.697 -156.693На GIS-Lab не так много геодезистов. Люди пользуются GIS, где X ≡ Easting, Y ≡ Northing.stout писал(а):во-первых, менять надо не местные X и Y, а то, что ТС называет "Пулково-42"
Та версия findkey рабатает с конформным преобразованием, при котором накладываются ограничения на значения параметров: A = E, B = −D. Получается такой набор:stout писал(а):Воспользовался вашей программкой findkey ... не совпадает с приведённым вами
Код: Выделить всё
0.94680720816955477, 0.036204771253028442, -9079696.4055638388, -0.036204771253028442, 0.94680720816955477, -5831653.4125548694Наверное плохо объяснил. Топикстартер привёл скан из каталога в местной системе, где икс назван иксом, а игрек — игреком. Думаю никто не будет возражать, что Пулково-42 это синоним СК-42. А в СК-42 X ≡ Northing, Y ≡ Easting. И упоминается 9 зона. Именно из этих соображений я и исходил.ErnieBoyd писал(а):На GIS-Lab не так много геодезистов. Люди пользуются GIS, где X ≡ Easting, Y ≡ Northing.
Предлагая аффинное преобразование, думаю не лишним будет напомнить, что в этом случае мы отказываемся не только от ортогональности осей, но и масштаб по осям будет разный.ErnieBoyd писал(а):Приведённые мной выше ключ и соответствующие невязки представляют аффинное преобразование без названных ограничений.
Это верно, однако в OP без обиняков сказано, что координаты сколоты с кривого скана:stout писал(а):Предлагая аффинное преобразование, думаю не лишним будет напомнить, что в этом случае мы отказываемся не только от ортогональности осей, но и масштаб по осям будет разный.
Skaywalker писал(а):Имеется не очень точно-привязанное графическое приложение … с него оцифрованы точки
catalog.zip