python gdal使用shp裁剪遥感影像

from osgeo import gdal
import time

t0 = time.time()

# dst = r"D:\temp\aus11.tif"
# shp = r"D:\temp\shp\shp.shp"
# src = r"D:\data\aus11.tif"

dst = r"D:\temp\sz3.tif"
shp = r"D:\data\gmy\shp\rect.shp"
src = r"D:\data\gmy\building\sz3\sz3.tif"

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

t1 = time.time()
print("gdal time is %s" % (t1 - t0))
# ————————————————
# 版权声明：本文为CSDN博主「独孤尚亮dugushangliang」的原创文章，遵循CC
# 4.0
# BY - SA版权协议，转载请附上原文出处链接及本声明。
# 原文链接：https: // blog.csdn.net / weixin_40450867 / article / details / 121258033