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是坐标参照

2.3、坐标转换

由于高德数据是火星坐标系，我们需要对坐标点进行坐标系转换，将火星坐标系转为CGCS2000坐标系，转换代码在这：

3、转换后效果