Python有关folium库在地图上批量标记地点的示例(基于高德地图)

文章正文
发布时间:2024-09-03 19:09

1. 根据地名查找对应的高德标准经纬度
2. 将百度标准经纬度转换成高德标准经纬度
3. 将excel中的经纬度等信息导入
4. 将高德标准经纬度对应的点位标注在Leaflet(网页)上
5. 将两个高德标准经纬度所连路段标注在Leaflet(网页)上

二、代码及使用方法 1. 本文所使用的python库 import folium from folium.plugins import FloatImage import pandas import requests import math

下载库:在cmd中输入以下代码或在pycham中下载。(File - Settings - Project:untitled - Poject Interpreter - "+" - 搜索库名)

pip install folium pip install pandas pip install requests 2. 由具体地名获得其对应高德标准经纬度 def gd_map(addr): para = {'key': 'xxxxxxxxxxxxxxxxxxx', 'address': addr} # key填入自己在高德开放平台上申请的key url = 'https://restapi.amap.com/v3/geocode/geo?' # 高德地图地理编码API服务地址 result = requests.get(url, para) # GET方式请求 result = result.json() lon_lat = result['geocodes'][0]['location'] # 获取经纬度 lon_lat = lon_lat.split(",") # 把经纬度以","为界分割 lon = float(lon_lat[0]) # 获取的经纬度格式为[经度,纬度] lat = float(lon_lat[1]) print(addr, "的高德标准经纬度为: (", lon, ",", lat, ")") return lon, lat # 返回经纬度(float) addr = '杭州市杭州东站' # 输入想要查询的文字地址(尽量详细) lon, lat = gd_map(addr) # 调用函数 函数名   gd_map  
传入数据   查询地址addr(str)  
返回数据   经度lon(float),纬度lat(float)  

首先要自己注册一个高德地图开放平台的开发者账号并进行认证(官方网址:高德开放平台 | 高德地图API) ,登陆后点击右上角头像旁边的“控制台”,进入高德控制台页面。点击左侧”应用管理“及下面的”我的应用“,点击右上角的”创建新应用“,设置好名称和类型。

 在我的应用中找到刚刚创建的应用,点击”添加key“。

之后在”服务平台“中勾选”Web服务“,我们主要需要的是其中的地理编码API。

然后把key下面的一串字符串复制下来,填入代码中第二行的'key'=' '内即可。

查询时在倒数第二行的addr = ' '内输入想查询的地点,该函数会返回其对应的经度lon(float)以及纬度lat(float)。

示例:addr = '杭州市杭州东站'

 

3. 将百度标准经纬度转换成高德标准经纬度 # 将百度标准经纬度转换成高德标准经纬度 def bd_gd(lon, lat): pi = 3.14159265358979324 * 3000.0 / 180.0 x = lon - 0.0065 y = lat - 0.006 z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * pi) theta = math.atan2(y, x) - 0.000003 * math.cos(x * pi) lon = z * math.cos(theta) lat = z * math.sin(theta) print("高德标准经纬度为:(", lon, ",", lat, ")") return lon, lat print("百度标准经纬度为:(", lon, ",", lat, ")") bd_gd(lon, lat) 函数名   bd_gd  
传入数据   百度经度lon(float),百度纬度lat(float)  
返回数据   高德经度lon(float),高德纬度lat(float)  

首先说一下该函数的使用背景。起初在尝试标记地点时发现经纬度标记后和目标地点有不小的误差,经过了解后发现当前主要有三种坐标系:WGS84、GCJ02(高德使用)和BD09(百度使用),相同的地点在三种坐标系下的经纬度坐标不同。

本文章使用的高德地图采用了GCJ02坐标系,因此如果你把在百度地图上获取的BD09经纬度标记在高德地图上就会出现偏差。本函数提供一个将BD09经纬度坐标转换到GCJ02经纬度坐标的方法。

使用时将百度坐标系下的经纬度传递给该函数即可。例:

lon_gd, lat_gd = bd_gd(120.212427, 30.291538)

本部分算法转载自:https://www.cnblogs.com/weihengblog/p/9995989.html

4. pandas库将excel中的经纬度等信息导入 def exl_load(): df = pandas.read_excel('经纬度表.xlsx') data_num = df['序号'].tolist() data_name = df['名称'].tolist() data_lon = df['经度'].tolist() data_lat = df['纬度'].tolist() data_times = df['次数'].tolist() return data_num, data_name, data_lon, data_lat, data_times num, name, lon, lat, times = exl_load() 函数名   exl_load  
传入数据    
返回数据   可根据实际情况自定义(list)  

因为要批量标记大量地点,所以一个一个标记过于麻烦。可以建立一个excel表格然后用pandas库中的方法调取表格中数据进行批量标记。表格中的第一行写上每一列所对应的数据名称(与代码中一致),示例如下图所示。

使用时可以将excel表格加入代码文件根目录中,这样就只需要在第二行df = pandas.read_excel('xxxx.xlsx')的单引号内写入文件名称即可。每列的信息会以列表(list)形式返回,可以用列表相关的方法进行数据处理。

5. folium库生成地图图层 # 设置地图 m = folium.Map( [35.033740, 110.873590], # 中心点经纬度(先纬后经) tiles='https://wprd01.is.autonavi.com/appmaptile?x={x}&y={y}&z={z}&lang=zh_cn&size=1&scl=1&style=7', # 地图风格 attr='高德-常规图', # 地图名称 zoom_start=15, # 默认比例尺 control_scale=True, # 是否在地图上添加比例尺 width='100%', # 缩放度 ) FloatImage('颜色.png', left=88, bottom=8, width=8).add_to(m) # 在地图上添加图片 FloatImage('标题.png', left=30, bottom=88, width=40).add_to(m) m.add_child(folium.LatLngPopup()) # 点击显示经纬度

 本段代码用于生成一个地图图层,并且能在地图上添加提示信息(以图片形式)。在FloatImage函数中,单引号内输入图片的路径,left和bottom分别在左侧和下侧确定图片位置,width则确定图像大小。

示例:

6. folium库将高德标准经纬度对应的点位标注在Leaflet(网页)上 # 设置标记点 for i in num: if times[i - 1] >= 76: color = 'red' elif 52 <= times[i - 1] < 76: color = 'orange' elif 35 <= times[i - 1] < 52: color = 'beige' else: color = 'green' folium.Marker( location=[lat[i - 1], lon[i - 1]], # 要标记的点经纬度(注意经纬度是反过来的) popup=name[i - 1], # 点击标记出现的文本 tooltip=times[i - 1], # 提示语(鼠标放到标记点上的提示语) icon=folium.Icon(color=color, icon='map-marker-alt'), # color设置颜色,icon设置标记点外形 ).add_to(m)

前面几行if条件语句是用来针对某一属性的值来设置不同的颜色的,比如在统计某些地点人流量大小的时候可以用红色表示爆满,橙色表示拥挤,黄色标识中等,绿色表示畅通。没有相关需求可以忽略,重点是后面的folium.Marker().add_to(m)。

使用时,在folium.Marker()中的location=[ , ]中填入想标记的高德标准经纬度就能实现在地图上标记点位。用for循环遍历4中得到的列表可以做到批量标记。

示例:

如果想用圆形标记,可以用folium.CircleMarker()实现。具体示例代码及效果图如下:

for i in num: if times[i - 1] >= 10: color = '#FF0000' # 颜色代码 elif 7 <= times[i - 1] < 10: color = '#FF8000' elif 4 <= times[i - 1] < 7: color = '#00FF00' elif 2 <= times[i - 1] < 4: color = '#00FFFF' else: color = '#0000FF' if lat[i - 1] == -1.0 or lon[i - 1] == -1.0: # 判断是否为空,为空则跳过本次循环 continue folium.CircleMarker( location=[lat[i - 1], lon[i - 1]], # 要标记的点经纬度(注意经纬度是反过来的) radius=10, # 圆形的半径 popup=folium.Popup(name[i - 1] + ":" + str(times[i - 1]) + "次", max_width=200), # 点击标记出现的文本 tooltip=times[i - 1], # 提示语(鼠标放到标记点上的提示语) color=color, # 外框颜色 fill=True, # 是否填充 fill_color=color # 填充颜色 ).add_to(m)

7. folium库将两个高德标准经纬度所连路段标注在Leaflet(网页)上 folium.PolyLine( locations=[[30.1944, 120.2347], [30.1940, 120.2666]], # 要标记的线段起始经纬度(注意经纬度是反过来的) popup="建设一路", # 点击标记出现的文本 color='red', # 颜色 weight=6, # 线段宽度 opacity=0.8, # 线段透明度(0~1) ).add_to(m)

作用是用线段来表示某一路段,同样也可以添加不同颜色来表示道路拥挤程度。本文没有编写实现相关功能的代码,但是可以参考上一条自行编写。

使用时,将需要标记的线段起始两组高德标准经纬度填入folium.PolyLine()中的locations=[ [ , ] , [ , ] ],即可得到一条标记线段。

示例:

locations=[[30.1944, 120.2347], [30.1940, 120.2666]]  # (建设一路)

8. 保存生成的地图为html # 保存地图 m.save('地图.html') print("地图已绘制完毕,请在左侧代码根目录中用浏览器打开生成的地图网址查看。")

在第一行可以更改html的保存路径以及名称。如果在pycharm中像这样直接保存,可以在左侧根目录中找到html。右键并选择在浏览器打开就能看到生成的带有标记的地图了。也可以选择将html存放在其他路径下。

三、总结及补充说明

1. folium库在标记经纬度时要先写纬度后写经度,这点在5.6.7三个部分上均有体现,一定要注意。

2. 本文各部分代码可以根据实际需求选择性参考,有些也可以连起来实现一些其他的功能。

例1:excel中获取到的地点名称列表可以再导入到2中进行批量查询经纬度,再用6批量标记。

例2:手中的excel数据集中的经纬度采用的是百度坐标系,可以将其导出为list后再导入进3中转为高德坐标系经纬度,最后再用6进行批量标记。

3. folium库官方文档:Folium — Folium unknown 文档

首页
评论
分享
Top