MeteoInfo脚本教程（十四）：轨迹数据

MeteoInfo

前面讲了格点和站点数据，除此以外轨迹数据也是气象领域常用的数据格式，比如气团轨迹、台风路径等。MeteoInfo目前直接支持MICAPS第7类格式（台风数据）和HYSPLIT模式输出的气团轨迹数据格式。其实数据格式有无数，使用脚本可以自己读各种格式中的轨迹数据。这里演示一个绘制HYSPLIT模式输出气团轨迹的例子。

用MeteoDataInfo类的openHYSPLITTrajData方法打开轨迹数据文件，并获取HYSPLITTrajDataInfo类的对象：
    mdi = MeteoDataInfo()
    mdi.openHYSPLITTrajData(fn)
    dataInfo = mdi.getDataInfo()

利用HYSPLITTrajDataInfo类的createTrajLayer生成气团轨迹线图层，createTrajPointLayer生成气团轨迹点图层，createTrajStartPointLayer生成轨迹起点点图层。

要绘制一维图，首先要构造相应的XY数据（XYDataset类），然后创建XYPlot类对象，再生成Chart类对象，最后在MapLayout中添加一个LayoutChart即可。

脚本代码如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-19
   
4. # Purpose: Plot HYSPLIT trajectory data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. #---- Import classes
   
8. print 'Import classes...'
   
9. from org.meteoinfo.layout import MapLayout, PaperSize
   
10. from org.meteoinfo.data.mapdata import MapDataManage
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo
    
12. from org.meteoinfo.global import Extent
    
13. from org.meteoinfo.drawing import PointStyle
    
14. from org.meteoinfo.chart import Chart
    
15. from org.meteoinfo.chart.plot import XYPlot, ChartPlotMethod
    
16. import os
    
17. from java.awt import Color
    
18. from javax.swing import JFrame
    
19. 
    
20. #---- Set directories
    
21. print 'Set directories...'
    
22. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
23. dataDir = os.path.join(baseDir, 'sample')
    
24. mapDir = os.path.join(baseDir, 'map')
    
25. figDir = 'D:/Temp/test'
    
26. 
    
27. #---- Create MapLayout object
    
28. mapLayout = MapLayout()
    
29. mapFrame = mapLayout.getActiveMapFrame()
    
30. ps = PaperSize('Custum', 650, 640)
    
31. mapLayout.setPaperSize(ps)
    
32. 
    
33. #---- Load country layer
    
34. print 'Load country layer...'
    
35. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
36. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
37. lb.setOutlineColor(Color.black)
    
38. lb.setOutlineSize(2)
    
39. mapFrame.addLayer(countryLayer)
    
40. 
    
41. #---- Create MeteoDataInfo object
    
42. mdi = MeteoDataInfo()
    
43. 
    
44. #---- Open a trajectory data file
    
45. print 'Open trajectory data file...'
    
46. fn = os.path.join(dataDir, 'HYSPLIT/tdump')
    
47. mdi.openHYSPLITTrajData(fn)
    
48. dataInfo = mdi.getDataInfo()
    
49. 
    
50. #---- Create trajectory layers
    
51. print 'Create trajectory layers...'
    
52. trajLayer = dataInfo.createTrajLineLayer()
    
53. startLayer = dataInfo.createTrajStartPointLayer()
    
54. lb = startLayer.getLegendScheme().getLegendBreaks().get(0)
    
55. lb.setStyle(PointStyle.Star)
    
56. 
    
57. #---- Add layer
    
58. mapFrame.addLayer(trajLayer)
    
59. mapFrame.addLayer(startLayer)
    
60. 
    
61. #---- Add pressure profile chart
    
62. print 'Add pressure profile chart...'
    
63. dataset = dataInfo.getXYDataset(12)    #Pressure is in 12th column
    
64. plot = XYPlot(True, dataset)
    
65. plot.setChartPlotMethod(ChartPlotMethod.LINE_POINT)
    
66. plot.getXAxis().setLabel('Time')
    
67. plot.getYAxis().setLabel('hPa')
    
68. plot.getYAxis().setInverse(True)
    
69. ls = trajLayer.getLegendScheme()
    
70. for i in range(0, dataset.getSeriesCount()):
    
71.   plot.setPolylineBreak(i, ls.getLegendBreaks().get(i))
    
72. chart = Chart(None, plot)
    
73. chart.setDrawLegend(False)
    
74. layoutChart = mapLayout.addChart(10, 410)
    
75. layoutChart.setWidth(590)
    
76. layoutChart.setHeight(200)
    
77. layoutChart.setChart(chart)
    
78. 
    
79. #---- Add title
    
80. title = mapLayout.addText('MeteoInfo script demo - Trajectory', 300, 30, 'Arial', 16)
    
81. 
    
82. #---- Zoom map
    
83. print 'Zoom map...'
    
84. mapLayout.getActiveLayoutMap().setWidth(550)
    
85. mapLayout.getActiveLayoutMap().setHeight(350)
    
86. mapFrame.getMapView().zoomToExtent(trajLayer.getExtent().extend(2, 2))
    
87. 
    
88. #---- Output figure
    
89. print 'Output figure...'
    
90. mapLayout.exportToPicture(os.path.join(figDir, 'traj_test.png'))
    
91. 
    
92. #---- Set mapframe
    
93. mapFrame.setGridXDelt(5)
    
94. mapFrame.setGridYDelt(5)
    
95. 
    
96. frame = JFrame('MeteoInfo Script Sample', size = (650, 680))
    
97. frame.add(mapLayout)
    
98. frame.visible = True
    
99. 
    
100. print 'Finish...'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-19<br /> # Purpose: Plot HYSPLIT trajectory data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.layout import MapLayout, PaperSize<br /> from org.meteoinfo.data.mapdata import MapDataManage<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.global import Extent<br /> from org.meteoinfo.drawing import PointStyle<br /> from org.meteoinfo.chart import Chart<br /> from org.meteoinfo.chart.plot import XYPlot, ChartPlotMethod<br /> import os<br /> from java.awt import Color<br /> from javax.swing import JFrame<br /> <br /> #---- Set directories<br /> print 'Set directories...'<br /> baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> dataDir = os.path.join(baseDir, 'sample')<br /> mapDir = os.path.join(baseDir, 'map')<br /> figDir = 'D:/Temp/test'<br /> <br /> #---- Create MapLayout object<br /> mapLayout = MapLayout()<br /> mapFrame = mapLayout.getActiveMapFrame()<br /> ps = PaperSize('Custum', 650, 640)<br /> mapLayout.setPaperSize(ps)<br /> <br /> #---- Load country layer<br /> print 'Load country layer...'<br /> countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))<br /> lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)<br /> lb.setOutlineColor(Color.black)<br /> lb.setOutlineSize(2)<br /> mapFrame.addLayer(countryLayer)<br /> <br /> #---- Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #---- Open a trajectory data file<br /> print 'Open trajectory data file...'<br /> fn = os.path.join(dataDir, 'HYSPLIT/tdump')<br /> mdi.openHYSPLITTrajData(fn)<br /> dataInfo = mdi.getDataInfo()<br /> <br /> #---- Create trajectory layers<br /> print 'Create trajectory layers...'<br /> trajLayer = dataInfo.createTrajLineLayer()<br /> startLayer = dataInfo.createTrajStartPointLayer()<br /> lb = startLayer.getLegendScheme().getLegendBreaks().get(0)<br /> lb.setStyle(PointStyle.Star)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(trajLayer)<br /> mapFrame.addLayer(startLayer)<br /> <br /> #---- Add pressure profile chart<br /> print 'Add pressure profile chart...'<br /> dataset = dataInfo.getXYDataset(12)    #Pressure is in 12th column<br /> plot = XYPlot(True, dataset)<br /> plot.setChartPlotMethod(ChartPlotMethod.LINE_POINT)<br /> plot.getXAxis().setLabel('Time')<br /> plot.getYAxis().setLabel('hPa')<br /> plot.getYAxis().setInverse(True)<br /> ls = trajLayer.getLegendScheme()<br /> for i in range(0, dataset.getSeriesCount()):<br />   plot.setPolylineBreak(i, ls.getLegendBreaks().get(i))<br /> chart = Chart(None, plot)<br /> chart.setDrawLegend(False)<br /> layoutChart = mapLayout.addChart(10, 410)<br /> layoutChart.setWidth(590)<br /> layoutChart.setHeight(200)<br /> layoutChart.setChart(chart)<br /> <br /> #---- Add title<br /> title = mapLayout.addText('MeteoInfo script demo - Trajectory', 300, 30, 'Arial', 16)<br /> <br /> #---- Zoom map<br /> print 'Zoom map...'<br /> mapLayout.getActiveLayoutMap().setWidth(550)<br /> mapLayout.getActiveLayoutMap().setHeight(350)<br /> mapFrame.getMapView().zoomToExtent(trajLayer.getExtent().extend(2, 2))<br /> <br /> #---- Output figure<br /> print 'Output figure...'<br /> mapLayout.exportToPicture(os.path.join(figDir, 'traj_test.png'))<br /> <br /> #---- Set mapframe<br /> mapFrame.setGridXDelt(5)<br /> mapFrame.setGridYDelt(5)<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (650, 680))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> <br /> print 'Finish...'

运行结果：

tmzly2008

正需要学习了，老师今天怎么有空发教程的啊

1008611

谢啦 哥们                                         

caflj

好东西，谢谢。

Nancy_fq

向王老师致敬！

juditxh901016

谢谢啦！{:eb502:}{:eb502:}

djselina

老师您好 我想问一下 您示范的例图，如果想要把气压廓线换成其他物理量比如比湿或者位温 可以忙？请问要怎么改？

MeteoInfo

djselina 发表于 2016-4-10 21:09
老师您好 我想问一下 您示范的例图，如果想要把气压廓线换成其他物理量比如比湿或者位温 可以忙？请问要怎 ...

参考这里：http://www.meteothinker.com/docs ... /data_tutorial.html

djselina

谢谢老师及时的回复 命令中dataset = dataInfo.getXYDataset(12)    #Pressure is in 12th column这一句我理解为气压的数据是在第12列，请问应该怎样查看其他列的数据分别代表什么物理意义呢？您给的链接我学习了 可是水平有限并没能找到答案。

MeteoInfo

djselina 发表于 2016-4-11 11:34
谢谢老师及时的回复 命令中dataset = dataInfo.getXYDataset(12)    #Pressure is in 12th column这一句我 ...

这个你需要看HYSPLIT模式的用户手册里轨迹节点文件格式的说明。