MeteoInfo脚本教程（七）：格点数据

MeteoInfo

MeteoInfo里有两个最重要的数据类：GridData（格点数据）和StationData（站点数据）。气象模式输出的数据通常是格点数据，对某一变量而言可以有4维数据（空间三维和时间维）。从MeteoDataInfo对象中获取GridData的方法是getGridData(String variableName)，可以通过setTimeIndex和setLevelIndex方法设置时间和高度维，从而获取水平的2维GridData数据。其实也可以通过其它设置获取其它的2维或1维格点数据，比如剖面格点数据，这个以后再讲。

有了GridData数据后可以用DrawMeteoData类的一些静态方法来创建相应的图层达到绘制数据的目的。从GridData数据创建图层的方法主要如下：
    DrawMeteoData.createContourLayer(GridData data,String layerName,String fieldName, Bollean isSmooth)    ->    生成等值线图层（VectorLayer，Polyline）
    DrawMeteoData.createShadedLayer(GridData data,String layerName,String fieldName, Bollean isSmooth)    ->    生成等值线填色图层（VectorLayer，Polygon）
    DrawMeteoData.createRasterLayer(GridData data,String layerName)    ->    生成栅格图层（RasterLayer）
    DrawMeteoData.createGridFillLayer(GridData data,String layerName,String fieldName)    ->    生成格点填色图层（VectorLayer，Polygon）
    DrawMeteoData.createGridPointLayer(GridData data,String layerName,String fieldName)    ->    生成格点点图层（VectorLayer，Point）


下面是一个Demo脚本：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-11-30
   
4. # Purpose: Get and plot 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.blue)
    
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 grid data
    
44. mdi.setTimeIndex(2)
    
45. mdi.setLevelIndex(3)
    
46. gdata = mdi.getGridData('Z')
    
47. gdata.extendToGlobal()
    
48. 
    
49. #---- Create shaded layer from the grid data
    
50. layer = DrawMeteoData.createContourLayer(gdata, 'Z_contour', 'Z', True)
    
51. #layer = DrawMeteoData.createShadedLayer(gdata, 'Z_shaded', 'Z', True)
    
52. #layer = DrawMeteoData.createRasterLayer(gdata, 'Z_raster')
    
53. #layer = DrawMeteoData.createGridFillLayer(gdata, 'Z_gridfill', 'Z')
    
54. #layer = DrawMeteoData.createGridPointLayer(gdata, 'Z_gridpoint', 'Z')
    
55. 
    
56. #---- Add layer
    
57. mapFrame.addLayer(layer)
    
58. 
    
59. #--- Move pressure 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. #---- Zoom layout map
    
66. print 'Zoom layout map...'
    
67. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(0, 360, -90, 90))
    
68. 
    
69. #---- Set mapframe
    
70. mapFrame.setGridXDelt(30)
    
71. mapFrame.setGridYDelt(30)
    
72. 
    
73. #---- Add legend
    
74. legend = mapLayout.addLegend(150, 440)
    
75. legend.setLegendStyle(LegendStyles.Bar_Horizontal)
    
76. legend.setLegendLayer(layer)
    
77. 
    
78. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
79. frame.add(mapLayout)
    
80. frame.visible = True
    
81. print 'Finished!'
    

</font>#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-11-30<br /> # Purpose: Get and plot 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.blue)<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 grid data<br /> mdi.setTimeIndex(2)<br /> mdi.setLevelIndex(3)<br /> gdata = mdi.getGridData('Z')<br /> gdata.extendToGlobal()<br /> <br /> #---- Create shaded layer from the grid data<br /> layer = DrawMeteoData.createContourLayer(gdata, 'Z_contour', 'Z', True)<br /> #layer = DrawMeteoData.createShadedLayer(gdata, 'Z_shaded', 'Z', True)<br /> #layer = DrawMeteoData.createRasterLayer(gdata, 'Z_raster')<br /> #layer = DrawMeteoData.createGridFillLayer(gdata, 'Z_gridfill', 'Z')<br /> #layer = DrawMeteoData.createGridPointLayer(gdata, 'Z_gridpoint', 'Z')<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(0, 360, -90, 90))<br /> <br /> #---- Set mapframe<br /> mapFrame.setGridXDelt(30)<br /> mapFrame.setGridYDelt(30)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(150, 440)<br /> legend.setLegendStyle(LegendStyles.Bar_Horizontal)<br /> legend.setLegendLayer(layer)<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (800, 600))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'<br /> <font style=\"font-size:14px\">

兰溪之水

色标会不会少一个标值？

MeteoInfo

兰溪之水 发表于 2014-12-1 19:09
色标会不会少一个标值？

不少，对于等值线来说，每条线是固定值而不是区间值。

兰溪之水

MeteoInfo 发表于 2014-12-1 19:20
不少，对于等值线来说，每条线是固定值而不是区间值。

我错了，等值线。。。

zhangbote

王老师，想请教您一下。java二次开发，画降水色斑图，调用clip方法，layer.clip(jjjMap);，layer是色斑图的图层，jjjMap是京津冀地图的图层，报了错TopologyException: found non-noded intersection between LINESTRING ( 116.82000012309558 42.679323430660034, 116.82000011879595 42.67932341079954 ) and LINESTRING ( 116.82 42.67932286206939, 116.82000020243557 42.67932379713915 ) [ (116.82000011908272, 42.67932341212417, NaN) ]。没有交点。这个是怎么回事呢。