1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
#include <octave/oct.h>
#include <octave/ov-struct.h>
#include "hdf5.h"
// this special treatment is necessary because Win32-Octave ships with a very old hdf5 version (1.6.10)
void CloseH5Object(hid_t obj)
{
H5Oclose(obj);
}
DEFUN_DLD (h5readatt_octave, args, nargout, "h5readatt_octave(<File_Name>,<DataSet_Name>,<Attribute_Name>)")
{
octave_value retval;
int nargin = args.length();
if (nargin != 3)
{
print_usage();
return retval;
}
if ((args(0).is_string()==false) || (args(1).is_string()==false) || (args(2).is_string()==false))
{
print_usage();
return retval;
}
//suppress hdf5 error output
H5Eset_auto2(H5E_DEFAULT, NULL, NULL);
hid_t file = H5Fopen( args(0).string_value().c_str(), H5F_ACC_RDONLY, H5P_DEFAULT );
if (file==-1)
{
error("h5readatt_octave: opening the given File failed");
return retval;
}
hid_t obj = H5Oopen(file, args(1).string_value().c_str(), H5P_DEFAULT);
if (obj==-1)
{
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: opening the given Object failed");
return retval;
}
hid_t attr = H5Aopen_by_name(obj, ".", args(2).string_value().c_str(), H5P_DEFAULT, H5P_DEFAULT);
if (attr==-1)
{
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: opening the given Attribute failed");
return retval;
}
hid_t type = H5Aget_type(attr);
if (type<0)
{
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: dataset type error");
return retval;
}
if (H5Tget_class(type)!=H5T_FLOAT)
{
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: attribute type not supported");
return retval;
}
size_t numVal = H5Aget_storage_size(attr)/H5Tget_size(type);
double value[numVal];
if (H5Tget_size(type)==sizeof(float))
{
float f_value[numVal];
if (H5Aread(attr, H5T_NATIVE_FLOAT, f_value)<0)
{
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: reading the given Attribute failed");
return retval;
}
for (size_t n=0;n<numVal;++n)
value[n] = f_value[n];
}
else if (H5Tget_size(type)==sizeof(double))
{
if (H5Aread(attr, H5T_NATIVE_DOUBLE, value)<0)
{
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: reading the given Attribute failed");
return retval;
}
}
else
{
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
error("h5readatt_octave: reading the given Attribute failed: unknown type");
return retval;
}
H5Aclose(attr);
CloseH5Object(obj);
H5Fclose(file);
Matrix mat(numVal,1);
for (size_t n=0;n<numVal;++n)
mat(n)=value[n];
retval = octave_value(mat);
return retval;
}
|
