MeteoInfo脚本教程（十一）：站点风场图

MeteoInfo

风场需要两个变量来表示：U/V分量或者风向/风速。MICAPS第1类格式数据中用的是后者。DrawMeteoData类中有生成站点风场图层的静态方法：
    DrawMeteoData.createSTBarbLayer(StationData windDir, StationData windSpeed, String layerName, Bollean isUV)    ->    生成站点风场Barb图层
    DrawMeteoData.createSTVectorLayer(StationData windDir, StationData windSpeed, String layerName, Bollean isUV)    ->    生成站点风场矢量图层

脚本示例程序：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-6
   
4. # Purpose: Get and plot station wind 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.black)
    
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 wind direction/speed station data
    
44. windDir = mdi.getStationData('WindDirection')
    
45. windSpeed = mdi.getStationData('WindSpeed')
    
46. 
    
47. #---- Create barb and vector wind layers
    
48. bLayer = DrawMeteoData.createSTBarbLayer(windDir, windSpeed, 'WindBarb_Point', False)
    
49. vLayer = DrawMeteoData.createSTVectorLayer(windDir, windSpeed, 'WindVector_Point', False)
    
50. 
    
51. #---- Add layers
    
52. mapFrame.addLayer(bLayer)
    
53. mapFrame.addLayer(vLayer)
    
54. 
    
55. #---- Add title
    
56. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
57. 
    
58. #---- Set layout map
    
59. print 'Set layout map...'
    
60. mapLayout.getActiveLayoutMap().setWidth(580)
    
61. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
62. 
    
63. #---- Set mapframe
    
64. mapFrame.setGridXDelt(10)
    
65. mapFrame.setGridYDelt(10)
    
66. 
    
67. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
68. frame.add(mapLayout)
    
69. frame.visible = True
    
70. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-6<br /> # Purpose: Get and plot station wind 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.black)<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 wind direction/speed station data<br /> windDir = mdi.getStationData('WindDirection')<br /> windSpeed = mdi.getStationData('WindSpeed')<br /> <br /> #---- Create barb and vector wind layers<br /> bLayer = DrawMeteoData.createSTBarbLayer(windDir, windSpeed, 'WindBarb_Point', False)<br /> vLayer = DrawMeteoData.createSTVectorLayer(windDir, windSpeed, 'WindVector_Point', False)<br /> <br /> #---- Add layers<br /> mapFrame.addLayer(bLayer)<br /> mapFrame.addLayer(vLayer)<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 /> frame = JFrame('MeteoInfo Script Sample', size = (800, 600))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

yuxingzhan

在最新版MeteoInfoLab 运行NameError: name'DrawMeteoData'is not defined，什么原因？

MeteoInfo

yuxingzhan 发表于 2020-3-9 20:52
在最新版MeteoInfoLab 运行NameError: name'DrawMeteoData'is not defined，什么原因？

这个帖子很旧了，建议你参考这里的脚本：http://bbs.06climate.com/forum.p ... &extra=page%3D1