版主帮帮忙，自己写的HDF文件不能读入

小寶

我自己写的HDF文件不能读入，但是在HDF Viewer是可以看到数据的。
这个是要符合一定的格式么？
如果自己写的HDF不行的话，请问还有什么方法可以解决

kongfeng0824

你用什么语言写的hdf？

MeteoInfo

需要HDF数据符合EOS规范才能自动读出数据的空间范围和格点坐标。
不符合EOS规范的HDF数据可以通过写脚本给定格点坐标来读取数据并生成图层，下面给一个Jython脚本的例子。


用你给的数据生成的图层：


脚本代码：

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 + 's66.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('ddd\dddd')
    
53.         dimList = [yDim, xDim]
    
54.         var.setDimensions(dimList)
    
55.         gData = mdi.getGridData(var.getName())
    
56.         gData = gData.multiply(0.01)
    
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!'

#-----------------------------------------------------<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 + 's66.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('ddd\dddd')<br />         dimList = [yDim, xDim]<br />         var.setDimensions(dimList)<br />         gData = mdi.getGridData(var.getName())<br />         gData = gData.multiply(0.01)<br />         aLS = LegendManage.createLegendSchemeFromGridData(gData, LegendType.GraduatedColor, ShapeTypes.Polygon)<br />         aLayer = DrawMeteoData.createRasterLayer(gData, "Test_HDF", aLS)<br />         aLayer.setProjInfo(dataInfo.getProjectionInfo())<br />         mf = miapp.getMapDocument().getActiveMapFrame()<br />         mf.addLayer(aLayer)<br />         mf.moveLayer(aLayer, 0)<br /> <br /> print 'Finished!'

小寶

kongfeng0824 发表于 2013-10-29 22:04

                               
登录/注册后可看大图

你用什么语言写的hdf？

IDL 一个遥感处理软件 不过不普及我觉得

小寶

MeteoInfo 发表于 2013-10-29 22:33

                               
登录/注册后可看大图

需要HDF数据符合EOS规范才能自动读出数据的空间范围和格点坐标。
不符合EOS规范的HDF数据可以通过写脚本给 ...

你好用脚本写入图层，但是不能选择是哪个map frame，怎么可以把脚本写入的图层加入其他frame？

MeteoInfo

小寶 发表于 2013-11-2 08:37

                               
登录/注册后可看大图

你好用脚本写入图层，但是不能选择是哪个map frame，怎么可以把脚本写入的图层加入其他frame？

代码里是将图层加入Active MapFarme（可以有多个MapFrame，但Active MapFrame只有一个）。可以获得MapFrame的列表,然后可以将图层加入某个MapFrame中，比如：

mfs = miapp.getMapDocument().getMapFrames()
mfs.get(1).addLayer(aLayer)

小寶

MeteoInfo 发表于 2013-11-2 10:47

                               
登录/注册后可看大图

代码里是将图层加入Active MapFarme（可以有多个MapFrame，但Active MapFrame只有一个）。可以获得MapFra ...

脚本里面有没有吧图像上下翻转和左右翻转的命令呢？

MeteoInfo

小寶 发表于 2014-1-4 16:36

                               
登录/注册后可看大图

脚本里面有没有吧图像上下翻转和左右翻转的命令呢？

没有这样的功能

小寶

MeteoInfo 发表于 2013-10-29 22:33

                               
登录/注册后可看大图

需要HDF数据符合EOS规范才能自动读出数据的空间范围和格点坐标。
不符合EOS规范的HDF数据可以通过写脚本给 ...

老师您好，这里有个问题就是如果绘制等值线
在代码中改？怎么改  哪个教程是对应这个的？
如果不用代码如何弄，搞半天没弄出来，打扰了

MeteoInfo

小寶 发表于 2014-3-7 12:58

                               
登录/注册后可看大图

老师您好，这里有个问题就是如果绘制等值线
在代码中改？怎么改  哪个教程是对应这个的？
如果不用代码 ...

可以参考此帖中的脚本：绘制闪电密度图-2
http://bbs.06climate.com/forum.p ... 824&fromuid=106