MeteoInfo脚本教程（九）：站点数据

MeteoInfo

之前提到过站点数据是StationData类，其中包含了站点号、站点坐标和站点观测值。还有一类站点数据是StationInfoData类，可以包含更多的站点数据信息，包括非数字信息。StationData和StationInfoData类的典型应用是创建站点分布图，分别由DrawMeteoData类的下面方法创建：

    DrawMeteoData.createSTPointLayer(StationData stdata, String layerName, String fieldName)
    DrawMeteoData.createSTInfoLayer(StationInfoData sidata, String layerName)

脚本示例如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-1
   
4. # Purpose: Get and plot station data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. #---- Import classes
   
8. print 'Import classes...'
   
9. from org.meteoinfo.layout import MapLayout
   
10. from org.meteoinfo.data.mapdata import MapDataManage
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
12. from org.meteoinfo.layout import LegendStyles
    
13. from org.meteoinfo.global import Extent
    
14. import os
    
15. from java.awt import Color
    
16. from javax.swing import JFrame
    
17. 
    
18. #---- Set directories
    
19. print 'Set directories...'
    
20. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
21. dataDir = os.path.join(baseDir, 'sample')
    
22. mapDir = os.path.join(baseDir, 'map')
    
23. 
    
24. #---- Create MapLayout object
    
25. mapLayout = MapLayout()
    
26. mapFrame = mapLayout.getActiveMapFrame()
    
27. 
    
28. #---- Load country layer
    
29. print 'Load country layer...'
    
30. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
31. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
32. lb.setDrawFill(False)
    
33. lb.setOutlineColor(Color.blue)
    
34. mapFrame.addLayer(countryLayer)
    
35. 
    
36. #---- Create MeteoDataInfo object
    
37. mdi = MeteoDataInfo()
    
38. 
    
39. #---- Open a MICAPS data file
    
40. fn = os.path.join(dataDir, 'MICAPS/10101414.000')
    
41. mdi.openMICAPSData(fn)
    
42. 
    
43. #---- Get station data
    
44. sdata = mdi.getStationData('Visibility')
    
45. sidata = mdi.getStationInfoData()
    
46. 
    
47. #---- Create station point layer from the grid data
    
48. layer = DrawMeteoData.createSTPointLayer(sdata, 'Visibility_STPoint', 'Visibility')
    
49. 
    
50. #---- Add layer
    
51. mapFrame.addLayer(layer)
    
52. 
    
53. #--- Move layer to bottom
    
54. mapFrame.moveLayer(layer, 0)
    
55. 
    
56. #---- Add title
    
57. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
58. 
    
59. #---- Set layout map
    
60. print 'Set layout map...'
    
61. mapLayout.getActiveLayoutMap().setWidth(580)
    
62. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
63. 
    
64. #---- Set mapframe
    
65. mapFrame.setGridXDelt(10)
    
66. mapFrame.setGridYDelt(10)
    
67. 
    
68. #---- Add legend
    
69. legend = mapLayout.addLegend(620, 150)
    
70. legend.setLegendStyle(LegendStyles.Normal)
    
71. legend.setLegendLayer(layer)
    
72. legend.setTitle('Visibility')
    
73. 
    
74. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
75. frame.add(mapLayout)
    
76. frame.visible = True
    
77. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-1<br /> # Purpose: Get and plot station data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.layout import MapLayout<br /> from org.meteoinfo.data.mapdata import MapDataManage<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.layout import LegendStyles<br /> from org.meteoinfo.global import Extent<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 /> <br /> #---- Create MapLayout object<br /> mapLayout = MapLayout()<br /> mapFrame = mapLayout.getActiveMapFrame()<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.setDrawFill(False)<br /> lb.setOutlineColor(Color.blue)<br /> mapFrame.addLayer(countryLayer)<br /> <br /> #---- Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #---- Open a MICAPS data file<br /> fn = os.path.join(dataDir, 'MICAPS/10101414.000')<br /> mdi.openMICAPSData(fn)<br /> <br /> #---- Get station data<br /> sdata = mdi.getStationData('Visibility')<br /> sidata = mdi.getStationInfoData()<br /> <br /> #---- Create station point layer from the grid data<br /> layer = DrawMeteoData.createSTPointLayer(sdata, 'Visibility_STPoint', 'Visibility')<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(layer)<br /> <br /> #--- Move layer to bottom<br /> mapFrame.moveLayer(layer, 0)<br /> <br /> #---- Add title<br /> title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)<br /> <br /> #---- Set layout map<br /> print 'Set layout map...'<br /> mapLayout.getActiveLayoutMap().setWidth(580)<br /> mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))<br /> <br /> #---- Set mapframe<br /> mapFrame.setGridXDelt(10)<br /> mapFrame.setGridYDelt(10)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(620, 150)<br /> legend.setLegendStyle(LegendStyles.Normal)<br /> legend.setLegendLayer(layer)<br /> legend.setTitle('Visibility')<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (800, 600))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-1
   
4. # Purpose: Get and plot station info data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. #---- Import classes
   
8. print 'Import classes...'
   
9. from org.meteoinfo.layout import MapLayout
   
10. from org.meteoinfo.data.mapdata import MapDataManage
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
12. from org.meteoinfo.layout import LegendStyles
    
13. from org.meteoinfo.global import Extent
    
14. import os
    
15. from java.awt import Color
    
16. from java.awt import Font
    
17. from javax.swing import JFrame
    
18. 
    
19. #---- Set directories
    
20. print 'Set directories...'
    
21. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
22. dataDir = os.path.join(baseDir, 'sample')
    
23. mapDir = os.path.join(baseDir, 'map')
    
24. 
    
25. #---- Create MapLayout object
    
26. mapLayout = MapLayout()
    
27. mapFrame = mapLayout.getActiveMapFrame()
    
28. 
    
29. #---- Load country layer
    
30. print 'Load country layer...'
    
31. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
32. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
33. lb.setDrawFill(False)
    
34. lb.setOutlineColor(Color.blue)
    
35. mapFrame.addLayer(countryLayer)
    
36. 
    
37. #---- Create MeteoDataInfo object
    
38. mdi = MeteoDataInfo()
    
39. 
    
40. #---- Open a MICAPS data file
    
41. fn = os.path.join(dataDir, 'MICAPS/10101414.000')
    
42. mdi.openMICAPSData(fn)
    
43. 
    
44. #---- Get station info data
    
45. sidata = mdi.getStationInfoData()
    
46. 
    
47. #---- Create station info layer from the station info data
    
48. layer = DrawMeteoData.createSTInfoLayer(sidata, 'StationInfo')
    
49. labset = layer.getLabelSet()
    
50. labset.setFieldName('CloudCover')
    
51. font = Font('Arial', Font.PLAIN, 14)
    
52. labset.setLabelFont(font)
    
53. labset.setYOffset(15)
    
54. layer.addLabels()
    
55. 
    
56. #---- Add layer
    
57. mapFrame.addLayer(layer)
    
58. 
    
59. #--- Move layer to bottom
    
60. mapFrame.moveLayer(layer, 0)
    
61. 
    
62. #---- Add title
    
63. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
64. 
    
65. #---- Set layout map
    
66. print 'Set layout map...'
    
67. mapLayout.getActiveLayoutMap().setWidth(580)
    
68. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
69. 
    
70. #---- Set mapframe
    
71. mapFrame.setGridXDelt(10)
    
72. mapFrame.setGridYDelt(10)
    
73. 
    
74. #---- Add legend
    
75. legend = mapLayout.addLegend(620, 150)
    
76. legend.setLegendStyle(LegendStyles.Normal)
    
77. legend.setLegendLayer(layer)
    
78. legend.setTitle('Station')
    
79. 
    
80. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
81. frame.add(mapLayout)
    
82. frame.visible = True
    
83. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-1<br /> # Purpose: Get and plot station info data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.layout import MapLayout<br /> from org.meteoinfo.data.mapdata import MapDataManage<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.layout import LegendStyles<br /> from org.meteoinfo.global import Extent<br /> import os<br /> from java.awt import Color<br /> from java.awt import Font<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 /> <br /> #---- Create MapLayout object<br /> mapLayout = MapLayout()<br /> mapFrame = mapLayout.getActiveMapFrame()<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.setDrawFill(False)<br /> lb.setOutlineColor(Color.blue)<br /> mapFrame.addLayer(countryLayer)<br /> <br /> #---- Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #---- Open a MICAPS data file<br /> fn = os.path.join(dataDir, 'MICAPS/10101414.000')<br /> mdi.openMICAPSData(fn)<br /> <br /> #---- Get station info data<br /> sidata = mdi.getStationInfoData()<br /> <br /> #---- Create station info layer from the station info data<br /> layer = DrawMeteoData.createSTInfoLayer(sidata, 'StationInfo')<br /> labset = layer.getLabelSet()<br /> labset.setFieldName('CloudCover')<br /> font = Font('Arial', Font.PLAIN, 14)<br /> labset.setLabelFont(font)<br /> labset.setYOffset(15)<br /> layer.addLabels()<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(layer)<br /> <br /> #--- Move layer to bottom<br /> mapFrame.moveLayer(layer, 0)<br /> <br /> #---- Add title<br /> title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)<br /> <br /> #---- Set layout map<br /> print 'Set layout map...'<br /> mapLayout.getActiveLayoutMap().setWidth(580)<br /> mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))<br /> <br /> #---- Set mapframe<br /> mapFrame.setGridXDelt(10)<br /> mapFrame.setGridYDelt(10)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(620, 150)<br /> legend.setLegendStyle(LegendStyles.Normal)<br /> legend.setLegendLayer(layer)<br /> legend.setTitle('Station')<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (800, 600))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

MeteoInfo

小蜡图图 发表于 2015-2-4 11:23
要画的数据导入应该修改的是哪一行呢？希望指教~

fn = os.path.join(dataDir, 'MICAPS/10101414.000')

总之需要将fn设为你计算机中micaps 1文件的全路径。

hugh39

~~~~~zhege keyi you a ~

rceclx

速度！快成一个系列了，赞王老师！

zihan6688

继续给我们带来新知识

sa13226139

感谢楼主！

小蜡图图

想问问老师，新手刚刚接触，直接复制的脚本出错：Import classes...
Set directories...
Load country layer...
Traceback (most recent call last):
  File "<string>", line 31, in <module>
AttributeError: 'NoneType' object has no attribute 'getLegendScheme'
想画第一幅图，想问问baseDir= 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'的路径里存放的是什么？问题略二，求解

MeteoInfo

小蜡图图 发表于 2015-2-4 11:01
想问问老师，新手刚刚接触，直接复制的脚本出错：Import classes...
Set directories...
Load country la ...

需要修改baseDir为你计算机上MeteoInfo的路径。

小蜡图图

MeteoInfo 发表于 2015-2-4 11:04
需要修改baseDir为你计算机上MeteoInfo的路径。

要画的数据导入应该修改的是哪一行呢？希望指教~

星雨

王老师，Visibility怎么会出现70这样的数值？是Micaps里数据读取有问题吗？