grads如何读取HDF-EOS文件

ziqiangbuxi

最近在学习grads处理HDF文件，其中一个示例文件数据为：AMSR_E_L3_SeaIce12km_B02_20020619.hdf，看到hdf官网说明可用sdfopen、xdfopen以及open打开，上面这个数据用sdfopen、xdfopen无法打开，只能用open处理，但是由于是ease-grid投影，需要编写ctl文件。在ctl文件中需要用到pdef参数，参考官网示例，ctl文件如下：

DSET AMSR_E_L3_SeaIce12km_B02_20020619.hdf
DTYPE hdfsds
UNDEF 0 _FillValue
PDEF 608 896 BILIN STREAM BINARY AE_SI12.bilin
XDEF 360 linear -179.5 1
YDEF 180 linear -89.5 1
ZDEF 1 levels 0
TDEF 1 linear 00Z19jun2002 1mo
VARS 1
SI_12km_NH_SNOWDEPTH_5DAY=>snow 0 y,x SI_12 Snow Depth
ENDVARS
现在对pdef选项有些疑问：
（1）pdef中的BILIN需要重新生成一个文件，此文件包含每个格点中的3个float型的数据：i values, j values, and wind rotation values
但是这个文件中没有wind这项啊，请问这个wind rotation是根据什么数据计算的，是怎么计算的？整个AE_SI12.bilin文件如何获得，是有程序调用接口还是需要自编程实现？
（2）这个AE_SI12.bilin文件用MATLAB读取后，显示为777600*1的矩阵，跟XDEF和YDEF是什么对应关系？
（3）参考grads官网对于pdef的说明，如果原始文件是polar stereographic投影，pdef可以定义如下：
        PDEF isize jsize NPS ipole jpole lonref gridinc
        那么ipole jpole lonref gridinc这几个参数是怎么获取？
       

ziqiangbuxi

请各位高手给予指导啊!

river

你说的那几个相关数据，应该在你下载hdf资料的网站有说明的吧，我觉得

ziqiangbuxi

river 发表于 2016-1-3 19:16
你说的那几个相关数据，应该在你下载hdf资料的网站有说明的吧，我觉得

这个数据是官网提供的学习示例数据，AE_SI12.bilin也提供了，但是里面只是说用程序可以实现，我看了半天，还是不知道怎么生成这个文件。

ziqiangbuxi

请教大家指点啊，不能沉啊！