有没有一种灵活的方法可以通过从shapefile中读取或使用一组顶点创建的多边形裁剪栅格？

Question

我正在尝试用python中的多边形剪辑栅格（ascii或tif文件）。可以从shapefile中读取多边形，也可以通过一组顶点group.getChildren();来创建多边形。

我知道如何使用POLYGON((0 0, 1 0, 1 1 , 0 1, 0 0))裁剪带有参数gdal.Warp的栅格，该参数代表shapefile名称。

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

Answer 1

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

{"type":"Polygon","coordinates":[[[0,0],[1,0],[1,1],[0,1],[0,0]]]}。

您可以使用http://kostas.krevatas.net/wkt-convertor/

之类的在线转换器

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

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