python gdal裁剪 最小外接矩形

Python的GDAL库提供了裁剪栅格影像的功能，并且可以根据需要通过最小外接矩形来裁剪影像。下面是用Python GDAL进行最小外接矩形裁剪的步骤：

首先，导入需要的库，包括GDAL库和OGR库。GDAL库用于处理栅格数据，OGR库用于处理矢量数据。

from osgeo import gdal, ogr

然后，打开需要裁剪的栅格影像，并获取其投影信息和地理转换信息。

# 打开栅格影像
ds = gdal.Open("input.tif")

# 获取投影信息和地理转换信息
proj = ds.GetProjection()
gt = ds.GetGeoTransform()

接下来，定义最小外接矩形的位置和大小。可以根据需要手动指定，也可以通过矢量数据计算得到。

# 手动指定最小外接矩形位置和大小
xmin, ymin, xmax, ymax = 10, 20, 30, 40

然后，创建最小外接矩形的矢量数据。

# 创建最小外接矩形的矢量数据
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(xmin, ymin)
ring.AddPoint(xmin, ymax)
ring.AddPoint(xmax, ymax)
ring.AddPoint(xmax, ymin)
ring.AddPoint(xmin, ymin)

poly = ogr.Geometry(ogr.wkbPolygon)
poly.AddGeometry(ring)

接下来，设置裁剪的参数，并执行裁剪操作。

# 设置裁剪参数
options = ["ALL_TOUCHED=TRUE"]

# 执行裁剪操作
gdal.Warp("output.tif", ds, format="GTiff", cutlineDSName=poly, dstSRS=proj, outputType=gdal.GDT_Float32, options=options)

最后，关闭数据集并输出裁剪后的影像。

# 关闭数据集
ds = None

# 输出裁剪后的影像
output = gdal.Open("output.tif")

以上是使用Python GDAL进行最小外接矩形裁剪的步骤。通过指定最小外接矩形的位置和大小，可以裁剪出指定区域的栅格影像。