arcpy 做批量裁切 (用面裁切点)

首先需要准备好两类数据 

一、被裁切的点数据

二、面数据

注意:因arcpy对中文不敏感,最好把两类数据分别保存到路径为英文的文件夹下,且被裁切的点数据和进行裁切的面数据最好不要放在同一个文件夹中,不然遍历时可能会出错,因为都是shp文件。

import arcpy
import glob
import sys
reload(sys)
sys.setdefaultencoding('utf-8')          #这里进行编码的强制转化,事实证明用处不大


path = r'E:\shp'                        #工作路径
arcpy.env.workspace = path 
lis = glob.glob(path+"\*.shp")          #glob 文件通配符 *.shp 查找该目录下所有shp文件
print(len(lis))                         #查看列表长度 看看是否正常索引到该目录
print(lis[1][7:-4])                     #字符串切片,以便给输出的文件加名字
in_features =r"\point\stable_Project.shp"  #输入文件,即被裁切的点文件

for i in range(0,len(lis)):
    clip_features = lis[i]      #输入文件,裁切的面文件,这里把列表中的字符串遍历进去,作为路径
    prefix=lis[i][7:-4]         #字符串切片好的名字做前缀
    suffix='.shp'                #文件类型做后缀
    out_feature_class =prefix+'res'+suffix #输出文件,即裁切好的点文件,名字是原文件名+res+.后缀
    arcpy.Clip_analysis(in_features, clip_features, out_feature_class) #用clip工具

 

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值