利用Python的folium包绘制城市道路图

2023-12-01 14:10

本文主要是介绍利用Python的folium包绘制城市道路图,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

写在前面

很长一段时间内,我都在研究在线地图的开发者文档,百度地图和高德地图的开发者中心提供了丰富的在线地图服务,虽然有一定的权限限制,但不得不说,还是给我的科研工作提供了特别方便的工具,在博客前面我先放上这两个在线地图开放平台的web API的地址链接:

  • 百度地图开放平台
  • 高德地图开放平台

基于这两个平台,博主进行了一系列的开发研究工作,本文介绍其中一项技术,如何用folium包绘制城市道路图,当然,也可绘制非城市道路图,只要提供正确的路名就行了。
开发工具:

  • Python3.7
  • Spyder编译器(也可以用pycharm,不过建议用Spyder,因为编译过程中产生的变量太多,基本上都是json数据,我都是一边看一边写,这里Spyder优势明显)
  • chrome浏览器

folium介绍及相关设置

folium基础功能

folium的开发包在这里
简单来说,它是一个地理信息可视化的包,目前除了pyecharts,我用的最多的就是这个包,支持在在线地图上添加点、线、面等要素,而且还支持画热力图,不过热力图的效果真心不咋地,我看中的是它添加点、线、面形状的功能,而且各种要素可以设置颜色、大小、文字标记等属性,可视化效果还是不错的。这篇博客,也是应用了它画线的功能,绘制道路轮廓线。
在这里插入图片描述
上图图片来源:https://www.jianshu.com/p/32ec6afcc7a6?utm_campaign=hugo
比较遗憾的是,目前folium支持的地图底图有限,像openstreetmap是支持的,高德也支持,但不支持百度地图。当然,不支持百度地图并不能限制咱们开发者的脚步,本博客中也有相关介绍。
pip install folium
使用之前,先安装一下这个包。

folium参数设置

先看两行代码:

import folium
line_road = folium.Map(location=[31.596730,120.233516],zoom_start=15,tiles = 'http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=7&x={x}&y={y}&z={z}',attr = 'default')
  • location参数,设置展示地图的中心坐标点,就是说,比如你想看无锡市,可以设置成无锡市市中心的经纬度坐标
  • zoom_start是地图缩放等级,最高差不多可以到19还是20,如果想看大场景,就设小一点,想看局部地图就设大一点
  • tiles这个参数很重要,设置的是你的地图格式,默认的是OpenStreetMap,我这里把它改成了’http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=7&x={x}&y={y}&z={z}’,表示我用的是高德地图作为底图。为啥要改呢,因为每个不同的地图公司,用的坐标系不一样,高德地图和google地图、soso地图、aliyun地图、mapabc地图所用坐标相同,都是国测局(GCJ02)坐标,和百度地图用的坐标系不一样,如果直接拿百度坐标系下的经纬度画在高德地图上,那就会整体偏移,使用之前必须进行坐标转换
  • 其他用默认参数

获取道路参数

本博客的目的是画道路的轮廓图,首先必须得有数据才能画图。博主知道,目前这些地图公司,都是自己把车在开在路上去采集路上的经纬度,只要我获取到了这些经纬度,那我不就能绘制道路了吗?本着这个想法,我就到处搜索资料,开发者的潜能是无限了,同样在CSDN上我找到了一篇博客,链接在这里 ,真的很棒,不过他是用js写的,无所谓,方法是通的,用这篇博文提供的接口,真的实现了在地图上绘制道路的功能。
但是,但是,但是,,,
用了一段时间后,这个功能被封了,为此,我特意联系了高德地图开发者中心,他们的解释如下:
在这里插入图片描述

意思就是,这个功能用不了了,花钱也别想用。
其实很正常,这个功能太牛逼了。
虽然用不了了,但我还是介绍一下怎么实现的,万一以后又能用了呢。

高德地图获取道路经纬度的API介绍

先看接口:

http://restapi.amap.com/v3/road/roadname?parameters

这个接口和高德地图其他功能的接口一样,后面的parameters是需要写的参数,每个参数之间用&隔开,其中keywords是道路名,这个参数必填,当然还有key也是必填的。现在来看看这一段的具体怎么写,比如我想获取的是无锡市钱荣路的经纬度:

# -*- 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 nproad = quote('钱荣路')
key = YourKey  # 换成你自己申请的key
url = 'http://restapi.amap.com/v3/road/roadname?city=0510&key=%s&keywords=%s' % (key, road)req = urlopen(url)
res = req.read().decode()
temp = json.loads(res)roads = temp['roads']
pos = []
# 由于道路可能分段,比如钱荣路会分成钱荣路普通段和钱荣路高架,这都属于钱荣路的路段,因此必须要都取出来
for p in range(len(roads)):pos.extend(roads[p]['polylines'])pos_cal = []
line_qrroad = folium.Map(location=[31.596730,120.233516],zoom_start=15,tiles = 'http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=7&x={x}&y={y}&z={z}',attr = 'default')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)folium.PolyLine(lat_lon,weight = 5, color = 'red',opacity = 0.8).add_to(line_qrroad)
line_qrroad.save('lineqrroad.html')    map_qrroad = folium.Map(location=[31.596730,120.233516],zoom_start=15,tiles = 'http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=7&x={x}&y={y}&z={z}',attr = 'default')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)map_qrroad.save('render.html')

现在这段程序,已经没法解析出经纬度了,运行的话,会报如下错误:

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 osdef 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?locations=%f,%f&coordsys=baidu'\'&key=YourKey'%(data[i][j][1], data[i][j][0])req = urlopen(url)res = req.read().decode()temp = json.loads(res)location = temp['locations'].split(',')location = list(map(float,location))location[0], location[1] = location[1], location[0] poly.append(location)polylines.append(poly)return polylines

常规的接口访问和数据分析代码,不做过多解释,不过给大家看一下函数的输入格式:
在这里插入图片描述
结合数据格式,大家应该能看明白这段代码
用相同的方法来画地图,看看结果
在这里插入图片描述
不错哦,再看看细节:
在这里插入图片描述
细节不如之前丰富,不过也很不错了。
再给大家看看,如果直接用百度坐标系下的经纬度点画到高德地图上是个啥效果:
在这里插入图片描述
看到没,整体偏了不少,所以坐标转换很重要。。

结语

写这篇博客,也是为了记录最近一段时间的工作。这篇博客实际上没解决任何问题,根本没有取到任何有用的数据。但是,其中坐标转换那个函数还是有点用的。里面的数据分析方法,也是博主花了一些心血才写出来的。
不断学习,不断进步吧。

这篇关于利用Python的folium包绘制城市道路图的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/441468

相关文章

python: 多模块(.py)中全局变量的导入

文章目录 global关键字可变类型和不可变类型数据的内存地址单模块(单个py文件)的全局变量示例总结 多模块(多个py文件)的全局变量from x import x导入全局变量示例 import x导入全局变量示例 总结 global关键字 global 的作用范围是模块(.py)级别: 当你在一个模块(文件)中使用 global 声明变量时,这个变量只在该模块的全局命名空

【Python编程】Linux创建虚拟环境并配置与notebook相连接

1.创建 使用 venv 创建虚拟环境。例如,在当前目录下创建一个名为 myenv 的虚拟环境: python3 -m venv myenv 2.激活 激活虚拟环境使其成为当前终端会话的活动环境。运行: source myenv/bin/activate 3.与notebook连接 在虚拟环境中,使用 pip 安装 Jupyter 和 ipykernel: pip instal

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

【WebGPU Unleashed】1.1 绘制三角形

一部2024新的WebGPU教程,作者Shi Yan。内容很好,翻译过来与大家共享,内容上会有改动,加上自己的理解。更多精彩内容尽在 dt.sim3d.cn ,关注公众号【sky的数孪技术】,技术交流、源码下载请添加微信号:digital_twin123 在 3D 渲染领域,三角形是最基本的绘制元素。在这里,我们将学习如何绘制单个三角形。接下来我们将制作一个简单的着色器来定义三角形内的像素

nudepy,一个有趣的 Python 库!

更多资料获取 📚 个人网站:ipengtao.com 大家好,今天为大家分享一个有趣的 Python 库 - nudepy。 Github地址:https://github.com/hhatto/nude.py 在图像处理和计算机视觉应用中,检测图像中的不适当内容(例如裸露图像)是一个重要的任务。nudepy 是一个基于 Python 的库,专门用于检测图像中的不适当内容。该

Flutter 进阶:绘制加载动画

绘制加载动画:由小圆组成的大圆 1. 定义 LoadingScreen 类2. 实现 _LoadingScreenState 类3. 定义 LoadingPainter 类4. 总结 实现加载动画 我们需要定义两个类:LoadingScreen 和 LoadingPainter。LoadingScreen 负责控制动画的状态,而 LoadingPainter 则负责绘制动画。

pip-tools:打造可重复、可控的 Python 开发环境,解决依赖关系,让代码更稳定

在 Python 开发中,管理依赖关系是一项繁琐且容易出错的任务。手动更新依赖版本、处理冲突、确保一致性等等,都可能让开发者感到头疼。而 pip-tools 为开发者提供了一套稳定可靠的解决方案。 什么是 pip-tools? pip-tools 是一组命令行工具,旨在简化 Python 依赖关系的管理,确保项目环境的稳定性和可重复性。它主要包含两个核心工具:pip-compile 和 pip

HTML提交表单给python

python 代码 from flask import Flask, request, render_template, redirect, url_forapp = Flask(__name__)@app.route('/')def form():# 渲染表单页面return render_template('./index.html')@app.route('/submit_form',

利用matlab bar函数绘制较为复杂的柱状图,并在图中进行适当标注

示例代码和结果如下:小疑问:如何自动选择合适的坐标位置对柱状图的数值大小进行标注?😂 clear; close all;x = 1:3;aa=[28.6321521955954 26.2453660695847 21.69102348512086.93747104431360 6.25442246899816 3.342835958564245.51365061796319 4.87