Python中shapefile轉(zhuǎn)換geojson的示例-創(chuàng)新互聯(lián)

shapefile轉(zhuǎn)換geojson

創(chuàng)新互聯(lián)建站于2013年開始，先為印臺等服務建站，印臺等地企業(yè)，進行企業(yè)商務咨詢服務。為印臺企業(yè)網(wǎng)站制作PC+手機+微官網(wǎng)三網(wǎng)同步一站式服務解決您的所有建站問題。

import shapefile
import codecs
from json import dumps
# read the shapefile
def shp2geo(file="line出產(chǎn).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')

當前名稱：Python中shapefile轉(zhuǎn)換geojson的示例-創(chuàng)新互聯(lián)
標題來源：http://muchs.cn/article6/cdicog.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供服務器托管、電子商務、企業(yè)網(wǎng)站制作、定制開發(fā)、外貿(mào)建站、用戶體驗

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)