Как обрезать карту в qgis
We will download some public domain raster data for Brazil and view it in QGIS. Next, we will merge these into a single mosaic and clip it using a country boundary to get a single seamless dataset for the country.
We need Brazil country boundary to clip our raster. You can get the Admin 0 — Countries shapefile from Natural Earth.
NASA/GSFC, Rapid Response site has a good collection of near real-time satellite imagery. A good regional product is USDA Foreign Agricultural Service (FAS) subsets.. We will use 2km resolution FAS subsets for Brazil for this tutorial.
Here is how to search and download the revelant data.
Open the South America region subsets. Find the Brazil subsets shown in the FAS Subsets section. Click on any one of them.
In the details page, click the 2km link under the product of your choice. Here we will download the NDVI product. Learn more about NDVI.
Click the Download GeoTIFF file link to download the raster image.
Repeat the process for all 7 FAS subsets for Brazil.
For convenience, you can directly download sample data used in this tutorial from links below.
Open QGIS and go to Layer ‣ Add Raster Layer.. .
Browse to the directory with the individual images. Hold down the Ctrl key and click on the image files to make a multiple selection. Click Open .
You will see the images load up in the Table of Content on the left panel. Now let us create a single Mosaic image from all these individual images. Click on Raster ‣ Miscellaneous ‣ Merge .
The Raster menu in QGIS comes from a core plugin called GdalTools . If you do not see the Raster menu, enable the GdalTools plugin from Plugins ‣ Manage and install plugins ‣ Installed . See Использование модулей расширения for more details.
In the Merge dialog, click Select… next to Input files and browse to the directory containing all the individual geotiffs. Keep holding Ctrl key and select all the. subsets. Now click Select… next to Output file and name the output as Brazil_mosaic.tif . At the bottom, check the box next to Load into canvas when finished . Click OK .
You will get a pop-up message saying Processing complete , once the mosaic is created and loaded to the QGIS Canvas. You will see that the individual images and now combined and mosaiced into a single layer. You can now turn off individual layers by un-checking the box next to them.
Another Raster operation you can do is to subset or crop an image. We can use a polygon from a vector layer to crop the raster to the exact shape. Let’s load the country polygons shapefile we downloaded from Natural Earth. Go to Layer ‣ Add Vector Layer .
Select the ne_10m_admin_0_countries.zip file and click Open . When prompted to select the layer within the zip file, select ne_10m_admin_0_countries.shp .
Once the vector layer is loaded, we want to select and extract the polygon for Brazil. Select the Select Single Feature tool from the toolbar.
Click anywhere on the Brazil polygon and it will be selected.
Right-click the ne_10m_admin_0_countries layer and select Save Selection As… .
Name your output as brazil_boundary.shp and make sure Add saved file to map box is checked. click OK .
You will see the Brazil boundary polygon now loaded in QGIS.
Now go to Raster ‣ Extraction ‣ Clipper .
Select the input file (raster) as Brazil_mosaic . Name the Output file as Brazil_mosaic_clipped . In the Clipping mode section, choose Mask layer . Select the newly created brazil_boundary as the mask layer. Check the box next to Load into canvas when finished . Click OK .
The new cropped layer will be loaded into QGIS. You will notice the black pixels surrounding the actual mosaic. Let’s remove that. Right-click on the Brazil_mosaic_clipped layer and select Properties .
Go to the Transparency tab, and add 0 as an Additional no data value .
Now you have a nice mosaic cropped to a country boundary that you can use in your project as a background layer or do further analysis on.
ГИС-натуралист
Исходный растровый слой был размером 7651 х 7431 пикселей (113,8 Мб), результирующий слой – 461 х 443 пикселя (409 Кб). Дальнейшая работа с результирующим растром требует существенно меньше машинных ресурсов, нежели работа с исходным изображением.
2) Oбрезка изображения по маске
Обрезка изображения по маске подразумевает, что нам нужно фигурно вырезать часть геопривязанного изображения, по контуру заданного полигона – маске.
Для этого нам нужно иметь векторный полигональный слой с объектом, по границам которого необходимо обрезать изображение. При этом нужно чтоб изображение и векторный слой имели одну систему координат. В нашем примере это – WGS 84 / UTM zone 36N (EPSG:32636) – система координат исходного изображения Landsat.
Для облегчения восприятия изображения, создадим композитное изображение из 7-го, 6-го и 4-го каналов, и добавим его в проект QGIS. Также добавим в проект векторный слой, по которому будем проводить обрезку. Меню «Слой» → «Добавить слой» → «Добавить векторный слой» или жмем <Ctrl>+<Shift>+<V>, выбираем нужный нам файл.
Запускаем инструмент «Обрезка». Меню «Растр» → «Извлечение» → «Обрезка». В открывшемся окне выбираем растр, который нужно обрезать (можно выбрать как растр загруженный в QGIS, так и указать файл на диске), целевой файл – файл созданный в результате обрезки, режим обрезки «Слой маски» и выбираем векторный слой маски (можно выбрать как слой загруженный в QGIS, так и указать файл на диске), так же можем задать создание альфа-канала.
В результате получаем файл изображения обрезанного по маске векторного слоя.
Настраиваем прозрачность слоя и размещаем его в нужном нам контексте.
Если статья для тебя полезна, подписывайся на мои каналы и группы в соцсетях:
Если Вы хотите в совершенстве овладеть функционалом QGIS, приглашаю посетить страницу моего онлайн курса «QGIS от новичка до профессионала» и ознакомится с его программой.
В курс включен концентрат моего 9-ти летнего опыта работы с QGIS для решения самых разных задач.
7.14. Работа с растровыми данными¶
В этом разделе описывается вопросы визуализиции и определения свойств растрового слоя. QGIS использует библиотеку GDAL для чтения и записи растровых форматов. Растр поддержка GRASS поставляется нативной QGIS плагин поставщика данных. Растровые данные также могут быть загружены в режиме чтения из почтового индекса и GZIP архивов в QGIS.
Перечисленные операции выполняются модулем gdalTools. По умолчанию он установлен. Если вы не находите этих пунктов в меню, установите и включите этот модуль. Модуль предоставляет интерфейс к консольным утилитам gdal.
7.14.1. Проекции¶
Растровая операция запускается через меню Растр — пункт меню Проекции.
Перепроецирование¶
Операция запускается через меню Растр — пункт меню Проекция — Перепроецирование:
Рис. 7.280. Инструмент для работы с растровыми данными — Перепроецирование ¶
Пересохраняет растр в другую поддерживаемую проекцию. Есть возможность пакетной обработки. Полное описание параметров — в описании утилиты http://www.gdal.org/gdalwarp.html
Назначить проекцию¶
Операция запускается через меню Растр — пункт меню Проекция — Назначить проекцию:
Рис. 7.281. Инструмент для работы с растровыми данными — Назначить проекцию ¶
Привязать проекцию к растру, если он уже геопривязан, но отсутствует информация о проекции. Также с помощью этого инструмента можно поменять текущую проекцию. Есть возможность пакетной обработки.
Извлечение проекции¶
Операция запускается через меню Растр — пункт меню Проекция — Извлечь проекцию:
Рис. 7.282. Инструмент для работы с растровыми данными — Извлечь проекцию ¶
Создаёт wld-файл или prj-файл для заданного растра. Есть возможность пакетной обработки.
7.14.2. Преобразование¶
Операция запускается через меню Растр — пункт меню Преобразование.
Растеризация (вектор в растр)¶
Операция запускается через меню Растр — пункт меню Преобразование — Растеризация (вектор в растр):
Рис. 7.283. Инструмент для работы с растровыми данными — Растеризация (вектор в растр). ¶
Используется для переноса векторных данных (точки, линии, полигоны) в растровый файл. Векторные данные берутся из OGR-совместимого источника данных, и должны быть в той же системе координат, что и растр. Перепроецирование «на лету» не предусмотрено (на базе gdal_rasterize). Значения растра берутся из заданого числового атрибута.
Для примера покажем генерацию индексированного растра с данными землепользования из Openstreetmap.
В Wizzard вводим запрос «landuse=*». Генерируется запрос, который выдаёт полигоны с тегом landuse
Выполняем запрос, сохраняем результат в GeoJSON.
Открываем geojson в NextGIS QGIS.
Генерируем в векторном слое числовое поле, значение которого запишется в растр. В калькуляторе полей выбираем создание нового поля, и вводим выражение
« CASE WHEN «landuse»=“residential“ THEN 10 WHEN «landuse»=“industrial“ THEN 20 WHEN «landuse»=“farmland“ THEN 30 END «
Запускаем инструмент Растр — Преобразование — Растеризация (вектор в растр).
Этот инструмент создаёт GeoTIFF.
Векторизация (растр в вектор)¶
Операция запускается через меню Растр — пункт меню Преобразование — Векторизация (растр в векторы):
Рис. 7.284. Инструмент для работы с растровыми данными — Векторизация (растр в векторы). ¶
Используется для создание векторных полигонов для всех соседних пикселей растра, имеющих один и тот же цвет (значение). Для каждого полигона в атрибутивную таблицу записывается соответствующее значение растра. Используя растровую маску, можно задать область для обработки (на базе gdal_polygonize).
Преобразовать формат¶
Операция запускается через меню Растр — пункт меню Преобразование — Преобразовать формат:
Рис. 7.285. Инструмент для работы с растровыми данными — Преобразовать формат. ¶
Конвертация растровых файлов между форматами. Возможно изменение системы координат, сжатие, задание значения NO_DATA.
Есть возможность пакетной обработки.
Эта операция подойдёт вам если нужно преобразовать GeoTIFF в JPEG (не в TIFF c сжатием JPEG, а в файл с расширением jpg, в диалоге экспорта растрового слоя можно сохранять только в GeoTIFF). Так же тут можно сжать GeoTIFF в JPEG.
RGB-изображение в PCT¶
Операция запускается через меню Растр — пункт меню Преобразование — Преобразовать RGB-изображение в индексированное:
Рис. 7.286. Инструмент для работы с растровыми данными — Преобразовать RGB-изображение в индексированное. ¶
Преобразование из цветового пространства RGB в индексированные цвета. Производится расчёт оптимального соответствия цветов, либо можно задать вручную количество вариантов цвета. Выполняется dithering. , работает пофайлово и в пакетном режиме. Полное описание — на https://www.gdal.org/rgb2pct.html
PCT в RGB-изображение¶
Операция запускается через меню Растр — пункт меню Преобразование — Преобразовать индексированное изображение в RGB-изображение:
Рис. 7.287. Инструмент для работы с растровыми данными — Преобразовать индексированное изображение в RGB-изображение. ¶
Преобразование растра с индексированными цветами в цветовое пространство RGB. Работает пофайлово и в пакетном режиме (использует pct2rgb). Полное описание — наhttps://www.gdal.org/pct2rgb.html
7.14.3. Извлечение¶
Запускается через меню Растр — пункт меню Извлечение.
Создать изолинии¶
Построение изолиний. Генерация изолиний по цифровой модели рельефа (DEM) (на базе gdal_contour, см https://www.gdal.org/gdal_contour.html).
Рис. 7.288. Исходный растр (одноканальное псевдоцветное изображение) ¶
Рис. 7.289. Полученные изолинии ¶
Операция запускается через меню Растр — пункт меню Извлечение — Создать изолинии:
Рис. 7.290. Параметры создания изолиний из растра ¶
Обрезать растр по охвату¶
Операция запускается через меню Растр — пункт меню Извлечение — Обрезка.
В полне «Охват обрезки» можно вручную задать прямоугольник обрезки при помощи координат.
Рядом с полем «Охват обрезки» находится выпадающее меню, в котором можно выбрать один из способов обрезки растра:
Текущий охват карты Обрезает растр по границам видимой в окне NG QGIS части карты. Этот же вариант включается кнопкой справа от поля.
Рис. 7.291. Обрезка по охвату. Стрелочкой указана кнопка «Текущий охват карты» ¶
Указать на карте Позволяет выделить на карте произвольный прямоугольник, по которому будет произведена обрезка. Выделение прямоугольником имеет смысл, только если в окне QGIS выставлена та же проекция, что у самого растра.
Обрезать растр по маске¶
Для обрезки по маске нужен полигональный слой в той же системе координат, что и растр. В полигоне могут быть дыры и части. Из слоя маски будет использован только первый объект. Если вы не уверены, нужно ли вам включать генерацию альфа-канала, то включите её. Полное описание параметров — в описании утилиты http://www.gdal.org/gdal_translate.html
Если вам нужно обрезать растр, например космоснимок по сложной границе, то нужно:
Узнать систему координат растрового файла: Свойства слоя –> Общие, посмотреть какая система координат подхватилась из файла.
Нарисовать слой: Слой –> Создать временный слой. Выберите ту же систему координат что и у растра, тип геометрии — мультиполигон. Нарисуйте маску.
Сохраните слой маски в формат ESRI Shapefile с той же системой координат, что и растр. (geopackage он не может использовать)
Запустите Растр –> Извлечение –> Обрезка, выберите растр, Слой обрезки — ваш новый слой, включите «Создать Альфа-канал» и «Охват целевого слоя по линии обрезки»
Файл после обрезки может получиться слишком большого размера на диске. В этом окне нельзя задавать опции сохранения и сжатия, поэтому в таком случае вам нужно или запустить обрезку из консоли с указанием -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR , либо потом обработать получившийся растр операцией «Преобразовать формат».
7.14.4. Анализ¶
Операция запускается через меню Растр — пункт меню Анализ.
Отсеивание¶
Операция запускается через меню Растр — пункт меню Анализ — Отсеивание:
Рис. 7.292. Инструмент для работы с растровыми данными — Отсеивание. ¶
Находит растровые полигоны с размером меньше заданного порога (в пикселях) и заменяет их на значение, соответствующее наибольшему соседнему растровому полигону. Результат отсеивания может быть записан как в существующий, так и в новый файл (использует gdal_sieve).
Сбросить в черный¶
Операция запускается через меню Растр — пункт меню Сбросить в черный:
Рис. 7.293. Инструмент для работы с растровыми данными — Сбросить в черный. ¶
Cканирование изображение и замена почти черных (или почти белых) пикселей возле рамки на пиксели с точным значением цвета. Часто используется для «коррекции» аэрофотоснимков, сжатых с потерями, чтобы можно было задать прозрачный цвет при создании мозаики (использует nearblack).
Заполнение пустот¶
Операция запускается через меню Растр — пункт меню Заполнение пустот:
Рис. 7.294. Инструмент для работы с растровыми данными — Заполнение пустот. ¶
Находит пустые растровые полигоны (обычно области NODATA) и заполняет их интерполяцией значений от краёв этого полигона.
Карта близости (расстояния в растре)¶
Операция запускается через меню Растр — пункт меню Карта близости (расстояния в растре):
Рис. 7.295. Инструмент для работы с растровыми данными — Карта близости (расстояния в растре). ¶
Строит карту близости растра, в которой указаны расстояния от центра каждого пикселя к центру ближайшего целевого пикселя. Целевыми пикселями будут все пиксели исходного растра, значения которых попадают в набор указанных величин (использует gdal_proximity).
Сетка (интерполяция)¶
Операция запускается через меню Растр — пункт меню Сетка (интерполяция):
Рис. 7.296. Инструмент для работы с растровыми данными — Сетка (интерполяция). ¶
Создает регулярную сетку (растр) на основе рассеяных данных, полученных из OGR-совместимого источника. Исходные даные будут интерполированы одним из доступных методов для получения значений узлов (на базе gdal_grid).
DEM (Анализ рельефа)¶
Операция запускается через меню Растр — пункт меню Анализ рельефа:
Рис. 7.297. Инструмент для работы с растровыми данными — Анализ рельефа. ¶
Создаёт новый растр на основе имеющегося растра с цифровой моделью рельефа (DEM). Может строить:
Индекс пересечённости (TRI).
Индекс превышения (TPI).
Обращайте внимание на системы координат ваших файлов с ЦМР при работе с этим инструментом. Файлы с SRTM или ASTER распространяются в EPSG:4326, с единицами измерения координат в градусах, а высота у них записана в метрах, и может быть что некоторые алгоритмы, например расчёта угла уклонов, выдадут вам неверные значения. Тогда нужно перепроецировать ЦМР во что-нибудь с метрами, например WGS 84/UTM Zone…
7.14.5. Прочее¶
Операция запускается через меню Растр — пункт меню Прочее.
Создать виртуальный растр (каталог)¶
Операция запускается через меню Растр — пункт меню Прочее — Создать виртуальный растр (каталог):
Рис. 7.298. Инструмент для работы с растровыми данными — Создать виртуальный растр (каталог). ¶
Создаёт файл VRT — в нём находятся ссылки на отдельные растровые файлы, а сам файл VRT используется как один растровый слой.
Объединение¶
Операция запускается через меню Растр — пункт меню Прочее — Объединение:
Рис. 7.299. Инструмент для работы с растровыми данными — Объединение. ¶
Склеивает несколько растровых файлов в один. Требует указания файлов. При настройке Склеить поканально, создает один растровый файл, где каждый исходный файл будет отдельным слоем.
Информация¶
Операция запускается через меню Растр — пункт меню Прочее — Информация:
Рис. 7.300. Инструмент для работы с растровыми данными — Информация. ¶
Выводит на экран вывод утилиты gdalinfo для заданного слоя. В этой информации пишется система координат и охват слоя.
Построить пирамиды¶
Операция запускается через меню Растр — пункт меню Прочее — Построить пирамиды:
Рис. 7.301. Инструмент для работы с растровыми данными — Построить пирамиды. ¶
Используется для создания или восстановления уменьшенных копий изображения (пирамид). Наличие пирамид несколько увеличивает скорость отрисовки растра. Может работать пофайлово и в пакетном режиме, использует gdaladdo. То же самое, что построение пирамид в настройках растрового слоя, но может работать пакетно.
Индекс мозаики растров¶
Операция запускается через меню Растр — пункт меню Прочее — Индекс мозаики растров:
Рис. 7.302. Инструмент для работы с растровыми данными — Индекс мозаики растров. ¶
GDAL extraction¶
Clips any GDAL-supported raster file to a given extent. The algorithm is derived from the GDAL grid utility .
Parameters¶
Defines a value that should be inserted for the nodata values in the output raster.
Clipping extent [extent]
Extent that should be used for the output raster. Selects only features within the specified bounding box.
Additional creation parameters [string]
<put parameter description here>
Outputs¶
Console usage¶
See also¶
Обрезать растр по маске¶
Description¶
Clips any GDAL-supported raster by a vector mask layer. The algorithm is derived from the GDAL grid utility .
Parameters¶
Defines a value for the no-data areas of the input raster.
Create and output alpha band [boolean]
Creates an alpha band for the result. The alpha band then includes the transparency values of the pixels.
Keep resolution of output raster [boolean]
The resolution of the output raster will not be changed.
Additional creation parameters [string]
<put parameter description here>
Outputs¶
Console usage¶
See also¶
Contour¶
Description¶
Extracts contour lines from any GDAL-supported elevation raster. The algorithm is derived from the GDAL contour utility .
Parameters¶
Defines the interval between the contour lines in the given units of the elevation raster.
Attribute name (if not set, no elevation attribute is attached) [string]
Defines the attribute name for the field containing the values of the contour lines. If nothing is defines the default name will be ‘ELEV’.
