无论安装以下哪种库,强烈建议在新环境下安装,之前在base环境下安装各种报错!!!如何在新环境下安装可参见我之前的总结或网上其他文章。
不推荐方法
1.pygrib
pygrib库基于python2.4+、numpy、pyproj、Jasper和GRIB_API,依赖库基于Linux或MacOS,GRIB_API将于2020年6月30日正式被ecCodes取代。因此,不推荐使用pygrib库。
2.GDAL
可于window下安装,安装简单,但是需借助wgrib2查看grib变量所在通道,因此不推荐。
推荐方法
xarray搭配cfgrib
优点:Linux windows MacOS均可用,安装较简单,强烈建议在新环境下安装各所需库
安装步骤:
1,新建新环境
conda creat -n cfgrib_env python=3.7.0
2.启动虚拟环境
source activate cfgrib_env
3.安装spyder xarray eccodes cfgrib及其他所需依赖库
ps:可用pip或conda安装,eccodes无需配置(有些地方写需要配置路径)
4.在新环境下启动spyder(base或其他环境下启动的spyder无效)
切换到安装路径并执行:~/anaconda3/envs/cfgrib_env/bin$ ./spyder
5.写入代码
import xarray as xr import cfgrib ds = xr.open_dataset('/1.grib2',engine='cfgrib')
提示错误:
DatasetBuildError: multiple values for unique key, try re-open the file with one of: filter_by_keys={'typeOfLevel': 'unknown'} filter_by_keys={'typeOfLevel': 'surface'} filter_by_keys={'typeOfLevel': 'isobaricInhPa'} filter_by_keys={'typeOfLevel': 'meanSea'} filter_by_keys={'typeOfLevel': 'depthBelowLandLayer'} filter_by_keys={'typeOfLevel': 'heightAboveGround'} filter_by_keys={'typeOfLevel': 'heightAboveGroundLayer'} filter_by_keys={'typeOfLevel': 'tropopause'} filter_by_keys={'typeOfLevel': 'maxWind'} filter_by_keys={'typeOfLevel': 'heightAboveSea'} filter_by_keys={'typeOfLevel': 'isothermZero'} filter_by_keys={'typeOfLevel': 'pressureFromGroundLayer'} filter_by_keys={'typeOfLevel': 'sigmaLayer'} filter_by_keys={'typeOfLevel': 'sigma'} filter_by_keys={'typeOfLevel': 'potentialVorticity'}
莫慌!!!
这是因为目前cfgrib库无法同时读取多个typeOfLevel,因此我们只需要根据提示筛选我们需要的数据就行了。原程序添加所需数据为:
ds = xr.open_dataset('/1.grib2',engine='cfgrib',backend_kwargs= {'filter_by_keys': { 'typeOfLevel': 'surface' } } )
所需变量为:ds.变量名
折腾了好几天,做其他事情都心不在焉,搞出来一身轻~~~
总结