GrADS数据转为NetCDF格式数据

MeteoInfo

增加了将GrADS二进制格点数据转为NetCDF格式数据的功能，使用GrADSDataInfo类里添加的ConvertToNCData(string ncFilePath)方法，用GrADS的经典示例数据（model.ctl, model.dat）进行了测试，应该还不错，不过也仅做了这个测试，不能保证一定能输出正确结果，大家可以用自己的数据测试一下。

之前有人反映一些nc文件的时间解码得不正确（比如slp.1978.nc），也做了修改。此类数据的时间变量的单位为：units = "hours since 1-1-1 00:00:0.0"，不清楚为什么解出来的时间会多48小时（解码程序应该没有问题），为了符合实际数据时间，简单粗暴地给所有时间减去了48小时，以后有了更好的理解再说。

懒得把这个功能加在菜单里，可以写个很简单的脚本程序来完成（脚本程序也更灵活），示例如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-7-12                                                
   
4. # Purpose: Convert GrADS binaray grid data to NetCDF data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------     
   
7. import clr
   
8. clr.AddReference("MeteoInfoC.dll")
   
9. from MeteoInfoC import *
   
10. from MeteoInfoC.Data import *
    
11. from MeteoInfoC.Data.MeteoData import *
    
12. 
    
13. #---- Set data files
    
14. infile = 'F:\\Temp\\GrADS\\model.ctl'
    
15. outfile = 'F:\\Temp\\nc\\model.nc'
    
16. 
    
17. #---- Convert GrADS data to netCDF data
    
18. aDataInfo = GrADSDataInfo()
    
19. aDataInfo.ReadDataInfo(infile)
    
20. aDataInfo.ConvertToNCData(outfile)
    
21. print 'Finished...'

复制代码

mofangbao

楼主体谅一下比你更懒的“菜鸟”们啊，哈哈  不过既然脚本文件都出来了，运行一下也是很方便的

残月

楼主  菜鸟请教这这个怎么用啊？谢谢了

MeteoInfo

残月 发表于 2012-10-3 11:53

                               
登录/注册后可看大图

楼主  菜鸟请教这这个怎么用啊？谢谢了

MeteoInfo的Tools -> Script菜单打开MeteoInfo的脚本编辑器，将一楼帖子里的脚本代码复制、粘帖到脚本编辑器中，修改数据路径，点击脚本编辑器中Run Script按钮运行该脚本脚本即可。

残月

MeteoInfo 发表于 2012-10-3 12:46

                               
登录/注册后可看大图

MeteoInfo的Tools -> Script菜单打开MeteoInfo的脚本编辑器，将一楼帖子里的脚本代码复制、粘帖到脚本编辑 ...

非常感谢楼主~

但是出现了个错误：17层转换后只剩一层了，是不是还需要什么设置呢？谢谢楼主
下面是ctl的内容
dset D:\Linux\home\kaomiao\newdata\air_jan1tofeb10.dat
undef -9.96921e+36
title mm
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef 164 linear 00Z01jan1960 6hr
vars  1
aveair 17 99  *air      
ENDVARS

MeteoInfo

残月 发表于 2012-10-3 22:21

                               
登录/注册后可看大图

非常感谢楼主~

但是出现了个错误：17层转换后只剩一层了，是不是还需要什么设置呢？谢谢楼主

可能是bug，这个功能没有经过多少测试。你把GrADS格式数据发上来我测试一下。

残月

MeteoInfo 发表于 2012-10-3 22:56

                               
登录/注册后可看大图

可能是bug，这个功能没有经过多少测试。你把GrADS格式数据发上来我测试一下。

数据和CTL文件

MeteoInfo

残月 发表于 2012-10-4 09:40

                               
登录/注册后可看大图

数据和CTL文件

刚试了试，没有问题呀，多层数据都在。

残月

MeteoInfo 发表于 2012-10-4 10:20

                               
登录/注册后可看大图

刚试了试，没有问题呀，多层数据都在。

用Meteonfo打开是正常的，用grads打开是不对的

MeteoInfo

残月 发表于 2012-10-4 10:34

                               
登录/注册后可看大图

用Meteonfo打开是正常的，用grads打开是不对的

grads中的sdfopen对NC文件的识别能力很有限，我也无可奈何。