gdal 多边形裁剪

在GDAL中，你可以使用多边形来裁剪栅格数据。以下是使用GDAL Python绑定库进行多边形裁剪的一般步骤：

1. 导入必要的库和模块：

import ogr
import gdal
import osr

2. 打开待裁剪的栅格数据和多边形数据：

# 打开栅格数据
raster_ds = gdal.Open("input_raster.tif")
# 打开多边形数据
polygon_ds = ogr.Open("input_polygon.shp")

3. 获取栅格数据的空间参考系统信息，并根据该信息创建输出栅格数据集：

# 获取栅格数据的空间参考系统
raster_srs = osr.SpatialReference()
raster_srs.ImportFromWkt(raster_ds.GetProjection())

# 创建输出栅格数据集
clipped_ds = gdal.GetDriverByName("GTiff").Create("output_raster.tif", x_size, y_size, bands, data_type)
clipped_ds.SetProjection(raster_srs.ExportToWkt())
clipped_ds.SetGeoTransform(geotransform)

4. 获取多边形图层和要素，以及多边形图层的空间参考系统信息：

# 获取多边形图层和要素
polygon_layer = polygon_ds.GetLayer()
polygon_feature = polygon_layer.GetNextFeature()

# 获取多边形图层的空间参考系统
polygon_srs = polygon_layer.GetSpatialRef()

5. 创建裁剪器，并将多边形的几何形状设置为裁剪器的裁剪区域：

# 创建裁剪器
clipper = gdal.Clipper()

# 设置裁剪区域为多边形的几何形状
clipper.SetClipGeometry(polygon_feature.GetGeometryRef())

6. 执行裁剪操作，并将结果写入输出栅格数据集：

# 执行裁剪操作
clipper.Clip(raster_ds, clipped_ds)

# 关闭数据集
raster_ds = None
clipped_ds = None

这只是一个简单的示例，实际的使用可能会根据你的数据和需求有所不同。你可以根据情况进行适当的调整和修改。希望这能帮助到你！