Как сделать карту qgis

Добавил пользователь Дмитрий К.
Обновлено: 29.08.2024

QGIS

Геопортал Публичной кадастровой карты Росреестра.
Растровая подложка (растр) из Sas.Планет.
XML файлы выписок Росреестр по объекту недвижимости или кадастровому плану территории.

Подключение геопортала ПКК

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

Подключение через сторонний сервис

Для подключения Публичной кадастровой карты возможно использование стороннего прокси сервера. Для этого необходимо просто добавьте два новых слоя в разделе XYZ Tiles:

Подключение Публичной кадастровой карты через сторонний прокси сервис

Или смотрите видео по подключению ПКК через сторонний прокси-сервис.

Обновленное видео по версии QGIS 3.16

Загрузка растра ПКК

Растры с параметрами привязки в QIS удобно подключить в формате ECW, с другими форматами и файлами привязки как то туго, необходимо предварительно конвертировать в удобоваримое для QGIS. Формировать растровую подложку в виде Публичной кадастровой карты можно через программу SAS.Планет.

Если растры имеют иную проекцию, отличную от текущей в проекте, то система координат преобразуется к проекции растра.

Формирование карты через обменные форматы

В качестве формирования векторной кадастровой карты могут выступать xml файлы кадастровых планов территорий Росреестра. Но рабочих вариантов подключения КПТ в QGIS на момент написания статьи нет. Для возможности использования данных кадастровых планов территории приходится использовать обменные форматы данных:

Подключение OGR-совместимых векторных слоев в QGIS

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

КПТ в форматах shp, mif/mid и dxf

Кроме этого имеются сторонние онлайн агрегаторы по преобразованию из XML файла в вышеуказанные форматы, но тем самым вы дарите им свои xml-ки. Так что пока приводит данные ресурсы в этой статье не буду.

Подключение местных систем координат

В зависимости от версии QGIS подключение параметров местной системы координат (МСК) производится через следующие форматы:

PROJ4 — для версий QGIS 3.8 и ниже;
WKT — для версий QGIS 3.9 и выше.

Формирование пользовательской системы координат в QGIS версии 3.14

В версии QGIS 3.14 прописывание в формате PROJ4 приведет его к преобразованию в WKT формат проекции.

Для общего примера возьмем универсальные параметры местных систем координат для MapInfo. Ниже приведен пример строки проекции МСК-12 зоны 1 в формате для прописывания в файле проекций MapInfo:

МСК в формате PROJ4

+proj=tmerc +lat_0=0 +lon_0=47.55 +k=1 +x_0=1250000 +y_0=-5914743.504 +ellps=krass +towgs84=23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22 +units=m +no_defs

МСК в формате WKT

Вот честно не могу сказать с какой версии произошел переход от proj4 к wkt по прописыванию параметров пользовательских проекций. Но тенденция к переходу к описанию проекций или систем координат через wkt-формат прослеживалась давно и рекомендовалась разработчиками данной программы. Покажем пример формирования такого файла в версии QGIS 3.14.

BOUNDCRS[
SOURCECRS[
PROJCRS["MSK-12 zone 1",
BASEGEOGCRS["MSK",
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",47.55,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",1250000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",-5914743.504,
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]]]]

Кадастровое деление территории

В новой версии публичной кадастровой карты была преобразована отдача границ кадастрового деления в виде векторных тайлов, что позволило получать эти границы в векторном формате. Как это делается представлено в следующем видео, где показывается 2-а варианта получения границ кадастрового деления по необходимой территории:

Мне сказали, на пикабу читать инструкции не любят, ведь можно просто показать, куда ткнуть. Поэтому я постараюсь уложиться в минимум слов, чтобы оставить максимум места для тыкания.

Будем разбирать на примере Швейцарии из моей карты (точнее, уже подаренной):

Мы будем работать с shape-файлами.

Что такое shape-файлы

Это слои трех видов: точки, кривые, полигоны. В основном вы будете скачивать архив, в котором много файлов *.shp. Все они будут названы интуитивно понятно:

Waterways – реки, buildings – постройки и т.д.

Где найти shape-файлы

1. Geofabrik - собраны отдельные карты стран и городов, а также можно найти карту всего мира. Не смог на сайте найти административные границы. Возможно, потому что там собраны только физические карты.

Пример: Швейцария. Нас интересует архив *.shp.zip. Но его качать не нужно, так как слои слишком детальны.

Подходит, если вы собираетесь создавать карту области (например, ЦФО) на большом формате.

2. Openstreetmap - можно экспортировать любую выбранную область (сверху кнопка ЭКСПОРТ). Выбирается область приближением к нужному месту на карте.

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

3. Mapzen - собраны карты крупных городов.

Подходит, если выбранная область города на сайте Openstreetmap слишком большая для экспорта, но весь ЦФО вам качать не нужно.

4. Gadm - границы стран.

Подходит для одной страны, если нужна граница одной страны.

5. Naturalearthdata - если английский слабо прокачан, смело переводите в хроме и качайте нужные слои: административные границы (всех стран), реки, здания, дороги и многое другое. Можно также найти и наложить на карту изображение гор и глубины морей и океанов. Для удобства карты поделены по уровню масштаба (10 – наиболее детальные)

Подходит для больших областей (часть света, карта всего мира).

6. Mapcruzin - карты не такие подробные, как в Openstreetmap, но подробнее, чем Naturalearthdata. Удивительно, но и здесь административных границ я не обнаружил.

Подходит для карты небольшой страны, на небольшом формате (а4).

7. Просто забивайте в поисковик “Нужная область” shape.

В чем открывать скачанные shape-файлы

Скачиваем Qgis. Не забываем про разрядность системы. Устанавливаем и запускаем ярлык:

Как работать в программе Qgis

Скачиваем с сайта Gadm границу Швейцарии:

В программе Qgis добавляем слой к проекту:

Жмем ОБЗОР. Выбираем файл CHE_adm0.shp из распакованного архива. Жмем ОТКРЫТЬ. Вы должны получить примерно следующее (цвет может отличаться):

Сразу обращаем внимание, что это полигон, и что вид Швейцарии не внушает доверия. Нам нужно изменить систему координат проекта. О координатах можно говорить много, но я скажу только необходимое:

1. В нашем проекте будем использовать систему Sphere Mercator.

Подходит для больших областей - весь мир или часть света. Вид будет привычен, как на интернет картах или картах на уроке географии.

2. На картах, где изображен весь мир, координаты сильно искажены. Поэтому наиболее правильно использовать координаты WGS 84. Но использовать их можно только на небольших областях (город, страна, область которой умещается в одной UTM zone).

Для небольшой области UTM zone ищите в поисковике.

Пример: для Парижа правильно использовать координаты WGS 84/ UTM zone 31N.

Примечание: чтобы понять почему только на небольших областях, загрузите границы России с системой координат для Парижа.

Меняем систему координат проекта:

Включаем галочку сверху Enable ‘on the fly’. В поиске ищем Sphere_Mercator. Выбираем. Жмем ОК. Любуемся результатом:

Добавим к проекту озера. Заходим на сайт Naturalearthdata. Скачиваем архив с озерами по верхней кнопке и дополнительные озера Европы по нижней. Добавляем к проекту файлы слоев с озерами (расширение *.shp).

Реки я хочу взять с сайта Mapcruzin. Одна Швейцария будет пусто смотреться, если скачать реки с Naturalearthdata. Нам нужны Switzerland Waterways. Добавляем к проекту файл слоя с реками (расширение *.shp).

Замечаем, что озера – это тоже полигоны. А нам нужен векторный файл из одних линий. Идем в меню:

И преобразуем по очереди два слоя озер (граница нам еще пригодится в качестве полигона). После преобразования слой автоматически добавится в проект, и будет назван "Линии из полигонов".

Совет: сразу переименовывайте автоматически добавленные слои после преобразований, а старые слои можно удалять из проекта.

Управление слоями расположено слева внизу:

Помимо удаления и переименования, можно перейти в свойства слоя, нажав правой кнопкой мыши по нему, и, например, изменить цвет слоя.

Теперь нужно удалить всё, что выходит за границы страны. Переходим в меню:

Исходный слой – слой, который нужно обрезать (например, реки). Слой обрезки – полигон, границы которого определят линию обрезки (полигон границы). По очереди обрезаем реки и озера, не забывая переименовывать слоя после автоматического добавления в проект и удалять ненужные слои, чтобы не путаться.

Осталось совсем немного! Преобразуем полигон границы в линию. Должно получиться нечто похожее:

Теперь необходимо пересохранить слои, чтобы создать единый слой (почему-то без пересохранения операция объединения вызывает стабильную ошибку). Правой кнопкой по слою/сохранить как.

Формат выбираем – AutoCAD DXF.

На строке File name нажимаем обзор, и выбираем место сохранения и имя файла.

И ОБЯЗАТЕЛЬНО не забываем выбрать систему координат нашего проекта - Sphere Mercator.

После сохранения слой автоматически добавится в проект с новым именем. Всплывающее окно будет спрашивать про систему координат. По умолчанию там стоит система, в который вы сохраните слой, поэтому просто нажимаем ОК. Проделываем данную процедуру для каждого слоя.

Осталось объединить слои. Заходим в меню:

Выбираем слои объединения, нажав многоточие в соответствующем пункте окна. Выбираем все необходимые слои. Нажимаем ОК. Нажимаем RUN. Все выскакивающие окна только уточняют систему координат. Везде нажимаем ОК. Сами знаете, что сделать со всем слоями, кроме результата объединения . Любуемся Швейцарией:

Сохраняем слой, как сохраняли предыдущие слои в формате AutoCAD DXF. Открываем в любом редакторе, который понимает данный формат (Corel Draw, Компас).

После открытия не забываем себе сказать, что вы великолепны.

Небольшой бонус до комментариев. Ответ на вопрос, который терзал мне душу около недели:

На этом все! Если возникнут вопросы - спрашивайте, стараюсь отвечать на все. Всем удачи и душевного спокойствия!

ArcGIS Desktop имеет возможность использовать базовые карты из ArcGIS online.

Есть ли у QGIS такие опции?

Обновление 2015: новый плагин с еще большим количеством опций фоновой карты - QuickMapServices

Оригинал: используйте плагин OpenLayers для получения фоновых карт Google Maps, Bing, OSM или Yahoo.

Смотрите другой / новый ответ на этой странице. Должен быть самый высокий голос, как я пишу это. Плагин не требуется.

Существует базовый функционал провайдера XYZ Tile Server, который был реализован с некоторыми другими приятными улучшениями UX для плиточных сервисов (доступно с QGIS 2.18). Это означает, что нет необходимости во внешнем плагине, хотя для простой настройки вы все равно можете использовать внешние плагины (см. Нижнюю часть этого поста), и он предлагает различные улучшения по сравнению с решениями на основе чистых плагинов .

На панели браузера найдите запись Tile Server и щелкните ее правой кнопкой мыши, чтобы добавить новую службу.

Введите URL-адрес службы, которую вы хотите использовать, заменив части x, y и z подстановкой фигурных скобок, как показано ниже.

Подготовка к печати ортофотоплана в нужном масштабе. Экспорт из Agisoft Photoscan в QGIS

Задача: подготовка к печати ортофотоплана в заданном масштабе, с масштабной линейкой, координатной сеткой и картушем компаса, указывающим направление на север. То, что у нас должно получиться, можно считать полноценной картой. Географи?ческая ка?рта — это построенное в картографической проекции, уменьшенное, обобщённое изображение поверхности Земли. Для решения я рекомендую использовать связку программ Agisoft Photoscan и QGIS. Почему именно эти программы? Фотограмметрический комплекс Photoscan (новое название программы — Metashape) от компании Agisoft — доступный и распространенный аэрогеодезический софт. Эта программа очень популярна и хорошо документирована в Сети. Геоинформационная система Quantum GIS — это развивающаяся, так же хорошо документированная, бесплатная программа с открытым исходным кодом. Обе программы существуют в версиях для Mac OS и Windows.

Экспорт растрового ортофотоплана из Agisoft Photoscan

Для получения растровой карты с масштабной линейкой и координатной сеткой нам нужно правильно экспортировать ортофотоплан. Процесс создания ортофотоплана описан в отдельной статье. Экспорт ортофотоплана нужно производить в прямоугольных координатах. Наиболее распространенная проекция для этого — Pulkovo 1942 Gaus Kruger. В нашем случае будет использована 6-я из 30 градусных зон. Это будет EPSG:28406 (последняя цифра датума — номер 30-градусной зоны). Узнать номер зоны, в которой вы будете экспортировать растр можно, например по широко распространенным генштабовским картам. Смотрим номенклатурный номер листа — в нашем случае это M-36-90. Вычитаем 30 из 36 и получаем 6 — искомый номер 30-градусной зоны.

Фрагмент карты Генерального штаба. Лист М-36-90

Соответственно, наш датум будет EPSG:28406. Теперь, когда мы знаем правильную проекцию, мы должны выбрать ее в списке доступных для экспорта ортофотоплана проекций. Выбираем нашу систему координат, вводя номер датума в строке фильтра.

Выбор системы координат для экспорта ортофотоплана в Agisoft Photoscan

Установка параметров ортофотоплана при экспорте из Agisoft Photoscan

Импорт ортофотоплана в QGIS

Теперь, запускаем QGIS. Напомню, что программа QGIS не работает с файлами и каталогами, в именах которых есть символы кириллицы. Проверьте чтобы таких символов не встречалось. Создаем новый проект. Проверяем в свойствах проекта систему координат. Если система координат не совпадает — установите правильную проекцию в настройках проекта QGIS. Как и в Photoscan, ищем ее по коду датума: EPGS:28406.

Выбор системы координат проекта QGIS

Окно менеджера источника данных

Окно перепроецирования данных QGIS

После импорта ортофотоплана в виде растрового слоя в QGIS, можно приступать к подготовке карты. Отрегулируйте размер импортированного ортофотоплана примерно так, как хотите его видеть на будущей карте. Далее перейдите по меню “Project” => “New Print Layout…” (Макет => Создать макет) (Command+P). Задайте имя будущего макета латиницей, после чего откроется отдельная программная оболочка для создания картографических материалов.

Создаем макет карты в QGIS

QGIS Layout менеджер с имортированным растром

Читайте также: