FY2E HDF格式数据处理绘图

liuzijing333

老师 这是啥意思啊 谢谢老师

MeteoInfo

liuzijing333 发表于 2017-8-21 15:23
老师 这是啥意思啊 谢谢老师

脚本程序出错了，第10行

liuzijing333

MeteoInfo 发表于 2017-8-21 15:27
脚本程序出错了，第10行

是啊 我就用的您的脚本，改了自己的数据路径，就报错。

liuzijing333

这就是第十行[
v.setdim是什么意思啊 老师

sysun

MeteoInfo 发表于 2015-9-11 14:32
有投影的数据用投影后的坐标更方便。

王老师，这样格式的云图数据该怎么读取？
Z_SATE_C_BAWX_20180901080053_P_FY2E_FDI_ALL_NOM_20180901_0730.hdf

MeteoInfo

sysun 发表于 2018-9-7 11:00
王老师，这样格式的云图数据该怎么读取？
Z_SATE_C_BAWX_20180901080053_P_FY2E_FDI_ALL_NOM_20180901_0 ...

你这个数据文件好像有问题，打不开。

ljqtibet

谢谢楼主，有问题再请教

孤蓝et

王老师，不知道MI是否支持FY4的HDF数据文件？

MeteoInfo

孤蓝et 发表于 2019-1-23 16:23
王老师，不知道MI是否支持FY4的HDF数据文件？

支持

孤蓝et

MeteoInfo 发表于 2019-1-23 19:46
支持

谢谢老师。按照老师的FY2E示例，刚已经测试出来了。

1. #Add data file
   
2. fn = 'D:/FY/FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20190123043000_20190123043416_4000M_V0001.HDF'
   
3. f = addfile(fn)
   
4. print(f)
   
5. 
   
6. #Get data
   
7. NOMChannel02 = f['NOMChannel02']
   
8. print(NOMChannel02)

复制代码

利用上述代码：
print(f)，我可以看到这个HDF的信息，包括卫星位置，高度，变量等信息；
print(NOMChannel02)可以查看到具体的变量信息，具体如下：
short NOMChannel02(null,null):
        NOMChannel02: valid_range = 0US, 4095US
        NOMChannel02: FillValue = -1US
        NOMChannel02: Intercept = 0.0f
        NOMChannel02: Slope = 1.0f
        NOMChannel02: units = "DN"
        NOMChannel02: center_wavelength = "0.65um"
        NOMChannel02: band_names = "band2(band number is range from 1 to 14)"
        NOMChannel02: long_name = "0.65um channel 1km image data layer"
        NOMChannel02: _Unsigned = "true"

但是，还有三个疑问：
第一个，FY4的数据是全圆盘图和区域图穿插生成的，我此处使用的是区域图，投影是否和全圆盘图有差异？

第二个，后续，x  、y这个怎么理解？麻烦老师指点下？

• #Set x/y
• 
                                 
  登录/注册后可看大图
  x = linspace(-5750000.0,5750000.0, 2288)
• 
                                 
  登录/注册后可看大图
  y = linspace(-5750000.0,5750000.0, 2288)
  

例如其中的5750000,2288这两个量如何理解，我是否需要配合FY4的数据描述信息来修改？

第三个，如果绘制某一个变量的图形，我直接利用demo中的NOMChannel02 = f['NOMChannel02']，就可以了吗？运行程序有提示错误，是不是我读取数据的方式有问题？
>>> run script...
Traceback (most recent call last):
  File "D:\FY\test.py", line 13, in <module>
    layer = imshowm(x, y, NOMChannel02, 20, proj=ax.proj)
  File "D:\MeteoInfo_1.6.10\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 2212, in imshowm
    r = gca.imshow(*args, **kwargs)
  File "D:\MeteoInfo_1.6.10\MeteoInfo\pylib\mipylib\plotlib\mapaxes.py", line 967, in imshow
    gdata = minum.asgridarray(a, x, y, fill_value)
  File "D:\MeteoInfo_1.6.10\MeteoInfo\pylib\mipylib\numeric\minum.py", line 1711, in asgridarray
    gdata = GridArray(data.asarray(), x.asarray(), y.asarray(), fill_value)
AttributeError: DimVariable instance has no attribute 'asarray'