溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務(wù)條款》

利用python繪制中國地圖(含省界、河流等)

發(fā)布時間:2020-10-07 09:40:31 來源:腳本之家 閱讀:338 作者:王琛NUIST 欄目:開發(fā)技術(shù)

我們可以使用Basemap這個工具包來實現(xiàn)中國地圖的繪制

首先需要加載一些包:

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

Basemap包就是氣象畫圖的利器,現(xiàn)在我們就可以愉快的畫圖了!

plt.figure(1)
map=Basemap()
map.drawcoastlines()
plt.title(r'$World\ Map$',fontsize=24)
plt.show()

第2行創(chuàng)建一個地圖,第3行添加海岸線,這樣一個世界地圖就出來了,怎么樣,很簡單吧。(plt.show()這行代碼是用來顯示圖片的)

利用python繪制中國地圖(含省界、河流等)

我們發(fā)現(xiàn)這只是海岸線圖,那么怎么將國界線添加上去呢?很簡單,只要添加一行代碼就可以了。

map.drawcountries()

利用python繪制中國地圖(含省界、河流等)

那么怎么添加河流呢?可能有些同學(xué)已經(jīng)猜到了,就是drawrivers()

map.drawrivers(color='blue',linewidth=0.3)

利用python繪制中國地圖(含省界、河流等)

好了,現(xiàn)在我們可以開始畫中國地圖了!
其實只要在創(chuàng)建地圖時指定一下范圍就可以了,查閱資料發(fā)現(xiàn),中國的經(jīng)緯度范圍是東經(jīng)135度2分30秒-東經(jīng)73度40分,北緯3度52分-北緯53度33分。

map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54)

利用python繪制中國地圖(含省界、河流等)

好了,一個中國地圖就出來了!但是我們發(fā)現(xiàn),好像少了點什么,沒錯就是省界。我們可以在https://gadm.org/download_country_v3.html下載中國大陸和臺灣省的行政區(qū)域的shape文件,下載后解壓,然后加入下面的代碼。

CHN='G:\python_material\MapOfChina'

CHN的值就是解壓后的地圖文件所在的地址。
下面我們就可以加入省界了!

map.readshapefile(CHN+'\gadm36_CHN_shp\gadm36_CHN_1',
         'states',drawbounds=True)

別忘了把臺灣省加上去

map.readshapefile(CHN+'\gadm36_TWN_shp\gadm36_TWN_1',
         'taiwan',drawbounds=True)

利用python繪制中國地圖(含省界、河流等)

還可以在地圖上加上經(jīng)緯度,比如我們要畫5條經(jīng)緯線,可以這么做:

parallels = np.linspace(3,55,5)
map.drawparallels(parallels,labels=[True,False,False,False])
meridians = np.linspace(70,140,5)
map.drawmeridians(meridians,labels=[False,False,False,True])

利用python繪制中國地圖(含省界、河流等)

大功告成!但是,emmm,我們發(fā)現(xiàn)好像有點歪?我們可以在創(chuàng)建地圖時選擇投影參數(shù)。

map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=137,urcrnrlat=54,
      projection = 'lcc', lat_1 = 33, lat_2 = 45, lon_0 = 100)

利用python繪制中國地圖(含省界、河流等)

這回正式的完成了!

附上所有代碼:

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

plt.figure(1)
map=Basemap(llcrnrlon=70,llcrnrlat=3,urcrnrlon=139,urcrnrlat=54,
      projection = 'lcc', lat_1 = 33, lat_2 = 45, lon_0 = 100)
map.drawcoastlines()
map.drawcountries()
map.drawrivers(color='blue',linewidth=0.3)
CHN='G:\python_material\MapOfChina'
map.readshapefile(CHN+'\gadm36_CHN_shp\gadm36_CHN_1',
         'states',drawbounds=True)
map.readshapefile(CHN+'\gadm36_TWN_shp\gadm36_TWN_1',
         'taiwan',drawbounds=True)
parallels = np.linspace(3,55,5)
map.drawparallels(parallels,labels=[True,False,False,False])
meridians = np.linspace(70,140,5)
map.drawmeridians(meridians,labels=[False,False,False,True])
plt.title(r'$China\ Map$',fontsize=24)

plt.show()

以上就是利用python繪制中國地圖的詳細(xì)內(nèi)容,更多關(guān)于python 繪制地圖的資料請關(guān)注億速云其它相關(guān)文章!

向AI問一下細(xì)節(jié)

免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。

AI