本文示例文件已上傳至我的Github 倉(cāng)庫(kù)https://github.com/CNFeffery/DataScienceStudyNotes ?
1 簡(jiǎn)介國(guó)慶期間,抽空給大家分享在geopandas 中疊加各種在線瓦片底圖的方法,來(lái)制作出更多樣式的地圖作品。話不多說(shuō),我們直接進(jìn)入正題。 圖12 在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供大家日常選擇使用: https://a.tile./mobile-atlas/{z}/{x}/{y}.png?apikey=41f4f936f1d148f69cbd100812875c88 圖3http://{s}.tile.stamen.com/terrain/{z}/{x}/{y}.jpg 圖4http://{s}.tiles.maps.sputnik.ru/{z}/{x}/{y}.png 圖5https://c.tile./cycle/{z}/{x}/{y}.png?apikey=41f4f936f1d148f69cbd100812875c88 圖6http://{s}.basemaps./light_all/{z}/{x}/{y}.png 圖7http://{s}.tile.stamen.com/toner/{z}/{x}/{y}.png 圖8https://{s}.tile.openstreetmap.fr/hot/{z}/{x}/{y}.png 圖9http://{s}.tile.stamen.com/watercolor/{z}/{x}/{y}.jpg 圖10https://{s}.tile.openstreetmap.fr/osmfr/{z}/{x}/{y}.png 圖11https://stamen-tiles-a.a.ssl./terrain-background/{z}/{x}/{y}.png
「即stamen地形圖的無(wú)地名標(biāo)注版本」 圖12https://d.basemaps./light_nolabels/{z}/{x}/{y}.png
「即carto淺色地圖的無(wú)地圖標(biāo)注版本」 圖13
以上就是本文的全部?jī)?nèi)容,歡迎在評(píng)論區(qū)與我進(jìn)行討論~ 加入我們的知識(shí)星球【Python大數(shù)據(jù)分析】 愛(ài)上數(shù)據(jù)分析!
|