石家庄GIS数据怎么转GeoJSON？Shapely与Fiona实战技巧（附：代码示例）

在石家庄的城市规划、环境监测或物流配送项目中，你是否遇到过这样的困境：手中的GIS数据是老旧的 Shapefile 或复杂的 CAD 格式，而现代的 Web 前端框架（如 Leaflet、Mapbox）却只认 GeoJSON？格式转换不仅耗时，还容易丢失坐标系或属性信息。本文将深入浅出地介绍如何利用 Python 的两大神器 Shapely 和 Fiona，将石家庄的 GIS 数据高效、精准地转换为 GeoJSON，并附带实战代码示例，助你打通数据流转的最后 1 公里。

Shapely 与 Fiona：核心工具解析

在开始转换之前，我们需要理解这两个库在 GIS 数据处理中的分工。简单来说，Shapey 负责处理几何形状的逻辑运算，而 Fiona 负责读写地理空间文件。

Shapely 是基于 GEOS（Geometry Engine, Open Source）的 Python 封装。它不关心数据存储在哪里，只关心内存中的几何对象。你可以用它来计算缓冲区、判断相交、求并集等。在转换过程中，我们常利用它来修正几何体的拓扑错误。

Fiona 则是 GDAL/OGR 的 Pythonic 包装器。它是读写矢量数据（Shapefile, GeoJSON 等）的利器。Fiona 能够读取文件的元数据（如坐标参考系统 CRS）、遍历特征（Features），并将几何对象序列化为 GeoJSON 格式。

实战步骤：从 Shapefile 到 GeoJSON

假设你有一份石家庄的行政区划 Shapefile（例如 shijiazhuang_districts.shp），以下是如何将其转换为标准 GeoJSON 的具体步骤。

1. 环境准备与库安装：确保你的 Python 环境中安装了 shapely 和 fiona。建议使用 Anaconda 或 pip 安装：

   pip install shapely fiona

2. 读取源数据：使用 Fiona 打开 Shapefile。注意，Fiona 使用上下文管理器（with 语句）来安全地处理文件 I/O。你需要指定读取模式为 'r'。
3. 几何处理（可选但推荐）：读取到的几何对象是 Shapely 对象。如果你的石家庄数据存在自相交或多边形不闭合的问题，可以使用 shapely.geometry.shape() 进行修复，例如使用 buffer(0) 来修复微小的拓扑错误。
4. 写入 GeoJSON：创建一个新的 GeoJSON 文件。Fiona 需要一个 schema（模式）来定义输出文件的属性和几何类型。最后，遍历源数据，将处理后的几何对象和属性写入新文件。

代码示例：完整的转换脚本

以下是一个完整的 Python 脚本，演示了如何读取石家庄的 Shapefile 并将其转换为 GeoJSON。

import fiona
from shapely.geometry import shape, mapping

# 定义输入和输出文件路径
input_shp = 'shijiazhuang_road_network.shp'
output_geojson = 'shijiazhuang_road_network.geojson'

# 定义 GeoJSON 的 schema（元数据）
# 注意：GeoJSON 标准通常使用 WGS84 (EPSG:4326) 坐标系
schema = {
    'geometry': 'LineString',  # 根据实际数据类型调整，如 Polygon, Point
    'properties': {'name': 'str', 'length': 'float'}
}

with fiona.open(input_shp, 'r') as source:
    # 获取源数据的坐标系
    crs = source.crs
    
    # 打开目标文件进行写入
    with fiona.open(output_geojson, 'w', driver='GeoJSON', crs=crs, schema=schema) as sink:
        for feature in source:
            # 获取几何对象
            geom = shape(feature['geometry'])
            
            # 【高级技巧】修复微小的拓扑错误（如自相交）
            if not geom.is_valid:
                geom = geom.buffer(0)
            
            # 更新特征的几何信息
            feature['geometry'] = mapping(geom)
            
            # 写入新文件
            sink.write(feature)

print("转换完成！")

扩展技巧：处理坐标系与属性过滤

在实际的石家庄 GIS 项目中，直接转换往往不够，还需要处理以下两个关键问题。

坐标参考系统 (CRS) 的转换

许多国内的 GIS 数据（如西安80或北京54坐标系）使用的是投影坐标（单位为米），而 GeoJSON 标准（RFC 7946）强制要求使用 WGS84 (EPSG:4326) 经纬度坐标。如果在 Web 地图上直接加载非 WGS84 的数据，位置会严重偏移。

使用 Fiona 和 pyproj 可以在读写过程中实时转换坐标。建议在写入 GeoJSON 前，显式检查并转换 CRS。

属性数据的清洗与筛选

Shapefile 中的属性表往往包含无用的字段或复杂的编码。在转换为 GeoJSON（体积较小的文本格式）时，建议只保留关键字段。

你可以通过修改循环中的 `feature['properties']` 字典来实现。例如，只保留 "NAME" 和 "POPULATION" 字段，并将中文字段名转换为英文，以兼容不同的前端解析库。

FAQ 常见问题解答

1. 为什么我的 GeoJSON 在浏览器中显示为空白？

最可能的原因是坐标系不匹配。请检查源数据的 CRS。如果源数据是投影坐标（如 EPSG:3857），而你未将其转换为 EPSG:4326，地图容器将无法正确解析经纬度。使用 Fiona 打开时，可以通过 `source.crs` 查看，并在写入时使用 `pyproj` 转换器处理。

2. Shapely 和 Fiona 与 GDAL/OGR 有什么关系？

Shapely 是纯粹的几何计算库，不依赖 GDAL；而 Fiona 是 GDAL/OGR 的 Python 接口，专门用于读写矢量数据。两者结合使用时，Fiona 负责读取文件流，Shapely 负责处理读取到的几何对象，这是 Python GIS 处理中最经典的组合。

3. 转换后的 GeoJSON 文件过大怎么办？

GeoJSON 是纯文本格式，文件体积通常比二进制的 Shapefile 大。如果数据量很大（如石家庄全市的路网），建议：

• 简化几何：使用 Shapely 的 simplify(tolerance) 方法减少节点数。
• 分层导出：按行政区划或要素类型拆分成多个小文件。
• 使用矢量切片：将 GeoJSON 转换为 MBTiles 或 PMTiles 格式，这是 Web 地图的大数据最佳实践。

总结

通过结合 Shapely 的几何处理能力和 Fiona 的文件 I/O 能力，石家庄的 GIS 数据转换不再是一项繁琐的重复劳动。掌握坐标系校正和属性清洗的技巧，能让你的 GeoJSON 数据在 Web 应用中表现得更加稳定和高效。希望这篇实战教程能为你的 GIS 项目提供实质性的帮助，立即动手试试吧！