Я установил GDAL с помощью установщика Osgeo. Как программно обрезать растровый слой векторным слоем? Есть ли какой-нибудь GDAL API, который может помочь мне в этом? Я использую Python.
26
Я не уверен , о библиотеке GDAL апи, есть void* GDALWarpOptions::hCutline
в Options Деформации ссылки из API учебника Деформации , но никаких явных примеров. Вы уверены, что вам нужен программный ответ? Утилиты командной строки могут сделать это из коробки:
ogrinfo
для определения степени отсечения шейп-файлаgdal_translate
для обрезки по форме экстентовgdalwarp
с -cutline
параметромШаги 2 и 3 для оптимизации, вы можете обойтись только с gdalwarp -cutline ...
.
См. Обрезка растров с помощью GDAL с использованием полигонов из Linfinity для решения на основе linux, заключенных в один скрипт. Еще один пример вырезки можно увидеть в учебнике Майкла Кори, посвященном созданию холмов для Mapnik .
Джоэл Лоухед из GeospatialPython имеет полный пример Python в растре Clip с использованием shapefile , хорошо написанного учебника. Вам нужно будет установить библиотеку изображений Python (PIL), которая не включена в Osgeo4W (для этого вам может понадобиться добавить o4w python в реестр Windows, чтобы заставить программу установки работать).
источник
Кажется, эта тема всегда возвращается. Я сам не знал, что GDAL> 1.8 настолько продвинут, что уже дает вам честную обработку командной строки для выполнения этой задачи.
Комментарий от Mike Toews довольно полезен, но вы можете просто сделать, например:
Вы можете обернуть эту команду внутри скрипта Python с помощью превосходного модуля подпроцесса .
Одна вещь, которая была для меня действительно проблематичной, заключается в том, что мне нужно было найти минимальное решение этой проблемы, то есть максимально простое и не требующее многих внешних зависимостей. Использование Python Imaging Library, как в учебном пособии Джоэла Лоухеда, очень удобно, но я нашел следующее решение: использование масок Numpy.
Я не знаю, лучше ли это, но это было то, что я знал назад, чем (3 года назад ...).
Первоначально я создал допустимую область данных внутри исходного растра (например, экстент выходного растра, где тот же), но мне понравилась идея сделать растр также меньшим (например, -crop_to_cutline), поэтому я принял
world2Pixel
от Джоэла Лоухеда. Вот мое собственное решение:для полного описания
class MaskRaster
и его методов, смотрите github моего проекта .Используя этот код, вам все равно нужно будет использовать GDAL. Однако в будущем планируется использовать чистый Python, где я могу, потому что предполагаемая аудитория моего программного обеспечения испытывает трудности с слишком большим количеством зависимостей (я использую Debian для разработки программного обеспечения, а клиенты используют Windows 7 ...).
источник