QGIS 矢量图层批量导出geojson格式的自动化脚本实现

1. 为什么需要批量导出geojson？

在日常GIS数据处理工作中，我们经常遇到需要将多个矢量图层统一导出为geojson格式的情况。比如我最近接手的一个城市基础设施项目，需要处理超过200个不同类别的矢量图层，如果一个个手动导出，不仅耗时费力，还容易出错。这时候，一个自动化脚本就能帮上大忙。

geojson作为一种轻量级的地理数据交换格式，在web地图开发、移动端应用和数据共享中有着广泛的应用。它基于JSON标准，人类可读性强，同时又能完整保留地理要素的几何信息和属性数据。实测下来，相比shapefile等传统格式，geojson在跨平台使用时更加稳定可靠。

在QGIS中，虽然可以通过右键菜单手动导出单个图层，但当面对几十上百个图层时，这种方法就显得力不从心了。我踩过的坑告诉我，手动操作不仅效率低下，还容易遗漏某些图层或者搞混导出路径。这就是为什么我们需要借助Python脚本来实现批量导出的自动化流程。

2. 环境准备与基础配置

2.1 确保QGIS Python环境正常

在开始编写脚本前，我们需要确认QGIS的Python环境已经配置妥当。打开QGIS后，点击顶部菜单栏的"插件"→"Python控制台"，这会打开一个交互式的Python命令行界面。在这里输入简单的Python命令如print("Hello QGIS")，如果能够正常输出，说明环境没有问题。

我建议在正式编写脚本前，先熟悉几个常用的QGIS Python模块：

• qgis.core：包含QGIS核心功能的接口
• qgis.gui：提供图形界面相关的类
• qgis.utils：包含实用工具函数

2.2 创建脚本编辑器

虽然可以在Python控制台直接输入代码，但对于较长的脚本，使用脚本编辑器会更方便。点击Python控制台工具栏上的"显示编辑器"按钮（图标看起来像个文本文件），这会打开一个功能更完善的代码编辑器窗口。

这里有个小技巧：我习惯把常用的脚本保存为.py文件，存放在固定的目录下。这样下次使用时可以直接打开，不用重新编写。编辑器支持语法高亮和基本的代码补全，对于Python开发来说已经足够用了。

3. 核心脚本编写详解

3.1 理解脚本的基本结构

让我们仔细分析这个批量导出脚本的每个部分。首先需要导入必要的模块：

from qgis.core import QgsProject, QgsVectorFileWriter, QgsCoordinateReferenceSystem

这三个类分别负责：

• QgsProject：管理当前QGIS项目中的所有图层
• QgsVectorFileWriter：处理矢量数据的写入操作
• QgsCoordinateReferenceSystem：定义坐标参考系统

接下来定义主函数：

def export_layers_to_geojson(output_folder): # 获取当前项目中的所有图层 layers = QgsProject.instance().mapLayers().values() # 遍历所有图层 for layer in layers: # 检查是否为矢量图层 if layer.type() == QgsMapLayer.VectorLayer: # 构造输出路径 output_file = f"{output_folder}/{layer.name()}.geojson" # 设置导出选项 options = QgsVectorFileWriter.SaveVectorOptions() options.driverName = "GeoJSON" options.fileEncoding = "UTF-8" # 执行导出 QgsVectorFileWriter.writeAsVectorFormat( layer, output_file, options ) print(f"Exported: {layer.name()} to {output_file}")

3.2 关键参数解析与优化

在实际使用中，我发现有几个参数值得特别注意：

1. 输出路径处理：Windows和Linux系统的路径分隔符不同，建议使用os.path.join来确保跨平台兼容性。比如：

import os output_file = os.path.join(output_folder, f"{layer.name()}.geojson")

2. 编码设置：虽然UTF-8是通用编码，但遇到特殊字符时可能会出现问题。可以添加错误处理：

try: QgsVectorFileWriter.writeAsVectorFormat(layer, output_file, options) except Exception as e: print(f"Failed to export {layer.name()}: {str(e)}")

3. 坐标系转换：如果需要统一输出坐标系，可以在options中添加：

options.ct = QgsCoordinateTransform( layer.crs(), QgsCoordinateReferenceSystem("EPSG:4326"), # WGS84 QgsProject.instance() )

4. 高级功能扩展

4.1 按条件筛选图层

有时候我们只需要导出特定类型的图层。比如只导出多边形图层，可以修改判断条件：

if (layer.type() == QgsMapLayer.VectorLayer and layer.geometryType() == QgsWkbTypes.PolygonGeometry): # 导出操作

或者根据图层名称过滤：

if "buildings" in layer.name().lower(): # 导出操作

4.2 批量设置属性选项

geojson支持多种导出选项，我们可以通过修改SaveVectorOptions来定制输出：

options.includeZ = True # 包含Z值 options.attributes = ["name", "type"] # 只导出指定属性 options.symbologyExport = QgsVectorFileWriter.SymbolLayerSymbology # 导出符号系统

4.3 进度反馈与日志记录

对于大量图层的导出，添加进度反馈很有必要：

total = len([l for l in layers if l.type() == QgsMapLayer.VectorLayer]) current = 0 for layer in layers: if layer.type() == QgsMapLayer.VectorLayer: current += 1 print(f"Processing {current}/{total}: {layer.name()}") # 导出操作

更专业的做法是写入日志文件：

import logging logging.basicConfig(filename='export.log', level=logging.INFO) # 在导出成功后添加 logging.info(f"Successfully exported {layer.name()}")

5. 实际应用案例分享

去年我在处理一个智慧城市项目时，需要将全市的POI数据（约150个图层）导出为geojson供前端使用。手动导出不仅耗时，还经常因为内存不足导致QGIS崩溃。改用这个脚本后，整个导出过程只需要3分钟，而且完全自动化。

另一个案例是为移动应用准备离线地图数据。我们每周都需要更新一次数据，使用这个脚本后，只需简单修改输出路径就能完成全部工作，节省了大量人力成本。

有个小技巧：我会在脚本开头添加一个时间戳，这样导出的文件会自动包含日期信息，方便版本管理：

from datetime import datetime timestamp = datetime.now().strftime("%Y%m%d") output_folder = f"D:/exports/{timestamp}"

6. 常见问题排查

6.1 导出失败的可能原因

根据我的经验，导出失败通常有以下几个原因：

1. 路径权限问题：确保输出文件夹有写入权限
2. 图层锁定：某些图层可能被其他程序占用
3. 内存不足：处理大型图层时可能需要增加QGIS内存限制
4. 特殊字符：图层名称或属性中包含特殊字符可能导致问题

6.2 性能优化建议

当处理大量数据时，可以尝试以下优化：

1. 关闭不必要的QGIS插件
2. 增加Python的递归限制：sys.setrecursionlimit(10000)
3. 分批处理：将图层分成若干组分别导出
4. 使用多线程（高级技巧，需要谨慎处理QGIS的线程安全）

6.3 调试技巧

遇到问题时，可以逐步缩小范围：

1. 先尝试导出单个图层
2. 检查中间变量值
3. 添加详细的错误日志
4. 在QGIS Python控制台交互测试代码片段

7. 与其他工具的集成

这个脚本可以很容易地集成到自动化工作流中。比如结合Windows任务计划或Linux的cron，可以实现定时自动导出。我通常会在导出完成后自动调用压缩工具打包文件：

import zipfile with zipfile.ZipFile('output.zip', 'w') as zipf: for file in os.listdir(output_folder): if file.endswith('.geojson'): zipf.write(os.path.join(output_folder, file))

对于更复杂的场景，可以考虑将脚本封装为QGIS插件，或者通过PyQGIS的API与其他Python程序交互。比如导出完成后自动上传到Web服务器：

import requests for file in os.listdir(output_folder): if file.endswith('.geojson'): with open(os.path.join(output_folder, file), 'rb') as f: requests.post('http://example.com/upload', files={'file': f})

8. 最佳实践建议

经过多次项目实践，我总结出以下几点经验：

1. 保持代码可读性：添加适当的注释和文档字符串
2. 参数化配置：将输出路径等可变参数提取为脚本开头变量
3. 错误处理：预料可能的问题并妥善处理
4. 版本控制：将脚本纳入Git等版本管理系统
5. 代码复用：将通用功能封装为函数，方便其他脚本调用

一个更健壮的脚本模板如下：

#!/usr/bin/env python """ QGIS矢量图层批量导出geojson脚本 作者：你的名字 日期：2023-08-20 """ import os from datetime import datetime from qgis.core import QgsProject, QgsVectorFileWriter, QgsMapLayer def export_layers(output_folder, filter_func=None): """导出所有符合条件的矢量图层到指定文件夹 Args: output_folder: 输出目录路径 filter_func: 可选，图层过滤函数 Returns: tuple: (成功数, 失败数) """ if not os.path.exists(output_folder): os.makedirs(output_folder) layers = QgsProject.instance().mapLayers().values() success = 0 fail = 0 for layer in layers: if layer.type() != QgsMapLayer.VectorLayer: continue if filter_func and not filter_func(layer): continue try: output_file = os.path.join(output_folder, f"{layer.name()}.geojson") options = QgsVectorFileWriter.SaveVectorOptions() options.driverName = "GeoJSON" QgsVectorFileWriter.writeAsVectorFormat(layer, output_file, options) success += 1 except Exception as e: print(f"Error exporting {layer.name()}: {str(e)}") fail += 1 return (success, fail) if __name__ == "__main__": # 使用示例 output_dir = os.path.expanduser("~/geojson_exports") timestamp = datetime.now().strftime("%Y%m%d_%H%M") final_output = os.path.join(output_dir, timestamp) print(f"开始导出到 {final_output}") success, fail = export_layers(final_output) print(f"导出完成: {success}个成功, {fail}个失败")