Упражнение 7 Пространственные и атрибутивные запросы

Архив с исходными данными

Контрольный лист

7.1 Введение

Цель задания — знакомство с редактированием векторных пространственных данных, элементами базовых технологий ГИС (оверлей, пространственный запрос, атрибутивный запрос).

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

Необходимая практическая подготовка: Знание основных компонент интерфейса QGIS (менеджер источников данных, панель слоёв, фрейм карты). Добавление источников пространственных данных в проект. Настройка символики и подписей объектов. Привязка растра. Создание макета, добавление карты и зарамочного оформления, экспорт макета.

Исходные данные: Схема районов города с сайта НЕАГЕНТ по состоянию на 2021 г.; набор пространственных данных о зданиях и сооружениях, созданный на основе базы OpenStreetMap

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

7.1.1 Контрольный лист

  • Привязать растровую схему районов города
  • Создать набор пространственных данных о районах путем цифрования растровой карты
  • Добавить семантическую информацию о районах
  • Выполнить серию пространственных и атрибутивных запросов для определения структуры застройки в пределах районов
  • Построить картодиаграммы по полученным значениям
  • Подготовить проект карты с компоновкой

7.2 Подготовка исходных данных

В начало упражнения ⇡

  1. Загрузите архив с исходными данными и распакуйте его в вашу рабочую директорию.

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

  3. Привяжите схему районов. Если вы не знаете, как привязать схему, сверьтесь с упражнением 5

Снимок экрана №1. Окно QGIS с привязанным растровым изображением.

7.3 Создание набора пространственных данных

В начало упражнения ⇡

На этом этапе мы создадим новый векторный набор пространственных данных, в котором будут храниться контура районов. QGIS поддерживает множество форматов пространственных данных, выбор конкретного формата в каждом случае определяется задачами, требованиями и ограничениями проекта. В качестве формата по умолчанию в QGIS 3 используется GeoPackage.

Файл GeoPackage представляет собой базу данных SQLite с дополнительными функциями для хранения наборов пространственных данных. Вся необходимая информация хранится в единственном файле *.gpkg; в процессе работы создаются временные файлы. В одном файле может быть сохранено несколько наборов пространственных данных (иногда их также называют «таблицами», tables, или «слоями», layers)

  1. Запустите процедуру создания нового набора данных GeoPackage. Для этого нажмите кнопку на панели менеджера источников данных или сочетание клавиш Ctrl+Shift+N.

  2. В открывшемся окне введите следующие параметры:

    • База данных: сохраните файл *.gpkg в директорию вашего проекта. Назовите его по шаблону Districts_%Фамилия%.gpkg;
    • Имя таблицы: примите значение, предлагаемое по умолчанию после указания пути к базе;
    • Тип геометрии: полигон (площадная)
    • Чекбоксы «Include Z dimension» и «Include M values» оставьте выключенными;
    • Система координат: такая же, как система координат проекта.
  3. Помимо геометрии объектов, нам нужно будет хранить семантическую информацию — номера районов. Создайте новое поле целочисленного типа с названием district_id, как показано на рисунке ниже.

  1. Проверьте правильность заполнения параметров. Если всё заполнено правильно, нажмите ОК. Новый слой будет добавлен в проект.

7.4 Векторизация районов

В начало упражнения ⇡

  1. Измените стиль слоя на pattern dot blue. Эта символика будет удобна для цифрования, поскольку позволит не «скрывать» растровый слой под создаваемыми объектами и в то же время будет контрастно выделяться на их фоне.

  2. Теперь всё готово к редактированию слоя. Режим редактирования запускается при помощи кнопки или из контекстного меню слоя. Запустите режим редактирования.

  3. Когда режим редактирования запущен, нажмите кнопку , чтобы начать добавление нового объекта

  4. Увеличьте изображение до района, обозначенного цифрой «1». Вы начнёте векторизацию с него. По нажатию левой кнопки мыши устанавливается положение первого узла создаваемого контура. «Обходите» контур по часовой стрелке или против часовой стрелки, устанавливая новые узлы нажатием левой кнопки мыши, как показано на изображении ниже

  1. Чтобы завершить редактирование объекта, нажмите правую кнопку мыши.

  2. Когда вы завершаете создание объекта в QGIS, появляется окно ввода значений атрибутов. Введите номер района, как показано ниже

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

  1. Найдите панель инструментов прилипания

  1. Включите режим прилипания, нажав на иконку с изображением магнита

  2. Нажимая на соседние иконки в панели, включите прилипание к объектам во всех слоях, но только к вершинам объектов

  3. Установите порог прилипания равным 12 пикселей экрана.

Теперь, если вы будете подводить курсор к существующим объектам, его положение будет автоматически «притягиваться» к ближайшим узлам

  1. Прилипание позволяет реализовать создание новых объектов путём трассировки по существующим объектам. Включите режим трассировки, нажав на кнопку

  2. Оцифруйте соседний объект с помощью трассировки. Для этого подведите курсор к углу существующего объекта и установите первый узел. Затем переместите курсор вдоль общей границы и установите следующий узел. Промежуточные узлы будут добавлены автоматически. После оцифровки общей границы продолжайте векторизацию обычным способом. После добавления нового объекта не забудьте внести его атрибуты во всплывающем окне.

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

Важно! Время от времени сохраняйте ваши правки при помощи кнопки на панели редактирования. В QGIS, как и во всех ГИС-пакетах, сохранение правок в данных выполняется отдельно от сохранения проекта.

Несколько советов по редактированию данных:

Чтобы отменить последний установленный узел, нажмите Backspace;

Чтобы полностью удалить последнюю оцифрованную фигуру, нажмите Ctrl+Z.

Чтобы изменить положение отдельного узла уже созданного объекта, можно воспользоваться инструментом редактирования узлов .

Трассировка в QGIS работает только тогда, когда весь путь трассировки отображается на экране. Если часть сегмента вылезла за пределы фрейма, трассировка не сработает.

Изучите панель «Дополнительные инструменты оцифровки» — она может вам пригодиться

Плавные кривые линии (берега, реки и т.п.) при цифровании обычно аппроксимируются ломаными линиями с относительно большим числом узлов — так, чтобы отличия не были визуально заметны. Хотя многие форматы данных позволяют создавать криволинейные объекты, этот подход не являются общепринятым.

  1. По окончании цифрования объектов отключите растровый слой и отобразите данные в полном охвате. Сделайте снимок экрана.

Снимок экрана №2. Окно QGIS с оцифрованными границами районов.

7.5 Пространственные и атрибутивные запросы

В начало упражнения ⇡

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

Для подсчёта числа зданий вы воспользуетесь следующим алгоритмом:

  • Выбрать район на карте

  • Выбрать здания, попадающие в его пределы (пространственный запрос).

  • Из полученной выборки оставить только здания, принадлежащие к определённому типу (атрибутивный запрос).

  • Записать число отобранных зданий в соответствующий атрибут текущего района.

Операции повторяются для каждого района.

  1. Добавьте на карту набор данных о зданиях, который вы загрузили в начале занятия.

  2. Откройте таблицу атрибутов слоя зданий. Для этого нажмите правой кнопкой мыши на слой в панели слоёв и выберите опцию «Открыть таблицу атрибутов» в контекстном меню.

Таблица атрибутов — это представление базы данных в виде таблицы. В большинстве случаев работа с семантической информацией в ГИС выполняется через таблицу атрибутов.

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

  1. Закройте таблицу атрибутов слоя зданий и откройте таблицу атрибутов слоя районов.

  2. Убедитесь, что для слоя районов включён режим редактирования.

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

  4. Добавьте поле с именем number_of_houses целочисленного типа.

На этом подготовительные операции закончены, мы переходим к запросам.

  1. Выберите любой из районов: выделите строку в таблице атрибутов или воспользуйтесь инструментом выборки в окне карты . Выделенный объект будет подсвечен жёлтым цветом на карте и синим цветом в таблице атрибутов.

  1. Теперь выберите те объекты из слоя зданий, которые находятся внутри выбранного района. Для этого запустите инструмент «Вектор» — «Выбор» — «Выбрать по расположению…». Используйте условие «находятся внутри» (are within). Включите опцию «Использовать только выделенные объекты» для выбирающего слоя.

Нажмите «Выполнить», чтобы выбрать объекты в слое зданий.

  1. Имея выборку в слое зданий и не закрывая окно пространственного запроса («Выбрать по расположению…»), осуществите выборку по атрибутам. Для этого выберите в таблице слоёв слой зданий и нажмите кнопку «Выделить объекты, удовлетворяющие условию» (или нажмите Ctrl+F3. Откроется форма ввода атрибутивного запроса.

Атрибутивные запросы в ГИС, как правило, создаются с использованием диалектов языка SQL. Само выражение представляет собой только условие (where clause) и часто использует значения атрибутов. В QGIS можно выполнить запрос и без составления выражения на языке SQL, но поскольку SQL более функционален, мы призываем вас сразу начать пользоваться им.

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

Атрибутивный запрос представляет собой логическое выражение, которое применяется к кортежу базы данных (строке таблицы). Результатом вычисления логического выражения является либо логическая 1 (ИСТИНА), либо логический 0 (ЛОЖЬ). В выборку попадают объекты, для которых при вычислении выражения получено значение 1.

В простейшем случае выражение состоит из проверки равенства. Например, выражение "type" = 'house' будет истинным для всех объектов, у которых в поле type записано значение house. Можно конструировать более сложные логические выражения с использованием операторов AND (И), OR (ИЛИ) и других. Например, выражение "type" = 'apartments' AND "name" IS NOT NULL будет истинным для объектов, у которых в поле type записано значение apartments (т.е. многоэтажные жилые дома) и одновременно в поле name записано какое-либо значение, отличное от пустого (NULL) (т.е. имеющие название). А запрос "type" = 'garage' OR "type" = 'garages' выберет из слоя зданий объекты, у которых в поле type записано либо значение garage, либо значение garages.

Примечание 1: в интерфейсе атрибутивного запроса и во всех остальных инструментах QGIS, использующих тот же интерфейс, названия полей обрамляются двойными кавычками (“), а строковые переменные — одинарными кавычками (’) Примечание 2: в QGIS и в большинстве систем управления базами данных (СУБД) для работы с пустыми значениями (NULL) используются отдельные операторы. Например, выражение "type" = NULL будет воспринято как некорректное, вместо него следует писать "type" IS NULL

  1. Составьте выражение, реализующее следующее условие:

Выбрать объекты, у которых атрибуту type присвоено значение apartments или house

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

  1. Мы сформулировали условие для выборки, однако нам нужно выбрать объекты не просто из слоя, а из уже существующей выборки. Для этого откройте выпадающий список кнопки «Выбрать объекты» и выберите функцию «Фильтровать текущее выделение»:

  1. После применения фильтра появится всплывающее сообщение с числом выбранных объектов. Эта информация будет продублирована внизу окна QGIS:

  1. Введите полученную цифру в таблицу атрибутов слоя районов в строке, соответствующей выбранному объекту (подсвечена синим).

  2. Теперь, не закрывая окна пространственного и атрибутивного запроса, выбирайте по очереди каждый следующий объект в слое границ районов и повторяйте для них шаги 3-13. Таким образом вы заполните весь столбец number_of_houses в таблице атрибутов. В процессе выполнения не забывайте периодически сохранять правки.

  3. Скопируйте таблицу атрибутов в любой табличный процессор (Microsoft Excel, Google Sheets, LibreOffice Calc). Для этого при помощи сочетания клавиш Ctrl+A выделите все записи в таблице, скопируйте при помощи Ctrl+C и вставьте записи без форматирования в табличный процессор при помощи Ctrl+Shift+V.

  4. Удалите столбец wkt_geometry, если он есть в вашей таблице. Скопируйте остальные столбцы и вставьте их в отчётный файл.

Таблица №1. Таблица атрибутов слоя районов.

  1. Завершите редактирование слоя районов.

  2. Отключите отображение слоя зданий.

7.6 Визуализация слоя при помощи картодиаграмм

В начало упражнения ⇡

Как вы уже знаете из курса картографии, не существует прямого взаимно-однозначного соответствия между способами изображения, принятыми в картографии, и способами визуализации данных в ГИС. Чтобы отобразить абсолютную величину в пределах единиц статистического учёта (или аналогичных контуров), в картографии применяют способ картодиаграмм. Для создания условных знаков, соответствующих способу картодиаграмм, в QGIS применяется вкладка «Диаграммы» в свойствах слоя.

  1. Откройте свойства слоя районов и перейдите на вкладку «Диаграммы». По умолчанию отображение диаграмм отключено. Включите отображение круговых диаграмм, выбрав опцию Pie Chart из выпадающего списка вверху.

  2. В настройках диаграмм перейдите на вкладку «Атрибуты». Добавьте атрибут number_of_houses к диаграмме. Задайте для него мягкий оттенок оранжевого или жёлтого цвета.

    Примечание: вкладка «Атрибуты» нужна для настройки секторов круговых диаграмм. Если вы добавляете несколько атрибутов к диаграмме, программа автоматически рассчитывает размеры секторов для каждого объекта. Сейчас перед нами не стоит такая задача, так что, фактически, можно было не добавлять атрибуты или добавить любой другой.

  3. На вкладке «Отрисовка» («Рендеринг») установите для диаграмм настройку прозрачности (70 %)

  4. На вкладке «Размер» измените способ задания размера с фиксированного на изменяющийся («масштабируемый»).

    Изменение размера диаграмм в QGIS работает следующим образом. Пользователь задаёт максимальное значение показателя, который будет управлять размером диаграммы, и соответствующий ему максимальный диаметр диаграммы. Размер круга масштабируется пропорционально величине показателя. Из картографических соображений следует всегда выбирать масштабирование площади (в новых версиях QGIS некорректно переведено как «к поверхности»), а не линейного размера (диаметра).

  5. В выпадающем списке «Атрибут» выберите нужный атрибут и нажмите кнопку «Найти, чтобы отобразить его максимальное значение. Оно будет отражено во вкладке «Макс. значение».

  1. Округлите полученное максимальное значение в большую сторону до величины, кратной 100.

  2. Установите максимальный размер диаграмм равным 14 (мм). Диаграмма такого диаметра будет соответстовать объекту с указанным максимальным значением атрибута.

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

  3. Перейдите на вкладку «Размещение» и проверьте, что для ваших диаграмм указана опция размещения над центроидом.

  4. Перейдите на вкладку «Легенда». Здесь вы настроите комбинацию условных знаков для отображения в легенде. Нажмите кнопку «Условные обозначения для размеров диаграмм» (Show Legend Entries for Diagram Size).

  5. Настройте отображение значков следующим образом:

    • «Коллапсируйте» значки легенды;
    • В качестве символа используйте белый маркер круглой формы с тёмно-серой обводкой;
    • Задайте заголовок («Число домов»);

Интерфейс настройки диаграммы должен принять приблизительно следующий вид:

Важное замечание: такой набор значков не является картографически корректным для легенды к абсолютной непрерывной шкале значков. По состоянию на февраль 2022 г. ни один широко используемый ГИС-пакет не может сделать легенду к размерам кругов картодиаграммы согласно принятым картографическим правилам. Легенды к таким картам следует составлять или исправлять вручную.

По окончании настройки символов окно QGIS будет выглядеть следующим образом:

  1. Самостоятельно измените настройки отображения слоёв так, чтобы базовая карта не «доминировала» над тематическим содержанием, и сделайте снимок экрана.

Снимок экрана №3. Окно QGIS с настроенной символикой.

7.7 Создание и экспорт макета компоновки

В начало упражнения ⇡

Изображение, которое вы видите во фрейме данных, можно экспортировать «как есть» (с помощью опции «Проекты» — «Импорт/экспорт» — «Экспортировать карту как изображение…»). Однако для картографических целей, как правило, формируется компоновка карты. На листе заданного формата размещается картографическое изображение, добавляется название, легенда, масштабная линейка и элементы зарамочного оформления.

Сейчас мы создадим макет компоновки с расчётом на то, что итоговая карта-схема будет вставлена в отчёт.

  1. Создайте новый макет компоновки («Проект» — «Создать Макет…») или Ctrl+P.

  2. В качестве названия макета введите свою фамилию.

В отличие от проектов QGIS и наборов пространственных данных, на названия макетов компоновки в QGIS не накладывается ограничений.

После ввода названия откроется окно компоновки (Layout)

Layout
Layout
  1. Измените ориентацию страницы с альбомной, предлагаемой по умолчанию, на портретную, соответствующую ориентации страниц отчёта. На изображении ниже показано, как это сделать.
Layout
Layout
  1. Добавьте картографическое изображение (фрейм карты) на макет
Add map
Add map
  1. Измените положение и размер фрейма карты таким образом, чтобы на листе оставались поля: левое: 3 см, правое: 1,5 см.

  2. Изменяя высоту фрейма, масштаб изображения, а также перемещая картографическое изображение внутри фрейма при помощи кнопки «Переместить содержимое элемента» , добейтесь вида фрейма, аналогичного рисунку ниже:

Map layout
Map layout
  1. Не закрывая окно макета, вернитесь в основное окно QGIS и переименуйте слой районов. Назовите его «Районы Новосибирска».

  2. Добавьте на макет легенду с помощью кнопки . Поместите легенду в наиболее свободное место в пределах карты.

  3. Отключите автообновление легенды и удалите ненужные записи, как показано на рисунке ниже:

  1. В свойствах элемента легенды найдите вкладку отступ и удалите (установите равным нулю) отступ под заголовком легенды.

  2. Добавьте масштабную линейку с помощью кнопки . Разместите масштабную линейку в юго-западном углу изображения. При необходимости измените обозначение единиц измерения.

  3. Добавьте название карты (текстовый элемент) с помощью кнопки . Самостоятельно выберите оптимальное место для размещения названия.

  4. Экспортируйте карту-схему («Макет» — «Экспорт в Изображение» или кнопка ) с разрешением 96 точек на дюйм (dpi). Используйте опцию «Кадрировать по содержимому» (Crop to Content).

    Примечание: разрешение 96 точек на дюйм считается довольно низким для картографических целей. Изображения с таким разрешением не годятся для печати, но иногда могут быть пригодны для размещения в Интернете. ВЫ снижаете разрешение для того, чтобы сохранить читаемость базовой карты.

  5. Вставьте полученное изображение в отчётный файл.

Итоговое картографическое изобаржение

  1. Ответьте на контрольные вопросы в конце отчётного файла