以前使用了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脚本工具),效果如下: