python讀取nc數(shù)據(jù)并繪圖的方法實(shí)例
獲取nc數(shù)據(jù)的相關(guān)信息
from netCDF4 import Dataset import numpy as np import pandas as pd import os import matplotlib.pyplot as plt path = "F:\\OCO2.SIF.all.daily.2001.nc" csv_path = "F:\\test.csv" dst = Dataset(path, mode='r', format="netCDF4") print(dst.variables.keys()) data = dst.variables['all_daily_sif'][:] print(data.shape) # 輸出結(jié)果如下: # dict_keys(['lat', 'lon', 'doy', 'all_daily_sif']) # (92, 360, 720) #可見有92個(gè)時(shí)間序列,經(jīng)度(lon)、緯度(lat)的取值有720,360個(gè) # # 查看數(shù)據(jù)經(jīng)緯度范圍,經(jīng)度-179.75~179.75,其中負(fù)值為西經(jīng),正值為東經(jīng);緯度正為北緯,負(fù)為南緯 # # 格點(diǎn)分辨率為0.5度 long = dst.variables['lon'][:] lati = dst.variables['lat'][:] print(long[0], long[-1], lati[0], lati[-1]) print(long.shape, lati.shape)
繪圖
用matplotlib繪圖
參考文獻(xiàn)1
# plt對某個(gè)doy的全球sif值作圖。左半部分為西半球,右邊是東半球 # 選了doy為10的sif數(shù)據(jù)作圖 plt.contourf(long, lati, data[10, :, :] ) plt.colorbar(label="Sif", orientation="horizontal") plt.show()
運(yùn)行結(jié)果:
用Basemap繪圖
參考文獻(xiàn)2
from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as plt lat = dst.variables['lat'][:] lon = dst.variables['lon'][:] data = dst.variables['all_daily_sif'][:] data[10] = data[10] # use .shape function to check that arrays have # the correct size. # e.g. lon.shape print(data[10].shape) lon0 = lon.mean() lat0 = lat.mean() # 設(shè)置投影方式:cyl為圓柱投影、還可設(shè)置merc為mercator投影 llcrnrlat為起始lat;urcrnrlat為終止lat # m = Basemap(projection='merc', llcrnrlat=lat[0], urcrnrlat=lat[-1], \ # llcrnrlon=lon[0], urcrnrlon=lon[-1], ax=ax1) # 參數(shù) "resolution" 用于控制地圖面積邊緣的精細(xì)程度,有'l'和'h'兩種取值 m = Basemap(lat_0=lat0, lon_0=lon0,projection='cyl',resolution='l') # 繪制等經(jīng)緯度線 緯度每隔20度畫一條線,且標(biāo)注經(jīng)緯度 m.drawparallels(np.arange(-90., 91., 20.), labels=[1, 0, 0, 0], fontsize=10) m.drawmeridians(np.arange(-180., 181., 40.), labels=[0, 0, 0, 1], fontsize=10) m.drawcoastlines()# 繪制海岸線 # m.drawcountries(linewidth=0.25) # 繪制國界線 # m.readshapefile('F:\E\data\grass_yield\shp\quhua\\省', 'states') # 讀取中國各省邊界,并繪圖 lon, lat = np.meshgrid(lon, lat) xi, yi = m(lon, lat) # cmap是顏色,還可選‘jet'、‘spring'、‘winter'、'summer'、'autumn' cs = m.contourf(xi, yi, data[10], cmap='summer') # pad指位置, cbar = m.colorbar(cs, location='bottom', pad="10%",format='%.1f') # cbar = m.colorbar(C, 'right', ticks=np.arange(-128, 128, 40), format='%.1f') font1 = {'family': 'DejaVu Sans', 'weight': 'normal', 'size': 16} plt.title('CSIF', font1) plt.show()
運(yùn)行效果:
用Cartopy繪圖
參考文獻(xiàn)3
此前 Python 最常用的地圖包是 Basemap,然而它將于 2020 年被棄用,官方推薦使用 Cartopy 包作為替代。Cartopy 是英國氣象局開發(fā)的地圖繪圖包,實(shí)現(xiàn)了 Basemap 的大部分功能,還可以通過 Matplotlib 的 API 實(shí)現(xiàn)豐富的自定義效果。
安裝Cartopy包
下載安裝OSGeo4W4
參考文獻(xiàn):https://zhuanlan.zhihu.com/p/129351199
參考文獻(xiàn):https://blog.csdn.net/weixin_39618339的plt畫圖像圖例的位置怎么寫代碼_用basemap畫氣象圖
參考文獻(xiàn): https://zhajiman.github.io/
https://trac.osgeo.org/osgeo4w/
總結(jié)
到此這篇關(guān)于python讀取nc數(shù)據(jù)并繪圖的文章就介紹到這了,更多相關(guān)python讀取nc數(shù)據(jù)繪圖內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
教你用Python代碼實(shí)現(xiàn)合并excel文件
近幾天一直因?yàn)閑xcel文件太多太雜的原因苦惱,今天特地整理了本篇文章,文章介紹的非常詳細(xì),對正在學(xué)習(xí)python的小伙伴們有很好地幫助,需要的朋友可以參考下2021-05-05用Python創(chuàng)建簡易網(wǎng)站圖文教程
今天給大家?guī)淼氖顷P(guān)于Python的相關(guān)知識,文章圍繞著用Python創(chuàng)建簡易網(wǎng)站展開,文中有非常詳細(xì)的介紹及圖文示例,需要的朋友可以參考下2021-06-06Python中使用json.load()和json.loads()加載json數(shù)據(jù)的方法實(shí)例
在python編程中,我們經(jīng)常要用到j(luò)son對象作為數(shù)據(jù)交換格式,下面這篇文章主要給大家介紹了關(guān)于Python中使用json.load()和json.loads()加載json數(shù)據(jù)的方法實(shí)例,文中通過實(shí)例代碼介紹的非常詳細(xì),需要的朋友可以參考下2022-08-08django+xadmin+djcelery實(shí)現(xiàn)后臺管理定時(shí)任務(wù)
這篇文章主要介紹了django+xadmin+djcelery實(shí)現(xiàn)后臺管理定時(shí)任務(wù),小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,也給大家做個(gè)參考。一起跟隨小編過來看看吧2018-08-08PyTorch一小時(shí)掌握之神經(jīng)網(wǎng)絡(luò)氣溫預(yù)測篇
這篇文章主要介紹了PyTorch一小時(shí)掌握之神經(jīng)網(wǎng)絡(luò)氣溫預(yù)測篇,本文給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2021-09-09