python中怎么插一個(gè)基礎(chǔ)的地圖分析類庫

本篇內(nèi)容主要講解“python中怎么插一個(gè)基礎(chǔ)的地圖分析類庫”，感興趣的朋友不妨來看看。本文介紹的方法操作簡(jiǎn)單快捷，實(shí)用性強(qiáng)。下面就讓小編來帶大家學(xué)習(xí)“python中怎么插一個(gè)基礎(chǔ)的地圖分析類庫”吧!

python中最常用的類庫可以使用mpl_toolkits.basemap中的Basemap，數(shù)據(jù)是北京的空氣質(zhì)量歷史數(shù)據(jù)，要使用現(xiàn)有的矢量數(shù)據(jù)，在地圖上畫一個(gè)北京政區(qū)圖。

矢量數(shù)據(jù)用shp或者geojson格式都可以，為了減少數(shù)據(jù)文件數(shù)量，轉(zhuǎn)了一個(gè)geojson數(shù)據(jù)，通過fiona進(jìn)行數(shù)據(jù)讀取，引用fiona的類庫，讀取geojson數(shù)據(jù)，可以參見python開發(fā)：特征數(shù)據(jù)處理。

import fiona
shpdata = fiona.open(path, 'r')

引用Basemap的類庫，循環(huán)遍歷geojson的數(shù)據(jù)，在地圖上進(jìn)行數(shù)據(jù)的繪制，數(shù)據(jù)根據(jù)其中的情況，獲得坐標(biāo)信息，使用plot函數(shù)進(jìn)行地圖的繪制線操作，設(shè)置線的顏色和寬度等屬性，使用的數(shù)據(jù)是polygon類型，要把其中的每條線都遍歷畫出來。

import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap as Basemap

fig, ax = plt.subplots()
m = Basemap(epsg=4326, resolution='l', llcrnrlat=bounds[1], urcrnrlat=bounds[3],
            llcrnrlon=bounds[0], urcrnrlon=bounds[2])

while True:
    try:
        evdata = shpdata.next()
        evshape = evdata['geometry']['coordinates']

for evls in evshape:
    xx, yy = zip(*evls)
    m.plot(xx, yy, linewidth=1, color='black')

有時(shí)需要在地圖上畫注記信息，就需要使用shapely計(jì)算一下區(qū)域的中心點(diǎn)，在地圖上畫注記。

from shapely.geometry import shape, point

ax.annotate(evdata['properties']['NAME99'], 
(X, Y), xytext=(-10, 5), textcoords='offset points')

效果圖：

到此，相信大家對(duì)“python中怎么插一個(gè)基礎(chǔ)的地圖分析類庫”有了更深的了解，不妨來實(shí)際操作一番吧！這里是億速云網(wǎng)站，更多相關(guān)內(nèi)容可以進(jìn)入相關(guān)頻道進(jìn)行查詢，關(guān)注我們，繼續(xù)學(xué)習(xí)！