关于如何使用MeteoInfo打开FY-2C云分类数据

墨墨卡托

我这里有2种格式的数据，第一种是附件1的hdf标准格式，第二种是AWX格式的数据；我想尝试的是使用MeteoInfo将这两种数据打开，经过一些处理，要有经纬度坐标的，然后保存为shp格式的文件或者是jpg格式的图片均可；但我处理时遇到的问题是：1、打不开我上传附件里云分类的hdf数据（图1）；2、可以打开awx的数据（图2），但不能加经纬度坐标，虽然可以看到在左下方有经纬度显示，而且无法保存成我需要的格式文件或图片；请哪位会的老师帮忙解答，谢谢！

MeteoInfo

出图的设置（包括经纬度线等）在“版面”中设置，这是MeteoInfo的基本操作，建议你仔细看看MeteoInfo的帮助文档。

建议你是用MeteoInfo最新的Java版本。FY的hdf文件C#版MeteoInfo无法处理，Java版可以，但也要写脚本，具体参见此帖：FY2E HDF格式数据处理绘图
http://bbs.06climate.com/forum.p ... 197&fromuid=106
(出处: 气象家园)

MeteoInfo

绘制你提供的HDF文件的脚本：

1. #Add data file
   
2. fn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'
   
3. f = addfile(fn)
   
4. #Get data variable
   
5. v = f['FY2C_Hourly_Cloud_Classification']
   
6. #Set x/y
   
7. x = linspace(-5750000.0,5750000.0, 2288)
   
8. y = linspace(-5750000.0,5750000.0, 2288)
   
9. #Set x/y dimensions
   
10. v.setdim('Y', y, 0)
    
11. v.setdim('X', x, 1)
    
12. #Get data array
    
13. cloud = v[:,:]
    
14. #Plot
    
15. axesm()
    
16. ax = axesm(proj='geos', lon_0=104.5, h=35785864, gridlabel=False, frameon=False)
    
17. mlayer = shaperead('D:/Temp/map/country1.shp')
    
18. geoshow(mlayer)
    
19. layer = imshowm(x, y, cloud, proj=ax.proj, cmap='grads_rainbow')
    
20. colorbar(layer)

#Add data file<br /> fn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'<br /> f = addfile(fn)<br /> #Get data variable<br /> v = f['FY2C_Hourly_Cloud_Classification']<br /> #Set x/y<br /> x = linspace(-5750000.0,5750000.0, 2288)<br /> y = linspace(-5750000.0,5750000.0, 2288)<br /> #Set x/y dimensions<br /> v.setdim('Y', y, 0)<br /> v.setdim('X', x, 1)<br /> #Get data array<br /> cloud = v[:,:]<br /> #Plot<br /> axesm()<br /> ax = axesm(proj='geos', lon_0=104.5, h=35785864, gridlabel=False, frameon=False)<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = imshowm(x, y, cloud, proj=ax.proj, cmap='grads_rainbow')<br /> colorbar(layer)

墨墨卡托

MeteoInfo 发表于 2015-6-30 12:36
绘制你提供的HDF文件的脚本：

王老师好，其他的都看懂了，但layer = imshowm(x, y, cloud, proj=proj, cmap='grads_rainbow')的rainbow比例尺您是怎么设的没看懂~，只能用grads吗，我没接触过这个软件，也不会用~

MeteoInfo

墨墨卡托 发表于 2015-6-30 14:19
王老师好，其他的都看懂了，但layer = imshowm(x, y, cloud, proj=proj, cmap='grads_rainbow')的rainbow ...

cmap参数是设定颜色的，这个参数也可以不给（有缺省的设置），你可以看看MeteoInfo软件目录中的colormaps目录，里面有很多颜色设置文件，cmap就是用这些文件来设置的（不用写文件后缀.rgb）。比如改为：cmap='hotcold_18lev'：

MeteoInfo

墨墨卡托 发表于 2015-6-30 14:19
王老师好，其他的都看懂了，但layer = imshowm(x, y, cloud, proj=proj, cmap='grads_rainbow')的rainbow ...

也可以把经纬线和标注加上：

1. #Add data file
   
2. fn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'
   
3. f = addfile(fn)
   
4. #Get data variable
   
5. v = f['FY2C_Hourly_Cloud_Classification']
   
6. #Set x/y
   
7. x = linspace(-5750000.0,5750000.0, 2288)
   
8. y = linspace(-5750000.0,5750000.0, 2288)
   
9. #Set x/y dimensions
   
10. v.setdim('Y', y, 0)
    
11. v.setdim('X', x, 1)
    
12. #Get data array
    
13. cloud = v[:,:]
    
14. #Plot
    
15. ax, proj = axesm(proj='geos', lon_0=104.5, h=35785864, gridlabel=True, gridline=True, frameon=False)
    
16. mlayer = shaperead('D:/Temp/map/country1.shp')
    
17. geoshow(mlayer, edgecolor=(0,0,255))
    
18. layer = imshowm(x, y, cloud, proj=proj, cmap='hotcold_18lev')
    
19. colorbar(layer)

#Add data file<br /> fn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'<br /> f = addfile(fn)<br /> #Get data variable<br /> v = f['FY2C_Hourly_Cloud_Classification']<br /> #Set x/y<br /> x = linspace(-5750000.0,5750000.0, 2288)<br /> y = linspace(-5750000.0,5750000.0, 2288)<br /> #Set x/y dimensions<br /> v.setdim('Y', y, 0)<br /> v.setdim('X', x, 1)<br /> #Get data array<br /> cloud = v[:,:]<br /> #Plot<br /> ax, proj = axesm(proj='geos', lon_0=104.5, h=35785864, gridlabel=True, gridline=True, frameon=False)<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor=(0,0,255))<br /> layer = imshowm(x, y, cloud, proj=proj, cmap='hotcold_18lev')<br /> colorbar(layer)

可以用鼠标在图形区域拉框放大，在图形区域点击鼠标右键，在右键菜单中点击‘Undo zoom’返回原图形范围。

墨墨卡托

MeteoInfo 发表于 2015-6-30 14:55
也可以把经纬线和标注加上：

老师辛苦了，我按照您说的重新下载了Java版本的，但还是打不开图像

MeteoInfo

墨墨卡托 发表于 2015-6-30 15:26
老师辛苦了，我按照您说的重新下载了Java版本的，但还是打不开图像

我这里没问题呀，你下载的什么版本？

墨墨卡托

MeteoInfo 发表于 2015-6-30 16:10
我这里没问题呀，你下载的什么版本？

嗯嗯，可见光的awx格式的可以打开了，但云分类的hdf就卡到这里了，还有老师您这的py是怎么调用的，我知道是java的py文件，但我不懂Java语言；

MeteoInfo

墨墨卡托 发表于 2015-6-30 16:33
嗯嗯，可见光的awx格式的可以打开了，但云分类的hdf就卡到这里了，还有老师您这的py是怎么调用的，我知道 ...

说过了hdf文件需要用脚本处理。在MeteoInfo软件目录下有MeteoInfoLab.exe文件，双击就能打开MeteoInfoLab软件。把代码粘贴到左上角Editor中，工具栏有保存文件按钮，工具栏中最右边的三角按钮是运行脚本程序的。你可能需要修改脚本程序中文件路径。