前言
以前用 grads 和 Fortran 时候,经常会用到 grd 二进制文件,并拿 ctl 文件画图
可以用 Python 的 xgrads 包来实现这些操作
下载 xgrads
pip install xgrads
官方网站:https://xgrads.readthedocs.io/en/latest/
注意事项
在下载 xgrads 之后,就可以使用了,使用之前,请确保你的 ctl 文件和 grd 文件正确。
示例
-
打开多个空间场相同,时间不同的 ctl 文件时,可以批量读取
类似于 xarray 的 open_mfdataset
from xgrads import open_mfdataset
ds=open_mfdataset('/../*.ctl')
ds -
转 nc 格式
from xgrads import open_mfdataset
ds=open_dataset('/../xxxx.ctl')
ds.attrs['pdef']='None' # 原来的 ds 的 attr 是没有 pdef 这个属性的,所以如果报错,可以加上这个,大致报错的原因都是因为属性缺少有些东西,缺啥补啥
ds.to_netcdf('/../xx.nc')
更多的使用方法,具体参考官网