2 #include <grass/raster.h>
4 void c_kurt(DCELL * result, DCELL * values,
int n,
const void *closure)
6 DCELL sum, ave, sumsq, sumqt, var;
13 for (i = 0; i < n; i++) {
14 if (Rast_is_d_null_value(&values[i]))
22 Rast_set_d_null_value(result, 1);
31 for (i = 0; i < n; i++) {
34 if (Rast_is_d_null_value(&values[i]))
39 sumqt += d * d * d * d;
44 *result = sumqt / (
count * var * var) - 3;
47 void w_kurt(DCELL * result, DCELL(*values)[2],
int n,
const void *closure)
49 DCELL sum, ave, sumsq, sumqt, var;
56 for (i = 0; i < n; i++) {
57 if (Rast_is_d_null_value(&values[i][0]))
60 sum += values[i][0] * values[i][1];
61 count += values[i][1];
65 Rast_set_d_null_value(result, 1);
74 for (i = 0; i < n; i++) {
77 if (Rast_is_d_null_value(&values[i][0]))
80 d = values[i][0] - ave;
81 sumsq += d * d * values[i][1];
82 sumqt += d * d * d * values[i][1];
87 *result = sumqt / (
count * var * var) - 3;
void c_kurt(DCELL *result, DCELL *values, int n, const void *closure)
void w_kurt(DCELL *result, DCELL(*values)[2], int n, const void *closure)