溫馨提示×

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

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

python如何讀取與處理netcdf數(shù)據(jù)方式

發(fā)布時(shí)間:2021-05-11 10:09:48 來(lái)源:億速云 閱讀:333 作者:小新 欄目:開(kāi)發(fā)技術(shù)

小編給大家分享一下python如何讀取與處理netcdf數(shù)據(jù)方式,希望大家閱讀完這篇文章之后都有所收獲,下面讓我們一起去探討吧!

python是什么意思

Python是一種跨平臺(tái)的、具有解釋性、編譯性、互動(dòng)性和面向?qū)ο蟮哪_本語(yǔ)言,其最初的設(shè)計(jì)是用于編寫(xiě)自動(dòng)化腳本,隨著版本的不斷更新和新功能的添加,常用于用于開(kāi)發(fā)獨(dú)立的項(xiàng)目和大型項(xiàng)目。

netcdf是氣候數(shù)據(jù)中的主流格式,當(dāng)涉及到大范圍的全球數(shù)萬(wàn)個(gè)格網(wǎng)點(diǎn)數(shù)據(jù)時(shí),使用python腳本可以較快地讀取與處理。

import netCDF4
from netCDF4 import Dataset
import numpy as np
import sys
import os

#計(jì)算日期數(shù)
import datetime 
d1=datetime.date(1900,1,1) 
d3 = d1 + datetime.timedelta(days =100) 
print (d3)

#查看nc數(shù)據(jù)基本信息
nc_obj=Dataset('precip.nc')
print(nc_obj)

#查看nc數(shù)據(jù)各個(gè)變量的信息
print(nc_obj.variables.keys())
for i in nc_obj.variables.keys():
  print('___________________________________________')
  print(i)
  print(nc_obj.variables[i])

precip=(nc_obj.variables['precip'][:])
lat=(nc_obj.variables['lat'][:])
lon=(nc_obj.variables['lon'][:])

舉例子

#39.76,116.25所在經(jīng)緯度網(wǎng)格大致位于北京市大興區(qū),查詢2012年年總降水量為523.8mm,大致吻合常規(guī)結(jié)果(百度百科大興區(qū)年平均降水量約為556mm)。
l1=[]
for i in range(23377,23742):     #2012年時(shí)間區(qū)間
  l1.append(data[i][0][100][592]) #100,592為經(jīng)緯度標(biāo)記 
print(np.sum(l1))

看完了這篇文章,相信你對(duì)“python如何讀取與處理netcdf數(shù)據(jù)方式”有了一定的了解,如果想了解更多相關(guān)知識(shí),歡迎關(guān)注億速云行業(yè)資訊頻道,感謝各位的閱讀!

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

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

AI