Python：Shapefile矢量转化为GeoJSON格式_python

概述在最近的项目中，完成了许多python处理矢量数据的算法程序，比如缓冲区分析、叠置分析、统计分析等，主要用到的是GDAL/OGR库，很多功能都参照了此链接中的示例：http://pcjericks.github.io/py-gdalogr-cookbook/# 由于完成的功能在上面链接中都能找到实现的代码，所以原来不打算写博客的，但是今天写的这个SHP转GeoJSON的算法把我坑得不轻，着实郁闷了

在最近的项目中，完成了许多python处理矢量数据的算法程序，比如缓冲区分析、叠置分析、统计分析等，主要用到的是GDAL/OGR库，很多功能都参照了此链接中的示例：http://pcjericks.github.io/py-gdalogr-cookbook/#

由于完成的功能在上面链接中都能找到实现的代码，所以原来不打算写博客的，但是今天写的这个SHP转GeoJsON的算法把我坑得不轻，着实郁闷了好久才解决，因此就写一下权当纪念吧。有时候写程序真的是一件充满玄学的事情，原理上怎么看怎么对，但是运行起来却偏偏不是那回事，就算调试也给你不报错直接闪退，特别是处理空间数据时属性的编码问题，更是搞得人头皮发麻。这种情况下，能不能调的出来全靠运气，搞不好一个小问题卡一个星期都是正常，这可能也是程序员普遍发量稀少的原因吧。更加玄学的是，解决问题的方式又常常充满戏剧性，折腾了一个早上加半个下午没找到原因，起来喝杯水却突然灵感乍现，三四行代码调整一下顺序，一顿骚 *** 作，点击运行，突然就能运行成功了，再看一下结果，卧槽! 还真是对的！山穷水尽疑无路，柳暗花明又一村，这种感觉真是刺激。

好了，将Shapefile转化为GeoJsON的完整代码如下：

 1 # -*- Coding: utf-8 -*- 2 from osgeo import ogr 3 import gdal 4 import sys 5 import os 6  7 def ChangetoJson(vector,output): 8     print("Starting........") 9     #打开矢量图层10     gdal.SetConfigOption("GDAL_filename_IS_UTF8","YES")11     gdal.SetConfigOption("SHAPE_ENCoding","GBK")12     shp_ds = ogr.Open(vector)13     shp_lyr = shp_ds.GetLayer(0)14 15     # 创建结果GeoJson16     basename = os.path.basename(output)17     out_driver = ogr.GetDriverByname(‘GeoJsON‘)18     out_ds = out_driver.CreateDataSource(output)19     if out_ds.GetLayer(basename):20         out_ds.DeleteLayer(basename)21     out_lyr = out_ds.CreateLayer(basename,shp_lyr.GetSpatialRef())22     out_lyr.CreateFIElds(shp_lyr.schema)23     out_feat = ogr.Feature(out_lyr.GetLayerDefn())24 25     #生成结果文件26     for feature in shp_lyr:27         out_feat.SetGeometry(feature.geometry())28         for j in range(feature.GetFIEldCount()):29             out_feat.SetFIEld(j,feature.GetFIEld(j))30         out_lyr.CreateFeature(out_feat)31 32     del out_ds33     del shp_ds34     print("Success........")35 36 if __name__ == ‘__main__‘:37     shapefile = ‘D:/python37Work/HainanAOI/affected_brIDge.shp‘38     out = ‘D:/python37Work/HainanAOI/affected_brIDge.Json‘39     ChangetoJson(shapefile,out)