发布MeteoInfo Java 0.9.8.5

MeteoInfo

主要是一些功能的完善和bug修正。

应HYSPLIT论坛上一些用户的需求，增加了下列功能：

给MeteoDataInfo类增加了getArrivalTimeData(String varName, double threshold)方法用了获取变量达到某个临界值所需的时间（距起始时间的小时数），两个参数分别是变量名和临界值。目前只能用于脚本，以后再看有没有必要加入GUI界面中。示例脚本程序：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-5-3
   
4. # Purpose: Get time of arrival data and plot
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. print 'Loading classes...'
   
8. from org.meteoinfo.plot import MapLayoutPlot
   
9. from org.meteoinfo.data import GridData
   
10. from org.meteoinfo.data.meteodata import MeteoDataInfo
    
11. from org.meteoinfo.data.meteodata import DrawMeteoData
    
12. from org.meteoinfo.legend import LegendScheme
    
13. from org.meteoinfo.legend import LegendManage
    
14. from org.meteoinfo.legend import LegendType
    
15. from org.meteoinfo.shape import ShapeTypes
    
16. import os.path
    
17. 
    
18. #Set directory
    
19. dataDir = 'D:/Temp/HYSPLIT'
    
20. infn = os.path.join(dataDir, 'HYSPLIT_cdump')
    
21. 
    
22. #Open data
    
23. print 'Open data...'
    
24. mdi = MeteoDataInfo()
    
25. mdi.openHYSPLITConcData(infn)
    
26. tData = mdi.getArrivalTimeData('DUST', 0.001)
    
27. #Plot data
    
28. print 'Plot data...'
    
29. aLS = LegendManage.createLegendSchemeFromGridData(tData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
30. #aLayer = DrawMeteoData.createShadedLayer(tData, aLS, 'TimeOfArrial_DUST', 'Hour', True)
    
31. aLayer = DrawMeteoData.createRasterLayer(tData, 'TimeOfArrial_DUST', aLS)
    
32. aLayer.setProjInfo(mdi.getProjectionInfo())
    
33. mf = miapp.getMapDocument().getActiveMapFrame()
    
34. mf.addLayer(aLayer)
    
35. mf.moveLayer(aLayer, 0)
    
36. miapp.getMapDocument().getMapLayout().paintGraphics()
    
37. 
    
38. print 'Finished'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-5-3<br /> # Purpose: Get time of arrival data and plot<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> print 'Loading classes...'<br /> from org.meteoinfo.plot import MapLayoutPlot<br /> from org.meteoinfo.data import GridData<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata import DrawMeteoData<br /> from org.meteoinfo.legend import LegendScheme<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.legend import LegendType<br /> from org.meteoinfo.shape import ShapeTypes<br /> import os.path<br /> <br /> #Set directory<br /> dataDir = 'D:/Temp/HYSPLIT'<br /> infn = os.path.join(dataDir, 'HYSPLIT_cdump')<br /> <br /> #Open data<br /> print 'Open data...'<br /> mdi = MeteoDataInfo()<br /> mdi.openHYSPLITConcData(infn)<br /> tData = mdi.getArrivalTimeData('DUST', 0.001)<br /> #Plot data<br /> print 'Plot data...'<br /> aLS = LegendManage.createLegendSchemeFromGridData(tData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> #aLayer = DrawMeteoData.createShadedLayer(tData, aLS, 'TimeOfArrial_DUST', 'Hour', True)<br /> aLayer = DrawMeteoData.createRasterLayer(tData, 'TimeOfArrial_DUST', aLS)<br /> aLayer.setProjInfo(mdi.getProjectionInfo())<br /> mf = miapp.getMapDocument().getActiveMapFrame()<br /> mf.addLayer(aLayer)<br /> mf.moveLayer(aLayer, 0)<br /> miapp.getMapDocument().getMapLayout().paintGraphics()<br /> <br /> print 'Finished'

给MeteoDataInfo类增加了toStation方法，有几个重载，用来获取任意时空轨迹上的气象数据（线性差值）。示例脚本程序：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-5-4
   
4. # Purpose: Get meteorological data along trajectory
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. import os
   
8. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
9. from org.meteoinfo.global.table import DataTypes
   
10. from java.text import SimpleDateFormat
    
11. 
    
12. # Set working directory
    
13. trajDir = 'D:/Trajectory/bj'
    
14. meteoDir = 'P:/MetData/MetData/2011'
    
15. 
    
16. # Open trjactory data file
    
17. print 'Open trajectory data file ...'
    
18. # Create a MeteoDataInfo object
    
19. mdi_traj = MeteoDataInfo()
    
20. # Open trajectory data file
    
21. trajFile = '11033120'
    
22. trajfn = os.path.join(trajDir, trajFile)
    
23. print 'Trajectory file: ' + trajfn
    
24. mdi_traj.openHYSPLITTrajData(trajfn)
    
25. # Create trajectory layer
    
26. trajLayer = mdi_traj.getDataInfo().createTrajPointLayer()
    
27. trajLayer.setLayerName('Trajectorie' + '_' + trajFile)
    
28. 
    
29. # Open meteorological data file
    
30. print 'Open meteorological data file...'
    
31. # Create a MeteoDataInfo object
    
32. mdi_meteo = MeteoDataInfo()
    
33. # Open meteorological data file
    
34. meteoFile = 'gdas1.mar11.w5'
    
35. meteofn = os.path.join(meteoDir, meteoFile)
    
36. print 'Meteorological file: ' + meteofn
    
37. mdi_meteo.openARLData(meteofn)
    
38. 
    
39. # Get meteorological data along trajectory
    
40. print 'Get meteorological data along trajectory...'
    
41. varName = 'UWND'
    
42. trajLayer.editAddField(varName, DataTypes.Float)
    
43. dformat = SimpleDateFormat('yyyyMMddHH')
    
44. idx = 0
    
45. for ps in trajLayer.getShapes():
    
46.   lon = ps.getPoint().X
    
47.   lat = ps.getPoint().Y
    
48.   pres = trajLayer.getCellValue('Pressure', idx)
    
49.   t = dformat.parse(trajLayer.getCellValue('Date', idx))
    
50.   v = mdi_meteo.toStation(varName, lon, lat, pres, t)
    
51.   trajLayer.editCellValue(varName, idx, v)
    
52.   idx += 1
    
53. 
    
54. # Add the layer into the MeteoInfo application (miapp)
    
55. print 'Add trajectory layer...'
    
56. mf = miapp.getMapDocument().getActiveMapFrame()
    
57. mf.addLayer(trajLayer)
    
58. 
    
59. print 'Finish...'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-5-4<br /> # Purpose: Get meteorological data along trajectory<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> import os<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.global.table import DataTypes<br /> from java.text import SimpleDateFormat<br /> <br /> # Set working directory<br /> trajDir = 'D:/Trajectory/bj'<br /> meteoDir = 'P:/MetData/MetData/2011'<br /> <br /> # Open trjactory data file<br /> print 'Open trajectory data file ...'<br /> # Create a MeteoDataInfo object<br /> mdi_traj = MeteoDataInfo()<br /> # Open trajectory data file<br /> trajFile = '11033120'<br /> trajfn = os.path.join(trajDir, trajFile)<br /> print 'Trajectory file: ' + trajfn<br /> mdi_traj.openHYSPLITTrajData(trajfn)<br /> # Create trajectory layer<br /> trajLayer = mdi_traj.getDataInfo().createTrajPointLayer()<br /> trajLayer.setLayerName('Trajectorie' + '_' + trajFile)<br /> <br /> # Open meteorological data file<br /> print 'Open meteorological data file...'<br /> # Create a MeteoDataInfo object<br /> mdi_meteo = MeteoDataInfo()<br /> # Open meteorological data file<br /> meteoFile = 'gdas1.mar11.w5'<br /> meteofn = os.path.join(meteoDir, meteoFile)<br /> print 'Meteorological file: ' + meteofn<br /> mdi_meteo.openARLData(meteofn)<br /> <br /> # Get meteorological data along trajectory<br /> print 'Get meteorological data along trajectory...'<br /> varName = 'UWND'<br /> trajLayer.editAddField(varName, DataTypes.Float)<br /> dformat = SimpleDateFormat('yyyyMMddHH')<br /> idx = 0<br /> for ps in trajLayer.getShapes():<br />   lon = ps.getPoint().X<br />   lat = ps.getPoint().Y<br />   pres = trajLayer.getCellValue('Pressure', idx)<br />   t = dformat.parse(trajLayer.getCellValue('Date', idx))<br />   v = mdi_meteo.toStation(varName, lon, lat, pres, t)<br />   trajLayer.editCellValue(varName, idx, v)<br />   idx += 1<br /> <br /> # Add the layer into the MeteoInfo application (miapp)<br /> print 'Add trajectory layer...'<br /> mf = miapp.getMapDocument().getActiveMapFrame()<br /> mf.addLayer(trajLayer)<br /> <br /> print 'Finish...'

veronicasyd

看上去很便捷好用的样子~~~
是不是比较适合用于对气候预测模型的分析呢？最近在分析CMIP5的模型的模拟效果，各种软件弄得头疼~~
不知道这个是不是会比较直观一些呢？

MeteoInfo

veronicasyd 发表于 2014-5-4 21:20
看上去很便捷好用的样子~~~
是不是比较适合用于对气候预测模型的分析呢？最近在分析CMIP5的模型的模拟效果 ...

这个要看你具体想做什么了

chensr1989

谢谢楼主，学习了

veronicasyd

MeteoInfo 发表于 2014-5-4 21:43
这个要看你具体想做什么了

我想要利用CMIP5的各个模式来预测东北的气候变化状况（主要是降水和温度）。

MeteoInfo

veronicasyd 发表于 2014-5-5 10:30
我想要利用CMIP5的各个模式来预测东北的气候变化状况（主要是降水和温度）。

这个主要是要转好模式，后期分析和绘图也要看具体想分析什么了

灵心流萤

学习了，多谢楼主

fzghui

要写脚本，刚装软件还不知道怎么用，在哪里输入。。。。

MeteoInfo

fzghui 发表于 2014-9-29 20:24
要写脚本，刚装软件还不知道怎么用，在哪里输入。。。。

“工具” -> “脚本”菜单

fzghui

MeteoInfo 发表于 2014-9-29 21:17
“工具” -> “脚本”菜单

嗯，会了，谢拉