【数据处理】高德地图路径规划数据转存shp以及json转geojson介绍

已于 2024-11-18 16:45:38 修改
阅读量1.2k 收藏 20
点赞数 32
分类专栏： GIS 高德地图 geopandas 文章标签： json
于 2024-11-17 08:30:30 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kowkhing/article/details/143827183
版权

 

大家好，我是小白鸽，上一期文章挖了一个坑，留了两个问题，分别是如何把json转为GeoJSON以及如何把高德路径规划结果转存成矢量数据，本期我们就来详细说说这两个问题。

其实要解决这两个问题的方法挺多的，但是进行数据分析首先要做的应该是分析数据结构，我们先来看一下路径规划返回的结果示例。

{
  "status" :"1",
  "info" :"ok",
  "infocode" :"10000",
  "count" :"1",
  "route" :{
    "origin" :"116.481028,39.989643",
    "destination" :"116.434446,39.90816",
    "paths" :[
      "0" :{
        "distance" :"11068",
        "duration" :"8854",
        "steps" :[
          "0" :{
            "instruction" :"向南步行16米左转",
            "orientation" :"南",
            "road" :[ ],
            "distance" :"16",
            "duration" :"13",
            "polyline" :"116.480885,39.989371;116.480907,39.989353;116.48089,39.989227",
            "action" :"左转",
            "assistant_action" :[ ],
            "walk_type" :"0"
          },
          "1" :{
            "instruction" :"步行292米右转",
            "orientation" :[ ],
            "road" :[ ],
            "distance" :"292",
            "duration" :"234",
            "polyline" :"116.480885,39.989223;116.480933,39.989201;116.480933,39.989201;116.480955,39.989184;116.481155,39.989058;116.481155,39.989058;116.481428,39.98888;116.481428,39.98888;116.481489,39.988841;116.481489,39.988841;116.482574,39.988125;116.48322,39.987665;116.48322,39.987665;116.48332,39.987565;116.48332,39.987565;116.48342,39.987435",
            "action" :"右转",
            "assistant_action" :[ ],
            "walk_type" :"0"
          }
        ]
      }
    ]
  }
}

其实，如果数据量不大，可以手动存成表格数据，再通过ArcGis转存成矢量数据。

e4b73249a89d42079377a12f6e368c68.png

那如果数据量很大呢？假设我们有20000个村采集的数据需要转存成矢量数据，那要一个个手动填吗？

 

6eba41ff2f2b1f89cb18dba02194ca03.jpeg

大可不必，可以看到路径规划返回的结果是一个json对象，我们本期的目标是把steps里面的内容解析出来，steps里面的其实是一个json对象数组，数组里面每一个元素都是json对象。json数据是以键值对的形式组织数据，每个键值对由一个键（key）和一个值（value）组成，它们之间用冒号 : 分隔。例如：

#name是键，John Doe是值{  "name": "John Doe"}

因此对于复杂的json对象也可以通过逐级索引的方式获取目标数据。以下是一个示例：

import json

# 假设这是你的JSON字符串
your_json_string = '''
{
  "person": {
    "name": "John Doe",
    "age": 30,
    "address": {
      "street": "123 Main St",
      "city": "Anytown",
      "zipCode": "12345"
    },
    "phoneNumbers": [
      {
        "type": "home",
        "number": "212 555-1234"
      },
      {
        "type": "office",
        "number": "212 555-4567"
      }
    ]
  }
}
'''

# 解析JSON字符串
json_data = json.loads(your_json_string)

# 逐层索引以获取办公室电话号码
office_number = json_data["person"]["phoneNumbers"][1]["number"]

print(office_number)  # 输出: 212 555-4567

获取到目标数据只是第一步，开头说了要把数据转换成GeoJSON，以及转存矢量数据，本文在此介绍3个python库（模块），分别是arcpy、geojson和geopandas，其中arcpy库和geopandas库都可以实现将结果数据直接转存成矢量数据，geojson库则可以实现将结果数据转换成GeoJSON格式，并且可以进一步结合geopandas库将GeoJSON格式数据转存成矢量数据。

 

b8d696f9de252688b6ada9a9a338f9a9.png

 

    利用arcpy转存矢量数据

arcpy其实是一个python库，通过引用（import）arcpy，可以让我们在python中使用arcgis的空间分析和数据管理等能力，相当于在python中使用ArcGis。本文在此只做简要介绍，以后将会开个专栏讲arcpy的使用以及如何利用arcpy去处理一些日常遇到的问题。

c8cbb4ef836a46028b2b2cd75d6d1dde.png

官方介绍请查阅：

https://desktop.arcgis.com/zh-cn/arcmap/latest/analyze/arcpy/what-is-arcpy-.htm

以下是利用arcpy将路径规划结果中的分段数据转存成矢量数据的示例：

import requests
import json
import arcpy

key = '你的高德地图api key'
#起点经纬度
origin = '116.481028,39.989643'
#终点经纬度
destination = '116.434446,39.90816'
#拼接url，上一期输出是xml格式，本期尝试输出json格式
url = 'https://restapi.amap.com/v3/direction/walking?origin='+ origin +'&destination='+ destination +'&key=' + key +'&output=JSON'
#请求url
responce = requests.get(url)
#获取返回内容
str = responce.content
#将返回的字符串转换为json对象
jstr = json.loads(str)
#定位到steps里面的内容
steps = jstr['route']['paths'][0]['steps']
#创建空list存数据
distance,duration,linestr,road = [],[],[],[]
#读取steps里面每一段路径的信息，包括distance、duration、road和polyline
for step in steps:
    distance.append(step['distance'])
    duration.append(step['duration'])
    linestr.append(step['polyline'])
    road.append(step['road'])

path="你的文件输出路径"
#输出可以是XXXX.shp
outputname="你的输出文件名称"
#设置空间参考，一般为4326，即wgs84
spatR = arcpy.SpatialReference(4326)
#创建一个空的线要素
line = arcpy.CreateFeatureclass_management(path, outputname,'POLYLINE',"","","",spatR)
#给shp添加字段属性，包括distance,duration,road
arcpy.AddField_management(line,'distance','Short Integer')
arcpy.AddField_management(line,'duration','Short Integer')
arcpy.AddField_management(line,"road",'Text')
fields =['SHAPE@',"distance", "duration", "road"]
#创建游标，arcpy对数据的操作通过游标进行
cur = arcpy.InsertCursor(line)
#通过循环对数据进行操作
for i in range(0,len(linestr)):
    #准备好属性数据，实际上如果熟练可以把该循环整合到上一个循环
    distance_ = distance[i]
    duration_ = duration[i]
    road_ = road[i]
    #获取单挑线段的坐标集合
    line = linestr[i]
    #获取该坐标集合中所有的点
    pnts = line.split(';')
    #创建空的array
    array = arcpy.Array()
    #对坐标集合中的点进行遍历，创建Point，并添加到array中，形成符合arcgis格式要求的点集
    for pnt in pnts:
        pnt_ = pnt.split(',')
        lon = pnt_[0]
        lat = pnt_[1]
        point = arcpy.Point(lon,lat)
        array.add(point)
    #根据点集创建Polyline线
    polyline = arcpy.Polyline(array,spatR)
    #清空array，为下一轮循环做准备
    array.removeAll()
    #根据上面的fields格式要求创建newFields并添加到原来的fields中
    newFields = [polyline,distance_,duration_,road_]
    cur.insertRow(newFields)
del cur

    利用geopandas转存矢量数据

geopandas也是一个python库，它是对pandas库的扩展，通过geopandas也可以进行空间分析和数据处理，与arcpy不同的是，不需要安装ArcGis就可以使用这些空间分析能力，它的易用性更好。

 

44919ba7c905b9494dd2e8511dc17d0c.png

官方介绍请查阅：

https://desktop.arcgis.com/zh-cn/arcmap/latest/analyze/arcpy/what-is-arcpy-.htm

以下是利用geopandas将路径规划结果中的分段数据转存成矢量数据的示例：

import requests
import json
from shapely.geometry import  LineString;
from geopandas import GeoSeries,GeoDataFrame
key = '你的高德地图api key'
#起点经纬度
origin = '116.481028,39.989643'
#终点经纬度
destination = '116.434446,39.90816'
#拼接url，上一期输出是xml格式，本期尝试输出json格式
url = 'https://restapi.amap.com/v3/direction/walking?origin='+ origin +'&destination='+ destination +'&key=' + key +'&output=JSON'
#请求url
responce = requests.get(url)
#获取返回内容
str = responce.content
#将返回的字符串转换为json对象
jstr = json.loads(str)
#定位到steps里面的内容
steps = jstr['route']['paths'][0]['steps']
#创建空list存数据
distance,duration,polyline,road = [],[],[],[]
#读取steps里面每一段路径的信息，包括distance、duration、road和polyline
for step in steps:
    distance.append(step['distance'])
    duration.append(step['duration'])
    polyline.append(step['polyline'])
    road.append(step['road'])
#创建一个空的GeoSeries
g = GeoSeries()
#通过循环向GeoSeries中添加数据
for i in range(0,len(polyline)):
    line = polyline[i]
    newstr = line.split(';')
    location_tuples = [tuple(map(float, loc.split(','))) for loc in newstr]
    pl = LineString(location_tuples)
    new_geoseries = GeoSeries(pl)
    g = pd.concat([g, new_geoseries])
#为避免出错，重置一下索引
g.reset_index(drop=True,inplace=True)
#将数据组装成GeoDataFrame，其中{'distance' :distance,'duration': duration,'road': road}代表属性信息，geometry=g代表空间信息，crs='epsg:4326'代表投影信息
gdf = GeoDataFrame({'distance' :distance,'duration': duration,'road': road}, geometry=g,crs='epsg:4326')
#保存成shp
gdf.to_file('line.shp')

    结合geojson和geopandas转存矢量数据

geojson是python的一个模块，geojson模块支持GeoJSON标准定义的所有数据结构，包括几何对象（点、线、多边形等）、特征、特征集合和坐标参考系统。该模块提供了函数来解析和序列化GeoJSON数据，包括从字符串和文件中解析GeoJSON数据，以及将Python对象序列化为GeoJSON字符串或写入文件。

422fdecd220643089892029d4a87d4b6.png

结合geojson和geopandas的处理方式分两步，第一步是将结果数据转换成GeoJSON格式，第二步是将GeoJSON格式数据转存成矢量数据。第一步时可以根据需求将多条线段转换成一个多段线要素（MultiLineString），也可以转换成多个线要素（LineString）。

MultiLineString和LineString的区别在于，MultiLineString一条记录可以对应多条线，LineString一条记录只对应一条线。

 

5aa528e0ab03c1a02edbcced28ade593.png

 

  • 转换成一个多段线要素示例：

import requests
import json
import geojson
key = '你的高德地图api key'
#起点经纬度
origin = '116.481028,39.989643'
#终点经纬度
destination = '116.434446,39.90816'
#拼接url，上一期输出是xml格式，本期尝试输出json格式
url = 'https://restapi.amap.com/v3/direction/walking?origin='+ origin +'&destination='+ destination +'&key=' + key +'&output=JSON'
#请求url
responce = requests.get(url)
#获取返回内容
str = responce.content
#将返回的字符串转换为json对象
jstr = json.loads(str)
#定位到steps里面的内容
steps = jstr['route']['paths'][0]['steps']
#创建空list存数据
distance,duration,polyline,road = [],[],[],[]
#读取steps里面每一段路径的信息，包括distance、duration、road和polyline
for step in steps:
    distance.append(step['distance'])
    duration.append(step['duration'])
    polyline.append(step['polyline'])
    road.append(step['road'])

line_list = []
#数据清洗，重新组织数据格式
for line in polyline:
    newstr = line.split(';')
    location_tuples = [tuple(map(float, loc.split(','))) for loc in newstr]
    line_list.append(location_tuples)
#创建MultiLineString
MultiLineString_geojson = geojson.MultiLineString(coordinates=line_list)
#创建FeatureCollection
Multi_feature_collection = geojson.FeatureCollection(features=[geojson.Feature(geometry=MultiLineString_geojson)])
#添加属性示例 属性名称name，属性值OD_name
Multi_feature_collection['features'][0]['properties']['name'] = 'OD_name'
#尝试打印
print(geojson.dumps(Multi_feature_collection, indent=2))

  • 转换成多个线要素示例

 

import requests
import json
import geojson
key = '你的高德地图api key'
#起点经纬度
origin = '116.481028,39.989643'
#终点经纬度
destination = '116.434446,39.90816'
#拼接url，上一期输出是xml格式，本期尝试输出json格式
url = 'https://restapi.amap.com/v3/direction/walking?origin='+ origin +'&destination='+ destination +'&key=' + key +'&output=JSON'
#请求url
responce = requests.get(url)
#获取返回内容
str = responce.content
#将返回的字符串转换为json对象
jstr = json.loads(str)
#定位到steps里面的内容
steps = jstr['route']['paths'][0]['steps']
#创建空list存数据
distance,duration,polyline,road = [],[],[],[]
#读取steps里面每一段路径的信息，包括distance、duration、road和polyline
for step in steps:
    distance.append(step['distance'])
    duration.append(step['duration'])
    polyline.append(step['polyline'])
    road.append(step['road'])

#创建一个空的LineString和FeatureCollection，用于存放后续数据
totalLineString_geojson = geojson.LineString()
total_feature_collection = geojson.FeatureCollection(features=[geojson.Feature(geometry=totalLineString_geojson)])
#通过循环将steps中每一段路径都转换成Feature，并添加进空的FeatureCollection
for i in range(0,len(polyline)):
    line = polyline[i]
    newstr = line.split(';')
    location_tuples = [tuple(map(float, loc.split(','))) for loc in newstr]
    one_LineString_geojson = geojson.LineString(coordinates=location_tuples)
    one_feature = geojson.Feature(geometry=one_LineString_geojson)
    #为每一段路径添加属性
    one_feature['properties']['distance'] = distance[i]
    one_feature['properties']['duration'] = duration[i]
    one_feature['properties']['road'] = road[i]
    #添加进FeatureCollection
    total_feature_collection['features'].append(one_feature)
#删除最初始的空feature
total_feature_collection['features'].remove(total_feature_collection['features'][0])
#尝试打印
print(geojson.dumps(total_feature_collection, indent=2))

  • 在上述代码基础上，只需如下代码即可将数据转存为矢量数据

#引用geopandas
import geopandas
# total_feature_collection 上一步生成的FeatureCollection
#将FeatureCollection转换成GeoDataFrame
#crs='epsg:4326'设置坐标为wgs84，最好提前进行坐标转换
total_gdf = geopandas.GeoDataFrame.from_features(total_feature_collection,crs='epsg:4326')
#重置一下索引
total_gdf.reset_index(drop=True,inplace=True)
#保存成line.shp
total_gdf.to_file('line.shp')

     结果展示

25bb119de7674ee09b231e5dd4e3f365.png

 

最后，其实数据处理的方法有很多，实现同一个目标也有多种不同的路径，最主要的是根据实际的需求选择适合自己的方法。本期分享就到此结束，我是小白鸽，欢迎大家点赞关注，欢迎评论交流。

 

663b7574c1da0100415f8f56be72e525.jpeg

 

高德地图中国省级边界线GEOJSON格式
07-13
完整的高德地图中国省级边界线，包含34个省/直辖市/自治区/行政区，坐标格式是WGS84,EPSG4326
shp转geojson工具.rar
06-26
"shp转geojson工具.rar" 提供了一个解决方案，帮助Java开发者处理与SHP文件相关的地理数据，将其转换为GeoJSON格式。这里我们将深入探讨SHP文件、GeoJSON以及如何在Java环境中进行这种转换。 1. **SHP文件**： - ...
基于高德API实现网络geoJSON功能（突出省份）
IIIIIAN
12-07 2199
【代码】基于高德API实现网络geoJSON功能（突出省份）
全国行政区划代码JSON数据包：Echarts地图可视化指南
最新发布
weixin_32661831的博客
03-16 1176
JSON（JavaScript Object Notation）作为一种轻量级的数据交换格式，由于其易于阅读和编写，同时易于机器解析和生成，已成为互联网应用中最为流行的数据交换语言之一。本章将带领读者初步了解JSON的基础知识和结构特性。Echarts，即Enterprise Charts，是一个使用 JavaScript 实现的开源可视化库，它由百度团队开发，并于2013年首次发布。Echarts 旨在提供简单易用的API，让开发者能够轻松地将图表集成到网页中。
利用高德api获取行政区边界geoJson，基于echart实现全国各省地磁分布统计图
杨丫丫的专栏
11-15 2404
利用高德api获取行政区边界geoJson，基于echart实现全国各省地磁分布统计图
【高德地图入门】--- 解析geojson
热门推荐
qq_39078783的博客
12-22 1万+
简介 GeoJSON是一种用于编码各种地理数据结构的数据。GeoJSON对象可以表示几何、特征或特征集合。GeoJSON支持以下几何类型：点（Point）、线（LineString）、面（Polygon）、多点（MultiPoint）、多线（MultiLineString）、多面（MultiPolygon）和几何集合（GeometryCollection）。GeoJSON中的功能包含几何对象和其他属性，特征集合表示一系列特性。 GeoJSON对象 GeoJSON总是由一个单独的对象组成。此对象（指的是下
JS 将Json转Geojson
Nanakris的博客
01-06 4674
JS 将Json转Geojson
基于高德API实现网络geoJSON功能（整体）
IIIIIAN
12-06 1053
【代码】基于高德API实现网络geoJSON功能。
python脚本实现shp文件转geojson文件
01-29
针对QGIS中shp转geojson数据容易造成中文乱码问题 提出针对python脚本的一键转换，在代码中只需要更改自己的文件地址，运行即可实现代码的转换。 脚本依靠GIS常用GDAL库实现，安装简单，熟悉基本的python语法即可...
geojson-geometry-shp格式相互转化总结
06-26
总结，GeoJSON和Geometry SHP格式的转换在Java开发中是一项常见的任务，涉及到数据类型转换、数据库操作以及地理空间数据的处理。理解并掌握这些转换方法，对于处理地理信息系统相关项目至关重要。
arcpy高德爬取路况信息数据json转shp
02-04
高德爬取路况信息数据可视化 利用arcpy将json文件转为shp数据，并且自动将火星坐标系转为CGCS2000坐标系 转换后的shp数据可直接在arcgis或者其他GIS平台种进行路况可视化
json2geojson:将 JSON 转换为 GeoJSON 的快速函数
06-01
JSON 到 GeoJSON 的转换 特征 将指定的 JSON 对象转换为具有指定初始属性的基本 GeoJSON 对象 用于向 GeoJSON 对象添加更多属性的补充函数 职能 这个库中有两个关键函数如下。 j2gConvert(jsonObject, 属性) jsonObject 是一个对象参数，它指的是您打算转换为 GeoJSON 的原始 JSON 对象。 property 是一个字符串参数，它指的是包含在 GeoJSON 对象的 'properties' 字段下的初始属性。 j2gAddproperty(geojson, jsonObject, add_property) geojson 是一个对象参数，指的是启动的 GeoJSON 对象。 jsonObject 是一个对象参数，它指的是您打算转换为 GeoJSON 的原始 JSON 对象。 add_property 是一个
最新成都geojson数据，含高新区、成都天府新区 适用于高德等地图数据
02-15
成都最新geojson数据，含高新东区、高新西区、高新南区、成都天府新区 适用于高德等地图数据
高德地图全国地址JSON文件（省、市、区、街道）
10-08
高德地图全国地址JSON文件（省、市、区、街道），前端可以直接一次性调用
高德地图城市json
02-08
用于遍历出各城市按钮，进行城市选择。如美团、饿了么等城市切换的城市数据源
高德地图，全国省级数据Json(包括：经纬度，省份Code，省份名称等)
01-08
高德地图，全国省级数据Json(包括：经纬度，省份Code，省份名称等)
利用Python实现Shp格式向GeoJSON的转换方法
09-19
本篇文章将详细介绍如何使用Python将Shapefile格式的数据转换为GeoJSON格式。 #### 二、GeoJSON格式说明 GeoJSON是一种用于存储地理空间数据的标准格式，它基于JSON格式。一个基本的GeoJSON文件由以下两部分组成：...
高德地址 AMap.GeoJSON解析geoJson并画出区域图 画出区域图标记出名称 获取地图的坐标
博观而约取，厚积而薄发~
01-08 1347
GeoJSON支持以下几何类型：点（Point）、线（LineString）、面（Polygon）、多点（MultiPoint）、多线（MultiLineString）、多面（MultiPolygon）和几何集合（GeometryCollection）。高德开放平台，矢量图形，使用GeoJSON根据标注GeoJSON配置在地图上对应覆盖物：https://lbs.amap.com/demo/javascript-api-v2/example/overlayers/geojson。
【高德地图进阶】--- 通过geoJson绘制城市版块
qq_39078783的博客
12-23 4422
我们已经学会了如何绘制点，线，面。今天我们就结合geoJson进行绘制 点，线，面。 通常的业务场景是希望绘制一个城市的版块，我们可以先拿到一个城市的geoJson（通常是只包含面（MultiPolygon）），然后通过高德的api绘制城市版块 如何获取快速城市的GeoJson数据 我们可以使用阿里云的DataV 可以通过点击选择想要的城市版块，这里以山东省为例。 复制链接，我们可以直接在项目中使用 全部代码展示 （其中xxxx 高德的秘钥，cccc是高德的key，需要自己申请） <!DOC.
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值