Matlab读写HDF-EOS文件

最近才发现Matlab集成了读写HDF文件（包括HDF-EOS文件）的功能。

HDF-EOS的三种数据类型及其对应的Matlab函数如下：
grid(GD) - hdfgd
point(PT) - hdfpt
swath(SW) - hdfsw

此外，Matlab还提供了GUI的HDF文件导入工具，可以通过 “hdftool”命令启动。

hdfgd的帮助需要用 “help hdfgd" 查看。
hdfgd可以认为是HDF-EOS库中所有Grid函数的一个共用接口，即通过hdfgd来调用不同的Grid函数。
以下是读取一个HDF-EOS文件的代码示例，可以看出完成相同的任务，Matlab比C/C++所需代码量少得多，当然执行效率会不如C/C++。

hdf = 'MOD15A2.A2006361.h10v04.005.2008134052818.hdf';

fileid = hdfgd('open', hdf, 'DFACC_READ');                 %对应GDopen

gridid = hdfgd('attach', fileid, 'MOD_Grid_MOD15A2');      %对应GDattach

[data, status] = hdfgd('readfield', gridid, 'Lai_1km', [], [], []);  %对应GDreadfield

hdfgd('detach', gridid);  %对应GDdetach

hdfgd('close', fileid);   %对应GDclose

全球LAI产品(Global LAI Products)

文献[1]对4种覆盖全球的LAI产品进行了综合评价：MODIS，CYCLOPES，ECOCLIMAP和GLOBCARBON，如下表所示。其它还有一些全球LAI产品，如POLDER和Land-SAF。还有一些区域性的LAI产品，如加拿大范围的LAI产品。这些产品中，既有基于物理模型，实现时采用查找表(LUT，Look-Up Table)或神经网络；也有基于LAI与植被指数(VI，Vegetation Index)统计关系的。

LAI产品名	时间分辨率	空间分辨率	年份	传感器	产品生成算法
MODIS	8天	1km	2000-	MODIS	主算法：3维辐射传输模型（LUT实现） 备用算法：基于NDVI的经验关系
CYCLOPES	10天	1km	1999-2007	VEGETATION	1维辐射传输模型（神经网络实现）
ECOCLIMAP	1月	1km	climatology	AVHRR	基于NDVI的经验关系
GLOBCARBON	1月	0.25°, 0.5°, 10km	1998-2007	AATSR, ATSR, MERIS, VEGETATION	基于VI的统计关系

(climatology的意思是：the mean state based on a historical set of observations)

以上四种产品的相互比较及与地面实测数据的比较[1]表明，MODIS和CYCLOPES对真实地表及其空间差异的表达，LAI时间序列特征及年份间的差异都优于ECOCLIMAP和GLOBCARBON。但另一方面，MODIS和CYCLOPES的反演成功率却低于ECOCLIMAP和GLOBCARBON，因此MODIS和CYCLOPES存在较多缺失值，时间序列不够连续，影响了产品的应用。

事实上，从上表也可以看出，目前只有MODIS的LAI产品在持续发布，其它产品处于产品评价和算法改进过程中。

[1] Garrigues, S., et al. (2008), Validation and intercomparison of global Leaf Area Index products derived from remote sensing data, J. Geophys. Res., 113, G02028, doi:10.1029/2007JG000635