（数据科学学习手札96）在geopandas中叠加在线地图

本文示例文件已上传至我的Github仓库https://github.com/CNFeffery/DataScienceStudyNotes

1 简介

　　国庆期间，抽空给大家分享在geopandas中叠加各种在线瓦片底图的方法，来制作出更多样式的地图作品。话不多说，我们直接进入正题。

图1

2 在geopandas中叠加在线地图

　　我们需要配合contextily这个第三方库来辅助geopandas叠加在线地图，在geopandas已经被正确安装的情况下，使用pip install contextily或conda install contextily安装contextily。

　　从下面这个小例子出发（使用“特殊”上网技巧的朋友请注意，需要将requests库降级到2.24.0才可以在使用“特殊”上网技巧的同时，正常使用叠加在线底图的功能）：

import geopandas as gpd
import contextily as ctx
import matplotlib.pyplot as plt

cq = gpd.read_file('重庆市.geojson').to_crs('EPSG:3857')

fig, ax = plt.subplots(figsize=(10, 10))
ax = cq.plot(ax=ax, alpha=0.1, edgecolor='k')


ax.axis('off')

ctx.add_basemap(ax, 
                source='https://{s}.tile.openstreetmap.fr/hot/{z}/{x}/{y}.png',
                zoom=8)

fig.savefig('图2 叠加在线地图示例.png', pad_inches=0, bbox_inches='tight', dpi=300)

图2 叠加在线地图示例

　　下面我们来划重点，在图2所示的例子中，我们前面正常读入矢量数据后一定要先变换投影为web墨卡托即EPSG:3857，接着正常绘图，在最后一步时将ax对象传入ctx.add_basemap中，并添加了参数source代表对应在线瓦片地图的url，参数zoom来控制地图缩放精度级别。

　　在稍事等待之后（如果没有“特殊”的上网技巧，可能要多等一段时间），我们的底图便自动获取拼接完毕，之后直接导出图像文件即可，是不是非常的方便~

　　在掌握了geopandas+contextily相互配合叠加在线底图之后，下面给大家推荐一些有意思的底图url供大家日常选择使用：