MeteoInfo脚本教程（八）：U/V格点数据

MeteoInfo

风场的表达需要两个变量：U/V或者风速/风向。模式通常会输出风场的U/V分量，可以绘制风场矢量图、风向杆图和流场图，分别用DrawMeteoData的下面方法：

    DrawMeteoData.createGridVectorLayer(GridData udata, GridData vdata,  String layerName, Boolean isUV)
    DrawMeteoData.createGridBarbLayer(GridData udata, GridData vdata,  String layerName, Boolean isUV)
    DrawMeteoData.createStreamlineLayer(GridData udata, GridData vdata,  String layerName, Boolean isUV)


示例脚本如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-1
   
4. # Purpose: Read and plot U/V grid 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 GrADS data file
    
40. fn = os.path.join(dataDir, 'GrADS/model.ctl')
    
41. mdi.openGrADSData(fn)
    
42. 
    
43. #---- Get U/V grid data
    
44. mdi.setTimeIndex(2)
    
45. mdi.setLevelIndex(3)
    
46. udata = mdi.getGridData('U')
    
47. vdata = mdi.getGridData('V')
    
48. 
    
49. #---- Create shaded layer from the grid data
    
50. layer = DrawMeteoData.createGridVectorLayer(udata, vdata,  'UV_Vector', True)
    
51. #layer = DrawMeteoData.createGridBarbLayer(udata, vdata,  'UV_Barb', True)
    
52. #layer = DrawMeteoData.createStreamlineLayer(udata, vdata, 'Z_Streamline', True)
    
53. 
    
54. #---- Add layer
    
55. mapFrame.addLayer(layer)
    
56. 
    
57. #--- Move pressure layer to bottom
    
58. mapFrame.moveLayer(layer, 0)
    
59. 
    
60. #---- Add title
    
61. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
62. 
    
63. #---- Zoom layout map
    
64. print 'Zoom layout map...'
    
65. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
66. 
    
67. #---- Set mapframe
    
68. mapFrame.setGridXDelt(10)
    
69. mapFrame.setGridYDelt(10)
    
70. 
    
71. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
72. frame.add(mapLayout)
    
73. frame.visible = True
    
74. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-1<br /> # Purpose: Read and plot U/V grid 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 GrADS data file<br /> fn = os.path.join(dataDir, 'GrADS/model.ctl')<br /> mdi.openGrADSData(fn)<br /> <br /> #---- Get U/V grid data<br /> mdi.setTimeIndex(2)<br /> mdi.setLevelIndex(3)<br /> udata = mdi.getGridData('U')<br /> vdata = mdi.getGridData('V')<br /> <br /> #---- Create shaded layer from the grid data<br /> layer = DrawMeteoData.createGridVectorLayer(udata, vdata,  'UV_Vector', True)<br /> #layer = DrawMeteoData.createGridBarbLayer(udata, vdata,  'UV_Barb', True)<br /> #layer = DrawMeteoData.createStreamlineLayer(udata, vdata, 'Z_Streamline', True)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(layer)<br /> <br /> #--- Move pressure 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 /> #---- Zoom layout map<br /> print 'Zoom layout map...'<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!'

MapLayout的addWindArrow(int left, int top)方法可以添加风矢量箭头图例：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-1
   
4. # Purpose: Read and plot U/V grid 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 GrADS data file
    
40. fn = os.path.join(dataDir, 'GrADS/model.ctl')
    
41. mdi.openGrADSData(fn)
    
42. 
    
43. #---- Get U/V grid data
    
44. mdi.setTimeIndex(2)
    
45. mdi.setLevelIndex(3)
    
46. udata = mdi.getGridData('U')
    
47. vdata = mdi.getGridData('V')
    
48. 
    
49. #---- Create wind vector layer from the U/V grid data
    
50. layer = DrawMeteoData.createGridVectorLayer(udata, vdata,  'UV_Vector', True)
    
51. 
    
52. #---- Add layer
    
53. mapFrame.addLayer(layer)
    
54. 
    
55. #--- Move pressure layer to bottom
    
56. mapFrame.moveLayer(layer, 0)
    
57. 
    
58. #---- Add title
    
59. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
60. 
    
61. #---- Add wind arrow
    
62. windArrow = mapLayout.addWindArrow(660, 420)
    
63. 
    
64. #---- Zoom layout map
    
65. print 'Zoom layout map...'
    
66. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
67. 
    
68. #---- Set mapframe
    
69. mapFrame.setGridXDelt(10)
    
70. mapFrame.setGridYDelt(10)
    
71. 
    
72. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
73. frame.add(mapLayout)
    
74. frame.visible = True
    
75. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-1<br /> # Purpose: Read and plot U/V grid 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 GrADS data file<br /> fn = os.path.join(dataDir, 'GrADS/model.ctl')<br /> mdi.openGrADSData(fn)<br /> <br /> #---- Get U/V grid data<br /> mdi.setTimeIndex(2)<br /> mdi.setLevelIndex(3)<br /> udata = mdi.getGridData('U')<br /> vdata = mdi.getGridData('V')<br /> <br /> #---- Create wind vector layer from the U/V grid data<br /> layer = DrawMeteoData.createGridVectorLayer(udata, vdata,  'UV_Vector', True)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(layer)<br /> <br /> #--- Move pressure 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 /> #---- Add wind arrow<br /> windArrow = mapLayout.addWindArrow(660, 420)<br /> <br /> #---- Zoom layout map<br /> print 'Zoom layout map...'<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!'