25 #include <grass/gis.h>
26 #include <grass/raster.h>
27 #include <grass/bitmap.h>
28 #include <grass/linkm.h>
29 #include <grass/interpf.h>
30 #include <grass/glocale.h>
42 int i, j, cfmask = -1, irev, MASKfd;
44 CELL *cellmask, *MASK;
47 if ((MASKfd = Rast_maskfd()) >= 0)
48 MASK = Rast_allocate_c_buf();
61 cellmask = Rast_allocate_c_buf();
62 cfmask = Rast_open_old(params->
maskmap, mapsetm);
67 for (i = 0; i < params->
nsizr; i++) {
68 irev = params->
nsizr - i - 1;
70 Rast_get_c_row(cfmask, cellmask, i);
72 Rast_get_c_row(MASKfd, MASK, i);
73 for (j = 0; j < params->
nsizc; j++) {
74 if ((cellmask && (cellmask[j] == 0 || Rast_is_c_null_value(&cellmask[j]))) ||
75 (MASK && (MASK[j] == 0 || Rast_is_c_null_value(&MASK[j]))))
76 BM_set(bitmask, j, irev, 0);
78 BM_set(bitmask, j, irev, 1);
94 double numbery,
double numberz,
int n_leafs)
108 for (ii = 0; ii < n_leafs; ii++)