點(diǎn)擊上方"藍(lán)字"關(guān)注我們

Python大數(shù)據(jù)分析

記錄? ?分享? ?成長(zhǎng)

添加微信號(hào)"CNFeffery"加入技術(shù)交流群

?
本文示例文件已上傳至我的Github倉(cāng)庫(kù)https://github.com/CNFeffery/DataScienceStudyNotes
?

1 簡(jiǎn)介

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

圖1

2 在geopandas中疊加在線地圖

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

從下面這個(gè)小例子出發(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 疊加在線地圖示例

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

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

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