7 #include <grass/glocale.h>
15 char buf[200], buf2[200], xname[512], xmapset[512];
20 G3d_error(_(
"G3d_openCellOldNoHeader: error in G3d_maskOpenOld"));
26 G3d_error(_(
"G3d_openCellOldNoHeader: error in G3d_malloc"));
31 sprintf(buf,
"%s/%s", G3D_DIRECTORY, xname);
32 sprintf(buf2,
"%s@%s", G3D_CELL_ELEMENT, xmapset);
36 sprintf(buf,
"%s/%s", G3D_DIRECTORY, name);
37 sprintf(buf2,
"%s", G3D_CELL_ELEMENT);
44 if (map->data_fd < 0) {
45 G3d_error(_(
"G3d_openCellOldNoHeader: error in G_open_old"));
86 G3D_Region *
window,
int typeIntern,
int cache)
91 int rows,
cols, depths, precision;
92 double ew_res, ns_res, tb_res;
93 int nofHeaderBytes, dataOffset, useXdr, hasIndex;
95 double north, south, east, west, top, bottom;
99 G3d_error(_(
"G3d_openCellOld: error in G3d_openCellOldNoHeader"));
103 if (lseek(map->data_fd, (
long)0, SEEK_SET) == -1) {
104 G3d_error(_(
"G3d_openCellOld: can't rewind file"));
110 &north, &south, &east, &west, &top, &bottom,
111 &rows, &cols, &depths,
112 &ew_res, &ns_res, &tb_res,
113 &tileX, &tileY, &tileZ,
114 &type, &compression, &useRle, &useLzw,
115 &precision, &dataOffset, &useXdr, &hasIndex, &unit)) {
116 G3d_error(_(
"G3d_openCellOld: error in G3d_readHeader"));
120 if (window == G3D_DEFAULT_WINDOW)
123 if (proj != window->proj) {
124 G3d_error(_(
"G3d_openCellOld: projection does not match window projection"));
127 if (zone != window->zone) {
128 G3d_error(_(
"G3d_openCellOld: zone does not match window zone"));
132 map->useXdr = useXdr;
137 &(map->indexLongNbytes), 1)) ||
139 &(map->indexNbytesUsed), 1))) {
140 G3d_error(_(
"G3d_openCellOld: can't read header"));
145 if (map->indexNbytesUsed >
sizeof(
long))
146 G3d_fatalError(_(
"G3d_openCellOld: index does not fit into long"));
150 G3d_error(_(
"G3d_openCellOld: error in G3d_malloc"));
155 if (read(map->data_fd, ltmp, map->indexLongNbytes) !=
156 map->indexLongNbytes) {
157 G3d_error(_(
"G3d_openCellOld: can't read header"));
160 G3d_longDecode(ltmp, &(map->indexOffset), 1, map->indexLongNbytes);
164 nofHeaderBytes = dataOffset;
166 if (typeIntern == G3D_TILE_SAME_AS_FILE)
170 type, precision, cache,
171 hasIndex, map->useXdr, typeIntern,
172 nofHeaderBytes, tileX, tileY, tileZ,
174 north, south, east, west, top, bottom,
175 rows, cols, depths, ew_res, ns_res, tb_res, unit)) {
176 G3d_error(_(
"G3d_openCellOld: error in G3d_fillHeader"));
217 int nofHeaderBytes, dummy = 0,
compression, precision;
219 char xname[512], xmapset[512];
223 G3d_error(_(
"G3d_openCellNew: error in G3d_maskOpenOld"));
232 G3d_error(_(
"G3d_openCellNew: error in G3d_malloc"));
237 map->fileName =
G_store(xname);
243 map->data_fd = open(map->tempName, O_RDWR | O_CREAT | O_TRUNC, 0666);
244 if (map->data_fd < 0) {
245 G3d_error(_(
"G3d_openCellNew: could not open file"));
256 else if (precision < -1)
259 else if (precision > 52)
261 else if (precision < -1)
265 if ((typeIntern == FCELL_TYPE) && (
g3d_file_type == DCELL_TYPE)) {
269 precision =
G3D_MIN(precision, 23);
273 precision = G3D_MAX_PRECISION;
278 map->indexLongNbytes =
sizeof(long);
288 &(map->indexLongNbytes), 1)) ||
290 G3d_error(_(
"G3d_openCellNew: can't write header"));
293 if (write(map->data_fd, &ldummy, map->indexLongNbytes) !=
294 map->indexLongNbytes) {
295 G3d_error(_(
"G3d_openCellNew: can't write header"));
301 nofHeaderBytes = lseek(map->data_fd, (
long)0, SEEK_CUR);
309 map->useXdr, typeIntern, nofHeaderBytes,
312 region->proj, region->zone,
313 region->north, region->south, region->east,
314 region->west, region->top, region->bottom,
315 region->rows, region->cols, region->depths,
316 region->ew_res, region->ns_res, region->tb_res,
318 G3d_error(_(
"G3d_openCellNew: error in G3d_fillHeader"));