Интерактивная карта в QGIS. Часть 1: обработка данных и автоматизация процессов

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

Ранее была опубликована интерактивная карта изменения береговой линии Симферопольского водохранилища, позволяющая отследить изменения внешнего вида водоема за 20 лет. Теперь рассмотрим процесс создания карты в QGIS. А процесс ее публикации в интернете будет описан во второй части.

Традиционно с сайта геологической службы США скачаны космические снимки рассматриваемой территории за последние 20 лет. Один снимок – один год. Период создания снимков апрель—октябрь. Ввиду сложившейся сегодня обстановки с водой в Крыму снимков 2020 года взято четыре. А снимки 2008 и 2012 года по техническим причинам включены не были.

Теперь предстоит обработка каждого из этих снимков. Добавляем в проект слои 4, 5 и 1 (в Landsat 8 слои 5, 6 и 2) и действуем по алгоритму:

Объединение слоев в указанной последовательности (Растр – Прочее – Результат/Merge (объединение)); располагаем каждый входной слой отдельно

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

Растр – Преобразование – RGB в PCT, с количеством цветов нужно экспериментировать, но зачастую нужный вариант находится в диапазоне 5–10

Вместо преобразования RGB в PCT можно использовать классификацию K-Means в модуле OrfeoToolBox. Результат идентичен.

Результат преобразования

Убираем в преобразованном слое все цвета, кроме того, который обозначает само водохранилище (водные объекты); видно, что контуры преобразованного слоя соответствуют контурам на снимке

Далее преобразованный слой переводим в векторный формат. Растр – Преобразование – Создание полигонов (растр в вектор).

Выделяем нужный объект и сохраняем в отдельный слой

Сравниваем новый слой с реальными контурами водохранилища

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

Результат упрощения, контуры пикселей исчезли

Как можно видеть, алгоритм длительный, поэтому решено было его автоматизировать хотя бы наполовину – к результату векторизации.

В нижней части расположенной слева панели будет две вкладки Inputs и Algorithms; в первой вкладке выбираются входные данные (слои), во второй – действия с ними

Во вкладке "Анализ данных" есть инструмент Graphical Modeler/Редактор моделей. Тогда откроется окно, в котором можно творить.

Промежуточные результаты в отдельный слой не выводим. Они не нужны. Слой true_raster (обрезанный снимок) вывели для сравнения результата со снимком.

Когда модель создана и сохранена, она появится во вкладке "Модели" в панели "Инструменты анализа".

Остается выбрать входные слои и нажать "выполнить"

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

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

Процесс создания html-страницы, содержащей карту, и публикация ее в интернете – во второй части статьи.

Читайте нас в Telegram и ВКонтакте