python讀取與處理netcdf數據

netcdf是氣候數據中的主流格式，當涉及到大范圍的全球數萬個格網點數據時，使用python腳本可以較快地讀取與處理。

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

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

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

#查看nc數據各個變量的信息
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所在經緯度網格大致位於北京市大興區，查詢2012年年總降水量為523.8mm，大致吻合常規結果（百度百科大興區年平均降水量約為556mm）。
l1=[]
for i in range(23377,23742):　　　　　#2012年時間區間
    l1.append(data[i][0][100][592]) #100,592為經緯度標記 
print(np.sum(l1))