裁切数据库
以前使用了GP工具来遍历数据库中的要素层,实现裁切。
现在尝试用Python,实现起来也比较方便,速度比写一个小的应用程序快得多
# --------------------------------------------------------------------------- # 功能: 裁切数据库 # 作者:gisweis # 时间:2015年1月27日 # --------------------------------------------------------------------------- #重载sys,设置编码 import sys reload(sys) sys.setdefaultencoding( "utf-8" ) #添加引用 import arcpy import string try: db = arcpy.GetParameterAsText(0) #裁剪数据库 clip_feat = arcpy.GetParameterAsText(1) #裁剪要素类 outworkspace = arcpy.GetParameterAsText(2) #输出数据库 #获取要素列表 arcpy.env.workspace=db fcs = arcpy.ListFeatureClasses() for fc in fcs: arcpy.Clip_analysis(fc,clip_feat,outworkspace+"\\"+str(fc)) except arcpy.ExecuteError: print arcpy.GetMessages()
把它导入ToolBox中(导入方法:ArcGIS使用Python脚本工具),效果如下:
作者:我也是个傻瓜
出处:http://www.cnblogs.com/liweis/
签名:成熟是一种明亮而不刺眼的光辉。