GIS-Lab - Вклад [ru]

Работа с публичной кадастровой картой на примере QGIS

2016-08-30T10:35:45Z

Institor: Изменился адрес доступа к сервисам

{{Статья|Черновик}}
Публичная кадастровая карта Росреестра (ПКК) создана на основе ArcGIS Server и может работать как [http://gis-lab.info/qa/ogc-intro.html '''стандартный WMS-сервис''']. 
Эта статья будет полезна тем, кто желает использовать данные '''Государственного кадастра недвижимости''' в QGIS.
 
'''ВНИМАНИЕ! Слои ПКК имеют настройки видимости (минимальный и максимальный масштабы) и вы не всегда сможете видеть все слои. Например, земельные участки видны только на крупных масштабах, на самом мелком масштабе выдно тольк кадастровое деление РФ - кадастровые округа, обычно соответствующие административно-территориальному делению РФ!''' 
 
Публичная кадастровая карта предназначена для публикации в интернете данных Государственного кадастра недвижимости (ГКН). К объектам недвижимости относятся земельные участки и их части, объекты капитального и незавершённого строительства, помещения.
 
Карта состоит из следующих слоёв, повторяющих единицы кадастрового деления:
* Граница РФ.
* Границы кадастровых округов. Кадастровый округ зачастую совпадает с субъектом РФ, но всегда. При ликвидации субъекта РФ, его кадастровый округ не ликвидируется.
* Границы кадастровых районов. Совпадают с административными районами субъектов РФ.
* Границы кадастровых кварталов. Совпадают с территоориальным зонированием территорий. Наименьшая единица кадастрового деления.
 
а так же из объектов недвижимости:
* земельных участков (вся территория РФ)
* объектов капитального строительства (Кемеровская обл, пилотная зона, январь 2013)
 
Самый нужный слой - групповой слой 3 "Земельные участки". Это земельные участки, но без всех остальных информация будет не полная, т.к. '''полные''' кадастровые номера состоят из совокупности номеров объектов кадастрового деления. Например, номер 70:21:023256:002. 002 - номер участка, 023256 кадастровый квартал в городе Томск (21) Томской области (70).
 
== QGIS 1.8 ==
 
Для того чтобы узнать что такое QGIS, как выбрать систему координат и как подключаются слои, есть[http://gis-lab.info/docs/qgis/ Руководство]. 

Для создания нового WMS-слоя нужно войти меню "Слой" и выберите инструмент "Новый WMS-слой". 
В появившимся диалоге нажать кнопку "Создать", появится диалог создания нового WMS-подключения (Рис.1).
В нём нужно ввести имя подключения, например, "ПКК" и его актуальный URL:

<pre>http://pkk5.rosreestr.ru/arcgis/services/Cadastre/CadastreWMS/MapServer/WMSServer</pre>

Неактуальные URL:

<pre>http://maps.rosreestr.ru/arcgis/services/Cadastre/CadastreWMS/MapServer/WMSServer?request=GetCapabilities&service=WMS
http://maps.rosreestr.ru/arcgis/rest/services/Cadastre/Cadastre
http://maps.rosreestr.ru/arcgis/services/CadastreNew/CadastreWMS/MapServer/WMSServer</pre>

Далее кнопку "ОК". 

[[Файл:NewWMS_QGIS.png|800px|thumb|center|<center>Рис.1</center>]] 
 
 
Далее, после закрытия диалога нового подключения, нужно нажать кнопку "Подключить". 
После этого будет доступен список доступных слоёв этого сервиса. Можно выбрать слои. (Рис.2). На первом уровне доступны групповые слои 1, 3, 8, 15, 25, образующие вместе весь кадастр. Их можно выбрать все сразу, поочерёдно нажимая на каждый из них мышью, а можно добавлять по отдельности, повторив предыдущие шаги столько раз, сколько вы хотите добавить слоёв. Второй подход полезен тогда, когда нужно иметь возможность отключать слои по отдельности. 
 

[[Файл:Qgs_pkk1.png|800px|thumb|center|<center>Рис.2</center>]] 
 
 
После этого нужно нажать кнопку "Добавить" и через некоторое время (время загрузки зависит от скорости интернет-соединения) WMS-слой добавится на карту. 
Далее нужно нажать кнопку "Close" для закрытия диалога добавления WMS-слоя. Увидеть добавленный слой скорее всего не получиться из-за настройки видимости в текущем масштабе. Необходимо добавить любые ваши данные как слой (на рисунке это слой "М_38") и увеличить охват карты до этого слоя.
 
На данном рисунке (Рис.3) показана часть территории РФ в месте расположения Омской (кадастровый округ 55) области. Кадастровые округа подписаны самым крупным шрифтом и их границы обозначены толстыми красными линиями. Более мелкие надписи и тонкие красные линии - кадастровые районы.
 
[[Файл:Qgs_pkk2.png|800px|thumb|center|<center>Рис.3</center>]]
 
 
Далее, можно посмотреть информацию на интересующий участок, например. Используя стандартный инструмент QGIS "Определить объеекты", указать интересующее место на карте и увидеть информацию. Информация предоставляется в не очень удобном виде, но основной атрибут - кадастровый номер можно найти в поле "CAD_NUM"(Рис.4).
 
 
[[Файл:View PCM in QGIS.png|800px|thumb|center|<center>Рис.4 Получения информации.</center>]]
 
 
[[Файл:View_ZU_PCM_in_QGIS.png|800px|thumb|center|<center>Рис.5 На рисунке видно, как кадастровый квартал "режет" участки. Это один из примеров некачественных кадастровых данных.</center>]] 
 
 
----

Как собрать и показать снимок Landsat из отдельных каналов в QGIS

2013-06-26T11:56:32Z

Institor: Добавил пример файла vrt

{{Статья|Опубликована|qgis-landsat-merge}}
{{Аннотация|Собираем Landsat из отдельных файлов-каналов.}}

Эта статья — пособие для новичков в ГИС, объясняющая как решить данную конкретную задачу. Подразумевается что QGIS у читателя уже установлена, если это не так — установите её используя [http://gis-lab.info/qa/qgis-osgeo4w.html инструкцию].

Также предполагается, что вы знаете как доставать данные Landsat из архива Glovis. Если это не так — ознакомьтесь со статьей «[http://gis-lab.info/qa/landsat-glovis.html Получение бесплатных космических снимков Landsat TM, ETM+ через Glovis]».

Если вы хотите потренироваться в создании снимка из отдельных каналов, то можно использовать [http://gis-lab.info/data/samples/lt51700282007141mor00.7z этот архив] (полная сцена 170/28, 7 каналов, формат TIF, UTM zone 38, Калмыкия, оригинал с Glovis).

=== Подготовка ===

Для работы нам необходимо расширение GDALTools, которое входит в состав QGIS начиная с версии 1.5.0. Проверить, установлено ли оно можно вызвав «Менеджер модулей» из меню «Модули → Управление модулями» и поискав расширение там. Если расширение есть и не активно, необходимо его включить.

[[Файл:qgis-landsat-merge-01.png|center]]

Если расширения нет в списке доступных, необходимо его установить. Для этого:

# Откройте «Установщик модулей» (в меню «Модули → Загрузить модули»)
# В окне установщика выберите вкладку «Репозитории» и проверьте подключен ли репозиторий «Faunalia Repository». Если репозиторий подключен переходим к пункту 5.
# Если нет, нажимаем кнопку «Добавить сторонние репозитории», читаем и принимаем к сведению предупреждение.
# Переходим на вкладку «Параметры» и в группе «Разрешённые модули» выставляем переключатель в положение «Показывать все модули, кроме помеченных как экспериментальные»
# Переходим на вкладку «Модули» и при помощи фильтра или используя полосу прокрутки находим в списке GDALTools, выделяем его и нажимаем «Установить модуль»
# После загрузки модуля он будет автоматически активирован.

Установленный и активированный модуль создает еще один пункт в главном меню QGIS (если он отсутствует) — «Растр». Или добавляет свои инструменты в существующее меню.

[[Файл:qgis-landsat-merge-02.png|center]]

=== Собираем растр из каналов ===

Итак, у нас имеется снимок Landsat в виде растров, представляющих отдельные каналы изображения, всего 7 файлов. Нам же необходимо получить один файл с 7 каналами.

Чтобы было удобнее работать, все 7 файлов складываем в одну папку. Набор файлов может выглядеть следующим образом, в зависимости от того, какой это Landsat, может добавляться еще файл _B80 значительно большего размера.

[[Файл:qgis-landsat-merge-09.png|center]]

Открываем меню «Растр» и выбираем пункт «Объединение». Начинаем заполнять поля диалога:

* в поле «Исходные файлы» через запятую указываются полные имена файлов-каналов. Так как все растры у нас лежат в одной папке, можно поступить проще: нажать «Выбрать...» и выделить все файлы, нажав Shift. '''Примечание''': Если указываются не все файлы (например исключается 6 канал, что обычно делается так как он тепловой), то необходимо выбрать файлы с первого по последний зажав Ctrl.
* в поле «Целевой файл» указываем полный путь к результирующему файлу. Как и в предыдущем случае, кнопка «Выбрать...» облегчает эту задачу. '''Примечание''': сохранить результат можно в любой, поддерживаемый GDAL, формат. Выбрать формат можно при помощи выпадающего списка «Типы файлов». Обычный формат для такого случая GeoTIFF. Имя файлу лучше давать осмысленное, например для случая выше: 170028_070521, первые цифры обозначают номенклатуру снимка в системе WRS-2, последние 6 — дату в формате ГГММДД.
* так как нам нужен растр с 7 каналами, ставим флажок «Склеить поканально». Если не поставить его, то получим одноканальный растр.
* если полученный растр нужно использовать в дальнейшей работе, можно установить флажок «Добавить результат на карту»

Отдельно стоит сказать о группе «Параметры создания». С ее помощью можно выполнить «тонкую» настройку создаваемого файла, например, включить использование сжатия, генерацию world-файла и др. У каждого формата свой набор возможных параметров, ознакомиться с ними можно на странице описания формата ([http://www.gdal.org/formats_list.html список форматов]).

Например, для формата TIFF можно установить следующие параметры: COMPRESS=LZW, TILED=YES, TFW=YES и BIGTIFF=IF_NEEDED. Этим мы включим сжатие результата алгоритмом Лемпеля — Зива — Велча (Lempel-Ziv-Welch), создадим world-файл, при необходимости будет задействован режим BIGTIFF, а сам растр будет создан с разбивкой на тайлы, что в некоторых случаях повышает скорость загрузки и отрисовки.

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

Если у вас нет особенных требований, не устанавливайте дополнительные параметры создания.

[[Файл:qgis-landsat-merge-03.png|center]]

После того как все необходимые поля заполнены нажимаем OK и ждем окончания процесса.

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

=== Собираем растр из каналов виртуально ===
Первый способ сборки растра физически объединяет все файлы. Если файлов много и они большого размера (например, наборы Landsat 8 имеют 10 каналов по 112 Мбайт каждый и панхроматический канал "весом" 450 Мбайт), результирующий файл может выйти за ограничения формата Tiff по размеру. В таком случае может быть полезен способ "виртуального" объединения каналов: создание виртуального растра ([http://www.gdal.org/ogr/drv_vrt.html GDAL Virtual Format]).

Итак, меню "Растр" - "Создать виртуальный растр (каталог)".

Диалог настройки очень похож на предыдущий:

[[Файл:GDAL VRT.png|center]]

Точно также выбираем исходные каналы, результирующий файл (это будет файл с расширением *.vrt и структурой обычного XML), значение "нет данных" (0), "сохранять в каналах" (ключ -separate, иначе получится одноканальный растр).

Пункт "Разрешение" отвечает за разрешение нашего виртуального растра: исходные файлы могут быть разного размера. Например, при выборе "максимальное" результирующий растр будет иметь разрешение, соответствующее максимальному разрешению из набора исходных файлов.

После нажатия "ОК" создается *.vrt файл, с которым можно работать точно так же, как с любым обычным растром. Операций с самими растрами не производится.
[https://dl.dropboxusercontent.com/u/997092/VRT.vrt Пример файла *.vrt (11 каналов Landsat 8)]

=== Настройка отображения растра ===

Результатом предыдущей операции был новый файл и теперь мы хотим на него посмотреть.

По умолчанию растр загружается с использованием стандартной комбинации каналов 1-2-3, что не всегда удобно. Изменить комбинацию каналов и некоторые другие параметры отображения растра можно в диалоговом окне «Свойства слоя» (вызывается из контекстного меню слоя, по правой кнопке мыши, пункт «Свойства»).

Комбинация каналов и алгоритм улучшения контраста настраиваются на вкладке «Стиль».

[[Файл:qgis-landsat-merge-04.png|center|thumb|400px]]

Ниже показан один и тот же растр, с разными комбинациями каналов и алгоритмом улучшения контраста.

[[Файл:qgis-landsat-merge-05.png|center|thumb|400px|<center>Комбинация каналов 1-2-3, улучшение контраста не используется</center>]]

[[Файл:qgis-landsat-merge-06.png|center|thumb|400px|<center>Комбинация каналов 5-4-3, растяжение до мин/макс + стандартное отклонение 2.0</center>]]

В этом же окне, на вкладке «Прозрачность» можно указать какое значение растра используется как NODATA. В нашем примере черные области вокруг растра имеют значение 0 (ноль), поэтому чтобы убрать их мы установили 0 в качестве значения «нет данных»

[[Файл:qgis-landsat-merge-07.png|center|thumb|400px]]

Результат применения этой настройки ниже.

[[Файл:qgis-landsat-merge-08.png|center|thumb|400px]]

Если вы хотите потренироваться в создании снимка из отдельных каналов, то можно использовать [http://gis-lab.info/data/samples/lt51700282007141mor00.7z этот архив] (полная сценаа 170/28, 7 каналов, формат TIF, UTM zone 38, Калмыкия, оригинал с Glovis).

=== Ссылки по теме ===

* [http://gis-lab.info/qa/landsat-bandcomb.html Интерпретация комбинаций каналов данных Landsat TM / ETM+]
* [http://linfiniti.com/2011/05/gdal-efficiency-of-various-compression-algorithms/ GDAL: efficiency of various compression algorithms]

Как собрать и показать снимок Landsat из отдельных каналов в QGIS

2013-06-26T11:22:46Z

Institor: Добавил создание витуального растра

{{Статья|Опубликована|qgis-landsat-merge}}
{{Аннотация|Собираем Landsat из отдельных файлов-каналов.}}

Эта статья — пособие для новичков в ГИС, объясняющая как решить данную конкретную задачу. Подразумевается что QGIS у читателя уже установлена, если это не так — установите её используя [http://gis-lab.info/qa/qgis-osgeo4w.html инструкцию].

Также предполагается, что вы знаете как доставать данные Landsat из архива Glovis. Если это не так — ознакомьтесь со статьей «[http://gis-lab.info/qa/landsat-glovis.html Получение бесплатных космических снимков Landsat TM, ETM+ через Glovis]».

Если вы хотите потренироваться в создании снимка из отдельных каналов, то можно использовать [http://gis-lab.info/data/samples/lt51700282007141mor00.7z этот архив] (полная сцена 170/28, 7 каналов, формат TIF, UTM zone 38, Калмыкия, оригинал с Glovis).

=== Подготовка ===

Для работы нам необходимо расширение GDALTools, которое входит в состав QGIS начиная с версии 1.5.0. Проверить, установлено ли оно можно вызвав «Менеджер модулей» из меню «Модули → Управление модулями» и поискав расширение там. Если расширение есть и не активно, необходимо его включить.

[[Файл:qgis-landsat-merge-01.png|center]]

Если расширения нет в списке доступных, необходимо его установить. Для этого:

# Откройте «Установщик модулей» (в меню «Модули → Загрузить модули»)
# В окне установщика выберите вкладку «Репозитории» и проверьте подключен ли репозиторий «Faunalia Repository». Если репозиторий подключен переходим к пункту 5.
# Если нет, нажимаем кнопку «Добавить сторонние репозитории», читаем и принимаем к сведению предупреждение.
# Переходим на вкладку «Параметры» и в группе «Разрешённые модули» выставляем переключатель в положение «Показывать все модули, кроме помеченных как экспериментальные»
# Переходим на вкладку «Модули» и при помощи фильтра или используя полосу прокрутки находим в списке GDALTools, выделяем его и нажимаем «Установить модуль»
# После загрузки модуля он будет автоматически активирован.

Установленный и активированный модуль создает еще один пункт в главном меню QGIS (если он отсутствует) — «Растр». Или добавляет свои инструменты в существующее меню.

[[Файл:qgis-landsat-merge-02.png|center]]

=== Собираем растр из каналов ===

Итак, у нас имеется снимок Landsat в виде растров, представляющих отдельные каналы изображения, всего 7 файлов. Нам же необходимо получить один файл с 7 каналами.

Чтобы было удобнее работать, все 7 файлов складываем в одну папку. Набор файлов может выглядеть следующим образом, в зависимости от того, какой это Landsat, может добавляться еще файл _B80 значительно большего размера.

[[Файл:qgis-landsat-merge-09.png|center]]

Открываем меню «Растр» и выбираем пункт «Объединение». Начинаем заполнять поля диалога:

* в поле «Исходные файлы» через запятую указываются полные имена файлов-каналов. Так как все растры у нас лежат в одной папке, можно поступить проще: нажать «Выбрать...» и выделить все файлы, нажав Shift. '''Примечание''': Если указываются не все файлы (например исключается 6 канал, что обычно делается так как он тепловой), то необходимо выбрать файлы с первого по последний зажав Ctrl.
* в поле «Целевой файл» указываем полный путь к результирующему файлу. Как и в предыдущем случае, кнопка «Выбрать...» облегчает эту задачу. '''Примечание''': сохранить результат можно в любой, поддерживаемый GDAL, формат. Выбрать формат можно при помощи выпадающего списка «Типы файлов». Обычный формат для такого случая GeoTIFF. Имя файлу лучше давать осмысленное, например для случая выше: 170028_070521, первые цифры обозначают номенклатуру снимка в системе WRS-2, последние 6 — дату в формате ГГММДД.
* так как нам нужен растр с 7 каналами, ставим флажок «Склеить поканально». Если не поставить его, то получим одноканальный растр.
* если полученный растр нужно использовать в дальнейшей работе, можно установить флажок «Добавить результат на карту»

Отдельно стоит сказать о группе «Параметры создания». С ее помощью можно выполнить «тонкую» настройку создаваемого файла, например, включить использование сжатия, генерацию world-файла и др. У каждого формата свой набор возможных параметров, ознакомиться с ними можно на странице описания формата ([http://www.gdal.org/formats_list.html список форматов]).

Например, для формата TIFF можно установить следующие параметры: COMPRESS=LZW, TILED=YES, TFW=YES и BIGTIFF=IF_NEEDED. Этим мы включим сжатие результата алгоритмом Лемпеля — Зива — Велча (Lempel-Ziv-Welch), создадим world-файл, при необходимости будет задействован режим BIGTIFF, а сам растр будет создан с разбивкой на тайлы, что в некоторых случаях повышает скорость загрузки и отрисовки.

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

Если у вас нет особенных требований, не устанавливайте дополнительные параметры создания.

[[Файл:qgis-landsat-merge-03.png|center]]

После того как все необходимые поля заполнены нажимаем OK и ждем окончания процесса.

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

=== Собираем растр из каналов (второй способ) ===
Первый способ сборки растра физически объединяет все файлы. Если файлов много и они большого размера (например, наборы Landsat 8 имеют 10 каналов по 112 Мбайт каждый и панхроматический канал "весом" 450 Мбайт), результирующий файл может выйти за ограничения формата Tiff. В таком случае может быть полезен способ "виртуального" объединения каналов: создание виртуального растра (GDAL Virtual Formats)
Итак, меню "растр" - "Создать виртуальный растр (каталог)".

Диалог настройки очень похож на предыдущий:

[[Файл:GDAL VRT.png|center]]

Точно также выбираем исходные каналы, результирующий файл (это будет файл с расширением *.vrt и структурой обычного XML), значение "нет данных" (0), "сохранять в каналах" (иначе получится одноканальный растр).
Пункт "Разрешение" отвечает за разрешение нашего виртуального растра: исходные файлы могут быть разного размера. Например, при выборе "максимальное" результирующий растр будет иметь разрешение, соответствующее максимальному разрешению из набора исходных файлов.

После нажатия "ОК" создается *.vrt файл, с которым можно работать точно так же, как с любым обычным растром.
=== Настройка отображения растра ===

Результатом предыдущей операции был новый файл и теперь мы хотим на него посмотреть.

По умолчанию растр загружается с использованием стандартной комбинации каналов 1-2-3, что не всегда удобно. Изменить комбинацию каналов и некоторые другие параметры отображения растра можно в диалоговом окне «Свойства слоя» (вызывается из контекстного меню слоя, по правой кнопке мыши, пункт «Свойства»).

Комбинация каналов и алгоритм улучшения контраста настраиваются на вкладке «Стиль».

[[Файл:qgis-landsat-merge-04.png|center|thumb|400px]]

Ниже показан один и тот же растр, с разными комбинациями каналов и алгоритмом улучшения контраста.

[[Файл:qgis-landsat-merge-05.png|center|thumb|400px|<center>Комбинация каналов 1-2-3, улучшение контраста не используется</center>]]

[[Файл:qgis-landsat-merge-06.png|center|thumb|400px|<center>Комбинация каналов 5-4-3, растяжение до мин/макс + стандартное отклонение 2.0</center>]]

В этом же окне, на вкладке «Прозрачность» можно указать какое значение растра используется как NODATA. В нашем примере черные области вокруг растра имеют значение 0 (ноль), поэтому чтобы убрать их мы установили 0 в качестве значения «нет данных»

[[Файл:qgis-landsat-merge-07.png|center|thumb|400px]]

Результат применения этой настройки ниже.

[[Файл:qgis-landsat-merge-08.png|center|thumb|400px]]

Если вы хотите потренироваться в создании снимка из отдельных каналов, то можно использовать [http://gis-lab.info/data/samples/lt51700282007141mor00.7z этот архив] (полная сценаа 170/28, 7 каналов, формат TIF, UTM zone 38, Калмыкия, оригинал с Glovis).

=== Ссылки по теме ===

* [http://gis-lab.info/qa/landsat-bandcomb.html Интерпретация комбинаций каналов данных Landsat TM / ETM+]
* [http://linfiniti.com/2011/05/gdal-efficiency-of-various-compression-algorithms/ GDAL: efficiency of various compression algorithms]

Файл:GDAL VRT.png

2013-06-26T11:19:08Z

Institor: загружена новая версия «Файл:GDAL VRT.png»: Диалог создания виртуального растра

Диалог создания виртуального растра

Файл:GDAL VRT.png

2013-06-26T11:17:55Z

Institor: Диалог создания виртуального растра

Диалог создания виртуального растра