文章目录

前言

一、遇到的问题？

1.中文编码问题

2.arcpy函数问题

二、实现步骤

1.代码思路

2.制作脚本工具

前言

使用的arcpy为arcgis10.2对应版本。（python27）

想遍历面要素的多条记录实现批量裁剪功能，所幸找到用shp批量裁剪栅格的代码。在此基础上进行改进，记录一下实现过程和问题的解决。

https://www.likecs.com/show-203727205.html

一、遇到的问题？

1.中文编码问题

遇到如下报错

UnicodeEncodeError: 'ascii' codec can't encode characters in position

解决方法：设置编码

import sys
reload(sys)
sys.setdefaultencoding('utf-8')

2.arcpy函数问题

起初使用Clip_management裁剪，报错参数无效。于是找了很久的对应函数，最终在官网文档找到完美解答。ps:对应的函数位置和arcgis中工具箱位置是一致的。(arcpy.Clip_analysis同样可以实现shp裁剪shp)

裁剪函数如下：

arcpy.analysis.Clip(inputshp,row[0],resultPath+"\\"+tempShpName)

二、实现步骤

1.代码思路

与参考案例类似，首先通过游标遍历面数据子要素，在每次遍历过程中实现clip，并添加目标字段，这里添加的是县域名称，最后实现合并。

字段添加并赋值代码如下：

arcpy.AddField_management(Clip, "qu", "TEXT", "", "", 20)
result=arcpy.CalculateField_management(Clip, "qu", '\"'+ str(row[1]) + '\"',"VB", "")

合并代码：

合并的函数一定要找对应版本的官方文档，arcgis pro使用的是arcpy.ia.Merge，在这里也踩了坑。

arcpy.Merge_management(inFeatures, resultPath+"\\"+'water_merge'+'.shp')

2.制作脚本工具

自定义脚本工具实现流程不再赘述，csdn也有很多很好的教程。推荐GIS荟

全部代码如下：

import arcpy
import sys
reload(sys)
sys.setdefaultencoding('utf-8')inputshp = arcpy.GetParameterAsText(0)
clipshp = arcpy.GetParameterAsText(1)
resultFolder = arcpy.GetParameterAsText(2)
inFeatures = []# 实现批量clip
def clip(inputshp, clipshp, resultPath):#遍历面要素，SHAPE@为几何对象with arcpy.da.SearchCursor(clipshp, ["SHAPE@",'NAME']) as cursor:for row in cursor:tempShpName = str(row[1]) + '.shp'Clip=arcpy.analysis.Clip(inputshp,row[0],resultPath+"\\"+tempShpName)#添加字段arcpy.AddField_management(Clip, "qu", "TEXT", "", "", 20)#字段赋值result=arcpy.CalculateField_management(Clip, "qu", '\"'+ str(row[1]) + '\"',"VB", "")arcpy.AddMessage(row[1]+"success")inFeatures.append(result)#合并arcpy.Merge_management(inFeatures, resultPath+"\\"+'water_merge'+'.shp')
if __name__ == '__main__':clip(inputshp, clipshp, resultFolder)

实现效果如下：

shp批量裁剪工具-基于Arcpy相关推荐

批量处理|基于ENVI的遥感影像批量裁剪工具
遥感影像裁剪的目的是将研究之外的区域去除,常用的方法是按照行政区划边界或者自然区划边界进行影像裁剪,在基础数据生产中,还经常要进行标准分幅裁剪. ENVI软件支持规则裁剪(矩形范围)和不规则裁剪(任意 ...
Python地理数据处理十五：基于arcpy的批量操作
目录 1. 批量格式转换(grd to tiff) 2. 批量定义投影(Albers) 3. 批量投影转换(WGS 1984 to WGS 1984 Albers) 3.1 转换前的栅格图像属性: 3 ...
基于arcpy开发arcgis工具
基于arcpy开发arcgis工具 1.添加工具 2.使用方法 2.1批量添加坐标系脚本工具简介工具源码 2.2批量添加字段脚本工具简介工具源码 2.3更新属性表脚本工具简介工具源码 2. ...
android源码下载方法批量下载基于windows os
安装win版的Gitbash, 在这里 http://msysgit.googlecode.com/files/Git-1.6.0.2-preview20080923.exe. 选择默认安装路径(否则 ...
Android 批量打包基于Walle的多渠道快速打包自动脚本
Android 批量打包基于Walle的多渠道快速打包自动脚本 Android Signature V2 Scheme签名下的新一代渠道包打包神器 github主页:https://github.c ...
Python基于Arcpy对多时间尺度NDVI进行最大值合成
Python基于Arcpy对多时间尺度NDVI进行最大值合成 1. 前期数据准备 2. 处理代码前言:在现实研究中很容易拿到日尺度.月尺度的NDVI栅格影像,繁多的栅格数据需要按照研究的目的进行重新 ...
ArcPy 使用Shp批量裁剪栅格 Clip
小记一下 Clip_management 的使用方法 inRaster = "d.tif" #输入栅格路径 outRaster = r"C:/Users/xhz36/Do ...
Arcpy批量导出shp文件属性表——使用arcpy.da.SearchCursor函数
1.问题描述我有10个文件夹,每一个文件夹有18个shp文件(温度数据),需要提取每个shp文件的温度数据,即下图属性表中RASTERVALU列. 2.上代码 #coding:utf-8 impo ...
基于FME实现的地理数据库批量建库的解决方案，支持gdb、mdb、shapefile等数据格式，gdb批量建库，mdb批量建库，shp批量建库，shapefile批量建库，地理数据批量建库
目录一.建库效果 1.图层及字段效果 2.图层坐标系效果二.实现过程 1.制作建库结构表 2.字段映射 3.映射几何类型 4.构建schema 5.输出设置 6.图层坐标系设置 7.FME模板的整 ...

shp批量裁剪工具-基于Arcpy

文章目录

前言

一、遇到的问题？

1.中文编码问题

2.arcpy函数问题

二、实现步骤

1.代码思路

2.制作脚本工具

shp批量裁剪工具-基于Arcpy相关推荐

最新文章

热门文章