GIS 数据转换：使用 GDAL 将 Shp 转换为 GeoJSON 数据

前言

❝

GeoJSON 作为一种通用的地理数据格式，可以很方便地用于共享交换。在 GIS 开发中，经常需要进行数据的转换处理，其中常见的便是将 Shp 转换为 GeoJSON 数据进行展示。

有关GeoJSON数据的详细介绍，请参考往期文章：GeoJSON 数据简介

在之前的文章中讲了如何使用GDAL或者ogr2ogr工具将txt以及csv文本数据转换为Shp格式，本篇教程在之前一系列文章的基础上讲解如何使用GDAL将Shp转换为GeoJSON数据。

• GDAL 简介
• GDAL 下载安装
• GDAL 开发起步

如果你还没有看过，建议从以上内容开始。

1. 开发环境

本文使用如下开发环境，以供参考。

2. 数据准备

如下是本文选取的世界边界范围的Shp数据结构：

3. 导入依赖

Shp作为一种矢量数据格式，可以使用矢量库OGR进行处理，以实现Shp数据转换为GeoJSON格式。还需要引入osr模块用于坐标定义以及json模块用于几何数据转换。

from osgeo import ogr,osr
import os
import json

4. 数据读取与转换

定义一个方法Shp2GeoJSON(shpPath,jsonPath)用于将Shp数据转换为GeoJSON数据。

"""
说明：将 GeoJSON 文件转换为 Shapfile 文件
参数：
-shpPath：Shp 文件路径
-jsonPath：GeoJSON 文件路径
"""
def Shp2GeoJSON(shpPath,jsonPath):

在进行GeoJSON数据格式转换之前，需要检查Shp数据路径是否存在。

# 检查文件是否存在
ifos.path.exists(shpPath):
print("shp 文件存在。")
else:
print("shp 文件不存在，请检查数据路径！")
return

打开Shp数据源。

# 读取Shp文件
shpDataSource = ogr.Open(shpPath)
shpLayer = shpDataSource.GetLayer()

首先构造一个GeoJSON空数据结构，用于填充Shp属性数据。

# 构造GeoJSON对象
geoJSON = {
"type":"FeatureCollection",
"features":[]
}

读取并遍历Shp图层，将几何对象以及属性对象复制到GeoJSON对象中。使用ExportToJson方法将几何数据转换为JSON格式，然后使用json.loads方法进行加载，属性读取完成之后，将要素添加到要素集合中。

# 遍历所有要素
forfeatureinshpLayer:
# 几何对象
geom = feature.GetGeometryRef()
# 构造GeoJSON Feature对象
featureJSON = {
"type":"Feature",
"geometry":json.loads(geom.ExportToJson()),
"property":{}
}

# 获取属性
fieldCount = feature.GetFieldCount()
foriinrange(fieldCount):
fieldName = feature.GetFieldDefnRef(i).GetName()
fieldValue = feature.GetField(i)
featureJSON["property"][fieldName] = fieldValue

# 添加要素
geoJSON["features"].append(featureJSON)

在Shp数据读取完成之后，将其保存到GeoJSON文件中，并关闭数据源。

# 写入文件
with open(jsonPath,"w",encoding="UTF-8") as f:
json.dump(geoJSON,f,ensure_ascii=False, indent=2)

# 关闭数据源
shpDataSource = None