1 簡介

國慶期間，抽空給大家分享在geopandas中疊加各種在線瓦片底圖的方法，來制作出更多樣式的地圖作品。話不多說，我們直接進(jìn)入正題。

圖1

2 在geopandas中疊加在線地圖

我們需要配合contextily這個第三方庫來輔助geopandas疊加在線地圖，在geopandas已經(jīng)被正確安裝的情況下，使用pip install contextily或conda install contextily安裝contextily。

從下面這個小例子出發(fā)：

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 疊加在線地圖示例

下面我們來「劃重點(diǎn)」，在圖2所示的例子中，我們前面正常讀入矢量數(shù)據(jù)后「一定要先變換投影為web墨卡托即EPSG:3857」，接著正常繪圖，在最后一步時將ax對象傳入ctx.add_basemap中，并添加了參數(shù)source代表對應(yīng)在線瓦片地圖的url，參數(shù)zoom來控制地圖縮放精度級別。

在稍事等待之后（如果沒有“特殊”的上網(wǎng)技巧，可能要多等一段時間），我們的底圖便自動獲取拼接完畢，之后直接導(dǎo)出圖像文件即可，是不是非常的方便~

在掌握了geopandas+contextily相互配合疊加在線底圖之后，下面給大家推薦一些有意思的底圖url供大家日常選擇使用：