v.colors (C program) added

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@47563 15284696-431f-4ddb-bdfa-cd5b030d7da7

vector/v.colors/Makefile

@@ -0,0 +1,12 @@
+MODULE_TOPDIR = ../..
+
+PGM = v.colors
+
+LIBES = $(VECTORLIB) $(DBMILIB) $(GISLIB) $(RASTERLIB) $(RASTER3DLIB)
+DEPENDENCIES = $(VECTORDEP) $(DBMIDEP) $(GISDEP) $(RASTERDEP) $(RASTER3DDEP)
+EXTRA_INC = $(VECT_INC)
+EXTRA_CFLAGS = $(VECT_CFLAGS)
+
+include $(MODULE_TOPDIR)/include/Make/Module.make
+
+default: cmd

vector/v.colors/local_proto.h

@@ -0,0 +1 @@
+void scan_cats(const struct Map_info *, int, int *, int*);

vector/v.colors/main.c

@@ -0,0 +1,279 @@
+
+/****************************************************************
+ *
+ * MODULE:       v.colors
+ * 
+ * AUTHOR(S):    Martin Landa <landa.martin gmail.com>
+ *               
+ * PURPOSE:      Manage color tables for vector maps
+ *               
+ * COPYRIGHT:    (C) 2011 by the GRASS Development Team
+ *
+ *               This program is free software under the GNU General
+ *               Public License (>=v2). Read the file COPYING that
+ *               comes with GRASS for details.
+ *
+ **************************************************************/
+
+#include <string.h>
+#include <stdlib.h>
+
+#include <grass/gis.h>
+#include <grass/raster.h>
+#include <grass/raster3d.h>
+#include <grass/vector.h>
+#include <grass/glocale.h>
+
+#include "local_proto.h"
+
+int main(int argc, char *argv[])
+{
+    struct GModule *module;
+    struct {
+        struct Flag *r, *w, *l, *g, *a, *n, *e;
+    } flag; 
+
+    struct {
+        struct Option *map, *field, *colr, *rast, *volume, *rules;
+    } opt;
+
+    int min, max, layer;
+    int have_stats;
+    int overwrite, remove, is_from_stdin, stat, have_colors;
+    const char *mapset, *cmapset;
+    const char *name, *style, *rules, *cmap;
+
+    struct Map_info Map;
+    struct Colors colors, colors_tmp;
+    /* struct Cell_stats statf; */
+    
+    G_gisinit(argv[0]);
+    
+    module = G_define_module();
+    G_add_keyword(_("vector"));
+    G_add_keyword(_("color table"));
+    module->description =
+	_("Creates/modifies the color table associated with a vector map.");
+
+    opt.map = G_define_standard_option(G_OPT_V_MAP);
+
+    opt.field = G_define_standard_option(G_OPT_V_FIELD);
+    
+    opt.colr = G_define_standard_option(G_OPT_M_COLR);
+
+    opt.rast = G_define_standard_option(G_OPT_R_INPUT);
+    opt.rast->key = "raster";
+    opt.rast->required = NO;
+    opt.rast->description =
+        _("Raster map from which to copy color table");
+    opt.rast->guisection = _("Define");
+
+    opt.volume = G_define_standard_option(G_OPT_R3_INPUT);
+    opt.volume->key = "volume";
+    opt.volume->required = NO;
+    opt.volume->description =
+        _("3D raster map from which to copy color table");
+    opt.volume->guisection = _("Define");
+
+    opt.rules = G_define_standard_option(G_OPT_F_INPUT);
+    opt.rules->key = "rules";
+    opt.rules->required = NO;
+    opt.rules->label = _("Path to rules file");
+    opt.rules->description = _("\"-\" to read rules from stdin");
+    opt.rules->guisection = _("Define");
+
+    flag.r = G_define_flag();
+    flag.r->key = 'r';
+    flag.r->description = _("Remove existing color table");
+    flag.r->guisection = _("Remove");
+
+    flag.w = G_define_flag();
+    flag.w->key = 'w';
+    flag.w->description =
+        _("Only write new color table if one doesn't already exist");
+
+    flag.l = G_define_flag();
+    flag.l->key = 'l';
+    flag.l->description = _("List available rules then exit");
+    flag.l->suppress_required = YES;
+    flag.l->guisection = _("Print");
+
+    flag.n = G_define_flag();
+    flag.n->key = 'n';
+    flag.n->description = _("Invert colors");
+    flag.n->guisection = _("Define");
+
+    flag.g = G_define_flag();
+    flag.g->key = 'g';
+    flag.g->description = _("Logarithmic scaling");
+    flag.g->guisection = _("Define");
+
+    flag.a = G_define_flag();
+    flag.a->key = 'a';
+    flag.a->description = _("Logarithmic-absolute scaling");
+    flag.a->guisection = _("Define");
+
+    flag.e = G_define_flag();
+    flag.e->key = 'e';
+    flag.e->description = _("Histogram equalization");
+    flag.e->guisection = _("Define");
+    
+    if (G_parser(argc, argv))
+	exit(EXIT_FAILURE);
+
+    if (flag.l->answer) {
+	G_list_color_rules(stdout);
+	return EXIT_SUCCESS;
+    }
+
+    overwrite = !flag.w->answer;
+    remove = flag.r->answer;
+    name = opt.map->answer;
+    style = opt.colr->answer;
+    rules = opt.rules->answer;
+    have_stats = FALSE;
+    
+    if (!name)
+        G_fatal_error(_("No vector map specified"));
+
+    if (opt.rast->answer && opt.volume->answer)
+        G_fatal_error(_("Options <%s> and <%s> are mutually exclusive"),
+		      opt.rast->key, opt.volume->key);
+
+    cmap = NULL;
+    if (opt.rast->answer)
+        cmap = opt.rast->answer;
+    if (opt.volume->answer)
+        cmap = opt.volume->answer;
+    
+    if (!cmap && !style && !rules && !remove)
+        G_fatal_error(_("One of -%c or options <%s>, <%s> or <%s> "
+			"must be specified"), flag.r->key, opt.colr->key,
+			opt.rast->key, opt.rules->key);
+
+    if (!!style + !!cmap + !!rules > 1)
+        G_fatal_error(_("Options <%s>, <%s>, and <%s> are mutually "
+			"exclusive"), opt.colr->key, opt.rules->key,
+		      opt.rast->key);
+
+    if (flag.g->answer && flag.a->answer)
+        G_fatal_error(_("Flags -%c and -%c flags are mutually exclusive"),
+		      flag.g->key, flag.a->key);
+
+    is_from_stdin = rules && strcmp(rules, "-") == 0;
+    if (is_from_stdin)
+        rules = NULL;
+
+    /* open map and get min/max values */
+    Vect_open_old2(&Map, name, "", opt.field->answer);
+    name   = Vect_get_name(&Map);
+    mapset = Vect_get_mapset(&Map);
+    
+    stat = -1;
+    if (remove) {
+	stat = Vect_remove_colors(name, mapset);
+        if (stat < 0)
+            G_fatal_error(_("Unable to remove color table of vector map <%s>"), name);
+        if (stat == 0)
+            G_warning(_("Color table of vector map <%s> not found"), name);
+        return EXIT_SUCCESS;
+    }
+
+    G_suppress_warnings(TRUE);
+    have_colors = Vect_read_colors(name, mapset, &colors);
+
+    if (have_colors > 0 && !overwrite) {
+        G_fatal_error(_("Color table exists. Exiting."));
+    }
+
+    G_suppress_warnings(FALSE);
+
+    layer = Vect_get_field_number(&Map, opt.field->answer);
+    if (layer < 1)
+	G_fatal_error(_("Layer <%s> not found"), opt.field->answer);
+    
+    scan_cats(&Map, layer, &min, &max);
+    
+    if (is_from_stdin) {
+	/*
+        if (!read_color_rules(stdin, &colors, min, max, fp))
+            exit(EXIT_FAILURE);
+	*/
+    } else if (style) {
+	if (strcmp(style, "random") == 0) {
+	    Rast_make_random_colors(&colors, (CELL) min, (CELL) max);
+	} else if (strcmp(style, "grey.eq") == 0) {
+	    G_fatal_error(_("Color table <%s> not supported"), "grey.eq");
+	    /*
+	      if (!have_stats)
+	      have_stats = get_stats(name, mapset, &statf);
+	      Rast_make_histogram_eq_colors(&colors, &statf);
+	    */
+        } else if (strcmp(style, "grey.log") == 0) {
+	    G_fatal_error(_("Color table <%s> not supported"), "grey.log");
+	    /*
+	      if (!have_stats)
+	      have_stats = get_stats(name, mapset, &statf);
+	      Rast_make_histogram_log_colors(&colors, &statf, (CELL) min,
+	      (CELL) max);
+	    */
+	} else if (G_find_color_rule(style))
+	    Rast_make_colors(&colors, style, (CELL) min, (CELL) max);
+        else
+            G_fatal_error(_("Unknown color request '%s'"), style);
+    }
+    else if (rules) {
+	if (!Rast_load_fp_colors(&colors, rules, min, max))
+	    G_fatal_error(_("Unable to load rules file <%s>"), rules);
+    }
+    else {
+	/* use color from another map (cmap) */
+	if (opt.rast->answer) {
+            cmapset = G_find_raster2(cmap, "");
+            if (!cmapset)
+                G_fatal_error(_("Raster map <%s> not found"), cmap);
+
+            if (Rast_read_colors(cmap, cmapset, &colors) < 0)
+                G_fatal_error(_("Unable to read color table for raster map <%s>"), cmap);
+        } else {
+            cmapset = G_find_grid3(cmap, "");
+            if (!cmapset)
+                G_fatal_error(_("3D raster map <%s> not found"), cmap);
+
+            if (Rast3d_readColors(cmap, cmapset, &colors) < 0)
+                G_fatal_error(_("Unable to read color table for 3D raster map <%s>"), cmap);
+        }
+    }
+
+    if (flag.n->answer)
+        Rast_invert_colors(&colors);
+
+    if (flag.e->answer) {
+	G_fatal_error(_("Flag -%c needs to be implemeneted"));
+	/*
+	if (!have_stats)
+	    have_stats = get_stats(name, mapset, &statf);
+	Rast_histogram_eq_colors(&colors_tmp, &colors, &statf);
+	colors = colors_tmp;
+	*/
+    }
+
+    if (flag.g->answer) {
+        Rast_log_colors(&colors_tmp, &colors, 100);
+        colors = colors_tmp;
+    }
+
+    if (flag.a->answer) {
+        Rast_abs_log_colors(&colors_tmp, &colors, 100);
+        colors = colors_tmp;
+    }
+
+    Vect_write_colors(name, mapset, &colors);
+
+    G_message(_("Color table for vector map <%s> set to '%s'"), 
+	      G_fully_qualified_name(name, mapset), 
+              is_from_stdin ? "rules" : style ? style : rules ? rules :
+              cmap);
+    
+    exit(EXIT_SUCCESS);
+}

vector/v.colors/scan_cats.c

@@ -0,0 +1,49 @@
+#include <grass/vector.h>
+#include <grass/glocale.h>
+
+#include "local_proto.h"
+
+static void scan_layer(int, const struct line_cats *, int *, int *);
+
+void scan_cats(const struct Map_info *Map, int field, int *min, int* max)
+{
+    int ltype, cmin, cmax;
+    struct line_cats *Cats;
+
+    *min = *max = -1;
+    
+    Cats = Vect_new_cats_struct();
+
+    while(TRUE) {
+	ltype = Vect_read_next_line(Map, NULL, Cats);
+	if (ltype == -1)
+	    G_fatal_error(_("Unable to read vector map"));
+	if (ltype == -2)
+	    break; /* EOF */
+
+	scan_layer(field, Cats, &cmin, &cmax);
+
+	if (*min == -1 || cmin <= *min)
+	    *min = cmin;
+	if (*max == -1 || cmax >= *max)
+	    *max = cmax;
+    }
+
+    Vect_destroy_cats_struct(Cats);
+}
+
+void scan_layer(int field, const struct line_cats *Cats, int *cmin, int *cmax)
+{
+    int n, cat;
+
+    *cmin = *cmax = -1;
+    for (n = 0; n < Cats->n_cats; n++) {
+        if (Cats->field[n] == field) {
+	    cat = Cats->cat[n];
+	    if (*cmin == -1 || cat <= *cmin)
+		*cmin = cat;
+	    if (*cmax == -1 || cat >= *cmax)
+		*cmax = cat;
+        }
+    }
+}

vector/v.colors/v.colors.html

@@ -0,0 +1,69 @@
+<h2>DESCRIPTION</h2>
+
+<em>v.colors</em> allows the user to create and/or modify the color
+table for a vector map. The vector map (specified on the command line
+by <b>map=</b>) must exist in the user's current mapset search path.
+
+<p>
+The <b>rast</b> option allows user to specify a raster map <i>name</i>
+from which to copy the color map.
+
+<p>
+The <b>volume</b> option allows user to specify a volume (3D raster)
+map <i>name</i> from which to copy the color map.
+
+<p>
+If the user specifies the <b>-w</b> flag, the current color table file
+for the input map will not be overwritten. This means that the color
+table is created only if the <i>map</i> does not already have a color
+table. If this option is not specified, the color table will be
+created if one does not exist, or modified if it does.
+
+<p>
+In general, tables which associate colors with percentages (aspect,
+bcyr, byg, byr, elevation, grey, gyr, rainbow, ramp, ryb, ryg and
+wave) can be applied to any data, while those which use absolute
+values (aspectcolr, curvature, etopo2, evi, ndvi, population, slope,
+srtm, and terrain) only make sense for data with certain ranges.
+
+One can get a rough idea of the applicability of a colour table by
+reading the corresponding rules file
+(<tt>$GISBASE/etc/colors/&lt;name&gt;</tt>). For example
+the <em>ryg</em> rule is defined as:
+
+<div class="code"><pre>
+0% red
+50% yellow
+100% green
+</pre></div>
+
+<p>
+The <b>rules</b> color table type will cause <em>r.colors</em> to read
+color table specifications from standard input and will build
+the color table accordingly.
+
+<p>
+For details see <em><a href="r.colors.html">r.colors</a></em> manual page.
+
+<h2>EXAMPLES</h2>
+
+<h2>SEE ALSO</h2>
+
+<em>
+  <a href="d.vect.html">d.vect</a>,
+  <a href="r.colors.html">r.colors</a>,
+  <a href="r3.colors.html">r3.colors</a>,
+</em>
+
+<p>
+See also wiki
+page <a href="http://grass.osgeo.org/wiki/Color_tables">Color
+tables</a> (from GRASS User Wiki)
+
+<h2>AUTHORS</h2>
+
+Adopted <a href="r.colors.html">r.colors</a> by Martin Landa, Czech
+Technical University in Prague, Czech Republic
+
+<p>
+<i>Last changed: $Date$</i>