← 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_sum.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_sum.c
1
1
#include <stdlib.h>
2
2
#include <unistd.h>
3
3
#include <grass/gis.h>
 
4
#include <grass/raster.h>
4
5
#include <grass/glocale.h>
5
6
#include "method.h"
6
7
 
7
 
#define STATS "r.stats"
8
 
#define RECLASS "r.reclass"
9
 
 
10
8
/* function prototypes */
11
9
static void sum_out(FILE *, long, double);
12
10
 
13
 
 
14
 
int
15
 
o_sum(char *basemap, char *covermap, char *outputmap, int usecats,
16
 
      struct Categories *cats)
 
11
int o_sum(const char *basemap, const char *covermap, const char *outputmap,
 
12
          int usecats, struct Categories *cats)
17
13
{
18
 
    char *me = "o_sum";
19
 
    char command[1024];
20
 
 
21
14
    long catb, basecat, covercat;
22
15
    double x, area, sum1;
23
16
    int stat;
24
 
    char *tempfile1, *tempfile2;
25
 
    FILE *fd1, *fd2;
26
 
 
27
 
    tempfile1 = G_tempfile();
28
 
    tempfile2 = G_tempfile();
29
 
 
30
 
    sprintf(command, "%s -cn input=\"%s,%s\" fs=space output=\"%s\"", STATS, basemap,
31
 
            covermap, tempfile1);
32
 
 
33
 
    if (stat = G_system(command)) {
34
 
        unlink(tempfile1);
35
 
        G_fatal_error(_("%s: running %s command"), me, STATS);
36
 
    }
37
 
 
38
 
    fd1 = fopen(tempfile1, "r");
39
 
    fd2 = fopen(tempfile2, "w");
40
 
    if (fd1 == NULL || fd2 == NULL) {
41
 
        unlink(tempfile1);
42
 
        unlink(tempfile2);
43
 
        G_fatal_error(_("%s: unable to open temporary file"), me);
44
 
    }
45
 
    sum_out(fd2, 0L, 0.0);      /* force at least one reclass rule */
 
17
    struct Popen stats_child, reclass_child;
 
18
    FILE *stats, *reclass;
 
19
 
 
20
    stats = run_stats(&stats_child, basemap, covermap, "-cn");
 
21
    reclass = run_reclass(&reclass_child, basemap, outputmap);
 
22
 
 
23
    sum_out(reclass, 0L, 0.0);  /* force at least one reclass rule */
46
24
 
47
25
    catb = 0;
48
26
    sum1 = 0.0;
49
27
 
50
 
 
51
 
    while (fscanf(fd1, "%ld %ld %lf", &basecat, &covercat, &area) == 3) {
 
28
    while (fscanf(stats, "%ld %ld %lf", &basecat, &covercat, &area) == 3) {
52
29
        if (catb != basecat) {
53
 
            sum_out(fd2, catb, sum1);
 
30
            sum_out(reclass, catb, sum1);
54
31
            sum1 = 0.0;
55
32
            catb = basecat;
56
33
        }
57
34
        if (usecats)
58
 
            sscanf(G_get_cat((CELL) covercat, cats), "%lf", &x);
 
35
            sscanf(Rast_get_c_cat((CELL *) &covercat, cats), "%lf", &x);
59
36
        else
60
37
            x = covercat;
61
38
        sum1 += x * area;
62
39
        /*        fprintf(stderr,"sum: %d\n",(int)sum1); */
63
40
 
64
41
    }
65
 
    sum_out(fd2, basecat, sum1);
66
 
    fclose(fd1);
67
 
    fclose(fd2);
68
 
    sprintf(command, "%s input=\"%s\" output=\"%s\" rules=\"%s\"",
69
 
            RECLASS, basemap, outputmap, tempfile2);
70
 
    stat = G_system(command);
71
 
    unlink(tempfile1);
72
 
    unlink(tempfile2);
73
 
 
74
 
    return (stat);
 
42
    sum_out(reclass, basecat, sum1);
 
43
 
 
44
    G_popen_close(&stats_child);
 
45
    G_popen_close(&reclass_child);
 
46
 
 
47
    return stat;
75
48
}
76
49
 
77
 
 
78
 
static void sum_out(FILE * fd, long cat, double sum1)
 
50
static void sum_out(FILE *fp, long cat, double sum1)
79
51
{
80
52
    char buf[64];
81
53
 
86
58
        G_trim_decimal(buf);
87
59
    }
88
60
 
89
 
    fprintf(fd, "%ld = %ld %s\n", cat, cat, buf);
 
61
    fprintf(fp, "%ld = %ld %s\n", cat, cat, buf);
90
62
}
 
63