координаты в 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]

проверка