需求决定了动力，此时近凌晨一点，忙里偷闲，终于忙完了今天的“这点儿”事儿。参考帮助文档，从hdf的读写，捉摸hdf5的读写，总算弄明白了。稍作总结，以备候查。
hdf5作为hdf数据的补充与升级，目前被新型的传感器所采用（VIIRS，风云三号等），但是如果没有hdf5的sav扩展，低于5.1版本的envi是无法打开的（扩展可以到网上搜）方便起见，花一点儿时间捉摸了下hdf5的数据读取，提取部分内容，写为hdf文件，以方便查看。
Hdf5数据重要的是几个概念：文件(File),组(Group)，数据集(Dataset)。常用读的数据过程是：文件读取，组的读取和数据集的读取。
对应使用的函数是h5f_,h5g_,h5d_这三种。
常用的函数是：
h5_list可以列出一个文件中的所有的数据集，然后可以通过这些数据集，直接利用相应的read来获取数据。
H5F_IS_HDF5判断是不是h5格式的数据；H5F_OPEN打开一个文件，返回文件的id；H5G_GET_NMEMBERS可以根据组别的名字获取组的成员个数；H5G_OPEN是根据组别的名字打开一个组，得到组的id；H5G_GET_MEMBER_NAME根据元素在文件中的索引位置（从0开始）获取数据集，返回的是数据集的name；H5D_OPEN根据数据集的name打开数据集，返回数据集的id；H5D_READ是根据数据集的id获取该数据集的数据内容。然后是对应的close函数。记得，打开的数据要注意关闭。
例子：读取hdf5的数据（seadas产生的rrc数据），提取出bands组别内的所有数据集，写入到一个hdf文件中，便于其他程序访问。
PRO HDF5
;
COMPILE_OPT idl2
ENVI,/RESTORE_BASE_SAVE_FILES
ENVI_BATCH_INIT
;
hdf = 'F:\temp\A2014119053500.L2_LAC_OC_reprojected.h5'
is_h5 = H5F_IS_HDF5(hdf)
IF is_h5 EQ 0 THEN RETURN
;h5_list,hdf
hdf_id = H5F_OPEN(hdf)
group= '/bands'
elements_num = H5G_GET_NMEMBERS(hdf_id,group)
band_g_id = H5G_OPEN(hdf_id,group)
;
outpath = 'F:\temp\hdf_test.hdf'
if file_test(outpath) then file_delete,outpath
sd_id=HDF_SD_START(outpath,/CREATE)
;进度条
str=['Outpath='+outpath,'Reading and writing']
ENVI_REPORT_INIT, str, title="Processing...",base=base ,/INTERRUPT
ENVI_REPORT_INC, base, elements_num
FOR index = 0L, elements_num-1 DO BEGIN
;更新进度条
ENVI_REPORT_STAT,base, index, elements_num
;
item_Name = H5G_GET_MEMBER_NAME(hdf_id,group,index)
item_id = H5D_OPEN(band_g_id,item_Name)
;item_type_code = H5D_GET_TYPE(item_id)
;item_type = H5T_IDLTYPE(item_type_code)
item_data = H5D_READ(item_id)
;print,item_Name,typename(item_data)
;help,item_data
WRITE_HDF,sd_id,item_Name,item_data,typename(item_data)
h5d_close,item_id
ENDFOR
;close
h5g_close,band_g_id
h5f_close,hdf_id
HDF_SD_END,sd_id
envi_report_init,base=base,/finish
void = dialog_message('Processing completed!',/info)
; print, members
;envi_batch_exit
END
writehdf的pro：
PRO WRITE_HDF,sd_id,data_set_name,data,datatype
;add dataset into hdf file
;sd_id: hdf identifer;data_set_name: the name you will add
;datatype:the data type
sds_id = 1
dims = SIZE(data,/dimensions)
; datatype = '/'+datatype
CASE (datatype) OF
'LONG': BEGIN
sds_id=HDF_SD_CREATE(sd_id,data_set_name,dims[*,*],/long)
END
'FLOAT':BEGIN
sds_id=HDF_SD_CREATE(sd_id,data_set_name,dims[*,*],/float)
END
ELSE: BEGIN
sds_id=HDF_SD_CREATE(sd_id,data_set_name,dims[*,*],/DFNT_UINT32)
END
ENDCASE
HDF_SD_ADDDATA, sds_id, REFORM(data[*,*])
HDF_SD_SETINFO,sds_id,Label = data_set_name
HDF_SD_ENDACCESS,sds_id
;
END;
效果显示:
