Python用SHP数据裁剪TIF数据

第一次用python来处理栅格数据，虽然ArcGIS也可以裁剪，但是如果数据比较大的话,ArcGIS就可能会崩掉，所以尝试用python来处理数据。从网上看了一下别人的代码，虽然代码很简单，但是很有用。

from osgeo import gdal  # 如果python中没有安装GDAL库，会报错，需要安装GDAL库
input_shape = r"E:/temp/Export_Output.shp"
output_raster=r'E:/code/temporary_clip/test1.tif'
# tif输入路径，打开文件
input_raster = r"E:/code/Hmp.tif"
# 矢量文件路径，打开矢量文件
input_raster=gdal.Open(input_raster)
# 开始裁剪
ds = gdal.Warp(output_raster,
              input_raster,
              format = 'GTiff',#设置格式
              cutlineDSName = input_shape,
               cropToCutline=True,  # 保证裁剪后影像大小跟矢量文件的范围大小一致（当设置为False，结果大小会跟待裁剪影像大小一样，会出现大量的空值区域）
              dstNodata = 0)