python 利用gdal库 shp文件转geojson

代码转载（我在最下方复制了代码）：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()