Форумы GIS-Lab.info

Добрый день! Подскажите как написать в proj4 или wkt данные с выписки:
Населенный пункт: Охинский район
Параметры от УСК к: МСК-65
Х1:7777,78
У1:-64742,447
Х2:1390795,40
У2: 1320840,403
@: -1°2' 27,16"
m: 0,99993625

С таким вариантом не получается
+proj=tmerc +lat_0=0.15 +lon_0=142.716666666667 +k=0.99993625 +x_0=1256098.956 +y_0=1398573.18 +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.219999999995224 +units=m +no_defs

Это типа аффинные, но к чему?

Запросили ключ перехода в МСК-65, они прислали вот это. А как его загнать в Qgis вообще не понимаю.
И в выписке подписано считайте по формулам Гельмерта.

Системы координат с дополнительным преобразованием

Думаю очевидно как получить параметры A1, A2, A0, B1, B2, B0 по вашим данным

mark142 писал(а): ↑18 май 2023, 09:44С таким вариантом не получается

Попробуйте такой:

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

DERIVEDPROJCRS["USK Okhtinsky",
    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",-1410715.53760689,
            LENGTHUNIT["metre",1],
            ID["EPSG",8623]],
        PARAMETER["A1",0.999898732782675,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8624]],
        PARAMETER["A2",0.0181669032260463,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8625]],
        PARAMETER["B0",-1358881.19824162,
            LENGTHUNIT["metre",1],
            ID["EPSG",8639]],
        PARAMETER["B1",-0.0181669032260463,
            SCALEUNIT["coefficient",1],
            ID["EPSG",8640]],
        PARAMETER["B2",0.999898732782675,
            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 - Okhtinsky district"],
        BBOX[51.89,141.63,54.43,143.33]]]

Если не получается, выложите данные, на которых можно проверить.

Координаты не те. Выписку прилагаю шейп

65_23_0000016.zip: (382.25 КБ) 435 скачиваний

mark142 писал(а): ↑19 май 2023, 01:01Координаты не те. Выписку прилагаю шейп

В архиве слой 65_23_0000016_2022-11-28_kpt11.xml в МСК-65:

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

PROJCRS["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]],
        ID["EPSG",4200]],
    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]]],
    CS[Cartesian,2],
        AXIS["northing (X)",north,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["easting (Y)",east,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["Russia - Sakhalin - MSK65 zone 1"],
        BBOX[45.89,141.63,54.43,144.22]]]

Для проверки нужны либо данные в обеих системах координат, либо хотя бы данные в УСК, соотнесённые с картами или снимками.

Есть такие файлы, только вот почему то не могу загрузить их сюда

закинь в архиве

любой файл или архив пишет:
ОШИБКА
Не удалось закачать вложение ./files/31431_47106e3d60a78f8af91a1fb85dbefa94.

попробуй ещё

неа. не получается. все тоже самое

положи в облако, а здесь ссылку

Скрин совмещения выписок https://drive.google.com/file/d/1hvz5Jo ... share_link
Сама выписка 2018 года со старыми координатами https://drive.google.com/file/d/1iKJ65q ... sp=sharing

Спасибо!

Присвоил слою 2018 года "USK Okhtinsky", слою 2022 года "Pulkovo 1942 / MSK65 zone 1". И расхождений не вижу, совмещение идеальное!

Форумы GIS-Lab.info

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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