本文的文字及图片来源于网络,仅供学习 、交流使用,不具有任何商业用途,如有问题请及时联系我们以作处理 。

以下文章来源于python数据分析之禅  ,作者小dull鸟

 

前言

今天给大家分享一篇非常实用的文章 ,用folium制作北京市交通电子眼分布地图,再也不怕被隐藏摄像头偷拍了

成果图如下:

 

主要分为以下2个步骤:

1.电子眼地理数据获取

网址如下:

https://www.icauto.com.cn/weizhang/wzd/110000/

 

 

此网站可以获取全国各地区的电子眼数据

我们以北京市为例,用python爬取北京市所有电子眼名称 ,代码如下:

url1='https://www.icauto.com.cn/weizhang/wzd/110000/list_1.html'
response=requests.get(url1,headers=header)
soup=BeautifulSoup(response.text,'html.parser')
results=soup.find('div',class_='cdz-ccotent').find_all('li')
for result in results:
    time.sleep(0.5)
    name=result.find('a').text
    index=result.find('span').text.split('')[1]

 

经过抓包分析,可以得到根据电子眼名称获取经纬度信息的接口:

https://api.map.baidu.com/?qt=gc&wd=北京西路中坝隧道路段&cn=北京&ie=utf-8&oue=1&fromproduct=jsapi&res=api&ak=s8sS5dBsZ7bLRi3bcVRAaYMAnqlXoyeo

 

返回结果为:

 

大家主要看coord参数,刚一看确实有点懵逼 ,不知道这是啥玩意,我查了半天,才知道这种是百度地图采用的墨卡托平面坐标 ,利用百度地图api可以转化成经纬度,转化代码如下:

url3='http://api.map.baidu.com/geoconv/v1/?coords={}&from=6&to=5&ak=换成你的ak'.format(str(coord['x'])+','+str(coord['y']))
response=requests.get(url3,headers=header)
result=json.loads(response.text)['result'][0]
lon=result['x']
lat=result['y']

 

最后将结果保存到csv表格中:

with open('dianziyan.csv', 'a+', newline='', encoding='gb18030') as f:
    f_csv = csv.writer(f)
    f_csv.writerow([name, index, lon,lat])

 

2.用folium标点

folium地图标点的方法之前已经有介绍了,可以参考:

python生成广州全市停车场分布地图

代码如下:

import pandas as pd
data=pd.read_csv('dianziyan.csv',encoding='gbk')
import folium
from folium import plugins
Camera_map = folium.Map(location=[data['纬度'].mean(), data['经度'].mean()], zoom_start=10,control_scale=True,)
incidents = folium.map.FeatureGroup()
for name,row in data.iterrows():
    incidents.add_child(
        folium.CircleMarker(            #CircleMarker表示花圆
            [row["纬度"], row["经度"]],   #每个停车场的坐标
            radius=7,                   #圆圈半径
            color='yellow',             #标志的外圈颜色
            fill=True,                  #是否填充
            fill_color='red',           #填充颜色
            fill_opacity=0.4,            #填充透明度
            
        )
    )

Camera_map.add_child(incidents)
Camera_map.save('Camera_map1.html')

 

 

经过放大后可以发现 ,坐标定位不是很准,有的严重偏离路线,有的甚至定位到湖里了 ,这是因为底图坐标和电子眼经纬度标准不统一

为了定位更准确 ,我们把底图和电子眼经纬度统一规范为高德地图

百度经纬度转换成高德经纬度的函数如下:

import math
def bdToGaoDe(lat,lon):
    """
    百度经纬度转高德经纬度
    :param lon:
    :param lat:
    :return:
    """
    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)
    return lat,lon

 

读取数据,进行转换:

import pandas as pd
data=pd.read_csv('dianziyan.csv',encoding='gbk')
for name,row in data.iterrows():
    print(bdToGaoDe(row["纬度"],row["经度"]))

 

将底图换成高德地图,然后将标点转换成电子眼图标:

import folium
from folium import plugins
Camera_map = folium.Map(location=[data['纬度'].mean(), data['经度'].mean()], zoom_start=10,zoom_control='False',
                     tiles='http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=7&x={x}&y={y}&z={z}',attr='AutoNavi')
incidents = folium.map.FeatureGroup()
tooltip ='请点击我查看该点信息'
for name,row in data.iterrows():
    incidents.add_child(
        folium.Marker(            #CircleMarker表示花圆
            [bdToGaoDe(row["纬度"],row["经度"])[0],bdToGaoDe(row["纬度"],row["经度"])[1]],   #每个停车场的坐标
            icon=folium.Icon(color='green', prefix='fa', icon='bullseye')
        )
    )

Camera_map.add_child(incidents)
Camera_map.save('Camera_map2.html')

 

 

 

通过tiles可以设置不同的地图瓦片在 ,这里设置为高德地图瓦片

Folium.Icon类可以设置color, icon_color, icon, angle, prefix这5个参数:

color的可选项包括:[‘red’, ‘blue’, ‘green’, ‘purple’, ‘orange’, ‘darkred’, ‘lightred’, ‘beige’, ‘darkblue’, ‘darkgreen’, ‘cadetblue’, ‘darkpurple’, ‘white’, ‘pink’, ‘lightblue’, ‘lightgreen’, ‘gray’, ‘black’, ‘lightgray’] ,或者HTML颜色代码
icon_color同上
icon可以在Font-Awesome网站中找到对应的名字,并设置prefix参数为’fa’
angle以度为单位设置

 

本文只显示了北京地区的电子眼位置 ,感兴趣的朋友可以深入研究一下,做一张全国的交通电子眼分布图 。

文章来源于网络,如有侵权请联系站长QQ61910465删除
本文版权归趣营销www.SEOgUrublog.com 所有,如有转发请注明来出,竞价开户托管,seo优化请联系QQ卍61910465