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
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
|
/**
* \file set_window.c
*
* \brief Set window
*
* (C) 2001-2008 by the GRASS Development Team
*
* This program is free software under the GNU General Public License
* (>=v2). Read the file COPYING that comes with GRASS for details.
*
* \author GRASS GIS Development Team
*
* \date 1999-2007
*/
#include <grass/gis.h>
#include <grass/glocale.h>
#include "G.h"
/**
* \brief Get the current working window.
*
* The current working window values are returned in the structure
* 'window'.
*
* \param[out] window window structure to be set
* \return 1
*/
int G_get_set_window(struct Cell_head *window)
{
G__init_window();
G_copy((char *)window, (char *)&G__.window, sizeof(*window));
return 1;
}
/**
* \brief Establishes 'window' as the current working window.
*
* Any opened cell files has its file-to-window mapping reworked.
*
* \param[in] window window to become operative window
*
* \return -1 on error
* \return 1 on success
*/
int G_set_window(struct Cell_head *window)
{
int i;
int maskfd;
char *err;
/* adjust window, check for valid window */
/* adjust the real one, not a copy
G_copy (&twindow, window, sizeof(struct Cell_head));
window = &twindow;
*/
if ((err = G_adjust_Cell_head(window, 0, 0))) {
G_warning("G_set_window(): %s", err);
return -1;
}
/* except for MASK, cell files open for read must have same projection
* and zone as new window
*/
maskfd = G__.auto_mask > 0 ? G__.mask_fd : -1;
for (i = 0; i < G__.fileinfo_count; i++) {
if (G__.fileinfo[i].open_mode == OPEN_OLD) {
if (G__.fileinfo[i].cellhd.zone == window->zone &&
G__.fileinfo[i].cellhd.proj == window->proj)
continue;
if (i != maskfd) {
G_warning(_("G_set_window(): projection/zone differs from that of "
"currently open raster maps"));
return -1;
}
}
}
/* close the mask */
if (G__.auto_mask > 0) {
G_close_cell(maskfd);
/* G_free (G__.mask_buf); */
G__.mask_fd = -1;
G__.auto_mask = -1; /* turn off masking */
}
/* copy the window to the current window */
G_copy((char *)&G__.window, (char *)window, sizeof(*window));
G__.window_set = 1;
/* now for each possible open cell file, recreate the window mapping */
/*
* also the memory for reading and writing must be reallocated for all opened
* cell files
*/
for (i = 0; i < G__.fileinfo_count; i++) {
if (G__.fileinfo[i].open_mode != OPEN_OLD &&
G__.fileinfo[i].open_mode != OPEN_NEW_UNCOMPRESSED &&
G__.fileinfo[i].open_mode != OPEN_NEW_COMPRESSED &&
G__.fileinfo[i].open_mode != OPEN_NEW_RANDOM)
continue;
if (G__.fileinfo[i].open_mode == OPEN_OLD)
G__create_window_mapping(i);
/* code commented 10/1999 due to problems */
/* else */
/* opened for writing */
/* {
G_free (G__.fileinfo[i].data);
G__.fileinfo[i].data = (unsigned char *) G_calloc (G__.window.cols,
G_raster_size(G__.fileinfo[i].map_type));
}
*/
/* allocate null bitstream buffers for reading/writing null rows */
/* for (j=0;j< NULL_ROWS_INMEM; j++)
{
G_free (G__.fileinfo[i].NULL_ROWS[j]);
G__.fileinfo[i].NULL_ROWS[j] = G__allocate_null_bits(G__.window.cols);
}
*/
/* initialize : no NULL rows in memory */
/* G__.fileinfo[i].min_null_row = (-1) * NULL_ROWS_INMEM;
if(G__.fileinfo[i].null_cur_row > 0)
{
G_warning(
"Calling G_set_window() in the middle of writing map %s",
G__.fileinfo[i].name);
G__.fileinfo[i].null_cur_row = 0;
}
*/
}
/* turn masking (back) on if necessary */
G__check_for_auto_masking();
/* reallocate/enlarge the G__. buffers for reading raster maps */
G__reallocate_null_buf();
G__reallocate_mask_buf();
G__reallocate_temp_buf();
G__reallocate_work_buf(sizeof(DCELL));
G__reallocate_work_buf(XDR_DOUBLE_NBYTES);
/* we want the number of bytes per cell to be maximum
so that there is enough memory for reading and writing rows */
return 1;
}
|
