int o_divr(const char *basemap, const char *covermap, const char *outputmap,
int usecats, struct Categories *cats)
{
struct Popen stats_child, reclass_child;
FILE *stats_fp, *reclass_fp;
int first;
long basecat, covercat, catb, catc;
double area;
stats_fp = run_stats(&stats_child, basemap, covermap, "-an");
reclass_fp = run_reclass(&reclass_child, basemap, outputmap);
first = 1;
while (read_stats(stats_fp, &basecat, &covercat, &area)) {
if (first) {
first = 0;
catb = basecat;
catc = 0;
}
if (basecat != catb) {
write_reclass(reclass_fp, catb, catc, Rast_get_c_cat((CELL *) &catc, cats),
usecats);
catb = basecat;
catc = 0;
}
catc++;
}
if (!first)
write_reclass(reclass_fp, catb, catc, Rast_get_c_cat((CELL *) &catc, cats), usecats);
G_popen_close(&stats_child);
G_popen_close(&reclass_child);
return 0;
}
int
o_mode(const char *basemap, const char *covermap, const char *outputmap, int usecats,
struct Categories *cats)
{
struct Popen stats_child, reclass_child;
FILE *stats, *reclass;
int first;
long basecat, covercat, catb, catc;
double value, max;
stats = run_stats(&stats_child, basemap, covermap, "-an");
reclass = run_reclass(&reclass_child, basemap, outputmap);
first = 1;
while (read_stats(stats, &basecat, &covercat, &value)) {
if (first) {
first = 0;
catb = basecat;
catc = covercat;
max = value;
}
if (basecat != catb) {
write_reclass(reclass, catb, catc, Rast_get_c_cat((CELL *) &catc, cats),
usecats);
catb = basecat;
catc = covercat;
max = value;
}
if (value > max) {
catc = covercat;
max = value;
}
}
if (first) {
catb = catc = 0;
}
write_reclass(reclass, catb, catc, Rast_get_c_cat((CELL *) &catc, cats), usecats);
G_popen_close(&stats_child);
G_popen_close(&reclass_child);
return 0;
}
int main(int argc, char *argv[])
{
struct GModule *module;
struct
{
struct Option *base, *cover, *output;
} parm;
char *basemap, *base_mapset;
char *covermap, *cover_mapset;
char *outmap;
char command[1024];
struct Categories cover_cats;
FILE *stats_fd, *reclass_fd;
int first;
long basecat, covercat, catb, catc;
double area;
struct stats stats;
G_gisinit(argv[0]);
module = G_define_module();
module->keywords = _("raster, statistics");
module->description =
_("Finds the median of values in a cover map within "
"areas assigned the same category value in a "
"user-specified base map.");
parm.base = G_define_standard_option(G_OPT_R_INPUT);
parm.base->key = "base";
parm.base->description = _("Name of base raster map");
parm.cover = G_define_standard_option(G_OPT_R_INPUT);
parm.cover->key = "cover";
parm.cover->description = _("Name of cover raster map");
parm.output = G_define_standard_option(G_OPT_R_OUTPUT);
if (G_parser(argc, argv))
exit(EXIT_FAILURE);
basemap = parm.base->answer;
covermap = parm.cover->answer;
outmap = parm.output->answer;
base_mapset = G_find_cell2(basemap, "");
if (base_mapset == NULL)
G_fatal_error(_("Base raster map <%s> not found"), basemap);
cover_mapset = G_find_cell2(covermap, "");
if (cover_mapset == NULL)
G_fatal_error(_("Raster map <%s> not found"), covermap);
if (G_legal_filename(outmap) < 0)
G_fatal_error(_("<%s> is an illegal file name"), outmap);
if (strcmp(G_mapset(), base_mapset) == 0 && strcmp(basemap, outmap) == 0)
G_fatal_error(_("Base map and output map <%s> must be different"),
outmap);
if (G_read_cats(covermap, cover_mapset, &cover_cats) < 0)
G_fatal_error(_("Unable to read category labels of raster map <%s>"),
covermap);
strcpy(command, "r.stats -an \"");
strcat(command, G_fully_qualified_name(basemap, base_mapset));
strcat(command, ",");
strcat(command, G_fully_qualified_name(covermap, cover_mapset));
strcat(command, "\"");
/* strcpy (command,"cat /tmp/t"); */
G_debug(3, "command: %s", command);
stats_fd = popen(command, "r");
G_debug(3, "r.reclass i=\"%s\" o=\"%s\"",
G_fully_qualified_name(basemap, base_mapset), outmap);
sprintf(command, "r.reclass i=\"%s\" o=\"%s\"",
G_fully_qualified_name(basemap, base_mapset), outmap);
reclass_fd = popen(command, "w");
first = 1;
while (read_stats(stats_fd, &basecat, &covercat, &area)) {
if (first) {
stats.n = 0;
stats.nalloc = 16;
stats.cat = (long *)
G_calloc(stats.nalloc, sizeof(long));
stats.area = (double *)
G_calloc(stats.nalloc, sizeof(double));
first = 0;
catb = basecat;
}
if (basecat != catb) {
catc = median(&stats);
write_reclass(reclass_fd, catb, catc,
G_get_cat(catc, &cover_cats));
catb = basecat;
stats.n = 0;
}
stats.n++;
if (stats.n > stats.nalloc) {
stats.nalloc *= 2;
stats.cat = (long *)
G_realloc(stats.cat, stats.nalloc * sizeof(long));
stats.area = (double *)
G_realloc(stats.area, stats.nalloc * sizeof(double));
}
stats.cat[stats.n - 1] = covercat;
stats.area[stats.n - 1] = area;
}
if (!first) {
catc = median(&stats);
write_reclass(reclass_fd, catb, catc, G_get_cat(catc, &cover_cats));
}
pclose(stats_fd);
pclose(reclass_fd);
exit(EXIT_SUCCESS);
}
int main(int argc, char *argv[])
{
struct GModule *module;
struct
{
struct Option *base, *cover, *output;
} parm;
char *basemap, *base_mapset;
char *covermap, *cover_mapset;
char *outmap;
char command[1024];
struct Categories cover_cats;
FILE *stats, *reclass;
int first;
long basecat, covercat, catb, catc;
double value, max;
G_gisinit(argv[0]);
module = G_define_module();
module->keywords = _("raster, statistics");
module->description =
_("Finds the mode of values in a cover map within "
"areas assigned the same category value in a "
"user-specified base map.");
parm.base = G_define_option();
parm.base->key = "base";
parm.base->description = _("Base map to be reclassified");
parm.base->required = YES;
parm.base->type = TYPE_STRING;
parm.base->gisprompt = "old,cell,raster";
parm.cover = G_define_option();
parm.cover->key = "cover";
parm.cover->description = _("Coverage map");
parm.cover->required = YES;
parm.cover->type = TYPE_STRING;
parm.cover->gisprompt = "old,cell,raster";
parm.output = G_define_option();
parm.output->key = "output";
parm.output->description = _("Output map");
parm.output->required = YES;
parm.output->type = TYPE_STRING;
parm.output->gisprompt = "new,cell,raster";
if (G_parser(argc, argv))
exit(1);
basemap = parm.base->answer;
covermap = parm.cover->answer;
outmap = parm.output->answer;
base_mapset = G_find_cell2(basemap, "");
if (base_mapset == NULL) {
G_fatal_error(_("%s: base raster map not found"), basemap);
}
cover_mapset = G_find_cell2(covermap, "");
if (cover_mapset == NULL) {
G_fatal_error(_("%s: cover raster map not found"), covermap);
}
if (G_legal_filename(outmap) < 0) {
G_fatal_error(_("<%s> is an illegal file name"), outmap);
}
if (strcmp(G_mapset(), base_mapset) == 0 && strcmp(basemap, outmap) == 0) {
G_fatal_error(_("%s: base map and output map must be different"),
outmap);
}
if (G_read_cats(covermap, cover_mapset, &cover_cats) < 0) {
G_fatal_error(_("%s: Unable to read category labels"), covermap);
}
strcpy(command, "r.stats -an \"");
strcat(command, G_fully_qualified_name(basemap, base_mapset));
strcat(command, ",");
strcat(command, G_fully_qualified_name(covermap, cover_mapset));
strcat(command, "\"");
/* printf(command); */
stats = popen(command, "r");
sprintf(command, "r.reclass i=\"%s\" o=\"%s\"",
G_fully_qualified_name(basemap, base_mapset), outmap);
/* printf(command); */
reclass = popen(command, "w");
first = 1;
while (read_stats(stats, &basecat, &covercat, &value)) {
if (first) {
first = 0;
catb = basecat;
catc = covercat;
max = value;
}
if (basecat != catb) {
write_reclass(reclass, catb, catc, G_get_cat(catc, &cover_cats));
catb = basecat;
catc = covercat;
max = value;
}
if (value > max) {
catc = covercat;
max = value;
}
}
if (first) {
catb = catc = 0;
}
write_reclass(reclass, catb, catc, G_get_cat(catc, &cover_cats));
pclose(stats);
pclose(reclass);
exit(0);
}
