netCDF数据转为GrADS数据

MeteoInfo

写了一个将netCDF数据转为GrADS数据的脚本程序，需要MeteoInfo最新文件的支持（见置顶帖子）。GrADS的二进制格式数据比较简单，容易用Fortran等程序读取，而netCDF数据格式比较复杂，用程序读取有些难度。这个脚本程序需要定义一些GrADSDataInfo的参数，然后就按照这些参数自动生成control文件。

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-5-28                                                
   
4. # Purpose: Convert NetCDF data to GrADS binary data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------     
   
7. import clr
   
8. from System import *
   
9. clr.AddReference("MeteoInfoC.dll")
   
10. from MeteoInfoC import *
    
11. from MeteoInfoC.Data import *
    
12. from MeteoInfoC.Data.MeteoData import *
    
13. 
    
14. #---- Set data folders
    
15. dataDir = 'D:\\Temp\\nc\\'
    
16. 
    
17. #---- Set input/output data files
    
18. inFile = dataDir + 'HadISST_SST_update.nc'
    
19. outCtlFile = dataDir + 'HadISST_SST_update.ctl'
    
20. outDataFile = dataDir + 'HadISST_SST_update.dat'
    
21. 
    
22. #---- Define data info
    
23. ncDataInfo = MeteoDataInfo()
    
24. gDataInfo = GrADSDataInfo()
    
25. 
    
26. #---- Open NetCDF data file
    
27. ncDataInfo.OpenNCData(inFile)
    
28. xArray = ncDataInfo.GetX()
    
29. yArray = ncDataInfo.GetY()
    
30. tNum = ncDataInfo.GetTimeNumber()
    
31. sTime = ncDataInfo.GetTime(0)
    
32. 
    
33. #---- Set output GrADS data info
    
34. gDataInfo.FileName = outCtlFile
    
35. gDataInfo.DSET = outDataFile
    
36. gDataInfo.TITLE = 'HadISST data'
    
37. gDataInfo.DTYPE = 'GRIDDED'
    
38. gDataInfo.XDEF.Type = 'LINEAR'
    
39. gDataInfo.XDEF.XNum = xArray.Length
    
40. gDataInfo.XDEF.XMin = xArray[0]
    
41. gDataInfo.XDEF.XDelt = xArray[1] - xArray[0]
    
42. gDataInfo.YDEF.Type = 'LINEAR'
    
43. gDataInfo.YDEF.YNum = yArray.Length
    
44. gDataInfo.YDEF.YMin = yArray[0]
    
45. gDataInfo.YDEF.YDelt = yArray[1] - yArray[0]
    
46. gDataInfo.ZDEF.Type = 'LINEAR'
    
47. gDataInfo.ZDEF.ZNum = 1
    
48. gDataInfo.ZDEF.SLevel = 1
    
49. gDataInfo.ZDEF.ZDelt = 1
    
50. gDataInfo.TDEF.Type = 'LINEAR'
    
51. gDataInfo.TDEF.TNum = tNum
    
52. gDataInfo.TDEF.STime = sTime
    
53. gDataInfo.TDEF.TDelt = '1mo'
    
54. aVar = Variable()
    
55. aVar.Name = 'sst'
    
56. aVar.Units = '99'
    
57. aVar.Description = 'Sea Surface Temporature'
    
58. gDataInfo.VARDEF.AddVar(aVar)
    
59. gDataInfo.MissingValue = ncDataInfo.MissingValue
    
60. 
    
61. #---- Write GrADS data file
    
62. gDataInfo.CreateDataFile(outDataFile)
    
63. for t in range(0, tNum):
    
64.         ncDataInfo.TimeIndex = t
    
65.         print t
    
66.         sst = ncDataInfo.GetGridData('sst')
    
67.         if t == 0:
    
68.                 gDataInfo.MissingValue = sst.MissingValue
    
69.                 print gDataInfo.MissingValue
    
70.         gDataInfo.WriteGridData(sst)
    
71. 
    
72. gDataInfo.CloseDataFile()
    
73. 
    
74. #---- Write GrADS control file
    
75. gDataInfo.WriteGrADSCTLFile()
    
76. 
    
77. print 'Finished...'

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2012-5-28                                                 <br /> # Purpose: Convert NetCDF data to GrADS binary data  <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------     <br /> import clr<br /> from System import *<br /> clr.AddReference("MeteoInfoC.dll")<br /> from MeteoInfoC import *<br /> from MeteoInfoC.Data import *<br /> from MeteoInfoC.Data.MeteoData import *<br /> <br /> #---- Set data folders<br /> dataDir = 'D:\\Temp\\nc\\'<br /> <br /> #---- Set input/output data files<br /> inFile = dataDir + 'HadISST_SST_update.nc'<br /> outCtlFile = dataDir + 'HadISST_SST_update.ctl'<br /> outDataFile = dataDir + 'HadISST_SST_update.dat'<br /> <br /> #---- Define data info<br /> ncDataInfo = MeteoDataInfo()<br /> gDataInfo = GrADSDataInfo()<br /> <br /> #---- Open NetCDF data file<br /> ncDataInfo.OpenNCData(inFile)<br /> xArray = ncDataInfo.GetX()<br /> yArray = ncDataInfo.GetY()<br /> tNum = ncDataInfo.GetTimeNumber()<br /> sTime = ncDataInfo.GetTime(0)<br /> <br /> #---- Set output GrADS data info<br /> gDataInfo.FileName = outCtlFile<br /> gDataInfo.DSET = outDataFile<br /> gDataInfo.TITLE = 'HadISST data'<br /> gDataInfo.DTYPE = 'GRIDDED'<br /> gDataInfo.XDEF.Type = 'LINEAR'<br /> gDataInfo.XDEF.XNum = xArray.Length<br /> gDataInfo.XDEF.XMin = xArray[0]<br /> gDataInfo.XDEF.XDelt = xArray[1] - xArray[0]<br /> gDataInfo.YDEF.Type = 'LINEAR'<br /> gDataInfo.YDEF.YNum = yArray.Length<br /> gDataInfo.YDEF.YMin = yArray[0]<br /> gDataInfo.YDEF.YDelt = yArray[1] - yArray[0]<br /> gDataInfo.ZDEF.Type = 'LINEAR'<br /> gDataInfo.ZDEF.ZNum = 1<br /> gDataInfo.ZDEF.SLevel = 1<br /> gDataInfo.ZDEF.ZDelt = 1<br /> gDataInfo.TDEF.Type = 'LINEAR'<br /> gDataInfo.TDEF.TNum = tNum<br /> gDataInfo.TDEF.STime = sTime<br /> gDataInfo.TDEF.TDelt = '1mo'<br /> aVar = Variable()<br /> aVar.Name = 'sst'<br /> aVar.Units = '99'<br /> aVar.Description = 'Sea Surface Temporature'<br /> gDataInfo.VARDEF.AddVar(aVar)<br /> gDataInfo.MissingValue = ncDataInfo.MissingValue<br /> <br /> #---- Write GrADS data file<br /> gDataInfo.CreateDataFile(outDataFile)<br /> for t in range(0, tNum):<br />         ncDataInfo.TimeIndex = t<br />         print t<br />         sst = ncDataInfo.GetGridData('sst')<br />         if t == 0:<br />                 gDataInfo.MissingValue = sst.MissingValue<br />                 print gDataInfo.MissingValue<br />         gDataInfo.WriteGridData(sst)<br /> <br /> gDataInfo.CloseDataFile()<br /> <br /> #---- Write GrADS control file<br /> gDataInfo.WriteGrADSCTLFile()<br /> <br /> print 'Finished...'

netCDF数据sst等值线图：


转为GrADS格点数据sst等值线图：

topmad

来顶王老师咯....身边用meteoinfo的人越来越多了

shawmelo

好东西，来学习下

sangsang

继续学习

huzengyun

求助： 利用ArcGis  处理netcdf 出问题了。怎么处理啊

MeteoInfo

huzengyun 发表于 2013-1-24 12:33

                               
登录/注册后可看大图

求助： 利用ArcGis  处理netcdf 出问题了。怎么处理啊

那就用MeteoInfo试试。

HZK

请问王老师，MeteoInfo怎么打开grd文件啊？

HZK

不是用MeteoInfo直接打开grd文件对应的ctl文件吗？打开什么都不显示，但是“数据信息”和“显示数据”都正常

MeteoInfo

HZK 发表于 2013-3-3 22:08

                               
登录/注册后可看大图

不是用MeteoInfo直接打开grd文件对应的ctl文件吗？打开什么都不显示，但是“数据信息”和“显示数据”都正常 ...

也许是因为你的数据范围太小了，你选择生成的图层名，在工具栏中点击“缩放至图层范围”按钮试试。

HZK

MeteoInfo 发表于 2013-3-4 08:35

                               
登录/注册后可看大图

也许是因为你的数据范围太小了，你选择生成的图层名，在工具栏中点击“缩放至图层范围”按钮试试。

出图了，不过没在设定的经纬度上，都跑到地图外了。。这个grd是grads生成的，ctl自己写的
dset e:\grads\surft.grd
undef -9.99e33
xdef 6 linear 27.5 2.5
ydef 6 linear 110 2.5
tdef 11689 linear jan1981 1dy
zdef 1 levels 0
vars 1
air 1 0 t surfaceair
endvars


应该没错啊