координаты в proj4

[mark142]

Да спасибо. Совмещение идеальное.
Подскажите еще в МапИнфо ключ перехода вот такой: "МСК Росреестр ОХА", 8, 1001, 7, 144, 0, 1, 0, -5900000. Я так понимаю он с тем что выложили первым "USK Okhtinsky" должен совпадать?

И сегодня пришла выписка по Ногликскому району (тоже самое можно):
Параметры от УСК к: МСК-65
Х1:50370,312
У1:2040,113
Х2:1262728,09
У2: 1327088,457
@: -1°0' 9,69"
m: 0,99997174

Архив с выписками https://drive.google.com/file/d/1sehdUt ... p=drivesdk

[Игорь Белов]

в МапИнфо ключ перехода вот такой … должен совпадать?

Нет, это совсем не то. Должно быть примерно такое (упрощённый вариант без охвата):

Код: Выделить всё

"USK Okhinsky", 1008, 1001, 7, 142.716666666667, 0.15, 1, 1300000, -4500000, 7, 0.99989873, 0.01816690, -1410715.538, -0.01816690, 0.99989873, -1358881.198

Вот WKT для Ногликского:

Код: Выделить всё

DERIVEDPROJCRS["USK Nogliksky",
    BASEPROJCRS["Pulkovo 1942 / MSK65 zone 1",
        BASEGEOGCRS["Pulkovo 1942",
            DATUM["Pulkovo 1942",
                ELLIPSOID["Krassowsky 1940",6378245,298.3,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["Degree",0.0174532925199433]]],
        CONVERSION["MSK65 zone 1",
            METHOD["Transverse Mercator",
                ID["EPSG",9807]],
            PARAMETER["Latitude of natural origin",0.15,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8801]],
            PARAMETER["Longitude of natural origin",142.716666666667,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8802]],
            PARAMETER["Scale factor at natural origin",1,
                SCALEUNIT["unity",1],
                ID["EPSG",8805]],
            PARAMETER["False easting",1300000,
                LENGTHUNIT["metre",1],
                ID["EPSG",8806]],
            PARAMETER["False northing",-4500000,
                LENGTHUNIT["metre",1],
                ID["EPSG",8807]]]],
    DERIVINGCONVERSION["Affine",
        METHOD["Affine parametric transformation",
            ID["EPSG",9624]],
        PARAMETER["A0",-1346980.21157804,
            LENGTHUNIT["metre",1],
            ID["EPSG",8623]],
        PARAMETER["A1",0.999875130637336,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8624]],
        PARAMETER["A2",0.0174998722550523,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8625]],
        PARAMETER["B0",-1188976.22347953,
            LENGTHUNIT["metre",1],
            ID["EPSG",8639]],
        PARAMETER["B1",-0.0174998722550523,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8640]],
        PARAMETER["B2",0.999875130637336,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8641]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["Russia - Sakhalin - Nogliksky district"],
        BBOX[50.73,142.23,52.96,143.67]]]

Совпадение между объектами 2019 и 2022 года получается идеальное. Но! Слой 2019 года запороли, узлы с отрицательными координатами схлопнуты на абсциссу и ординату.

[mark142]

Спасибо! Но Ногликский нужно в MSK65 zone 1.

[Игорь Белов]

Ногликский нужно в MSK65 zone 1

Как и в предыдущем наборе данных на Охинский, старый слой в УСК, его просто трансформируем в МСК-65.
А вот со слоем 2022 года тяжёлый случай: первые 195 объектов в МСК-65, оставшиеся 119 в УСК. Нужно этот слой разделить на два, второй трансформировать в МСК, после чего воссоединить две половинки.
И да, та часть слоя 2022 года, что в УСК, так же изуродована, как и слой 2019 года.

[mark142]

Спасибо огромное!

[trir]

проверка

[mark142]

Добрый день!
Наткнулся на такую статью по по МСК-65 https://terraingis.ru/msk-65.html
И смотрю, что те параметры которые здесь давались немного различаются , если смотреть ГИС QGIS МСК-65 зона 1.
Какие параметры все таки использовать, можете подсказать?

[Игорь Белов]

Какие параметры все таки использовать, можете подсказать?

• Это в рекомендуемом формате WKT 2, набор параметров исходный.
• То в устаревшем формате PROJ 4, набор параметров адаптирован под древние навигаторы.

В плане координатных преобразований оба варианта эквивалентны, можете использовать любой.

[mark142]

Добрый день! А еще можно вопрос. Почему не совмещается СК-42 и МСК65. И как это исправить?

[AlexRomantsov]

Обе системы - кривые как сабли. Их параметры не универсальны даже для региона, идеально подобранные описания для Южно-Сахалинска (с полевыми измерениями, калибровками, локализациями) не сработают в Углегорске. Но, конечно, если Вам достаточно точности в 2-3 метра, то ок.
В вашем случае надо понять, что лучше лежит, зеленая или сиреневая. И поправить proj зеленой или сиреневой (например другие сдвиги x_0 y_0 задать)

[mark142]

Да только сдвиги не совсем помогут он еще и крутится по своей оси (на 270 км крутит почти на 360 градусов) приложил файлик с примерным углом и сдвигом

[AlexRomantsov]

При том, что на одной долготе, практически... Каталоги (пунктов) в СК42 более деформированы, чем в МСК (а могут быть "встречные" деформации в МСК). И "архивные" материалы в СК42 наследуют эти деформации. Для "танковых" совковых карт М1:50000 5 метров - и не ошибка, СК42 - для таких карт создавалась. Поскольку СК42 давно отменена, каталоги не обновляются. Совместить можно, если сезон убить на полевые измерения, локализовав обе системы для Охи, Ногликов и еще 5-6 мест. Совсем правильно создавать сетки трансформаций типа NTv2, но это еще дороже (по полевым)

[mark142]

Спасибо за информацию!