python 利用gdal库 shp文件转geojson

最新推荐文章于 2024-08-08 17:34:50 发布

只喝粥!不弹吉他

最新推荐文章于 2024-08-08 17:34:50 发布

阅读量2.3k

点赞数 2

分类专栏： Python相关问题文章标签： python 开发语言

本文链接：https://blog.csdn.net/weixin_43844000/article/details/123370412

版权

Python相关问题专栏收录该内容

2 篇文章 0 订阅

订阅专栏

本文介绍了如何使用Python的GDAL库将Shp文件转换为Geojson格式，包括关键代码示例。在转换过程中，可能会遇到如GDAL安装、proj.db配置和Shp文件放置等问题，文中提供了解决方案。对于新手，需要注意从官方源下载GDAL whl文件安装，设置环境变量以解决proj.db找不到的问题，并确保所有Shp文件组件在同一目录下。

摘要由CSDN通过智能技术生成

代码转载（我在最下方复制了代码）：https://www.cnblogs.com/MatthewHome/p/11383808.html

gdal在线文档：http://pcjericks.github.io/py-gdalogr-cookbook/#

上面的程序可以跑，但是对于新手有几个注意点（我遇到、解决并可以跑出来的！！）

注意点：

1、安装gdal需要在官网下载，不能直接pip安装

下载网址：https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal

下载如图的对应版本whl文件，再使用pip install GDAL-3.4.1-cp37-cp37m-win_amd64.whl 安装对应python版本的whl文件即可。

2、报错：ERROR 1: PROJ: proj_identify: Cannot find proj.db

原因：未配置环境(由于安装了postgis或者其他gis相关软件引起的)

解决方法：找到你的gdal包对应的文件夹，再找出proj.db文件，然后把proj.db所在目录如图进行配置，变量名为PROJ_LIB，变量值为所在路径，然后关闭重启相关编程软件（我是vscode，关闭重启vscode）即可。

3、shp文件放置问题

shp文件由很多子文件构成，执行python算法时需要把相关文件放在同一目录下，才能识别，要不然就会报getLayer为0的错误

代码：

# -*- coding: utf-8 -*-

from osgeo import ogr
from osgeo import gdal
import sys
import os

def ChangeToJson():
    print("Starting........")
    vector = "E:/P_new11.shp"
    output = 'E:/creatJson.json'
    #打开矢量图层
    gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES")
    gdal.SetConfigOption("SHAPE_ENCODING", "GBK")
    shp_ds = ogr.Open(vector)
    print("sadasd",shp_ds)
    shp_lyr = shp_ds.GetLayer(0)

    # 创建结果Geojson
    baseName = os.path.basename(output)
    out_driver = ogr.GetDriverByName('GeoJSON')
    out_ds = out_driver.CreateDataSource(output)
    if out_ds.GetLayer(baseName):
        out_ds.DeleteLayer(baseName)
    out_lyr = out_ds.CreateLayer(baseName, shp_lyr.GetSpatialRef())
    out_lyr.CreateFields(shp_lyr.schema)
    out_feat = ogr.Feature(out_lyr.GetLayerDefn())

    #生成结果文件
    for feature in shp_lyr:
        out_feat.SetGeometry(feature.geometry())
        for j in range(feature.GetFieldCount()):
            out_feat.SetField(j, feature.GetField(j))
        out_lyr.CreateFeature(out_feat)

    del out_ds
    del shp_ds
    print("Success........")

if __name__ == '__main__':
    ChangeToJson()