求助，如何利用NECP FNL资料、MetoInfo里的插件TrajStat，绘制气流轨迹，聚类分析？

驴小杏

已经下载了NECP FNL资料，grib2格式的，安装好了MeteoInfo Java 版和TrajStat插件，已经看了TrajStat的帮助文档，没找着该怎么生成第一步的气流轨迹，困在第一步，恳请知道的老师、同学帮忙解答一下，不胜感激啊~刚接触，确实啥也不懂，看了半天也没明白

MeteoInfo

HYSPLIT模式需要专门格式的气象数据（ARL格式），grib2格式的气象数据需要先转换。建议你直接下载ARL格式的气象数据来做轨迹，详见此贴：http://bbs.06climate.com/forum.p ... mp;page=3#pid307656

如果一定要用其它格式气象数据，也可以通过写MeteoInfo脚本程序进行格式转换。

驴小杏

老师，格式转换麻烦吗？刚才试了一下下载ARL格式的，想分析一年的，下载速度好像太慢了

驴小杏

MeteoInfo 发表于 2014-12-23 10:44
HYSPLIT模式需要专门格式的气象数据（ARL格式），grib2格式的气象数据需要先转换。建议你直接下载ARL格式的 ...

老师，转换格式麻烦吗? 刚才试了一下ARL格式的，我想分析全年的，下载速度好像太慢了点，

MeteoInfo

驴小杏 发表于 2014-12-23 11:56
老师，转换格式麻烦吗? 刚才试了一下ARL格式的，我想分析全年的，下载速度好像太慢了点，

之前写过的一个脚本程序，根据你自己的需要修改：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2013-10-324                                               
   
4. # Purpose: Convert GRIB data to ARL data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata.arl import ARLDataInfo
   
9. from org.meteoinfo.data.meteodata.arl import DataLabel
   
10. import os
    
11. 
    
12. #---- Set directories
    
13. dataDir = "D:/Temp"
    
14. 
    
15. #---- Set output data file
    
16. outFile = os.path.join(dataDir, 'arl/test2.arl')
    
17. 
    
18. #---- Read a GRIB data file
    
19. mydata = MeteoDataInfo()
    
20. infile = os.path.join(dataDir, 'grib/201001011800.pgbh06.gdas.20100101-20100105.grb2')
    
21. print infile
    
22. mydata.openNetCDFData(infile)
    
23. print 'GRIB file has been opened...'
    
24. 
    
25. #---- Set output ARL data info
    
26. arlDI = ARLDataInfo()
    
27. 
    
28. #---- Set variable and level list
    
29. gvars = ['Pressure_surface','Temperature_height_above_ground',\
    
30.   'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',\
    
31.   'Geopotential_height_isobaric','Temperature_isobaric',\
    
32.   'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',\
    
33.   'Relative_humidity_isobaric']
    
34. avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH']
    
35. levels = [0,1000,975,950,925,900,875,850,825,800,775,750,700,\
    
36.   650,600,550,500,450,400,350,300,250,225,200,175,150,\
    
37.   125,100,70,50,30,20,10,7,5,3,2,1]
    
38. for l in levels:
    
39.   arlDI.levels.add(l)
    
40.   if l == 0:
    
41.     arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M'])
    
42.   else:
    
43.     arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])
    
44. 
    
45. #---- Write ARL data file
    
46. dataInfo = mydata.getDataInfo()
    
47. arlDI.createDataFile(outFile)
    
48. arlDI.X = dataInfo.getXDimension().getValues()
    
49. arlDI.Y = dataInfo.getYDimension().getValues()
    
50. variables = dataInfo.getVariables()
    
51. tNum = dataInfo.getTimeNum()
    
52. for t in range(0, tNum):
    
53.   mydata.setTimeIndex(t)
    
54.   atime = dataInfo.getTimes().get(t)
    
55.   aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'FNL1', 2)
    
56.   arlDI.writeIndexRecord(atime, aDH)
    
57.   lidx = 0
    
58.   for l in arlDI.levels:
    
59.     print l
    
60.     for v in arlDI.LevelVarList[lidx]:
    
61.       vName = gvars[avars.index(v)]
    
62.       print vName
    
63.       if lidx == 0:
    
64.         mydata.setLevelIndex(lidx)
    
65.       else:
    
66.         variable = dataInfo.getVariable(vName)
    
67.         nlidx = variable.getZDimension().getDimValue().indexOf(l*100.0)
    
68.         mydata.setLevelIndex(nlidx)
    
69.       gData = mydata.getGridData(vName)
    
70.       if v == 'PRSS' or v == 'WWND':
    
71.         gData = gData.divide(100)      
    
72.       aDL = DataLabel(atime)
    
73.       aDL.setLevel(lidx)
    
74.       aDL.setVarName(v)
    
75.       aDL.setGrid(99)
    
76.       aDL.setForecast(0)
    
77.       arlDI.writeGridData(aDL, gData)
    
78.     lidx += 1
    
79. 
    
80. arlDI.closeDataFile()
    
81. 
    
82. print 'Finished!'

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2013-10-324                                               <br /> # Purpose: Convert GRIB data to ARL data  <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata.arl import ARLDataInfo<br /> from org.meteoinfo.data.meteodata.arl import DataLabel<br /> import os<br /> <br /> #---- Set directories<br /> dataDir = "D:/Temp"<br /> <br /> #---- Set output data file<br /> outFile = os.path.join(dataDir, 'arl/test2.arl')<br /> <br /> #---- Read a GRIB data file<br /> mydata = MeteoDataInfo()<br /> infile = os.path.join(dataDir, 'grib/201001011800.pgbh06.gdas.20100101-20100105.grb2')<br /> print infile<br /> mydata.openNetCDFData(infile)<br /> print 'GRIB file has been opened...'<br /> <br /> #---- Set output ARL data info<br /> arlDI = ARLDataInfo()<br /> <br /> #---- Set variable and level list<br /> gvars = ['Pressure_surface','Temperature_height_above_ground',\<br />   'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',\<br />   'Geopotential_height_isobaric','Temperature_isobaric',\<br />   'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',\<br />   'Relative_humidity_isobaric']<br /> avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH']<br /> levels = [0,1000,975,950,925,900,875,850,825,800,775,750,700,\<br />   650,600,550,500,450,400,350,300,250,225,200,175,150,\<br />   125,100,70,50,30,20,10,7,5,3,2,1]<br /> for l in levels:<br />   arlDI.levels.add(l)<br />   if l == 0:<br />     arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M'])<br />   else:<br />     arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])<br /> <br /> #---- Write ARL data file<br /> dataInfo = mydata.getDataInfo()<br /> arlDI.createDataFile(outFile)<br /> arlDI.X = dataInfo.getXDimension().getValues()<br /> arlDI.Y = dataInfo.getYDimension().getValues()<br /> variables = dataInfo.getVariables()<br /> tNum = dataInfo.getTimeNum()<br /> for t in range(0, tNum):<br />   mydata.setTimeIndex(t)<br />   atime = dataInfo.getTimes().get(t)<br />   aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'FNL1', 2)<br />   arlDI.writeIndexRecord(atime, aDH)<br />   lidx = 0<br />   for l in arlDI.levels:<br />     print l<br />     for v in arlDI.LevelVarList[lidx]:<br />       vName = gvars[avars.index(v)]<br />       print vName<br />       if lidx == 0:<br />         mydata.setLevelIndex(lidx)<br />       else:<br />         variable = dataInfo.getVariable(vName)<br />         nlidx = variable.getZDimension().getDimValue().indexOf(l*100.0)<br />         mydata.setLevelIndex(nlidx)<br />       gData = mydata.getGridData(vName)<br />       if v == 'PRSS' or v == 'WWND':<br />         gData = gData.divide(100)      <br />       aDL = DataLabel(atime)<br />       aDL.setLevel(lidx)<br />       aDL.setVarName(v)<br />       aDL.setGrid(99)<br />       aDL.setForecast(0)<br />       arlDI.writeGridData(aDL, gData)<br />     lidx += 1<br /> <br /> arlDI.closeDataFile()<br /> <br /> print 'Finished!'

驴小杏

MeteoInfo 发表于 2014-12-23 13:46
之前写过的一个脚本程序，根据你自己的需要修改：

好的，谢谢老师！

驴小杏

MeteoInfo 发表于 2014-12-23 13:46
之前写过的一个脚本程序，根据你自己的需要修改：

老师，看不明白，

MeteoInfo

驴小杏 发表于 2014-12-23 14:24
老师，看不明白，

如果你一点Python的基础也没有，看不明白很正常。

驴小杏

MeteoInfo 发表于 2014-12-23 14:49
如果你一点Python的基础也没有，看不明白很正常。

真没接触过Python，先下载数据看看了

枫之山峰

好高端，帮顶