GIS实验之文件夹中的批量图层合并(Python实现)

1. 实验任务

数据描述 :提供含有多个数据文件的 20 个文件夹,每中的数据名称都一致。例如每个文件夹中有为“ 名称都一致。例如每个文件夹中有为“ 名称都一致。例如每个文件夹中有为“ city_area” city_area” city_area” city_area” 的文件。

处理要求 :将 20 个文件夹中名称相同的数据分别合并为一。

2. 实验数据分析

在这里插入图片描述

实现思路

  1. 根据data2的路径(path),打开data2文件,遍历出所有的文件名(name)
  2. data2文件与文件名通过字符串连接,形成一个新的路径(path_name),以此路径遍历出所有的shp文件
  3. 通过条件选择,如过shp文件名为"全国铁路.shp"(shpname),根据path_name 与shpname连接成为最终的shp文件的完整路径
  4. 把所有shpname相同的完整路径用一个列表(list)来存储起来
  5. 最后通过arcpy.Merge_management(list,“保存的名字.shp”)即可将列表中的shp文件合并为同一图层的数据

3. Python实现图层合并

# 导入os,arcpy包
import os,arcpy
# data2文件的路径,一定要转为Unicode字符串,否则后面路径连接会出错
path = u'C:\\Users\\17520\\Desktop\\work1\\data2'
# 获取data2文件中的所有文件名
fileList = os.listdir(path)
# 每一个省中有7个图层数据,因此创建7个空列表用于存储数据
a = []
b = []
c = []
d = []
e = []
f = []
g = []
# 遍历data2文件中的文件名
for i in fileList:
	# 生成子文件夹的路径
    path_name =  path + u"\\" + unicode(i)
    # 设置工作空间
    arcpy.env.workspace = path_name
    # 获取所有的图层名称,结果为"名称.shp",结果是Unicode字符串
    fsc = arcpy.ListFeatureClasses()
    # 遍历工作空间中的所有shp文件
    for j in fsc:
    	# 通过判断shp文件名选择用哪一个列表来存储数据
        if j == u"全国铁路.shp":
            a.append(path_name + u"\\" + j)
        elif j == u'县级行政区(面).shp':
            b.append(path_name + u"\\" + j)
        elif j == u'县级行政区(点).shp':
            c.append(path_name + u"\\" + j)
        elif j == u'公路.shp':
            d.append(path_name + u"\\" + j)
        elif j == u'河流及水域四级.shp':
            e.append(path_name + u"\\" + j)
        elif j == u'地级市(点).shp':
            f.append(path_name + u"\\" + j)
        elif j == u'省份.shp':
            g.append(path_name + u"\\" + j)
# 设置图层合并保存的文件地理数据库位置           
arcpy.env.workspace = r'C:\Users\17520\Desktop\work1\model2\data.gdb'
# 设置如果有相同的名称则覆盖
arcpy.env.overwriteoutput = True
# 合并列表中的图层数据
arcpy.Merge_management(a,"全国铁路.shp")
arcpy.Merge_management(b,"县级行政区(面).shp")
arcpy.Merge_management(c,"县级行政区(点).shp")
arcpy.Merge_management(d,"公路.shp")
arcpy.Merge_management(e,"河流及水域四级.shp")
arcpy.Merge_management(f,"地级市(点).shp")
arcpy.Merge_management(g,"省份.shp")

运行过程与结果展示

在这里插入图片描述
在这里插入图片描述

  • 4
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
### 回答1: 代码如下:import arcpymxd = arcpy.mapping.MapDocument("CURRENT")for lyr in arcpy.mapping.ListLayers(mxd):if lyr.supports("DATASOURCE"):arcpy.mapping.ExportToPDF(mxd, lyr.name + ".pdf")del mxd ### 回答2: 为了通过GIS编写一个批量图层导图纸的代码,我们可以使用Python语言,并且利用ArcPy模块的一些工具和函数。下面是一个简单的示例代码: ```python import arcpy from arcpy import mapping # 设置工作空间和输出目录 arcpy.env.workspace = r"C:\data\workspace" output_folder = r"C:\data\output" # 获取地图文档的列表 mxd_list = arcpy.ListFiles("*.mxd") # 循环遍历地图文档列表 for mxd_file in mxd_list: # 加载地图文档 mxd = mapping.MapDocument(mxd_file) # 获取地图文档的数据框 df_list = mapping.ListDataFrames(mxd) # 循环遍历数据框列表 for df in df_list: # 获取数据框图层 layer_list = mapping.ListLayers(mxd, "", df) # 循环遍历图层列表 for layer in layer_list: # 设置导出路径和文件名 output_path = output_folder + "\\" + layer.name + ".pdf" # 导出图层为PDF文件 mapping.ExportToPDF(mxd, output_path, df) # 释放对地图文档的引用 del mxd ``` 在这个代码,我们首先设置工作空间和输出目录。然后,使用`ListFiles()`函数获取所有以`.mxd`为扩展名的地图文档文件列表。 接下来,我们循环遍历地图文档列表。对于每个地图文档,我们使用`mapping.MapDocument()`函数加载它,并使用`ListDataFrames()`函数获取数据框列表。 然后,我们循环遍历数据框列表,并使用`ListLayers()`函数获取每个数据框图层列表。 接下来,我们循环遍历图层列表,并为每个图层设置输出路径和文件名。 最后,我们使用`ExportToPDF()`函数将每个图层导出为PDF文件,并设置地图文档、输出路径和数据框作为参数。 在代码的最后,我们使用`del`语句释放对地图文档的引用。 这是一个简单的示例代码,你可以根据自己的需求来进行修改和扩展。 ### 回答3: 要编写一个批量图层导图纸的GIS代码,首先需要明确输入和输出的数据格式。假设输入数据为一个包含多个图层GIS文件,输出为多个图层的图纸文件。以下是一种可能的代码设计: 1. 首先,我们需要定义输入和输出路径: 输入路径 input_path = "input.gis" 输出路径 output_path = "output" 2. 接下来,需要读取输入的GIS文件gis_data = read_gis(input_path) 可以使用适当的GIS库函数来实现读取GIS文件的功能。根据具体的GIS库和文件格式,读取的方式可能有所不同。 3. 然后,需要对每个图层进行处理: for layer in gis_data.layers: layer_name = layer.name layer_data = layer.data # 你可以在这里添加其他图层处理步骤 # 根据图层名称创建相应的图纸文件 create_drawing(layer_data, output_path, layer_name) 对于每个图层,可以按照需要进行进一步的处理,比如数据清理、地理操作或空间分析等。 4. 创建图纸文件的函数create_drawing的具体实现取决于使用的GIS库和目标文件格式,一般流程如下: - 创建一个新的图纸对象,并设置适当的图层属性(比如图层名称、符号样式)。 - 将图层数据添加到图纸对象。 - 将图纸对象保存为图纸文件。 要实现批量导出多个图纸文件,可以使用不同的图层名称作为文件名称或编号,以保持文件的唯一性。 5. 最后,保存所有图纸文件并结束程序。 这是一个简单的批量图层导图纸的GIS代码示例。具体的实现细节会根据使用的GIS库和文件格式而有所不同,可以根据实际需求进行相应的调整和扩展。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

而又何羡乎

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值