Python如何将GrADs常用文件转换为NetCDF格式？

之前有写过文章使用Ruby和NCL读取转换grd文件，现在有国人开发的GrADs的Python接口xgrads可用于文件格式转换。(点击可跳转！)

首先需要确保xgrads库的安装： pip install xgrads Install from github 或者 git clone https://github.com/miniufo/xgrads.git cd xgrads python setup.py install 链接https://github.com/miniufo/xgrads ， 有提供示例ctl和dat文件，下面我们是使用的ctl和grd文件转换的，方法类似： #import sys #sys.path.append('/home/gavin/miniconda3/envs/atmpy/lib/python3.8/site-packages') #sys.path from xgrads import CtlDescriptor, open_CtlDataset ds = open_CtlDataset('lst.ctl') ctl = CtlDescriptor(file='lst.ctl') ds.attrs['pdef' ] = 'None' ds.to_netcdf('lst.nc') data = ds.ro1 data.where(data!=ctl.undef).plot(figsize=(9,5), cmap='jet') 以上需要注意两点： 1.如果在jupyter-lab中无法加载xgrads需要手动添加其路径，使用到的是：import sys 2. xgrads存在bug，如果不添加语句ds.attrs['pdef' ] = 'None'会一直报错，无法生成nc文件！