GDAL遥感影像裁剪

主要使用gdal.Warp这个函数来进行影像裁剪，需要有矢量文件

from osgeo import gdal
def read_tiff(filename):
    # 读取影像，获取影像位置
    dataset = gdal.Open(filename)
    # 获取影像波段数,获取影像长宽
    im_band = dataset.GetRasterBand  # 波段数
    im_width = dataset.RasterXSize  # 宽，列数
    im_height = dataset.RasterYSize  # 高,行数
    # 仿射矩阵
    im_GeoTransform = dataset.GetGeoTransform()
    # 地图投影信息
    img_proj = dataset.GetProjection()

    im_data = dataset.ReadAsArray(0, 0, im_width, im_height)
    return img_proj, im_GeoTransform, im_data

porj, geo,data = read_tiff(workspace2)
        input_raster = gdal.Open(workspace2)
        output_raster = workspace2.replace('.tif', '.cut.tif')
        ds = gdal.Warp(output_raster,
                       input_raster,
                       format='GTiff',
                       cropToCutline=True,  # True,
                       cutlineDSName=workspace1,  # or any other file format
                       # cutlineWhere="FID_1 = 'whatever'",#根据shp属性表字段类型进行裁切
                       # optionally you can filter your cutline (shapefile) based on attribute values
                       dstNodata=-9999)

gdal.Warp是GDAL（Geospatial Data Abstraction Library）中的一个函数，用于执行栅格数据的重投影、裁剪、变换等操作。它可以用来处理各种格式的栅格数据，如地理影像、DEM（数字高程模型）等。

这个函数的基本语法是：

gdal.Warp(destNameOrDestDS, srcDSOrSrcDSTab, options=None, **kwargs)
其中参数含义如下：

destNameOrDestDS: 输出数据集的文件名或已存在的GDAL数据集对象。
srcDSOrSrcDSTab: 输入数据集的文件名、已存在的GDAL数据集对象、或包含多个输入数据集的列表。
options: 一些可选参数，用于控制重投影、裁剪、变换等操作。
**kwargs: 其他关键字参数，用于设置其他选项。
gdal.Warp函数可用于许多常见的栅格数据处理任务，如：

投影变换：将数据集从一个投影转换到另一个投影。
裁剪：根据指定的范围或掩膜文件对数据集进行裁剪。
重采样：根据需要调整栅格数据的分辨率。
格式转换：将一个数据集从一种格式转换为另一种格式。
总之，gdal.Warp函数提供了强大的功能，可用于处理和转换各种栅格数据，是空间数据处理和分析中的重要工具之一。