简单例子:
import arcpy
shuju = r"E:/HJShape/h49B001001"
arcpy.env.workspace = shuju
fcs = arcpy.ListFeatureClasses()
lspt = []
for fc in fcs:
lspt.append(fc)
arcpy.Merge_management(lspt,"ceshi.shp")/
这里引入arcgis自带的arcpy工具包,这个只有在arcgis自带的python环境中才有。arcpy.env.workspace表示为运行空间,即数据目录。arcpy.Merge_management(A,“B”),其中A表示数据目录下的图层要素组成的数组,B为新生成的shp名称。 注意:这里只支持同类型的矢量数据合并,即只能是点数据、或者线数据或者面数据,不能两两混合。 如有必要,需要识别数据类型,代码如下:
arcpy.Describe(fc).shapeType
打印上面的值,结果为Polygon,Polyline,Point,MultiPoint,MultiPatch值。通过if判断语句相应的进行合并。 识别数据类型后,我们需要获取同类型名称前缀的数据将其合并,如HYDNT_arc_21011、HYDNT_arc_21012等,统一合并成HYDNT,其他数据依此类推。因此,我们需要遍历文件夹。代码如下:
def readFilename(path, allfile): # 读取原始文件夹下文件目录
filelist = os.listdir(path) # 获取目录下文件列表
for filename in filelist: # 循环取出列表中的元素
filepath = os.path.join(path, filename) # 将元素名与路径名合并
if os.path.isdir(filepath):
readFilename(filepath, allfile) # 子目录下递归
else:
allfile.ap

本文介绍了如何使用arcgis的arcpy工具包进行矢量图层的批量合并。首先设置工作空间,然后通过识别图层数据类型(点、线、面),将相同类型的数据复制到新文件夹,最后对每个文件夹的内容进行合并。代码示例展示了遍历文件夹、复制和合并图层的过程,实现高效的数据整理。
最低0.47元/天 解锁文章
406





