python处理FNL数据的grib文件和nc文件(纬度存储的问题)
python处理FNL数据的grib文件和nc文件(纬度存储的问题)
在使用python处理FNL数据时,2007年及之前的数据存储为grib文件,2007年以后为grib2文件,在数据处理时,将grib2转化为nc文件处理Windows上python读取grib2文件(不用Linux),对grib文件则安装了Anaconda环境处理Windows下xarray+cfgrib读取grib文件。但在根据经纬度读取变量的时候,发现不同文件的存储不同,而变量是根据经纬度存储的下标去获取值,所以需要注意一下不同的文件的存储内容。
1.python获取FNL的nc文件对应的变量
nc_obj =Dataset('F:/Zhu/download/2000-2015/NC/NC2008/fnl_20080113_06_00.nc')
nc_obj.variables[name][0][lat][lon] 根据变量名,经纬度提取变量值
从FNL的nc文件中获取变量的方式如上,需要根据变量名,第一个是time,这里是一维存储,下标只有0即可。然后是纬度,经度获取变量的值。
先看一下,FNL的nc文件中纬度的存储:
print(nc_obj.variables['latitude'][:])
从结果可以看到