← Back to branch summary

~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to raster/r.statistics/o_kurt.c

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
raster/r.statistics/o_kurt.c
3
3
#include <stdlib.h>
4
4
#include <unistd.h>
5
5
#include <grass/gis.h>
 
6
#include <grass/raster.h>
6
7
#include <grass/glocale.h>
7
8
#include "method.h"
8
9
 
12
13
static int kurt(double *, int, double *);
13
14
 
14
15
 
15
 
int
16
 
o_kurt(char *basemap, char *covermap, char *outputmap, int usecats,
17
 
       struct Categories *cats)
 
16
int o_kurt(const char *basemap, const char *covermap, const char *outputmap,
 
17
           int usecats, struct Categories *cats)
18
18
{
19
 
    char command[1024];
 
19
    struct Popen stats_child, reclass_child;
20
20
    FILE *stats, *reclass;
21
21
    int first, mem, i, count;
22
22
    long basecat, covercat, catb, catc;
26
26
    mem = MEM * sizeof(double);
27
27
    tab = (double *)G_malloc(mem);
28
28
 
29
 
    sprintf(command, "r.stats -cn input=\"%s,%s\" fs=space", basemap,
30
 
            covermap);
31
 
    stats = popen(command, "r");
32
 
 
33
 
    sprintf(command, "r.reclass i=\"%s\" o=\"%s\"", basemap, outputmap);
34
 
    reclass = popen(command, "w");
 
29
    stats = run_stats(&stats_child, basemap, covermap, "-cn");
 
30
    reclass = run_reclass(&reclass_child, basemap, outputmap);
35
31
 
36
32
    first = 1;
37
33
    while (read_stats(stats, &basecat, &covercat, &value)) {
52
48
        }
53
49
 
54
50
        if (usecats)
55
 
            sscanf(G_get_cat((CELL) covercat, cats), "%lf", &x);
 
51
            sscanf(Rast_get_c_cat((CELL *) &covercat, cats), "%lf", &x);
56
52
        else
57
53
            x = covercat;
58
54
 
73
69
    kurt(tab, count, &var);
74
70
    fprintf(reclass, "%ld = %ld %f\n", catb, catb, var);
75
71
 
76
 
 
77
 
    pclose(stats);
78
 
    pclose(reclass);
79
 
 
80
 
    return (0);
 
72
    G_popen_close(&stats_child);
 
73
    G_popen_close(&reclass_child);
 
74
 
 
75
    return 0;
81
76
}
82
77
 
83
78