Python中shapefile转换geojson的示例
shapefile转换geojson
import shapefile import codecs from json import dumps # read the shapefile def shp2geo(file="line出产.shp"): reader = shapefile.Reader(file) fields = reader.fields[1:] field_names = [field[0] for field in fields] buffer = [] for sr in reader.shapeRecords(): record = sr.record record = [r.decode('gb2312', 'ignore') if isinstance(r, bytes) else r for r in record] atr = dict(zip(field_names, record)) geom = sr.shape.__geo_interface__ buffer.append(dict(type="Feature", geometry=geom, properties=atr)) # write the GeoJSON file geojson = codecs.open(file.split('.')[0] + "-geo.json", "w", encoding="gb2312") geojson.write(dumps({"type": "FeatureCollection", "features": buffer}, indent=2) + "\n") geojson.close() if __name__ == '__main__': # import os # for z,x,c in os.walk('.'): # for zz in c: # if zz.endswith(".shp"): # shp2geo(zz) # shp2geo(file='D.shp') shp2geo(file='ttttttttttt.shp')
以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,谢谢大家对安科网的支持。如果你想了解更多相关内容请查看下面相关链接
相关推荐
sailxu00 2020-06-18
somebodyoneday 2020-06-06
newthon 2020-02-21
geojson需要先制作shp,然后导入下面网站生成geojson。geojson,最好放后台,前台通过异步请求去加载json,然后显示。this.map.flyTo;将地图的视野范围添加到相应区域。
somebodyoneday 2020-02-17
baijinswpu 2020-01-28
fengchao000 2019-10-31
cfh00 2019-07-01
fengchao000 2017-06-28
qianqianxiao 2019-06-29
mingyangwang 2018-12-27
cfh00 2019-06-27
yigeng 2014-09-05
lnn 2019-05-29
kkpiece 2019-03-15
yangguangdblu 2019-03-01
LHpython 2019-01-03