调用arcpy包批量进行矢量掩膜提取

时间:2019-09-03 01:27:32   收藏:0   阅读:333

使用一个polygon矢量提取某个文件夹中所有的tif格式栅格数据

(要确保先安装好arcpy包)

import arcpy
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput=1
arcpy.env.workspace = "E:\\data"  #所有栅格影像所在文件夹
rasters = arcpy.ListRasters("*", "tif")
mask= "E:\\data\\bound.shp"  #用于提取的矢量掩膜
for raster in rasters:
    print(raster)
    out= "E:\\data\\images\\"+raster
    arcpy.gp.ExtractByMask_sa(raster, mask, out)
print("OK")

2019-09-02 23:53:36

原文:https://www.cnblogs.com/icydengyw/p/11450600.html

评论(0
© 2014 bubuko.com 版权所有 - 联系我们:wmxa8@hotmail.com
打开技术之扣,分享程序人生!