Привязка растра с неизвестной проекцией в QGIS

[UlvHare]

Здравствуйте!
Привязываю териогеографическую карту Алтайского края:

atlas_altay.jpg (280.68 КБ) 1124 просмотра

Проекция неизвестна, предположительно Каврайского. Добавил её в QGIS с таким определением (валидно):

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

PROJCRS["kavr",
    BASEGEOGCRS["pulk_42",
        DATUM["Pulkovo 1942",
            ELLIPSOID["Krassowsky 1940",6378245,298.3,
                LENGTHUNIT["metre",1]],
            ID["EPSG",6284]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["Degree",0.0174532925199433]]],
    CONVERSION["unnamed",
        METHOD["Equidistant Conic",
            ID["EPSG",1119]],
        PARAMETER["Latitude of false origin",54,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8821]],
        PARAMETER["Longitude of false origin",100,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8822]],
        PARAMETER["Latitude of 1st standard parallel",47,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Latitude of 2nd standard parallel",62,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8824]],
        PARAMETER["Easting at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8826]],
        PARAMETER["Northing at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8827]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]]

Однако, когда я в этой проекции отображаю референсную карту (ESRI topo), то она сильно поворачивается по часовой стрелке:

Screenshot.jpg (639.48 КБ) 1124 просмотра

Видно, что западная граница и Телецкое озеро встали почти вертикально.

В чём ошибка? Неправильно предположил проекцию растра? Неправильно её задал? Задача -- привязать карту с минимальными искажениями. По ней сделаны кальки маршрутов экспедиций с точками, которые тоже надо будет привязывать.

[AlexRomantsov]

Longitude of false origin",100 - это для всей РФ осевой, наблюдаемый поворот связан с этой долготой, так как эта проекция установлена для карты. Переключите карту в другую проекцию. Что важнее - проверяйте, насколько привязанная териогеографическая карта дружит с другими, сделайте ее полупрозрачной и сравните.

[tikhpetr]

Простите, внешне, это геологическая карта. Только СК-42 с соответствующей зоной

[UlvHare]

Спасибо! СК-42 зона 14 (Pulkovo 1942 / Gauss-Kruger zone 14, EPSG 28414) легла наиболее прилично. На Polynom2 чуть выгибался восточный край, поэтому таки сплайн. Собственно, карта попадает и в 15-ю зону, но привязывать сразу в двух я не представляю как

Screenshot_20240321_071139.jpg (613.93 КБ) 988 просмотров

[tikhpetr]

Не самая сложная процедура. Аккуратно, я подчеркиваю. Разбиваете область интереса на 2 территории по меридиану. Дальше создаете карту не в GK,а в Тrancwers Mercator с центральным меридианом ровно посредине Вашей территории. Все.

[UlvHare]

Спасибо!

[UlvHare]

Попробовал такое:

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

+proj=tmerc +lat_0=0 +lon_0=84 +k=1 +x_0=15500000 +y_0=0 +ellps=krass +units=m +no_defs

84 градуса - это точно посередине карты и как раз граница 2 зон СК-42. Пробовал также восточный сдвиг 15000000 (уносит привязанный растр на запад). В любом случае, даже со сплайнами появляются более сильные искажения по рекам и озёрам (до 3км) по сравнению с СК-42 зона 14. Похоже, для идеального совпадения надо резать карту по границе зон и привязывать по отдельности 14 и 15 зоны

[tikhpetr]

Я не работал в QGIS. Возможно, есть проблемы с параметрами перехода to WGS. Попробуйте создать пользовательскую СК(только восточное смещение без зоны) параметрами, к-рые в файле.

[tikhpetr]

И, кстати, масштаб тематической карты какой? На скрине не видно.

[UlvHare]

Спасибо! QGIS не поддерживает Coordinate Frame Rotation, только Position Vector, но на идею вы меня натолкнули. Масштаб 1:1600000

[tikhpetr]

Знак в "Scale_Difference" замените на+. Получится формат для Position Vector. А масштаб более чем странный.

[UlvHare]

В итоге лучше всего получилась компиляция из 14 и 15 зон:

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

PROJCRS["Pulkovo 1942 / Altay",
    BASEGEOGCRS["Pulkovo 1942",
        DATUM["Pulkovo 1942",
            ELLIPSOID["Krassowsky 1940",6378245,298.3,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4284]],
    CONVERSION["Altay Pulkovo-42",
        METHOD["Transverse Mercator",
            ID["EPSG",9807]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",84,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",15000000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            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["Topographic mapping (medium scale)."],
        AREA["Altay - 78?E to 90?E."],
        BBOX[41.04,78,81.27,90]]]

Интересно, что QGIS, предлагая оптимальные параметры трансформации в WGS84, выдала именно такие цифры, как у вас в файле.

На границе с Китаем всё равно оставались искажения, подозреваю, что либо не хватало данных у картографов, либо намеренно, чтобы испортить жизнь нарушителям границ. Подтянул сплайнами.

Спасибо, на этом задачу считаю решённой.

[tikhpetr]

Хорошо, что получилось. Но принципиально, ИМХО, что-то не совсем верно.

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

лучше всего получилась компиляция из 14 и 15 зон

Зачем? Из коробки на выбор EPSG:2604 "Pulkovo 1942 / 3-degree Gauss-Kruger CM 84E" и EPSG:2544 "Pulkovo 1942 / 3-degree Gauss-Kruger zone 28".

[tikhpetr]

Насколько я помню, все геологические карты выполнялись в GK6 СК-42. Но дело не в этом. Я работал на ГГС50 и на ГДП200. Папаша - на ГГС500. В смысле масштабов без трех 0. Миллионка, говорят, была визуальной при Царе-освободителе). А ТС указал, ваще, какой-то изуверский масштаб. Вязать такую карту в каноническом смысле - нонсенс. На скрине я, вроде, увидел сетку. ИМХО, правильно указав СК, этого за глаза достаточно