5#include <grass/raster.h>
6#include <grass/raster.h>
14static int icmp(
const void *aa,
const void *bb)
22static int fcmp(
const void *aa,
const void *bb)
34static int dcmp(
const void *aa,
const void *bb)
46int f_nmedian(
int argc,
const int *argt,
void **args)
50 int size = argc * Rast_cell_size(argt[0]);
56 for (i = 1; i <= argc; i++)
57 if (argt[i] != argt[0])
62 array = G_realloc(array, size);
69 CELL **argv = (CELL **) &args[1];
77 for (j = 0; j < argc; j++) {
78 if (IS_NULL_C(&argv[j][i]))
88 qsort(a, n,
sizeof(CELL), icmp);
93 *resc = (*resc + a1) / 2;
102 FCELL *res = args[0];
103 FCELL **argv = (FCELL **) &args[1];
108 for (i = 0; i <
columns; i++) {
111 for (j = 0; j < argc; j++) {
112 if (IS_NULL_F(&argv[j][i]))
122 qsort(a, n,
sizeof(FCELL), fcmp);
127 *resc = (*resc + a1) / 2;
136 DCELL *res = args[0];
137 DCELL **argv = (DCELL **) &args[1];
142 for (i = 0; i <
columns; i++) {
145 for (j = 0; j < argc; j++) {
146 if (IS_NULL_D(&argv[j][i]))
156 qsort(a, n,
sizeof(DCELL), dcmp);
161 *resc = (*resc + a1) / 2;