MeteoInfo脚本示例 - 格点数据计算

MeteoInfo

格点数据（GridData）是MeteoInfo数据结构中很重要的一个类，这个Jython脚本程序（需要用MeteoInfo Java版）演示了两个文件中同一变量格点数据的加、减、乘、除几个基本运算。需要在MeteoInfo Java版桌面程序中运行此脚本（脚本代码中有miapp变量的都需要在桌面程序中运行，因为miapp就是当前运行的桌面程序主界面对象）。Java不支持运算符重载，所以只能用add, sub, mul, div方法来实现相应的计算功能。

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-9-12
   
4. # Purpose: Add, sub, mul, div functions
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata import DrawMeteoData
   
9. from org.meteoinfo.legend import LegendManage
   
10. from org.meteoinfo.legend import LegendType
    
11. from org.meteoinfo.shape import ShapeTypes
    
12. import os.path
    
13. import jarray
    
14. 
    
15. #Set data directory
    
16. dataDir = 'D:/Temp/nc'
    
17. 
    
18. #Create MeteoDataInfo object
    
19. mdi1 = MeteoDataInfo()
    
20. mdi2 = MeteoDataInfo()
    
21. 
    
22. #Read data files
    
23. fn1 = os.path.join(dataDir, 'wrfout_d01_2008-01-01_12_00_00.nc')
    
24. fn2 = os.path.join(dataDir, 'wrfout_d01_2008-01-02_12_00_00.nc')
    
25. print 'File 1: ' + fn1
    
26. print 'File 2: ' + fn2
    
27. mdi1.openNetCDFData(fn1)
    
28. mdi2.openNetCDFData(fn2)
    
29. #Get grid data
    
30. print 'Get grid data...'
    
31. gData1 = mdi1.getGridData('Q2')
    
32. gData2 = mdi2.getGridData('Q2')
    
33. #Add two grid data
    
34. print 'Add grid data...'
    
35. addData = gData1.add(gData2)
    
36. #Subtract two grid data
    
37. print 'Subtract grid data...'
    
38. subData = gData1.sub(gData2)
    
39. #Multiple two grid data
    
40. print 'Multiple grid data...'
    
41. mulData = gData1.mul(gData2)
    
42. #Divide two grid data
    
43. print 'Divide grid data...'
    
44. divData = gData1.div(gData2)
    
45. #Create legend scheme
    
46. addLS = LegendManage.createLegendSchemeFromGridData(addData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
47. subLS = LegendManage.createLegendSchemeFromGridData(subData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
48. mulLS = LegendManage.createLegendSchemeFromGridData(mulData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
49. addLS = LegendManage.createLegendSchemeFromGridData(divData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
50. #Create shaded layers
    
51. print 'Create shaded layers...'
    
52. addLayer = DrawMeteoData.createShadedLayer(addData, addLS, 'Test_Add_Shaded', 'Add', True)
    
53. subLayer = DrawMeteoData.createShadedLayer(subData, subLS, 'Test_Sub_Shaded', 'Sub', True)
    
54. mulLayer = DrawMeteoData.createShadedLayer(mulData, addLS, 'Test_Mul_Shaded', 'Mul', True)
    
55. divLayer = DrawMeteoData.createShadedLayer(divData, addLS, 'Test_Div_Shaded', 'Div', True)
    
56. addLayer.setProjInfo(mdi1.getProjectionInfo())
    
57. subLayer.setProjInfo(mdi1.getProjectionInfo())
    
58. mulLayer.setProjInfo(mdi1.getProjectionInfo())
    
59. divLayer.setProjInfo(mdi1.getProjectionInfo())
    
60. #Add layers
    
61. print 'Add layers...'
    
62. mf = miapp.getMapDocument().getActiveMapFrame()
    
63. mf.getMapView().projectLayers(mdi1.getProjectionInfo())
    
64. mf.addLayer(addLayer)
    
65. mf.moveLayer(addLayer, 0)
    
66. mf.addLayer(subLayer)
    
67. mf.moveLayer(subLayer, 0)
    
68. mf.addLayer(mulLayer)
    
69. mf.moveLayer(mulLayer, 0)
    
70. mf.addLayer(divLayer)
    
71. mf.moveLayer(divLayer, 0)
    
72. 
    
73. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-9-12<br /> # Purpose: Add, sub, mul, div functions<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata import DrawMeteoData<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.legend import LegendType<br /> from org.meteoinfo.shape import ShapeTypes<br /> import os.path<br /> import jarray<br /> <br /> #Set data directory<br /> dataDir = 'D:/Temp/nc'<br /> <br /> #Create MeteoDataInfo object<br /> mdi1 = MeteoDataInfo()<br /> mdi2 = MeteoDataInfo()<br /> <br /> #Read data files<br /> fn1 = os.path.join(dataDir, 'wrfout_d01_2008-01-01_12_00_00.nc')<br /> fn2 = os.path.join(dataDir, 'wrfout_d01_2008-01-02_12_00_00.nc')<br /> print 'File 1: ' + fn1<br /> print 'File 2: ' + fn2<br /> mdi1.openNetCDFData(fn1)<br /> mdi2.openNetCDFData(fn2)<br /> #Get grid data<br /> print 'Get grid data...'<br /> gData1 = mdi1.getGridData('Q2')<br /> gData2 = mdi2.getGridData('Q2')<br /> #Add two grid data<br /> print 'Add grid data...'<br /> addData = gData1.add(gData2)<br /> #Subtract two grid data<br /> print 'Subtract grid data...'<br /> subData = gData1.sub(gData2)<br /> #Multiple two grid data<br /> print 'Multiple grid data...'<br /> mulData = gData1.mul(gData2)<br /> #Divide two grid data<br /> print 'Divide grid data...'<br /> divData = gData1.div(gData2)<br /> #Create legend scheme<br /> addLS = LegendManage.createLegendSchemeFromGridData(addData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> subLS = LegendManage.createLegendSchemeFromGridData(subData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> mulLS = LegendManage.createLegendSchemeFromGridData(mulData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> addLS = LegendManage.createLegendSchemeFromGridData(divData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> #Create shaded layers<br /> print 'Create shaded layers...'<br /> addLayer = DrawMeteoData.createShadedLayer(addData, addLS, 'Test_Add_Shaded', 'Add', True)<br /> subLayer = DrawMeteoData.createShadedLayer(subData, subLS, 'Test_Sub_Shaded', 'Sub', True)<br /> mulLayer = DrawMeteoData.createShadedLayer(mulData, addLS, 'Test_Mul_Shaded', 'Mul', True)<br /> divLayer = DrawMeteoData.createShadedLayer(divData, addLS, 'Test_Div_Shaded', 'Div', True)<br /> addLayer.setProjInfo(mdi1.getProjectionInfo())<br /> subLayer.setProjInfo(mdi1.getProjectionInfo())<br /> mulLayer.setProjInfo(mdi1.getProjectionInfo())<br /> divLayer.setProjInfo(mdi1.getProjectionInfo())<br /> #Add layers<br /> print 'Add layers...'<br /> mf = miapp.getMapDocument().getActiveMapFrame()<br /> mf.getMapView().projectLayers(mdi1.getProjectionInfo())<br /> mf.addLayer(addLayer)<br /> mf.moveLayer(addLayer, 0)<br /> mf.addLayer(subLayer)<br /> mf.moveLayer(subLayer, 0)<br /> mf.addLayer(mulLayer)<br /> mf.moveLayer(mulLayer, 0)<br /> mf.addLayer(divLayer)<br /> mf.moveLayer(divLayer, 0)<br /> <br /> print 'Finished!'

river

前来支持王老师！

rceclx

王老师的脚本示例很赞啊，学习后也要争取推出自己的脚本示例。

evita

谢谢王老师分享。好好学习了。

eroson

王老师太牛了，学习了~~

motonger

老师 ，您好 我想请教您一个问题 流域内只有降水格点数据，也没有实测站点信息，怎么利用格点数据计算流域面降雨量啊