arcpy高德爬取路况信息数据json转shp

本文主要是介绍arcpy高德爬取路况信息数据json转shp，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

最近工作上遇到爬取的高德路况信息数据需要在地图上展示出来，由于json数据不具备直接可视化的能力，又联想到前两个月学习了一点点arcpy的知识，就花了一些时间去写了个代码，毕竟手动处理要了老命了。

1、json文件解读

json文件显示如下：

通过api和json组织结构解读，我们只需要道路的polyline坐标点空间信息，道路名称，道路方向、道路限速和道路拥堵程度等属性信息。

所以我们需要用python代码去获取roads数组内的信息

2、python代码编写

2.1、json文件读取

读取json文件需要线导入json模块，根据需要，我们还导入了arcpy等模块

import json
import os
import pandas as pd
import math
import csv
import arcpy
from arcpy import da
from arcpy import env

###打开json文件###
with open(jsonfile) as f:fileread=f.readline()jsondata=json.loads(fileread)print(jsondata)trafficinfo = jsondata['trafficinfo']for row in trafficinfo['roads']:name = row['name']status = row['status']direction = row['direction']angle = row['angle']if "speed" in trafficinfo:speed = row['speed']lcodes = row['lcodes']print ("--------------这里是坐标点的信息------------------------")points_data = row['polyline']points_data = points_data.split(";")#将polyline的坐标点循环遍历取出来for i in range(1,len(points_data)):print (i)print (points_data[i])oldLon = float(points_data[i].split(",")[0])oldLat = float(points_data[i].split(",")[1])lon = gcj02_to_wgs84(oldLon,oldLat)[0]lat = gcj02_to_wgs84(oldLon,oldLat)[1]pnt = arcpy.Point(lon,lat)array.add(pnt)polyline = arcpy.Polyline(array,cgcs2000project)array.removeAll()newFields = [polyline,name,status,direction,angle,speed,lcodes,polyline.length]cur.insertRow(newFields)del cur    f.close

2.2、创建shp数据

###创建shp文件###
print ("创建shp")
#创建一个空白的线要素
arcpy.CreateFeatureclass_management(os.path.dirname(shp_output),os.path.basename(shp_output),'POLYLINE',"","","",cgcs2000project)
#给shp添加字段属性
arcpy.AddField_management(shp_output,'name','Text')
arcpy.AddField_management(shp_output,'status','Text')
arcpy.AddField_management(shp_output,"direction",'Text')
arcpy.AddField_management(shp_output,'angle','Text')
arcpy.AddField_management(shp_output,"speed",'Text')
arcpy.AddField_management(shp_output,"lcodes",'Text')
arcpy.AddField_management(shp_output,"Length",'Double')
fields =['SHAPE@',"name", "status", "direction", "angle", "speed", "lcodes",'Length']
#插入行
cur = da.InsertCursor(shp_output,fields)

其中outputpath 是输出文件夹路径，
jsonfile是json文件路径
shp_output是输出的文件名称
cgcs2000project是坐标参照