Python shapefile转GeoJson的2种方式实例

作者:小李AI飞刀^_^ 时间:2023-02-20 01:07:52 

GeoJson的简要介绍

GeoJson是用json的语法表达和存储地理数据,可以说是json的子集。

GeoJson以键值对的形式保存原有对象的信息,具有轻量化、易解析等优点。

GeoJson包括的地理要素有Point(点)、 MultiPoint(多点)、 LineString(线)、MultiLineString(多线)、 Polygon(面)、 MultiPolygon(多面)、 GeometryCollection(几何集合)

这些地理要素包括在geometry的type属性中,并且不同的type具有不同的coordinates值。更多的GeoJson相关内容可参考RFC7946标准。

{
        "type": "MultiPoint",
        "coordinates": [
            [100.0, 0.0],
            [101.0, 1.0]
        ]
    }

{
        "type": "MultiPolygon",
        "coordinates": [
            [
                [
                    [102.0, 2.0],
                    [103.0, 2.0],
                    [103.0, 3.0],
                    [102.0, 3.0],
                    [102.0, 2.0]
                ]
            ],
            [
                [
                    [100.0, 0.0],
                    [101.0, 0.0],
                    [101.0, 1.0],
                    [100.0, 1.0],
                    [100.0, 0.0]
                ],
                [
                    [100.2, 0.2],
                    [100.2, 0.8],
                    [100.8, 0.8],
                    [100.8, 0.2],
                    [100.2, 0.2]
                ]
            ]
        ]
    }

两种将shapefile文件转换为GeoJson的方式

1. 使用geopandas

核心代码:geopandas.GeoSeries 和out_data.to_file

import geopandas as gpd

def shp2geojson_gpd(shp_file, geojson_file):
   """
   将shapefile格式的文件转化为geojson
   :param shp_file: 需要转换的shapefile文件名,投影信息可以缺失,也可以指定
   :param geojson_file: 转换输出的geojson文件名
   """

if os.path.exists(geojson_file):
       os.remove(geojson_file)

out_data = gpd.read_file(shp_file)
   crs = out_data.crs
   out_data = gpd.GeoSeries(out_data.geometry, crs=crs)
   out_data.to_file(geojson_file, driver='GeoJSON', encoding="utf-8")
   print("successfully convert shapefile to geojson")

使用geopandas转换的时候两行核心代码即可搞定,简单粗暴。但是在实践过程中发现,采用geopandas转换后的GeoJson文件并没有保留shapefile中的属性properities信息,如area, name等,如下图所示:

Python shapefile转GeoJson的2种方式实例

2. 使用gdal

import gdal
import ogr
import os

def shp2geojson_gdal(shp_file, geojson_file):
   gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES")
   gdal.SetConfigOption("SHAPE_ENCODING", "GBK")
   src_ds = ogr.Open(shp_file)
   src_layer = src_ds.GetLayer(0)

# 创建结果Geojson
   baseName = os.path.basename(geojson_file)
   dst_driver = ogr.GetDriverByName('GeoJSON')
   dst_ds = dst_driver.CreateDataSource(geojson_file)
   if dst_ds.GetLayer(baseName):
       dst_ds.DeleteLayer(baseName)
   dst_layer = dst_ds.CreateLayer(baseName, src_layer.GetSpatialRef())
   dst_layer.CreateFields(src_layer.schema)
   dst_feat = ogr.Feature(dst_layer.GetLayerDefn())

# 生成结果文件
   for feature in src_layer:
       dst_feat.SetGeometry(feature.geometry())
       for j in range(feature.GetFieldCount()):
           dst_feat.SetField(j, feature.GetField(j))
       dst_layer.CreateFeature(dst_feat)

del dst_ds
   del src_ds
   print("successfully convert shapefile to geojson")

结果包含原始shapefile文件中的属性信息:

Python shapefile转GeoJson的2种方式实例

来源:https://blog.csdn.net/MLH7M/article/details/121024141

标签:python,shapefile,geojson
0
投稿

猜你喜欢

  • 10款最好的Web开发的 Python 框架

    2023-12-21 11:26:04
  • python not关键字实例用法

    2023-12-05 10:25:50
  • Mysql 5.6添加修改用户名和密码的方法

    2024-01-20 09:02:57
  • python解决12306登录验证码的实现

    2023-05-29 10:28:35
  • 前端模板引擎

    2010-07-27 12:33:00
  • go学习笔记读取consul配置文件详解

    2024-05-09 10:08:14
  • Python 函数返回值的示例代码

    2021-10-12 19:30:23
  • Python-split()函数实例用法讲解

    2023-12-12 07:13:20
  • pytorch进行上采样的种类实例

    2023-11-28 14:02:03
  • Web 标准设计实践:Google 的首页

    2008-10-12 12:14:00
  • python实现集中式的病毒扫描功能详解

    2022-03-04 03:16:00
  • javascript获取wx.config内部字段解决微信分享

    2024-04-23 09:15:50
  • SQL Server修改表所有者

    2010-02-04 08:33:00
  • ORM框架之Dapper简介和性能测试

    2024-05-03 15:30:44
  • Python中的集合介绍

    2022-09-10 03:49:52
  • 内容适应形式

    2010-03-18 16:09:00
  • python求众数问题实例

    2022-02-06 22:25:40
  • Python 统计字数的思路详解

    2023-01-29 00:17:44
  • CSS技巧及常见问题列表

    2008-04-06 14:00:00
  • mysql字符集和数据库引擎修改方法分享

    2024-01-21 22:02:49
  • asp之家 网络编程 m.aspxhome.com