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
|
#include <stdlib.h>
#include "display.h"
#include "raster.h"
#include "gis.h"
#include "D.h"
static int cell_draw( char *, char *, struct Colors *, int);
/*!
* \brief
*
* If the map is a floating-point map, read the map using
* <tt>G_get_d_raster_row()</tt> and plot using <tt>D_draw_d_cell()</tt>. If the
* map is an integer map, read the map using <tt>G_get_c_raster_row()</tt> and
* plot using <tt>D_draw_cell()</tt>.
*
* \param name
* \param mapset
* \param overlay
* \return int
*/
int Dcell(char *name, char *mapset, int overlay)
{
struct Cell_head wind ;
struct Colors colors ;
char buff[128] ;
int offset ;
G_get_set_window (&wind) ;
if (D_check_map_window(&wind))
G_fatal_error("Setting map window") ;
if (G_set_window(&wind) == -1)
G_fatal_error("Current window not settable") ;
/* Get existing map window for this graphics window, or save window */
/* cell maps wipe out a picture, so we clear info on the window too */
if (!overlay && D_clear_window())
G_fatal_error("Can't clear current graphics window") ;
/* Save the current map window with the graphics window */
D_check_map_window(&wind) ;
G_set_window (&wind);
/* Get color offset value for current window and pass to driver */
D_offset_is(&offset) ;
R_color_offset(offset) ;
/* Set the colors for the display */
if (G_read_colors(name, mapset, &colors) == -1)
{
sprintf(buff,"Color file for [%s] not available", name) ;
G_fatal_error(buff) ;
}
D_set_colors (&colors);
/* Go draw the cell file */
cell_draw(name, mapset, &colors, overlay) ;
/* release the colors now */
G_free_colors (&colors);
/* record the cell file */
/* If overlay add it to the list instead of setting the cell name */
/* added 6/91 DBS @ CWU */
if (overlay) {
sprintf(buff,"d.rast -o map=%s", G_fully_qualified_name(name,mapset));
D_add_to_list(buff) ;
}
else {
D_set_cell_name(G_fully_qualified_name(name, mapset));
}
return 0;
}
/* I modified this function to read and draw raster cell as doubles */
static int cell_draw(char *name, char *mapset,
struct Colors *colors, int overlay)
{
char buff[128] ;
int cellfile ;
DCELL *xarray ;
int cur_A_row ;
int t, b, l, r ;
/* Set up the screen, conversions, and graphics */
D_get_screen_window(&t, &b, &l, &r) ;
if (D_cell_draw_setup(t, b, l, r))
{
sprintf(buff,"Cannot use current window") ;
G_fatal_error(buff) ;
}
D_set_overlay_mode(overlay);
/* Make sure map is available */
if ((cellfile = G_open_cell_old(name, mapset)) == -1)
{
sprintf(buff,"Not able to open cellfile for [%s]", name);
G_fatal_error(buff) ;
}
/* Allocate space for cell buffer */
xarray = G_allocate_d_raster_buf() ;
/* loop for array rows */
for (cur_A_row = 0; cur_A_row != -1; )
{
/* Get window (array) row currently required */
G_get_d_raster_row(cellfile, xarray, cur_A_row) ;
/* Draw the cell row, and get the next row number */
cur_A_row = D_draw_d_raster(cur_A_row, xarray, colors) ;
}
R_flush() ;
/* Wrap up and return */
G_close_cell(cellfile) ;
free (xarray);
return(0) ;
}
|
