GRIB格式转为ARL格式

dolces · 发表于 2013-1-16 10:30:01

MeteoInfo 发表于 2013-1-15 09:42

登录/注册后可看大图

我下载一个0.5度的数据试试吧。

太谢谢了王老师

MeteoInfo · 发表于 2013-1-16 14:01:59

本帖最后由 MeteoInfo 于 2013-1-16 14:03 编辑

dolces 发表于 2013-1-16 10:30

登录/注册后可看大图

太谢谢了王老师

下了个0.5度的数据试了试，现在应该可以了，需要MeteoInfo最新文件（见置顶帖子，或者在www.meteothinker.com上下载更新文件）。示例脚本程序如下，脚本中读取0.5度ARL数据，挑出一些高度层数据输出。

#-----------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2013-1-16
# Purpose: Read and write ARL data
# Note: Sample
#-----------------------------------------------------------
import clr
from System import *
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *
#---- Set directories
dataDir = "D:\\SampleData"
#---- Set output data file
outFile = dataDir + '\\test.arl'
#---- Read a GRIB data file
mydata = MeteoDataInfo()
infile = dataDir + '\\20130113_gdas0p5'
mydata.OpenARLData(infile)
inarlDI = mydata.DataInfo
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
arlDI.X = mydata.GetX()
arlDI.Y = mydata.GetY()
arlDI.levels.Add(inarlDI.levels[0])
arlDI.LevelVarList.Add(inarlDI.LevelVarList[0])
for lidx in range(1, inarlDI.levels.Count, 2):
arlDI.levels.Add(inarlDI.levels[lidx])
arlDI.LevelVarList.Add(inarlDI.LevelVarList[lidx])
#---- Write ARL data file
arlDI.CreateDataFile(outFile)
tNum = mydata.GetTimeNumber()
for t in range(0, tNum):
mydata.TimeIndex = t
atime = mydata.GetTime(t)
print atime.ToString("yyyy-MM-dd HH")
aDH = arlDI.GetDataHead(mydata.ProjInfo, 'GDAS', 4)
arlDI.WriteIndexRecord(atime, aDH)
lidx = 0
for l in arlDI.levels:
print l
nlidx = inarlDI.levels.IndexOf(l)
if nlidx > 0:
nlidx -= 1
print nlidx
mydata.LevelIndex = nlidx
for v in arlDI.LevelVarList[lidx]:
print v
gData = mydata.GetGridData(v)
aDL = DataLabel(atime)
aDL.Level = lidx
aDL.Variable = v
aDL.Grid = 99
aDL.Forecast = 0
arlDI.WriteGridData(aDL, gData)
lidx += 1
arlDI.CloseDataFile()
print 'Finished!'

#-----------------------------------------------------------         
# Author: Yaqiang Wang                                           
# Date: 2013-1-16                                                
# Purpose: Read and write ARL data   
# Note: Sample                                                    
#----------------------------------------------------------- 
import clr 
from System import * 
from System.Collections.Generic import * 
clr.AddReference("MeteoInfoC.dll") 
from MeteoInfoC import * 
from MeteoInfoC.Data import * 
from MeteoInfoC.Data.MeteoData import * 
 
#---- Set directories 
dataDir = "D:\\SampleData" 
 
#---- Set output data file 
outFile = dataDir + '\\test.arl' 
 
#---- Read a GRIB data file 
mydata = MeteoDataInfo() 
infile = dataDir + '\\20130113_gdas0p5' 
mydata.OpenARLData(infile) 
inarlDI = mydata.DataInfo 
 
#---- Set output ARL data info 
arlDI = ARLDataInfo() 
 
#---- Set variable and level list 
arlDI.X = mydata.GetX() 
arlDI.Y = mydata.GetY() 
arlDI.levels.Add(inarlDI.levels[0]) 
arlDI.LevelVarList.Add(inarlDI.LevelVarList[0]) 
for lidx in range(1, inarlDI.levels.Count, 2): 
   arlDI.levels.Add(inarlDI.levels[lidx]) 
   arlDI.LevelVarList.Add(inarlDI.LevelVarList[lidx]) 
 
#---- Write ARL data file 
arlDI.CreateDataFile(outFile) 
tNum = mydata.GetTimeNumber() 
for t in range(0, tNum): 
    mydata.TimeIndex = t 
    atime = mydata.GetTime(t) 
    print atime.ToString("yyyy-MM-dd HH") 
    aDH = arlDI.GetDataHead(mydata.ProjInfo, 'GDAS', 4) 
    arlDI.WriteIndexRecord(atime, aDH) 
    lidx = 0 
    for l in arlDI.levels: 
        print l 
        nlidx = inarlDI.levels.IndexOf(l) 
        if nlidx > 0: 
            nlidx -= 1 
        print nlidx 
        mydata.LevelIndex = nlidx 
        for v in arlDI.LevelVarList[lidx]:    
            print v 
            gData = mydata.GetGridData(v) 
            aDL = DataLabel(atime) 
            aDL.Level = lidx 
            aDL.Variable = v 
            aDL.Grid = 99 
            aDL.Forecast = 0 
            arlDI.WriteGridData(aDL, gData) 
        lidx += 1 
 
arlDI.CloseDataFile() 
 
print 'Finished!'

dolces · 发表于 2013-1-17 13:30:49

本帖最后由 dolces 于 2013-1-17 13:42 编辑

MeteoInfo 发表于 2013-1-16 14:01

登录/注册后可看大图

下了个0.5度的数据试了试，现在应该可以了，需要MeteoInfo最新文件（见置顶帖子，或者在www.meteothinker ...

王老师，我用这个script转0.5度的数据，不间隔取层数，就是说56层全部重写到新文件里的话，出来的数据是1.7G，但是原数据只有479M，你测试的时候有没有出现这样的情况？

另外，把新写的数据放到HYSPLIT里算轨迹出错了，出错信息如下：
*ERROR* metset: 2nd time period INDX record missing
Calculation Started ... please be patient
*ERROR* metpos: start point not within (x,y,t) any data file
- start time before start of meteorology data

dolces · 发表于 2013-1-17 13:54:55

本帖最后由 dolces 于 2013-1-17 14:59 编辑

dolces 发表于 2013-1-17 13:30

登录/注册后可看大图

王老师，我用这个script转0.5度的数据，不间隔取层数，就是说56层全部重写到新文件里的话，出来的数据是1 ...

我查看了原数据和新数据的信息，有几个参数的值不同，不知道对你有没有用：
原数据：
File Name: E:\test\20110101_gdas0p5
File Start Time: 2011-01-01 00:00
File End Time: 2011-01-01 21:00
Record Length Bytes: 259970
Meteo Data Model: GHDA
Xsize = 720  Ysize = 361  Zsize = 56  Tsize = 8
Record Per Time: 236
Pole pnt lat/lon: 90  359.5
Reference pnt lat/lon: 0.5  0.5
Grid Size: 0
Orientation: 0
Tan lat/cone: 0
Syn pnt x/y: 1  1
Syn pnt lat/lon: -90  0

新数据：
File Name: E:\test\test.arl
File Start Time: 2011-01-01 00:00
File End Time: 2011-01-01 21:00
Record Length Bytes: 259970
Meteo Data Model: GDAS
Xsize = 720  Ysize = 361  Zsize = 56  Tsize = 8
Record Per Time: 841
Pole pnt lat/lon: 90  0
Reference pnt lat/lon: 0.5  0.5
Grid Size: 0
Orientation: 0
Tan lat/cone: 0
Syn pnt x/y: 1  1
Syn pnt lat/lon: -90  0

dolces · 发表于 2013-1-17 13:57:17

Record Per Time: 841 这个参数应该是新数据大小是原数据大概3.6倍的原因

MeteoInfo · 发表于 2013-1-17 16:27:20

dolces 发表于 2013-1-17 13:57

登录/注册后可看大图

Record Per Time: 841 这个参数应该是新数据大小是原数据大概3.6倍的原因

我试了试，所有层都输出和源文件大小是一样的，估计你的script程序改得有问题。你依据下面的script程序重新转一下，然后用HYSPLIT测试一下吧。

#-----------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2013-1-16
# Purpose: Read and write ARL data
# Note: Sample
#-----------------------------------------------------------
import clr
from System import *
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *
#---- Set directories
dataDir = "D:\\SampleData"
#---- Set output data file
outFile = dataDir + '\\test.arl'
#---- Read a GRIB data file
mydata = MeteoDataInfo()
infile = dataDir + '\\20130113_gdas0p5'
mydata.OpenARLData(infile)
inarlDI = mydata.DataInfo
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
arlDI.X = mydata.GetX()
arlDI.Y = mydata.GetY()
for lidx in range(0, inarlDI.levels.Count):
arlDI.levels.Add(inarlDI.levels[lidx])
arlDI.LevelVarList.Add(inarlDI.LevelVarList[lidx])
#---- Write ARL data file
arlDI.CreateDataFile(outFile)
tNum = mydata.GetTimeNumber()
for t in range(0, tNum):
mydata.TimeIndex = t
atime = mydata.GetTime(t)
print atime.ToString("yyyy-MM-dd HH")
aDH = arlDI.GetDataHead(mydata.ProjInfo, 'GDAS', 4)
arlDI.WriteIndexRecord(atime, aDH)
lidx = 0
for l in arlDI.levels:
print l
nlidx = inarlDI.levels.IndexOf(l)
if nlidx > 0:
nlidx -= 1
print nlidx
mydata.LevelIndex = nlidx
for v in arlDI.LevelVarList[lidx]:
print v
gData = mydata.GetGridData(v)
aDL = DataLabel(atime)
aDL.Level = lidx
aDL.Variable = v
aDL.Grid = 99
aDL.Forecast = 0
arlDI.WriteGridData(aDL, gData)
lidx += 1
arlDI.CloseDataFile()
print 'Finished!'

#-----------------------------------------------------------         
# Author: Yaqiang Wang                                           
# Date: 2013-1-16                                                
# Purpose: Read and write ARL data   
# Note: Sample                                                    
#----------------------------------------------------------- 
import clr 
from System import * 
from System.Collections.Generic import * 
clr.AddReference("MeteoInfoC.dll") 
from MeteoInfoC import * 
from MeteoInfoC.Data import * 
from MeteoInfoC.Data.MeteoData import * 
 
#---- Set directories 
dataDir = "D:\\SampleData" 
 
#---- Set output data file 
outFile = dataDir + '\\test.arl' 
 
#---- Read a GRIB data file 
mydata = MeteoDataInfo() 
infile = dataDir + '\\20130113_gdas0p5' 
mydata.OpenARLData(infile) 
inarlDI = mydata.DataInfo 
 
#---- Set output ARL data info 
arlDI = ARLDataInfo() 
 
#---- Set variable and level list 
arlDI.X = mydata.GetX() 
arlDI.Y = mydata.GetY() 
for lidx in range(0, inarlDI.levels.Count): 
   arlDI.levels.Add(inarlDI.levels[lidx]) 
   arlDI.LevelVarList.Add(inarlDI.LevelVarList[lidx]) 
 
#---- Write ARL data file 
arlDI.CreateDataFile(outFile) 
tNum = mydata.GetTimeNumber() 
for t in range(0, tNum): 
    mydata.TimeIndex = t 
    atime = mydata.GetTime(t) 
    print atime.ToString("yyyy-MM-dd HH") 
    aDH = arlDI.GetDataHead(mydata.ProjInfo, 'GDAS', 4) 
    arlDI.WriteIndexRecord(atime, aDH) 
    lidx = 0 
    for l in arlDI.levels: 
        print l 
        nlidx = inarlDI.levels.IndexOf(l) 
        if nlidx > 0: 
            nlidx -= 1 
        print nlidx 
        mydata.LevelIndex = nlidx 
        for v in arlDI.LevelVarList[lidx]:    
            print v 
            gData = mydata.GetGridData(v) 
            aDL = DataLabel(atime) 
            aDL.Level = lidx 
            aDL.Variable = v 
            aDL.Grid = 99 
            aDL.Forecast = 0 
            arlDI.WriteGridData(aDL, gData) 
        lidx += 1 
 
arlDI.CloseDataFile() 
 
print 'Finished!'

dolces · 发表于 2013-1-17 16:56:58

本帖最后由 dolces 于 2013-1-17 17:05 编辑

MeteoInfo 发表于 2013-1-17 16:27

登录/注册后可看大图

我试了试，所有层都输出和源文件大小是一样的，估计你的script程序改得有问题。你依据下面的script程序重 ...

这次数据大小没问题，轨迹也画出来了，辛苦了王老师。

dolces · 发表于 2013-1-17 17:14:00

MeteoInfo 发表于 2013-1-17 16:27

登录/注册后可看大图

我试了试，所有层都输出和源文件大小是一样的，估计你的script程序改得有问题。你依据下面的script程序重 ...

还有个小问题，level 0是地面，层数越大高度越高对吗？

MeteoInfo · 发表于 2013-1-17 17:36:08

dolces 发表于 2013-1-17 17:14

登录/注册后可看大图

还有个小问题，level 0是地面，层数越大高度越高对吗？

level 0是地面，hybrid坐标我也不是很清楚，你查查资料吧。

qxjy123 · 发表于 2014-4-3 07:53:36

感谢楼主。学习了

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