C语言对HDF文件数据集和属性的读写.doc

以三维数据读写为例打开，读取空间属性、维度，数据集属性，文件属性，数据集的数据；/ 写入txt文件（动态定义数组用指针），写入DAT文件（可以用ENVI打开查看图像），写入HDF文件写入HDF文件属性、数据属性、空间属性，数据压缩。hid_t file_id, space_id,dset_id;hsize_t dims3=0,0,0;int ndims=0;int status;int col = 0;int row =0;int high=0;char FileName200;char OutFileName200;int N=0;int i,j,k; memset(FileName,0,sizeof(FileName);memset(OutFileName,0,sizeof(OutFileName);strcpy(FileName,C:UsershanssDesktopFY3B_MERSI_GBAL_L1_20170301_0000_1000M_MS.HDF);strcpy(OutFileName,C:UsershanssDesktopFY3B.HDF);unsigned int* dset=new unsigned introw * col*high;memset(dset,0,sizeof(unsigned int)* row * col*high);一、 打开HDF文件并读取相关数据信息1、打开HDF文件file_id = H5Fopen(FileName, H5F_ACC_RDONLY, H5P_DEFAULT);2、打开数据集dset_id = H5Dopen(file_id,EV_1KM_RefSB,H5P_DEFAULT);3、获取空间地址space_id = H5Dget_space(dset_id);4、获取空间维度信息ndims=H5Sget_simple_extent_dims(space_id,dims,NULL);5、读取数据集status=H5LTread_dataset(file_id,EV_1KM_RefSB,H5T_NATIVE_UINT,dset);6、读取数据属性float Slope30;int _FillValue1;char band_name41;H5LTget_attribute_float(file_id,EV_1KM_RefSB,Slope,Slope);for(int i=0;i15;i+) printf(Slope:%fn,Slopei);H5LTget_attribute_int(file_id,EV_1KM_RefSB,_FillValue,_FillValue); printf(_FillValue:%dn,_FillValue0);H5LTget_attribute_string(file_id,EV_1KM_RefSB,band_name,band_name); printf(%sn,band_name);或者用H5LTget_attribute_float(dset_id,/EV_1KM_RefSB,Slope,Slope);（推荐）H5LTget_attribute_float(file_id,/EV_1KM_RefSB,Slope,Slope);或者char dset_name20=EV_1KM_RefSB;H5LTget_attribute_float(file_id,dset_name,Intercept,Intercept);【注释】1、 如果是file_id的话，不加斜杠也是可以的，但是如果时dse_id则一定要加斜杠2、 如果是file_id,加斜杠也可以3、 输出数组时遍历每个元素才行，直接输出数组名，其实是输出了数组地址，需要加*，那就是输出数组第一个元素值；也可以定义指针，首先分配空间float *Intercept;Intercept=(float*)malloc(sizeof(float)*15);H5LTget_attribute_float(dset_id,/EV_1KM_RefSB,Intercept,Intercept); printf(Intercept:n);for(int i=0;i15;i+) printf(%fn,*(Intercept+i);4、 不要把IDL语法放到C中5、 C+中的unsigned char就是byte型，0-255值 unsigned char DQ1;H5LTget_attribute_uchar(file_id,/,Data Quality,DQ); /读取文件属性 printf(Data Quality:%dn,*DQ);H5LTset_attribute_uchar(ofile_id,/,Data Quality,DQ,1); /写入文件属性最后显示的数据类型是 8-bit unsigned integer7、读取文件属性float BB200;H5LTget_attribute_float(file_id,/,BB start time,BB); for(int i=0;i200;i+) printf(BB:%fn,BBi);二、将数据和信息写入HDF文件hid_t file_id;hid_t dataset_id;hid_t filespace_id;1、创建HDF文件file_id = H5Fcreate(OutFileName,H5F_ACC_TRUNC, H5P_DEFAULT,H5P_DEFAULT);2、设置空间维度属性filespace_id = H5Screate_simple(3,hdims,NULL);3、创建数据集（或组）dataset_id = H5Dcreate(file_id,datasetName,H5T_NATIVE_UINT,filespace_id,H5P_DEFAULT,H5P_DEFAULT,H5P_DEFAULT); 4、写入数据H5Dwrite(dataset_id,H5T_NATIVE_UINT,H5S_ALL,H5S_ALL,H5P_DEFAULT,pData);5、写文件属性H5LTset_attribute_string(ofile_id,/,AdditionalAnotation,AA);H5LTset_attribute_double(ofile_id,/,BB start time,BB,200);6、写数据集属性H5LTset_attribute_float(ofile_id,EV_1KM_RefSB,Intercept,Intercept,15);H5LTset_attribute_string(odset_id,/EV_1KM_RefSB,band_name,band_name);7、 数据压缩hid_t plist_id; / 创建数据集生成的性质列表plist_id = H5Pcreate(H5P_DATASET_CREATE); hsize_t cdims3; cdims0 = 1; cdims1 = 256;cdims2 = 2048; / / 创建分块的性质status = H5Pset_chunk(plist_id,3,cdims); / 设置压缩方法, zlib同gzipstatus = H5Pset_deflate(plist_id,5);压缩后，创建数据集时要用odset_id = H5Dcreate(ofile_id,EV_1K