为了我的项目,我想从一个netCDF文件中提取数据,该文件表示海洋表面温度与纬度/经度的关系。我想选择一个具有精确坐标的点,然后对每个月份的多年平均值进行计算。
import netCDF4 as nc
import numpy as np
import cdo
import os
import pandas as pd
from scipy.io import netcdf
from netCDF4 import Dataset
from netCDF4 import num2date
import matplotlib as mpl
import matplotlib.pylab as plt
from cdo import Cdo
cdo = Cdo()
data = "/Users/bf/Desktop/Stage
tos_Omon_EC-Earth3_historical_r5i1p1f1_gn_185001-201412.nc"
ds = nc.Dataset(data)
infile = "/Users/bf/Desktop/Stage/tos_Omon_EC-Earth3_historical_r5i1p1f1_gn_185001-201412.nc"
#to select years : 1979-2014
#new name for the file with the selected years : tos_1979-2014.nc
cdo.seldate('1979-01-01,2014-12-31', input=infile, output="tos_1979-2014.nc")
#to select longitute and latitute -> no important now
#new name for the file with the selected longitudes/latitudes : tos_1979-2014_2.nc
cdo.sellonlatbox('-26.45,65.03', input="tos_1979-2014.nc", output="tos_1979-2014_2.nc")
#ymonavg : multi-year monthly statistical values average
cdo.ymonavg(input="tos_1979-2014_2.nc", output="tos_1979-2014_mean_2.nc")
我想提取我的温度数据(自1979年至2014年每个月的平均温度),并将它们放入一个数据框中,以后用mathplolib进行图形处理。
但是使用CDO
工具时,我无法选择一个点。也许我需要使用另一个工具,但我不知道该怎么做。
谢谢您的帮助!