剪切、提取到掩码和栅格子集是常见的 GIS 操作。使用 Python 的 gdal 软件包可以非常容易地实现这些过程的自动化。本教程将演示如何使用 gdal Python API 中的 Warp() 函数将栅格剪切到指定范围并将栅格剪切到多边形图层。这些操作通常也被称为掩码提取或提取到掩码。

安装 GDAL

您需要gdal在 Python 环境中安装才能完成本教程。不幸的是，gdal这并不是最容易安装的包，特别是如果您使用非 anaconda Python 发行版。在我的其他博客中也介绍了如何安装gdal，请自行参阅

输入数据

现在我们需要一些输入数据来处理。我将使用来自 USGS 的 1 m 激光雷达衍生 DEM。我还将使用我创建的与光栅图像的一部分重叠的多边形。您可以自行下载并创建这些数据。