调用arcpy包批量进行矢量掩膜提取栅格数据和矢量数据裁剪

使用一个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")
复制代码

批量使用一个矢量数据裁剪另一个矢量数据(点、线、面矢量)

复制代码
#批量Clip矢量文件
import arcpy
arcpy.env.workspace = "E:\\02WORK\\DATA"#所有矢量数据在这个文件夹下面
fcs=arcpy.ListFeatureClasses()
clipmask= "E:\\02WORK\\XiongAnBoundary.shp"
for fc in fcs:
    print(fc)
    out= "Clip\\"+fc[:-5] #Clip是这个工作空间下面的文件夹名称
    arcpy.Clip_analysis(fc, clipmask ,out)
复制代码

 

 

2019-09-02 23:53:36

获取掩膜提取到的栅格的属性:

https://pro.arcgis.com/zh-cn/pro-app/tool-reference/data-management/get-raster-properties.htm

参考:

https://pro.arcgis.com/zh-cn/pro-app/tool-reference/spatial-analyst/extract-by-mask.htm

 
posted @   icydengyw  阅读(2541)  评论(2编辑  收藏  举报
编辑推荐:
· 一个奇形怪状的面试题:Bean中的CHM要不要加volatile?
· [.NET]调用本地 Deepseek 模型
· 一个费力不讨好的项目,让我损失了近一半的绩效!
· .NET Core 托管堆内存泄露/CPU异常的常见思路
· PostgreSQL 和 SQL Server 在统计信息维护中的关键差异
阅读排行:
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库
· 上周热点回顾(2.17-2.23)
点击右上角即可分享
微信分享提示