本示例要解决什么问题？

有河北省东北部矢量数据，包含53个图层，例如城市名、立交桥、功率、水系、县乡道等，需要将北戴河行政区划内所有的相关图层数据都提取出来，生成53个新的图层，文件名称保持一致，只包含北戴河数据。

ArcGIS 可以通过 ArcToolbox - Analysis Tools - Extract - Clip 实现对矢量数据的裁剪功能。但是一次只能处理生成一个文件。

但该功能显然不适用于本示例，每次只能处理一个，要累死人。

使用 ArcPy (Python in ArcGIS) 编写脚本，即可批处理运行。

ArcPy入口

1. 【开始菜单】-【ArcGIS】-【Python (command line)】
2. C:\Python27\ArcGIS10.2\python.exe （与实际安装路径相关）

ArcPy官方

1. ArcPy快速浏览

单独编辑县界图层，只保留北戴河，保存为 北戴河.shp

具体方法就不说了。

ArcPy 进行Clip处理，主要代码如下

arcpy.Clip_analysis("clipfrom.shp","base.shp","C:/clipto.shp")

本次使用python遍历工程文件夹，将所有shp进行clip处理，所需代码如下

clippatch.py.txt

#conding=utf8  
 
import sys
reload(sys)
sys.setdefaultencoding('utf8')
 
import os 
 
import arcpy
 
from arcpy import env
 
path_in="C:/home/00-library/download/SHP"
path_out="C:/home/00-library/download/SHP_for"
path_clipfile="北戴河.shp"
 
env.workspace = path_in
g = os.walk(path_in) 
num=0
for path,dir_list,file_list in g:  
    for file_name in file_list:  
        if os.path.splitext(file_name)[-1]=='.shp' and file_name!=path_clipfile:
            num=num+1
            shpfileto=path_out+“/”+file_name
            print 'Clipping', file_name
            try:   
                #print 'test', shpfileto
                arcpy.Clip_analysis(file_name.decode("gb2312"),path_clipfile.decode("gb2312"),shpfileto.decode("gb2312"))
            except:
                print '-----error-----', num, '----------------------'
                #print '********* Sth is happened, but still go on.'
            print '-----over-----', num, '----------------------'

打开 ArcGIS Python (command line)，将代码粘贴进去，回车即可运行

运行成功后，处理的成果如下：

1. 使用“裁剪（clip）”工具出现了“000117 警告: 生成的输出为空”拓扑无效 [空多边形不完整] 原因和解决办法：

（1）拓扑无效 [空多边形不完整]错误:

• 原因：两个图册的空间坐标系不同
• 解决：“Data Management Tools”－“投影与变换”－“要素”－“投影”工具，统一图层的坐标系

（2）000117错误：

• 原因：失量图层有问题
• 解决：“Data Management Tools”－“要素”－“修复几何”工具先修复，再裁剪

2. Python 处理文件中文路径问题

1. python 编码问题：’ascii’ codec can’t encode characters in position 的解决方案