netCDF数据转为GrADS数据

MeteoInfo · 发表于 2012-5-29 15:53:29

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 MeteoInfo 于 2014-4-6 17:57 编辑

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

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

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

netCDF数据sst等值线图：

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

topmad · 发表于 2012-5-29 16:46:16

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

shawmelo · 发表于 2012-5-29 20:48:26

好东西，来学习下

sangsang · 发表于 2012-6-11 14:36:13

继续学习

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

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

MeteoInfo · 发表于 2013-1-24 15:09:38

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

登录/注册后可看大图

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

那就用MeteoInfo试试。

HZK · 发表于 2013-3-3 21:53:47

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

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

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

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

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

登录/注册后可看大图

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

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

HZK · 发表于 2013-3-4 11:29:38

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

应该没错啊

		自动登录	找回密码
密码			立即注册

netCDF数据转为GrADS数据

登录后查看更多精彩内容~

评分