裁切数据库

以前使用了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脚本工具），效果如下：