前言
❝GeoJSON 作为一种通用的地理数据格式，可以很方便地用于共享交换。在 GIS 开发中，经常需要进行数据的转换处理，其中常见的便是将 Shp 转换为 GeoJSON 数据进行展示。
有关GeoJSON数据的详细介绍，请参考往期文章：GeoJSON 数据简介
在之前的文章中讲了如何使用GDAL或者ogr2ogr工具将txt以及csv文本数据转换为Shp格式，本篇教程在之前一系列文章的基础上讲解如何使用GDAL将Shp转换为GeoJSON数据。
GDAL 简介 GDAL 下载安装 GDAL 开发起步
如果你还没有看过，建议从以上内容开始。
1. 开发环境
本文使用如下开发环境，以供参考。
时间：2025年
系统：Windows 11
Python：3.11.7
GDAL：3.11.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
