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
121
122
123
124
125
126
127
|
/*
----------------------------------------------------------------------------
| Copyright (C) 1999 Emergent IT Inc. and Raytheon Systems Company |
----------------------------------------------------------------------------
*/
#include <HE5_config.h>
#include <HE5_HdfEosDef.h>
/*
* In this example we will (1) open the "Grid.h5" HDF-EOS file, (2) attach to
* the "Grid1" grid, and (3) define the grid fields.
*/
int main()
{
herr_t status = FAIL;
hid_t gdfid = FAIL;
hid_t GDid1 = FAIL;
hid_t GDid2 = FAIL;
hid_t GDid3 = FAIL;
float fillval1 = -7.;
float fillval2 = -9999.;
int tilerank, compcode = 16;
int compparm[5] = {0, 0, 0, 0, 0};
hsize_t tiledims[8];
/*
* We first open the HDF-EOS grid file, "Grid.h5". Because this file
* already exist and we wish to write to it, we use the H5F_ACC_RDWR access
* code in the open statement. The HE5_GDopen routine returns the grid file
* id, gdfid, which is used to identify the file in subsequent routines.
*/
gdfid = HE5_GDopen("Grid.h5", H5F_ACC_RDWR);
/*
* If the grid file cannot be found, HE5_GDopen will return -1 for the file
* handle (gdfid). We there check that this is not the case before
* proceeding with the other routines.
*
* The HE5_GDattach routine returns the handle to the existing grid "Grid1",
* GDid. If the grid is not found, HE5_GDattach returns -1 for the handle.
*/
if (gdfid != FAIL)
{
GDid1 = HE5_GDattach(gdfid, "UTMGrid");
status = HE5_GDsetfillvalue(GDid1, "Pollution", H5T_NATIVE_FLOAT, &fillval1);
printf("Status returned by HE5_GDsetfillvalue(...,\"Pollusion\",...) : %d \n", status);
status = HE5_GDdeffield(GDid1, "Pollution", "Time,YDim,XDim", NULL, H5T_NATIVE_FLOAT, 0);
printf("Status returned by HE5_GDdeffield(..., \"Pollution\",...) : %d \n", status);
/* Define appendable field "Vegetation" */
tilerank = 2;
tiledims[0] = 100;
tiledims[1] = 60;
compparm[0] = 8;
status = HE5_GDdeftile(GDid1,HE5_HDFE_TILE,tilerank,tiledims);
printf("Status returned by HE5_GDdeftile(): %d\n", status);
#ifdef HAVE_FILTER_SZIP_ENCODER
status = HE5_GDdefcomp(GDid1,compcode, compparm);
printf("Status returned by HE5_GDdefcomp(): %d\n", status);
#else
printf("SZIP encoder is not available. The data compression test is skipped.\n");
#endif
status = HE5_GDdeffield(GDid1, "Vegetation", "YDim,XDim",NULL,H5T_NATIVE_FLOAT, 0);
printf("Status returned by HE5_GDdeffield(\"Vegetation\"): %d\n", status);
GDid2 = HE5_GDattach(gdfid, "PolarGrid");
status = HE5_GDdeffield(GDid2, "Temperature", "YDim,XDim", NULL, H5T_NATIVE_FLOAT, 0);
printf("Status returned by HE5_GDdeffield(..., \"Temperature\",...) : %d \n", status);
status = HE5_GDsetfillvalue(GDid2, "Pressure", H5T_NATIVE_FLOAT, &fillval2);
printf("Status returned by HE5_GDsetfillvalue(...,\"Pressure\",...) : %d \n", status);
status = HE5_GDdeffield(GDid2, "Pressure", "YDim,XDim", NULL, H5T_NATIVE_INT, 0);
printf("Status returned by HE5_GDdeffield(..., \"Pressure\",...) : %d \n", status);
status = HE5_GDdeffield(GDid2, "Soil Dryness", "YDim,XDim", NULL, H5T_NATIVE_FLOAT, 0);
printf("Status returned by HE5_GDdeffield(..., \"Soil Dryness\",...) : %d \n", status);
status = HE5_GDdeffield(GDid2, "Spectra", "Bands,YDim,XDim", NULL, H5T_NATIVE_DOUBLE, 0);
printf("Status returned by HE5_GDdeffield(..., \"Spectra\",...) : %d \n", status);
status = HE5_GDdetach(GDid1);
status = HE5_GDdetach(GDid2);
GDid3 = HE5_GDattach(gdfid, "GEOGrid");
if (GDid3 == -1)
{
printf("\t\tError: Cannot attach to grid \"GEOGrid\"\n");
return -1;
}
status = HE5_GDsetfillvalue(GDid3, "GeoSpectra", H5T_NATIVE_FLOAT, &fillval2);
if (status == -1)
{
printf("\t\tError: Cannot define fill value for field \"Pressure\"\n");
return -1;
}
status = HE5_GDdeffield(GDid3, "GeoSpectra", "YDim,XDim", NULL, H5T_NATIVE_FLOAT, 0);
if (status == -1)
{
printf("\t\tError: Cannot define field \"GeoSpectra\"\n");
return -1;
}
HE5_GDdetach(GDid3);
}
status = HE5_GDclose(gdfid);
return 0;
}
|
