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
python gdal使用shp裁剪遥感影像
最新推荐文章于 2023-10-25 20:18:46 发布