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_median(
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];
71 CELL *a1 = &a[(argc - 1) / 2];
72 CELL *a2 = &a[argc / 2];
77 for (j = 0; j < argc && !nv; j++) {
78 if (IS_NULL_C(&argv[j][i]))
87 qsort(a, argc,
sizeof(CELL), icmp);
88 res[i] = (*a1 + *a2) / 2;
97 FCELL **argv = (FCELL **) &args[1];
99 FCELL *a1 = &a[(argc - 1) / 2];
100 FCELL *a2 = &a[argc / 2];
102 for (i = 0; i <
columns; i++) {
105 for (j = 0; j < argc && !nv; j++) {
106 if (IS_NULL_F(&argv[j][i]))
115 qsort(a, argc,
sizeof(FCELL), fcmp);
116 res[i] = (*a1 + *a2) / 2;
124 DCELL *res = args[0];
125 DCELL **argv = (DCELL **) &args[1];
127 DCELL *a1 = &a[(argc - 1) / 2];
128 DCELL *a2 = &a[argc / 2];
130 for (i = 0; i <
columns; i++) {
133 for (j = 0; j < argc && !nv; j++) {
134 if (IS_NULL_D(&argv[j][i]))
143 qsort(a, argc,
sizeof(DCELL), dcmp);
144 res[i] = (*a1 + *a2) / 2;