GDAL栅格裁剪（Python）

def clip_raster(in_raster, out_raster, mask_shp):
    """
    :param in_raster: 输入栅格
    :param out_raster: 输出栅格
    :param mask_shp: 裁剪矢量
    :param wkid: wkid
    :return:
    """
    gdal.Warp(out_raster,
              in_raster,
              format='GTiff',
              dstSRS='EPSG:4326',
              cutlineDSName=mask_shp,
              cropToCutline=True,  # 按掩膜图层范围裁剪
              dstNodata=-9999,
              outputType=gdal.GDT_Float64)