MeteoInfo脚本教程（十二）：Maskout

MeteoInfo

如果希望专业内容的图形部分只出现在某个区域内，即屏蔽区域外的图形，可以用MeteoInfo中的Maskout功能。图形最终都是显示在MapView中，MapView类中有个变量是MaskOut类的对象，用MapView类的getMaskOut()方法获取。MaskOut类的setMask(Boolean isMask)方法设置了是否启用MaskOut功能，SetMaskLayer(String layerName)方法设置了作为Maskout使用的图层，图层的图元类型必须Polygon，最好不要太复杂（会影响图形绘制速度）。每一个图层都有一个setMaskout(Boolean isMask)方法来设置该图层是否被Maskout，这样可以实现只Maskout某些特定图层。

示例脚本如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-12-6
   
4. # Purpose: Maskout layer
   
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.legend import LegendManage, LegendType, GridLabelPosition
    
13. from org.meteoinfo.layout import LegendStyles
    
14. from org.meteoinfo.shape import ShapeTypes
    
15. from org.meteoinfo.global import Extent
    
16. import os.path
    
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/GrADS')
    
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. mapView = mapFrame.getMapView()
    
31. 
    
32. #---- Load country layer
    
33. print 'Load country layer...'
    
34. bou2Layer = MapDataManage.loadLayer(os.path.join(mapDir, 'bou2_4p.shp'))
    
35. lb = bou2Layer.getLegendScheme().getLegendBreaks().get(0)
    
36. lb.setDrawFill(False)
    
37. lb.setOutlineColor(Color.gray)
    
38. mapFrame.addLayer(bou2Layer)
    
39. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
40. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
41. lb.setDrawFill(False)
    
42. lb.setOutlineColor(Color.blue)
    
43. mapFrame.addLayer(countryLayer)
    
44. chinaLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'china.shp'))
    
45. chinaLayer.setVisible(False)
    
46. mapFrame.addLayer(chinaLayer)
    
47. 
    
48. #---- Open GrADS data
    
49. print 'Open GrADS data...'
    
50. mdi = MeteoDataInfo()
    
51. mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))
    
52. 
    
53. #---- Set time index
    
54. mdi.setTimeIndex(2)
    
55. 
    
56. #---- Get pressure grid data
    
57. gdata = mdi.getGridData('PS')
    
58. gdata.extendToGlobal()
    
59. 
    
60. #---- Create pressure shaded layer
    
61. print 'Create pressure shaded layer...'
    
62. pressLayer = DrawMeteoData.createShadedLayer(gdata, 'Pressure', 'PS', True)
    
63. pressLayer.setMaskout(True)
    
64. 
    
65. #---- Add layer
    
66. mapFrame.addLayer(pressLayer)
    
67. 
    
68. #--- Move pressure layer to bottom
    
69. mapFrame.moveLayer(pressLayer, 0)
    
70. 
    
71. #---- Add title
    
72. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
73. 
    
74. #---- Zoom layout map
    
75. print 'Zoom layout map...'
    
76. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(70, 140, 15, 55))
    
77. 
    
78. #---- Set maskout
    
79. mapView.getMaskOut().setMask(True)
    
80. mapView.getMaskOut().setMaskLayer(chinaLayer.getLayerName())
    
81. 
    
82. #---- Set mapframe
    
83. mapFrame.setGridXDelt(10)
    
84. mapFrame.setGridYDelt(10)
    
85. 
    
86. #---- Add legend
    
87. legend = mapLayout.addLegend(660, 120)
    
88. legend.setLegendStyle(LegendStyles.Bar_Vertical)
    
89. legend.setLegendLayer(pressLayer)
    
90. 
    
91. #---- Output figure
    
92. print 'Output figure...'
    
93. mapLayout.paintGraphics()
    
94. #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.ps'))
    
95. 
    
96. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
97. frame.add(mapLayout)
    
98. frame.visible = True
    
99. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-6<br /> # Purpose: Maskout layer<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.legend import LegendManage, LegendType, GridLabelPosition<br /> from org.meteoinfo.layout import LegendStyles<br /> from org.meteoinfo.shape import ShapeTypes<br /> from org.meteoinfo.global import Extent<br /> import os.path<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/GrADS')<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 /> mapView = mapFrame.getMapView()<br /> <br /> #---- Load country layer<br /> print 'Load country layer...'<br /> bou2Layer = MapDataManage.loadLayer(os.path.join(mapDir, 'bou2_4p.shp'))<br /> lb = bou2Layer.getLegendScheme().getLegendBreaks().get(0)<br /> lb.setDrawFill(False)<br /> lb.setOutlineColor(Color.gray)<br /> mapFrame.addLayer(bou2Layer)<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 /> chinaLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'china.shp'))<br /> chinaLayer.setVisible(False)<br /> mapFrame.addLayer(chinaLayer)<br /> <br /> #---- Open GrADS data<br /> print 'Open GrADS data...'<br /> mdi = MeteoDataInfo()<br /> mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))<br /> <br /> #---- Set time index<br /> mdi.setTimeIndex(2)<br /> <br /> #---- Get pressure grid data<br /> gdata = mdi.getGridData('PS')<br /> gdata.extendToGlobal()<br /> <br /> #---- Create pressure shaded layer<br /> print 'Create pressure shaded layer...'<br /> pressLayer = DrawMeteoData.createShadedLayer(gdata, 'Pressure', 'PS', True)<br /> pressLayer.setMaskout(True)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(pressLayer)<br /> <br /> #--- Move pressure layer to bottom<br /> mapFrame.moveLayer(pressLayer, 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 maskout<br /> mapView.getMaskOut().setMask(True)<br /> mapView.getMaskOut().setMaskLayer(chinaLayer.getLayerName())<br /> <br /> #---- Set mapframe<br /> mapFrame.setGridXDelt(10)<br /> mapFrame.setGridYDelt(10)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(660, 120)<br /> legend.setLegendStyle(LegendStyles.Bar_Vertical)<br /> legend.setLegendLayer(pressLayer)<br /> <br /> #---- Output figure<br /> print 'Output figure...'<br /> mapLayout.paintGraphics()<br /> #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.ps'))<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (800, 600))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

未启用Maskout：


启用Maskout:

mofangbao

楼主要出书的节奏啊

风儿飘飘

真希望王老师出一本书，这样就能全面学习了

wangxue910402

王老师，我想问一下，我有全中国的数据，想只做云南省的等值线图，假设用maskout选取出来云南省的区域，那能在云南省成一个闭合的等值线么？还是只是截取了云南省区域的的。因为我原来用grads做的时候只是直接截取了这个区域的。不知道我表达清楚了没？新手，有点表达不清。

MeteoInfo

wangxue910402 发表于 2015-11-16 11:35
王老师，我想问一下，我有全中国的数据，想只做云南省的等值线图，假设用maskout选取出来云南省的区域，那 ...

MeteoInfo里maskout分两种方式，maskout数据（和GrADS一样）和maskout图形。

maskout数据的例子：http://bbs.06climate.com/forum.p ... &extra=page%3D1

maskout图形的例子：http://bbs.06climate.com/forum.p ... &extra=page%3D1

冬青老先生

老师您好，请问如果我只想要海区的内容，应当如何去掉陆地部分呢？

zq01

请问he5做的图也可以屏蔽吗

MeteoInfo

zq01 发表于 2016-7-28 11:27
请问he5做的图也可以屏蔽吗

跟什么数据做的图没关系

黑色青春

脚本不错

突破的渴望

老师您好，脚本第25行这个test的格式是怎么样的呢？ 谢谢老师~~