IDLHDF5读取与转换

需求决定了动力，此时近凌晨一点，忙里偷闲，终于忙完了今天的“这点儿”事儿。参考帮助文档，从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;

效果显示: