Местная система координат в Qgis

[QgisUser]

Доброго времени суток.
У меня следующая проблема. Заказчик попросил предоставить графические материалы в местной системе координат МСК 86, и интегрироваться в структуру информационной системы обеспечения градостроительной деятельности (ИСОГД) Нижневартовского района.
Подскажите, как это можно сделать в Qgis? И задавать местные системы координат в Qgis?
Заранее спасибо

[trir]

Что именно?
Добавить МСК 86 можно

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

Создайте пользовательскую систему координат:
[Установки] → [Ввод системы координат…]
Нажмите зелёный плюс и введите имя "МСК-86 зона 4" и параметры:

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

+proj=tmerc +lat_0=0 +lon_0=78.05 +k=1 +x_0=4500000 +y_0=-5811057.628 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs

Сохраните копии своих слоёв в новой СК.

[dab]

Подскажите, как это можно сделать в Qgis? И задавать местные системы координат в Qgis?

http://www.mapbasic.ru/store

[Evgeniya]

Создайте пользовательскую систему координат:
[Установки] → [Ввод системы координат…]
Нажмите зелёный плюс и введите имя "МСК-86 зона 4" и параметры:

Игорь, добрый день! Скажите пожалуйста, а как понять, какую именно зону МСК применять для qgis?

[Evgeniya]

Создайте пользовательскую систему координат:
[Установки] → [Ввод системы координат…]
Нажмите зелёный плюс и введите имя "МСК-86 зона 4" и параметры:

Игорь, добрый день! Скажите пожалуйста, а как понять, какую именно зону МСК применять для qgis?

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

Если имеются координаты точек в МСК, номер зоны виден как миллионы метров в восточной координате.
Если же координат нет, но известно географическое положение объекта, смотрим на Геодезисте тему МСК-** разбивка на зоны.

Если же неясно, что означает число 86 в строке "МСК-86 зона 4", то это обычно номер кадастрового округа. Впрочем, существует великое множество исключений.

[cpu-ekb]

Добрый день. Дело было так. Создал местную СК для Тюмени (MSK-72/ zone 1). Слои легли корректно на местности, перепроецировал их в созданную CRS. А потом странное. Стоит перезагрузить проект - все настройки CRS падают и в проекте и по слоям. И там и там она значится как Unknown, хотя слои на своем месте, никуда не улетели. До кучи при попытке создать для проекта SpatialLite БД я не вижу своей пользовательской CRS в предлагаемых вариантах. В чем может быть ошибка?

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

Моё дилетантское видение ситуации таково. Мы живём в смутное время перемен.

1. QGIS вслед за GDAL в описании СК переходит на нотацию WKT. В этой нотации описание проекций как таковых довольно ограничено.
2. В описаниях трансформаций датумов QGIS может использовать реализованные в PROJ пайпы, создавая последовательные комбинации операций.

Резюмируя, если QGIS позволит в качестве СК использовать комбинации, основанные на пайпах, можно будет дополнять проекцию аффинным или конформным преобразованием, как в MapInfo или Global Mapper. Жить станет лучше, жить станет веселее.

С другой стороны, мы становимся свидетелями процесса, когда городские системы в РФ одна за другой переопределяются так, чтобы дополнительных конформных преобразований не было. Может быть, когда-нибудь в QGIS появится комбинирование СК, но оно будет ненужно.

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

Создал местную СК для Тюмени (MSK-72/ zone 1)

Это именно МСК-72 зона 1?

[cpu-ekb]

Игорь, здравствуйте.
Я впервые делал пользовательскую CRS. Параметры для первой зоны МСК-72 брал отсюда https://mapbasic.ru/msk72. Найденные мной примеры и имевшаяся на руках МСК-66 (для Екатеринбурга) описывали ввод параметров именно в формате Proj String. Слегка напрягши голову, сообразил аналогичную строку для Тюмени. Но с этой проекцией начались описанные выше непонятки. Вчера, танцуя с бубном, несколько по другим причинам стал упаковывать слои проекта в базу SpatiaLite. Откровенно говоря, уже не помню порядок действий, но для проекта, созданного из слоев этой базы мне удалось создать новую CRS, проекцию которой я задал в WKT. И после этого все заработало. Я смог её назначить для слоев и после перезагрузки новая CRS никуда не подевалась... Проблема, какбэ, решилась, но теперь меня терзает вопрос, как же я преобразовал описание параметров из Proj String в WKT?

[lam]

Как у Вас выглядит описание в WKT?

[cpu-ekb]

Так вот...


BOUNDCRS[
SOURCECRS[
PROJCRS["unknown",
BASEGEOGCRS["unknown",
DATUM["Unknown based on Krassovsky, 1942 ellipsoid",
ELLIPSOID["Krassovsky, 1942",6378245,298.3,
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8901]]],
CONVERSION["unknown",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",66.08333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",1500000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",-6000000,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
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]]]]],
TARGETCRS[
GEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["geodetic latitude (Lat)",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["geodetic longitude (Lon)",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4326]]],
ABRIDGEDTRANSFORMATION["Transformation from unknown to WGS84",
METHOD["Position Vector transformation (geog2D domain)",
ID["EPSG",9606]],
PARAMETER["X-axis translation",23.57,
ID["EPSG",8605]],
PARAMETER["Y-axis translation",-140.95,
ID["EPSG",8606]],
PARAMETER["Z-axis translation",-79.8,
ID["EPSG",8607]],
PARAMETER["X-axis rotation",0,
ID["EPSG",8608]],
PARAMETER["Y-axis rotation",0.35,
ID["EPSG",8609]],
PARAMETER["Z-axis rotation",0.79,
ID["EPSG",8610]],
PARAMETER["Scale difference",0.99999978,
ID["EPSG",8611]]]]

Немного разобрался. QGIS, когда показывает неизвестную пользовательскую проекцию, заданную в Proj String, тут же приводит её параметры в WKT. Я оттуда скопировал и создал уже новую, которая держится и не пропадает после перезагрузки.

[lam]

если QGIS позволит в качестве СК использовать комбинации, основанные на пайпах, можно будет дополнять проекцию аффинным или конформным преобразованием, как в MapInfo или Global Mapper. Жить станет лучше, жить станет веселее.

https://gis.stackexchange.com/questions ... transforma

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

Спасибо за ссылку! Это просто праздник какой-то!