？：风场区域的选择

solitaire

选择出风速大于某一阈值的区域  用脚本怎么实现
程序的数据用的micaps数据的UV
效果如图

MeteoInfo

想要什么样的最终结果？

solitaire

MeteoInfo 发表于 2015-4-13 19:25
想要什么样的最终结果？

就图片的样子 比如选择出风速大于20m/s的区域 再把这个区域的流线画出来。现在的困惑就是不知道怎么根据uv数据选择出某一阈值的区域

MeteoInfo

solitaire 发表于 2015-4-13 19:32
就图片的样子 比如选择出风速大于20m/s的区域 再把这个区域的流线画出来。现在的困惑就是不知道怎么根据u ...

示例脚本程序，注意看里面的注释：

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 import DataMath
    
11. from org.meteoinfo.data.mapdata import MapDataManage
    
12. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
13. from org.meteoinfo.layout import LegendStyles
    
14. from org.meteoinfo.global import Extent
    
15. import os
    
16. from java.awt import Color
    
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.black)
    
35. mapFrame.addLayer(countryLayer)
    
36. 
    
37. #---- Create MeteoDataInfo object
    
38. mdi = MeteoDataInfo()
    
39. 
    
40. #---- Open a GrADS data file
    
41. fn = os.path.join(dataDir, 'GrADS/model.ctl')
    
42. mdi.openGrADSData(fn)
    
43. 
    
44. #---- Get U/V grid data
    
45. mdi.setTimeIndex(2)
    
46. mdi.setLevelIndex(3)
    
47. udata = mdi.getGridData('U')
    
48. vdata = mdi.getGridData('V')
    
49. #---- Calculate wind speed from U/V data
    
50. speed = DataMath.magnitude(udata, vdata)
    
51. #---- Set the grid value to -1 if the speed is less than 10
    
52. speed.replaceValue(10, -1, False)
    
53. #---- Maskout U data, Only the grid with speed bigger than 10 has valid value
    
54. udata = udata.maskout(speed)
    
55. 
    
56. #---- Create wind layer from the U/V grid data
    
57. #layer = DrawMeteoData.createGridVectorLayer(udata, vdata,  'UV_Vector', True)
    
58. layer = DrawMeteoData.createGridBarbLayer(udata, vdata,  'UV_Barb', True)
    
59. #layer = DrawMeteoData.createStreamlineLayer(udata, vdata, 'Z_Streamline', True)
    
60. 
    
61. #---- Add layer
    
62. mapFrame.addLayer(layer)
    
63. 
    
64. #--- Move pressure layer to bottom
    
65. mapFrame.moveLayer(layer, 0)
    
66. 
    
67. #---- Add title
    
68. title = mapLayout.addText('MeteoInfo script demo', 350, 30, 'Arial', 16)
    
69. 
    
70. #---- Zoom layout map
    
71. print 'Zoom layout map...'
    
72. mapLayout.getActiveLayoutMap().zoomToExtentLonLatEx(Extent(0, 360, -90, 90))
    
73. 
    
74. #---- Set mapframe
    
75. mapFrame.setGridXDelt(10)
    
76. mapFrame.setGridYDelt(10)
    
77. 
    
78. frame = JFrame('MeteoInfo Script Sample', size = (800, 600))
    
79. frame.add(mapLayout)
    
80. frame.visible = True
    
81. 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 import DataMath<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 /> #---- Calculate wind speed from U/V data<br /> speed = DataMath.magnitude(udata, vdata)<br /> #---- Set the grid value to -1 if the speed is less than 10<br /> speed.replaceValue(10, -1, False)<br /> #---- Maskout U data, Only the grid with speed bigger than 10 has valid value<br /> udata = udata.maskout(speed)<br /> <br /> #---- Create wind layer from the U/V 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(0, 360, -90, 90))<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!'

结果：

solitaire

MeteoInfo 发表于 2015-4-14 23:34
示例脚本程序，注意看里面的注释：

多谢老师  看懂了  有没有个地方或者书籍可以系统的学习一下这些？

MeteoInfo

solitaire 发表于 2015-4-15 15:34
多谢老师  看懂了  有没有个地方或者书籍可以系统的学习一下这些？

MeteoInfo的网站