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
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
|
/*!
* \file lib/raster/set_window.c
*
* \brief Raster Library - Set window (map region)
*
* (C) 2001-2009 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 Original author CERL
*/
#include <grass/gis.h>
#include <grass/raster.h>
#include <grass/glocale.h>
#include "../gis/G.h"
#include "R.h"
static void update_window_mappings(void);
static void check_write_window(void);
void Rast__init_window(void)
{
if (G_is_initialized(&R__.window_set))
return;
G__init_window();
R__.rd_window = G__.window;
R__.wr_window = G__.window;
R__.split_window = 0;
G_initialize_done(&R__.window_set);
}
/*!
* \brief Establishes 'window' as the current working window.
*
* \param window window to become operative window
*/
void Rast_set_window(struct Cell_head *window)
{
Rast__init();
if (R__.split_window)
G_warning(_("Rast_set_window() called while window split"));
check_write_window();
G_adjust_Cell_head(window, 0, 0);
R__.wr_window = *window;
R__.rd_window = *window;
R__.split_window = 0;
update_window_mappings();
}
/*!
\brief Unset current window
*/
void Rast_unset_window(void)
{
G_debug(4, "Rast_unset_window()");
R__.window_set = 0;
}
/*!
* \brief Establishes 'window' as the current working window for output.
*
* \param window window to become operative window
*/
void Rast_set_output_window(struct Cell_head *window)
{
Rast__init();
check_write_window();
G_adjust_Cell_head(window, 0, 0);
R__.wr_window = *window;
R__.split_window = 1;
G_set_window(window);
}
/*!
* \brief Establishes 'window' as the current working window for input.
*
* Any opened cell files has its file-to-window mapping reworked.
*
* \param window window to become operative window
*/
void Rast_set_input_window(struct Cell_head *window)
{
Rast__init();
G_adjust_Cell_head(window, 0, 0);
R__.rd_window = *window;
R__.split_window = 1;
update_window_mappings();
}
static void update_window_mappings(void)
{
int i;
int maskfd;
/* adjust window, check for valid window */
/* adjust the real one, not a copy
G_copy (&twindow, window, sizeof(struct Cell_head));
window = &twindow;
*/
/* except for MASK, cell files open for read must have same projection
* and zone as new window
*/
maskfd = R__.auto_mask > 0 ? R__.mask_fd : -1;
for (i = 0; i < R__.fileinfo_count; i++) {
struct fileinfo *fcb = &R__.fileinfo[i];
if (fcb->open_mode == OPEN_OLD) {
if (fcb->cellhd.zone == R__.rd_window.zone &&
fcb->cellhd.proj == R__.rd_window.proj)
continue;
if (i != maskfd)
G_fatal_error(_("Rast_set_read_window(): projection/zone "
"differs from that of "
"currently open raster maps"));
}
}
/* close the mask */
if (R__.auto_mask > 0) {
Rast_close(maskfd);
/* G_free (R__.mask_buf); */
R__.mask_fd = -1;
R__.auto_mask = -1; /* turn off masking */
}
/* 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 < R__.fileinfo_count; i++) {
struct fileinfo *fcb = &R__.fileinfo[i];
if (fcb->open_mode != OPEN_OLD &&
fcb->open_mode != OPEN_NEW_UNCOMPRESSED &&
fcb->open_mode != OPEN_NEW_COMPRESSED)
continue;
if (fcb->open_mode == OPEN_OLD)
G_fatal_error(_("Input window changed while maps are open for "
"read. Map name <%s>"),
fcb->name);
}
/* turn masking (back) on if necessary */
Rast__check_for_auto_masking();
}
static void check_write_window(void)
{
int i;
for (i = 0; i < R__.fileinfo_count; i++) {
struct fileinfo *fcb = &R__.fileinfo[i];
if (fcb->open_mode == OPEN_NEW_UNCOMPRESSED ||
fcb->open_mode == OPEN_NEW_COMPRESSED)
G_fatal_error(_("Output window changed while maps are open for "
"write. Map name <%s>"),
fcb->name);
}
}
|
