下面的例子展示了如何新建一个HDF文件并写入数据，及打开已关闭的HDF文件再向其追加数据。
本例测试于Linux(Ubuntu 10.04)系统，MATLAB版本为2009b.Unix。

%------------------------------------------
% the program example demonstrates the process that create a HDF file and wirte some data set into it, furthermore, the steps that how to reopen the created HDF file and add some data set into it are also illustrated.
% written by WANG Zijun
% Jilin University, Email: wzj08@mails.jlu.edu.cn
%------------------------------------------

% process 1
%% create a hdf file and write two data set

% create a new HDF file named sat_atmo.hdf in current directory
sdId = hdfsd('start','sat_atmo.hdf','create');
% in my example the lat_sat is a 14*22 array
% lat_sat = ...
% define a scientific data set
dsName = 'lat_sat';
dsType = 'double';
dsRank = ndims(lat_sat);
dsDims = fliplr(size(lat_sat));
% and create it
sdsId = hdfsd('create', sdId, dsName, dsType,dsRank, dsDims);
% define the data set range
dsStart = zeros(1, ndims(lat_sat));
dsStride =[];
dsEdges = fliplr(size(lat_sat));
% and write the data set lat_sat into the HDF file
stat = hdfsd('writedata', sdsId, dsStart, dsStride, dsEdges, lat_sat);

% Add another data set lon_sat into the HDF file
dsName = 'lon_sat';
dsType = 'double';
dsRank = ndims(lon_sat);
dsDims = fliplr(size(lon_sat));
sdsId = hdfsd('create', sdId, dsName, dsType,dsRank, dsDims);

dsStart = zeros(1, ndims(lon_sat));
dsStride =[];
dsEdges = fliplr(size(lon_sat));
stat = hdfsd('writedata', sdsId, dsStart, dsStride, dsEdges, lon_sat);

% terminate the file access and close the file
stat = hdfsd('endaccess', sdsId);
stat = hdfsd('end', sdId)
hdfml('closeall');


% process 2
%% reopen the hdf file and add a data set
% in our example the wavelength is 14 *3 array
% wl = ...
sdId = hdfsd('start', 'sat_atmo.hdf', 'rdwr');
dsName = 'wavelegnth';
dsType = 'double';
dsRank = ndims(wl);
dsDims = fliplr(size(wl));
sdsId = hdfsd('create', sdId, dsName, dsType, dsRank, dsDims);
dsStart = zeros(1, ndims(wl));
dsStride = [];
dsEdges = fliplr(size(wl));
stat = hdfsd('writedata', sdsId, dsStart, dsStride, dsEdges, wl);

stat = hdfsd('endaccess', sdsId);
stat = hdfsd('end', sdId)

hdfml('closeall');


转载本文请联系原作者获取授权，同时请注明本文来自汪自军科学网博客。
链接地址：https://wap.sciencenet.cn/blog-43777-364662.html

上一篇：CWzjWAData类：c++数组与MATLAB数组相互转换
下一篇：基于MATLAB创建HDF之Vgroup及向存在Vgroup中追加数据