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
|
/*
*******************************************************************
* G_set_window (window)
* Cell_head *window window to become operative window
*
* Establishes 'window' as the current working window. Any opened
* cell files has its file-to-window mapping reworked.
*
*******************************************************************
* G_get_set_window (window)
* Cell_head *window
*
* The current working window values are returned in the structure
* 'window'.
*******************************************************************/
#include "gis.h"
#include "glocale.h"
#include "G.h"
int G_get_set_window (struct Cell_head *window)
{
G__init_window() ;
G_copy((char *) window, (char *) &G__.window, sizeof(*window) ) ;
return 1;
}
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 < MAXFILES; 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 files"));
return -1;
}
}
}
/* close the mask */
if (G__.auto_mask > 0)
{
G_close_cell (maskfd);
/* 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 < MAXFILES; 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 */
/* {
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++)
{
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;
}
|
