python实现把一个包括多个面的shp文件，导出为多个只含一个面的shp文件，并按照属性存储在对应文件夹

一、环境配置（在Pycharm中使用arcpy）

        由于需要使用到arcpy模块，而python3与arcpy是不兼容的，所以要用到python2.7。而我们安装的Arcgis10.2内置python2.7。所以如果想在Pycharm中使用arcpy，可以采取的办法是：将pycharm的python解释器替换为arcgis内置的解释器。具体操作如下：

        打开Pycharm，新建一个Project，将Interpreter选择为arcgis内置的2.7版本python.exe，再点击Create:

         至此，基本的环境配置工作已经完成！

二、数据介绍

        我用到的数据名为test_5000.shp，共包含5000个面，按照形状分为10个类型（分别为"E","F","H","I","L","O","T","U","Y","Z"），每种类型包含500个面，其属性表如下：

 三、处理内容

（1）把这一个test_5000.shp文件，导出为5000个单独的shp文件；

（2）每个shp文件，以字段“name”+“fid”命名；

（3）每种类型的500个文件，存放在一个文件夹内，例如500个E类文件，存在文件夹E_shps中；

（4）如果想要将导出的shp文件批量转为栅格，并且按照属性存储在不同的文件夹里，可以看我这篇文章的评论区。

四、程序实现

实现思路：

        （1）先将含有5000个面的test_5000.shp，按照字段“name”导出为10个shp文件,其中每个shp文件包含500个面，并以“name”字段命名；

        （ 2）再将每个shp导出为单独的shp，这样做的好处是，可以实现最终单独shp文件的顺序命名。

# - * - coding: utf-8 -
import arcpy
import os
#工作空间
arcpy.env.workspace = "E:/"
#输入要素
inFc = "E:/data_GISAlgorithm/test_5000r.shp"
#这里写的是我要输出的种类
dirs_name=["E","F","H","I","L","O","T","U","Y","Z"]
############（1）分别把十种类型的要素保存在单独的文件夹
with arcpy.da.SearchCursor(inFc, [ 'name']) as cursor:
    for row in cursor: #遍历字段
        name = str(row[0]) #获取name字段信息
        path="E:/data_GISAlgorithm/Class/folder_{}".format(name)
        if not os.path.exists(path):
            os.makedirs(path)
        outFc = path+'/'+ name  +".shp" #导出的要素路径和名称
        #where_clause = '"NAME" = \'%s\''%(city)
        where_clause = '"name"' + " = " + "'" + name + "'" #导出要素的条件，即根据不同name导出
        if (arcpy.Exists(outFc) == True):  #判断是否已有相同名称的导出要素
            pass
        else:
            arcpy.Select_analysis(inFc, outFc, where_clause) #根据条件导出要素
            print (u"成功导出：" + outFc)

############（2）把导出的每种类型文件分割为单独的shp文件
#创建保存结果的文件夹
for a in dirs_name:
    savepath = 'E:/data_GISAlgorithm/Class_shps/{}_shps'.format(a)
    if not os.path.exists(savepath):
        os.mkdir(savepath)

#第（2）步的正戏部分
for i in dirs_name:
    dir=r"E:/data_GISAlgorithm/Class/folder_{}".format(i)
    files = os.listdir(dir)
    for f in files:
        if os.path.splitext(f)[1] == ".shp":
            input=dir + os.sep + f

            with arcpy.da.SearchCursor(input ,["SHAPE@", 'fid', 'name']) as cursor:
                for row in cursor:
                     # 单个导出shp
                    out_name=str(row[2])+str(row[1])+'.shp'
                    print(out_name)

                    arcpy.FeatureClassToFeatureClass_conversion (row[0],'E:/data_GISAlgorithm/Class_shps/{}_shps'.format(i),out_name)

五、结果展示

第四部分中，第（1）步的结果是：

 

 其中，每一个文件夹内是这样的：

这里面的一个shp文件包含500个面；

第（2）步的结果是：

其中每一个文件夹内是这样的（包含500*4个文件）：

 

 

六、总结

 本文实现了将一个含有5000个面的shp文件按着类型导出为单个的shp文件，并将导出的单个shp文件按着类型存储在对应的文件夹中。