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

首先需要确保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文件！