用matlab读取cloudsat卫星数据

kongfeng0824

转自：http://hi.baidu.com/sarali/item/6d76380ebb41be183b53ee79

简单的把自己学到的东西总结一下，可能会比较乱

Height=hdfread('E:\cloudproperties\cloudsat\cloud\data\2008269004730_12825_CS_2B-TAU_GRANULE_P_R04_E02.hdf', '/2B-TAU/Geolocation Fields/Height', 'Index', {[1  1],[1  1],[37081    125]});

在这里面，Height是一个37081*125的矩阵。全部为整形。读取时的书写就是

Para = hdftead('fpath\fname','子集的存放目录','index',{[1 1],[1 1]},[37081 125]}),其中index中包含的是行列的起始位置以及步长的设置信息。

Profile_time = hdfread('E:\cloud properties\cloudsat\cloud\data\2008269004730_12825_CS_2B-TAU_GRANULE_P_R04_E02.hdf', '/2B-TAU/Geolocation Fields/Profile_time', 'Fields', 'Profile_time', 'FirstRecord',1 ,'NumRecords',37081);

显示的Profile_time 的数据结构式相同的。如果利用hdftool打开可以看到是结构体。因此这种结构的打开与Height的打开是不同的。基本的读取格式为

Para = hdftead('fpath\fname','子集的存放目录','Fields/para','Fields','para','FirstRecord',1,'NumRwcords',1);

其中的FirstRecord与NumRwcords代表的是结构体的行列数。

ID_SITE = hdfread('E:\cloud properties\cloudsat\cloud\data\2008269004730_12825_CS_2B-TAU_GRANULE_P_R04_E02.hdf', '/2B-TAU/Swath Attributes/ID_SITE', 'Fields', 'AttrValues', 'FirstRecord',1 ,'NumRecords',1);

如果是属性值就要用上述方法。

不过一般在matlab中输入hdftool之后，加载数据时会有相应的语句书写示例，可以把语句直接复制出来。

---

我个人的读取方法：
hdfname='E:\IDL\pf\JJA_3year_level3.HDF'
A=hdftool(hdfname)
LON = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/LON', 'Index', {[1  1],[1  1],[170  360]});
LAT = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/LAT', 'Index', {[1  1],[1  1],[170  360]});
Z = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/Z', 'Index', {[1],[1],[500]});
PIX_TOT = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_TOT', 'Index', {[1  1  1  1],[1  1  1  1],[2  500  170  360]});
PIX_CLD = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_CLD', 'Index', {[1  1  1  1],[1  1  1  1],[2  500  170  360]});
以上的最后两行：
PIX_TOT = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_TOT', 'Index', {[1 1 1 1],[1 1 1 1],[2 500 170 360]});
PIX_CLD = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_CLD', 'Index', {[1 1 1 1],[1 1 1 1],[2 500 170 360]});
是什么含义？我还是没有完全弄清楚。
应该如何提取出PIX_TOT和PIX_CLD的数据信息并画图呢？
PIX_TOT和PIX_CLD是4-Dsingle的数据。

wlzhongouc

PIX_TOT = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_TOT', 'Index', {[1 1 1 1],[1 1 1 1],[2 500 170 360]});
PIX_CLD = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_CLD', 'Index', {[1 1 1 1],[1 1 1 1],[2 500 170 360]});

这两行就是读取四维数据的语句，[2 500 170 360]分别对应着每一维的大小

kongfeng0824

wlzhongouc 发表于 2013-3-15 08:51

                               
登录/注册后可看大图

PIX_TOT = hdfread('E:\IDL\pf\JJA_3year_level3.HDF', '/PIX_TOT', 'Index', {[1 1 1 1],[1 1 1 1],[2 500 ...

嗯。我知道前两个中括号的元素代表的是起始数据和步长。但是我的所有维度有5个变量，我就不太清楚如何访问读取出这个结构体中的数据，并画出图来。因为我想用的就死以上两列代码读出来的结构体。TOT和CLD。所以现在接下来要做的就是如何访问这个结构体的内容。它是四维的，应该如何访问每一维，并计算和作图呢？

wlzhongouc

kongfeng0824 发表于 2013-3-15 21:47

                               
登录/注册后可看大图

嗯。我知道前两个中括号的元素代表的是起始数据和步长。但是我的所有维度有5个变量，我就不太清楚如何访问 ...

四维的没法作图出来显示吧 你可以用squeeze降维 比如说先一个固定时间和一层数然后就成为二维的 再画图

kongfeng0824

wlzhongouc 发表于 2013-3-16 12:53

                               
登录/注册后可看大图

四维的没法作图出来显示吧 你可以用squeeze降维 比如说先一个固定时间和一层数然后就成为二维的 再画图

这个应该属于结构体的知识。看来我得去再好好看看教材了。

qxtlyf

多维数据绘图，需要转化为两维进行绘图，不过也可以画3维的图

kongfeng0824

qxtlyf 发表于 2013-3-17 10:13

                               
登录/注册后可看大图

多维数据绘图，需要转化为两维进行绘图，不过也可以画3维的图

嗯。现在是如何把结构体的数据以相同的向量读取出来，这样作图。

幻想的春

这个先mark 很重要！！！

Jing123

楼主解决了吗。。我现在要画AOD的图，读取了hdf的AOD信息之后不会画图啊，求教！