关于第六讲的ShadedActionPerformed求解

小寶

第六讲的ShadedActionPerformed（）这个方法生成上图
想用来改自己的图片
我的图片这样

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2013-10-29
   
4. # Purpose: Read hdf4 data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata import Dimension
   
9. from org.meteoinfo.data.meteodata import DimensionType
   
10. from org.meteoinfo.data.meteodata import DrawMeteoData
    
11. from org.meteoinfo.projection import ProjectionInfo
    
12. from org.meteoinfo.legend import LegendManage
    
13. from org.meteoinfo.legend import LegendType
    
14. from org.meteoinfo.shape import ShapeTypes
    
15. import os.path
    
16. import jarray
    
17. 
    
18. #Set data directory
    
19. dataDir = 'e:/'
    
20. 
    
21. #Create MeteoDataInfo object
    
22. mdi = MeteoDataInfo()
    
23. 
    
24. #Read hdf data file
    
25. fn = dataDir + 'tend.HDF'
    
26. if os.path.isfile(fn):
    
27.         print fn
    
28.         mdi.openNetCDFData(fn)
    
29.         dataInfo = mdi.getDataInfo()
    
30.         print mdi.getInfoText()
    
31.         xmin = -180.0
    
32.         ymin = -90.0
    
33.         xnum = 360
    
34.         ynum = 180
    
35.         xdelt = 1.0
    
36.         ydelt = 1.0
    
37.         xlist = []
    
38.         ylist = []
    
39.         for i in range(0,xnum):
    
40.                 xlist.append(xmin + xdelt * i)
    
41.         for i in range(0,ynum):
    
42.                 ylist.append(ymin + ydelt * i)
    
43. 
    
44.         X = jarray.array(xlist, 'd')
    
45.         Y = jarray.array(ylist, 'd')
    
46.         xDim = Dimension(DimensionType.X)
    
47.         xDim.setValues(X)
    
48.         dataInfo.setXDimension(xDim)
    
49.         yDim = Dimension(DimensionType.Y)
    
50.         yDim.setValues(Y)
    
51.         dataInfo.setYDimension(yDim)
    
52.         var = dataInfo.getVariable('AOD_Mean_Tend')
    
53.         dimList = [yDim, xDim]
    
54.         var.setDimensions(dimList)
    
55.         gData = mdi.getGridData(var.getName())
    
56.         gData = gData.multiply(1)
    
57.         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
58.         aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)
    
59.         aLayer.setProjInfo(dataInfo.getProjectionInfo())
    
60.         mf = miapp.getMapDocument().getActiveMapFrame()
    
61.         mf.addLayer(aLayer)
    
62.         mf.moveLayer(aLayer, 0)
    
63. 
    
64. print 'Finished!'

复制代码

用之前老师给的例子，但是这个jython代码 好像和java不一样，改动的时候不知道如何下手，改了几次都出错
我用的是java版 ，本版内部是不是很多用的是c#呀  为什么一些按键我都没有

小寶

MeteoInfo@MeteoInfo打扰了有空帮看下？
二次开发一定要安装第一讲那个软件？

MeteoInfo

要生成等值线填充图层可以用DrawMeteoData类的createShadedLayer方法，由于要生成等值线并填充，此方法会比较慢，如果只是关注中国区域可以先将GridData的范围缩小（用GridData类的extract方法），然后再生成shaded图层，脚本代码如下：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2013-10-29
   
4. # Purpose: Read hdf4 data
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata import Dimension
   
9. from org.meteoinfo.data.meteodata import DimensionType
   
10. from org.meteoinfo.data.meteodata import DrawMeteoData
    
11. from org.meteoinfo.projection import ProjectionInfo
    
12. from org.meteoinfo.legend import LegendManage
    
13. from org.meteoinfo.legend import LegendType
    
14. from org.meteoinfo.shape import ShapeTypes
    
15. import os.path
    
16. import jarray
    
17. 
    
18. #Set data directory
    
19. dataDir = 'D:/Temp/hdf/'
    
20. 
    
21. #Create MeteoDataInfo object
    
22. mdi = MeteoDataInfo()
    
23. 
    
24. #Read hdf data file
    
25. fn = dataDir + 'tend.HDF'
    
26. if os.path.isfile(fn):
    
27.         print fn
    
28.         mdi.openNetCDFData(fn)
    
29.         dataInfo = mdi.getDataInfo()
    
30.         print mdi.getInfoText()
    
31.         xmin = -180.0
    
32.         ymin = -90.0
    
33.         xnum = 360
    
34.         ynum = 180
    
35.         xdelt = 1.0
    
36.         ydelt = 1.0
    
37.         xlist = []
    
38.         ylist = []
    
39.         for i in range(0,xnum):
    
40.                 xlist.append(xmin + xdelt * i)
    
41.         for i in range(0,ynum):
    
42.                 ylist.append(ymin + ydelt * i)
    
43. 
    
44.         X = jarray.array(xlist, 'd')
    
45.         Y = jarray.array(ylist, 'd')
    
46.         xDim = Dimension(DimensionType.X)
    
47.         xDim.setValues(X)
    
48.         dataInfo.setXDimension(xDim)
    
49.         yDim = Dimension(DimensionType.Y)
    
50.         yDim.setValues(Y)
    
51.         dataInfo.setYDimension(yDim)
    
52.         var = dataInfo.getVariable('AOD_Mean_Tend')
    
53.         dimList = [yDim, xDim]
    
54.         var.setDimensions(dimList)
    
55.         gData = mdi.getGridData(var.getName())
    
56.         #gData = gData.multiply(0.01)
    
57.         gData = gData.extract(70.0, 140.0, 15.0, 60.0)
    
58.         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
59.         #aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)
    
60.         aLayer = DrawMeteoData.createShadedLayer(gData, aLS, 'Test_HDF_Shaded', 'AOD', True)
    
61.         aLayer.setProjInfo(dataInfo.getProjectionInfo())
    
62.         mf = miapp.getMapDocument().getActiveMapFrame()
    
63.         mf.addLayer(aLayer)
    
64.         mf.moveLayer(aLayer, 0)
    
65. 
    
66. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2013-10-29<br /> # Purpose: Read hdf4 data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata import Dimension<br /> from org.meteoinfo.data.meteodata import DimensionType<br /> from org.meteoinfo.data.meteodata import DrawMeteoData<br /> from org.meteoinfo.projection import ProjectionInfo<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/hdf/'<br /> <br /> #Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #Read hdf data file<br /> fn = dataDir + 'tend.HDF'<br /> if os.path.isfile(fn):<br />         print fn<br />         mdi.openNetCDFData(fn)<br />         dataInfo = mdi.getDataInfo()<br />         print mdi.getInfoText()<br />         xmin = -180.0<br />         ymin = -90.0<br />         xnum = 360<br />         ynum = 180<br />         xdelt = 1.0<br />         ydelt = 1.0<br />         xlist = []<br />         ylist = []<br />         for i in range(0,xnum):<br />                 xlist.append(xmin + xdelt * i)<br />         for i in range(0,ynum):<br />                 ylist.append(ymin + ydelt * i)<br /> <br />         X = jarray.array(xlist, 'd')<br />         Y = jarray.array(ylist, 'd')<br />         xDim = Dimension(DimensionType.X)<br />         xDim.setValues(X)<br />         dataInfo.setXDimension(xDim)<br />         yDim = Dimension(DimensionType.Y)<br />         yDim.setValues(Y)<br />         dataInfo.setYDimension(yDim)<br />         var = dataInfo.getVariable('AOD_Mean_Tend')<br />         dimList = [yDim, xDim]<br />         var.setDimensions(dimList)<br />         gData = mdi.getGridData(var.getName())<br />         #gData = gData.multiply(0.01)<br />         gData = gData.extract(70.0, 140.0, 15.0, 60.0)<br />         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br />         #aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)<br />         aLayer = DrawMeteoData.createShadedLayer(gData, aLS, 'Test_HDF_Shaded', 'AOD', True)<br />         aLayer.setProjInfo(dataInfo.getProjectionInfo())<br />         mf = miapp.getMapDocument().getActiveMapFrame()<br />         mf.addLayer(aLayer)<br />         mf.moveLayer(aLayer, 0)<br /> <br /> print 'Finished!'

运行结果：

小寶

顺便问下 这些帮助文档哪里有？ java的?

MeteoInfo

小寶 发表于 2014-3-10 22:39

                               
登录/注册后可看大图

顺便问下 这些帮助文档哪里有？ java的?

1、MeteoInfo网站：http://www.meteothinker.com/
2、本论坛
3、博客：http://hi.baidu.com/meteoinfo

就这些了，主要是太忙了，难以抽出时间写帮助文档。

小寶

MeteoInfo 发表于 2014-3-10 22:50

                               
登录/注册后可看大图

1、MeteoInfo网站：http://www.meteothinker.com/
2、本论坛
3、博客：http://hi.baidu.com/meteoinfo
...

图一是原来的，标注的值也都可以每0.2个点标记一次


图2 是用等直线，发现虽然趋势不错，但是标注不对粉红色过多没有梯度的变化


图3 我是改后的值，之前改过好多，就不贴了，这张主要是想表示第一图东北的高值点0.8-0.1
等值线改过好多标注值，可是大多都不怎么变化 就是一片粉色值，不知道这是为什么？

MeteoInfo

小寶 发表于 2014-3-11 13:10

                               
登录/注册后可看大图

图一是原来的，标注的值也都可以每0.2个点标记一次

栅格图层（RasterLayer）和矢量图层（VectorLayer）是不同的。DrawMeteoData.createRasterLayer方法生成的是栅格图层，原始的格点数据保存在栅格图层中，更改图例会起到作用。DrawMeteoData.createShadedLayer方法生成的是矢量图层，保存的是多边形（Polygon）的坐标和属性信息，可以修改图例但图元形状并不能随图例改变。也就是说对于矢量图层需要在等值线分析的时候就把图例设置好。你可以先利用之前生成的图层编辑好图例，将图例保存保存为图例文件（.lgs），然后可以在脚本中加载图例文件，然后在生成等值线填充图层。

从文件加载图例可以参照下面的语句：
aLS = LegendScheme(ShapeTypes.Polygon)
aLS.importFromXMLFile('D:\\Temp\\rain1.lgs')