（數據科學學習手札96）在geopandas中疊加在線地圖

本文示例文件已上傳至個人Github倉庫https://github.com/CNFeffery/DataScienceStudyNotes

1 簡介

　　國慶期間，抽空給你們分享在geopandas中疊加各類在線瓦片底圖的方法，來製做出更多樣式的地圖做品。話很少說，咱們直接進入正題。

圖1

2 在geopandas中疊加在線地圖

　　咱們須要配合contextily這個第三方庫來輔助geopandas疊加在線地圖，在geopandas已經被正確安裝的狀況下，使用pip install contextily或conda install contextily安裝contextily。

　　從下面這個小例子出發：

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供你們平常選擇使用：

• https://a.tile.thunderforest.com/mobile-atlas/{z}/{x}/{y}.png?apikey=41f4f936f1d148f69cbd100812875c88

圖3

• http://{s}.tile.stamen.com/terrain/{z}/{x}/{y}.jpg

圖4

• http://{s}.tiles.maps.sputnik.ru/{z}/{x}/{y}.png

圖5

• https://c.tile.thunderforest.com/cycle/{z}/{x}/{y}.png?apikey=41f4f936f1d148f69cbd100812875c88

圖6

• http://{s}.basemaps.cartocdn.com/light_all/{z}/{x}/{y}.png

圖7

• http://{s}.tile.stamen.com/toner/{z}/{x}/{y}.png

圖8

• https://{s}.tile.openstreetmap.fr/hot/{z}/{x}/{y}.png

圖9

• http://{s}.tile.stamen.com/watercolor/{z}/{x}/{y}.jpg

圖10

• https://{s}.tile.openstreetmap.fr/osmfr/{z}/{x}/{y}.png

圖11

• https://stamen-tiles-a.a.ssl.fastly.net/terrain-background/{z}/{x}/{y}.png

　　即stamen地形圖的無地名標註版本

圖12

• https://d.basemaps.cartocdn.com/light_nolabels/{z}/{x}/{y}.png

　　即carto淺色地圖的無地圖標註版本

圖13

---

　　以上就是本文的所有內容，歡迎在評論區與我進行討論~