python：矢量数据裁剪遥感影像

ArcGIS里面的按掩膜提取功能，在python里面的代码如下：

裁剪前：

裁剪后：

from osgeo import gdal

# **********************************根据shp裁剪tif********************************** #

dst = r"E:\\RemoteSensing\\tif_clip.tif"  # 裁剪图像保存完整路径（包括文件名）
src = r'E:\\RemoteSensing\\tif.tif'  # 待裁剪的影像完整路径（包括文件名）
shp = r"E:\\RemoteSensing\\roi_84_UTM.shp"  # 矢量文件的完整路径

ds = gdal.Warp(dst,  # 裁剪图像保存完整路径（包括文件名）
               src,  # 待裁剪的影像
               # warpMemoryLimit=500 内存大小M
               format = 'GTiff',  # 保存图像的格式
               cutlineDSName = shp,  # 矢量文件的完整路径
               cropToCutline = True,
               copyMetadata = True,
               creationOptions = ['COMPRESS=LZW', "TILED=True"],
               dstNodata = -9999)

print("Clip Finish!")