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
|
/****************************************************************************
*
* MODULE: r.buffer
*
* AUTHOR(S): Michael Shapiro - CERL
*
* PURPOSE: This program creates distance zones from non-zero
* cells in a grid layer. Distances are specified in
* meters (on the command-line). Window does not have to
* have square cells. Works both for planimetric
* (UTM, State Plane) and lat-long.
*
* COPYRIGHT: (C) 2005 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.
*
****************************************************************************/
#include <stdlib.h>
#include "distance.h"
#include <grass/raster.h>
#include <grass/glocale.h>
/* read the input map. convert non-nulls to 1 */
int read_input_map(char *input, char *mapset, int ZEROFLAG)
{
int fd;
int row;
int hit;
register int col;
register CELL *cell;
register MAPTYPE *ptr;
map = (MAPTYPE *)G_malloc((size_t)window.rows * window.cols *
sizeof(MAPTYPE));
fd = Rast_open_old(input, mapset);
cell = Rast_allocate_c_buf();
ptr = map;
minrow = -1;
maxrow = -1;
mincol = window.cols;
maxcol = 0;
G_message(_("Reading input raster map <%s>..."),
G_fully_qualified_name(input, mapset));
count_rows_with_data = 0;
for (row = 0; row < window.rows; row++) {
hit = 0;
G_percent(row, window.rows, 2);
Rast_get_c_row(fd, cell, row);
for (col = 0; col < window.cols; col++) {
if (ZEROFLAG) {
if ((*ptr++ = (*cell++ != 0))) {
if (minrow < 0)
minrow = row;
maxrow = row;
if (col < mincol)
mincol = col;
if (col > maxcol)
maxcol = col;
if (!hit) {
count_rows_with_data++;
hit = 1;
}
}
}
else { /* use NULL */
if ((*ptr++ = !Rast_is_c_null_value(cell++))) {
if (minrow < 0)
minrow = row;
maxrow = row;
if (col < mincol)
mincol = col;
if (col > maxcol)
maxcol = col;
if (!hit) {
count_rows_with_data++;
hit = 1;
}
}
}
}
cell -= window.cols;
}
G_percent(row, window.rows, 2);
Rast_close(fd);
G_free(cell);
return 0;
}
|
