pip install folium

import folium
line_road = folium.Map(location=[31.596730,120.233516],zoom_start=15,
      tiles = 'http://webrd02.is.autonavi.com/appmaptile"external nofollow" target="_blank" href="https://www.jb51.net/article/193956.htm">链接在这里 ，真的很棒，不过他是用js写的，无所谓，方法是通的，用这篇博文提供的接口，真的实现了在地图上绘制道路的功能。


但是，但是，但是，，，


用了一段时间后，这个功能被封了，为此，我特意联系了高德地图开发者中心，他们的解释如下：



意思就是，这个功能用不了了，花钱也别想用。

其实很正常，这个功能太牛逼了。

虽然用不了了，但我还是介绍一下怎么实现的，万一以后又能用了呢。
高德地图获取道路经纬度的API介绍
先看接口：


http://restapi.amap.com/v3/road/roadname"htmlcode">

# -*- coding: utf-8 -*-
"""
Created on Mon Mar 30 16:54:32 2020

@author: HP
"""

import json
import pandas as pd
from urllib.request import urlopen, quote
import folium
import numpy as np

road = quote('钱荣路')
key = YourKey # 换成你自己申请的key
url = 'http://restapi.amap.com/v3/road/roadname"htmlcode">

runfile('D:/python/folium/qianrongroad.py', wdir='D:/python/folium')
Traceback (most recent call last):

 File "D:\python\folium\qianrongroad.py", line 26, in <module>
 roads = temp['roads']

KeyError: 'roads'



意思就是说，没有‘road'这个key，我试图把请求串输入浏览器，返回的结果如下：

{"info":"INSUFFICIENT_PRIVILEGES","infocode":"10012","status":"0","sec_code_debug":"d41d8cd98f00b204e9800998ecf8427e","key":"ea12ed719e4ed13862dd0876384c6512","sec_code":"d41d8cd98f00b204e9800998ecf8427e"}



说我没有足够的权限。

好了，暂且不说了，看看代码的意思

前面是常规的json数据解析，没啥好说的，只要接口正常，就能取出数据来。


# 由于道路可能分段，比如钱荣路会分成钱荣路普通段和钱荣路高架，这都属于钱荣路的路段，因此必须要都取出来
for p in range(len(roads)):
 pos.extend(roads[p]['polylines'])



上面这个循环，注释已经解释清楚了，一条路可能会被高德分成好几部分，当然这是科学的，比如完整的钱荣路是分成了钱荣路普通路段和高架路段的。也就是说解析出来的roads的长度是2，分别是roads[0]和roads[1]，而经纬度数据则在roads[p][‘polylines']里面。


for i in range(len(pos)):
 m = pos[i].split(';')
 lat_lon = []
 for j in range(len(m)):
  n = m[j].split(',')
  n = list(map(float, n))
  n[0],n[1] = n[1],n[0]
  lat_lon.append(n)
  pos_cal.append(n)



这里是数据的分析，看起来写的很简单，其实很复杂，可惜没有数据来配套解释了。首先，经度和纬度之间用的是',‘分割，每一小段路之间用的是';‘来分割，这个';'分割我理解为高德对数据的一种加密方式，完整的一条路被高德划分成了很多小段，我把数据取出来后，自己用matplotlib演示了一下完整的路的绘制过程，看下面几张动图（没法插视频）：








为什么我这里要这么做，因为只有这样，我才是真正的理解了这些解析出来的数据是怎么连成一条完整的道路的，这样才好到folium中去绘制道路，实际上就是循环绘制，每一小段一小段的画，最后会连成一条完整的道路，过程就是下面这段代码：


for point in range(len(pos_cal)):
 folium.CircleMarker(location=[pos_cal[point][0],pos_cal[point][1]],
      radius=4,popup='popup',
      color='red',fill=True,
      fill_color='red').add_to(map_qrroad)



循环可以简化，博主习惯了写range(len)这种方式

渲染成网页，就可以打开了，看下结果：





忽略图中的圆圈标记，是我添加的其他信息。

放大看细节：





很良心有木有，双向车道、辅道、支路全部都有了，可惜当初没有把数据保存下来，只保存了这么个图。

这样就完事儿了。


自从高德把这个接口封了之后，博主神伤了好久，想了各种办法，连付费使用都想出来了，但是高德一个字，不给用、没权限、有钱也不行。没办法，项目还要继续，功能还要继续实现。想到之前百度地图事业部某年轻有为的负责人来咱们单位交流过，一番交涉，发现百度地图API也没有公开这个功能，但是离线地图可以。于是，，，博主又开始忙活了。

因涉及相关隐私，博主不具体介绍。总之，一番操作，获取到了百度地图坐标系下的道路经纬度数据，但前面说了，folium不支持百度地图，强行用百度地图坐标系下的经纬度坐标数据是会出乱子的，但这点小问题难不倒博主，高德地图API有坐标转换的接口呢。
百度坐标系下的坐标点转换成高德坐标系下的坐标点
接口在这里


接着上代码


import json
from urllib.request import urlopen, quote
import folium
import os

 
def BaiduMap2AMap(data):
 polylines = []
 for i in range(len(data)):
  poly = []
  for j in range(len(data[i])):
   url = 'https://restapi.amap.com/v3/assistant'    '/coordinate/convert"text-align: center">


结合数据格式，大家应该能看明白这段代码

用相同的方法来画地图，看看结果





不错哦，再看看细节：





细节不如之前丰富，不过也很不错了。

再给大家看看，如果直接用百度坐标系下的经纬度点画到高德地图上是个啥效果：





看到没，整体偏了不少，所以坐标转换很重要。。
结语

                                
                                

                                    Python,folium城市道路图,Python,folium绘制 
                                
DDR爱好者之家 Design By 杰米