NCL怎么读取Cloudsat 2B-GEOPROF数据

denglin19901030

官网上的脚本里的个例都是E03的，现在已经变成E06了，用老的脚本会显示不存在变量，应该是数据格式发生了变化，但是要用新的数据怎么画图呢？

上楼无言

这个问题我之前也遇到过，困扰了很久，最后找到了一种很诡异的方法：把文件的后缀名由.hdf直接改为.hefeos，经纬度等信息就都能读取了。虽然我也不懂原理是什么。。。。。

longlivehj

denglin19901030 发表于 2015-3-12 16:05
这个试过的，我刚才还试了下filegrpdef这个命令想确定下我的文件里面到底有几个GROUP，可是出来的结果是0 ...

好吧！我下载了一个名为“2013119200420_37263_CS_2B-GEOPROF_GRANULE_P_R04_E06.hdf”的文件。ncl版本为6.1.2。

（1）用ncl_filedump查看文件信息。

ncl_filedump 2013119200420_37263_CS_2B-GEOPROF_GRANULE_P_R04_E06.hdf.he2

结果如下：


（2）在脚本中读取经纬度

f = addfile("2013119200420_37263_CS_2B-GEOPROF_GRANULE_P_R04_E06.hdf.he2", "r")
lon = f->Longitude_2B_GEOPROF
printVarSummary(lon)

结果如下：

flying鹰

这说明文件中不存在这个变量，有可能变量名有变化，你先看看文件中都有些什么变量~~~

denglin19901030

flying鹰 发表于 2015-3-11 19:55
这说明文件中不存在这个变量，有可能变量名有变化，你先看看文件中都有些什么变量~~~

查看过文件的，里面现在可以读的变量就只有Height;cpr_cloud_mask;gaseous_attenuation;radar_reflectivity这四个了，像老脚本里的Longitude都变成GEOFIELD这个结构体里的东西了，读不出来啊？

flying鹰

denglin19901030 发表于 2015-3-11 20:06
查看过文件的，里面现在可以读的变量就只有Height;cpr_cloud_mask;gaseous_attenuation;radar_reflectivi ...

你读radar_reflectivity就行了，这是个二维数组~其余类似

denglin19901030

flying鹰 发表于 2015-3-11 20:11
你读radar_reflectivity就行了，这是个二维数组~其余类似

但是我想要他的轨迹图，就是全球的那一种，没有经纬度的画怎么办呢？

longlivehj

http://www.ncl.ucar.edu/Applications/hdf5.shtml
用“=>”访问group，具体看链接中的例子。

denglin19901030

longlivehj 发表于 2015-3-11 21:46
http://www.ncl.ucar.edu/Applications/hdf5.shtml
用“=>”访问group，具体看链接中的例子。

多谢版主！！！！！！！！

denglin19901030

denglin19901030 发表于 2015-3-12 08:43
多谢版主！！！！！！！！

再请问下我看里面是结构体套着结构体的，但是按着一层一层往里写，最后会显示不存在变量？？

denglin19901030

longlivehj 发表于 2015-3-11 21:46
http://www.ncl.ucar.edu/Applications/hdf5.shtml
用“=>”访问group，具体看链接中的例子。

再请问下我看里面是结构体套着结构体的，但是按着一层一层往里写，最后会显示不存在变量？？求指点

longlivehj

刚仔细看了一下你发的图片，感觉问题复杂化了。楼主试一下，在addfile或者ncl_filedump时，文件名后加入he2后缀，即“…….hdf.he2”，看看坐标变量能不能出来。