python - 是否有一种灵活的方法可以通过从形状文件读取或使用一组顶点创建的多边形来裁剪光栅？

Question

我正在尝试在 python 中用多边形剪辑光栅(ascii 或 tif 文件)。多边形可以从 shapefile 中读取或由一组顶点创建，例如 POLYGON((0 0, 1 0, 1 1 , 0 1, 0 0)) .

我知道如何使用 gdal.Warp使用参数 cutlineDSName 裁剪栅格，代表 shapefile 名称。

但是如果我想使用由一组点定义的简单多边形来切割栅格该怎么办？

Answer 1

一种解决方案是将多边形 POLYGON((0 0, 1 0, 1 1 , 0 1, 0 0)) 转换为 geojson:

{"类型":"多边形","坐标":[[[0,0],[1,0],[1,1],[0,1],[0,0] ]]}。

您可以使用在线转换器，例如 http://kostas.krevatas.net/wkt-convertor/

然后，将 cutlineDSName 设置为 geojson 文件形式的多边形，例如 polygon.json

...因为 geojson 是一种有效的输入格式，与 shapefile 相同。