Geografic-Information-System/display/d.vect/zcoor.c

/* Print z coordinate value for each node */

#include <grass/gis.h>
#include <grass/vector.h>
#include <grass/display.h>
#include <grass/glocale.h>

#include "local_proto.h"
#include "plot.h"

int display_zcoor(struct Map_info *Map, int type, LATTR *lattr)
{
    int num, el, ltype;
    double xl, yl, zl;
    struct line_pnts *Points;
    char text[50];

    if (!Vect_is_3d(Map)) {
	G_warning(_("Vector map is not 3D. Unable to display z-coordinates."));
	return 1;
    }
    
    G_debug(1, "display zcoor:");
    Points = Vect_new_line_struct();
    
    D_RGB_color(lattr->color.R, lattr->color.G, lattr->color.B);
    D_text_size(lattr->size, lattr->size);
    if (lattr->font)
	D_font(lattr->font);
    if (lattr->enc)
	D_encoding(lattr->enc);

    Vect_rewind(Map);
    /* Points - no nodes registered */
    while (TRUE) {
        ltype = Vect_read_next_line(Map, Points, NULL);
        if (ltype == -1) {
            G_warning(_("Unable to read vector map"));
            return 1;
        }
        else if (ltype == -2) {
            break;
        }
        if ((ltype != GV_POINT) && (ltype & type))
            continue;
        
        sprintf(text, "%.2f", Points->z[0]);
        show_label(&Points->x[0], &Points->y[0], lattr, text);
    }
    Vect_destroy_line_struct(Points);
    
    if (Vect_level(Map) < 2) {
        /* no topology */
        return 0;
    }
    
    num = Vect_get_num_nodes(Map);
    G_debug(1, "n_nodes = %d", num);

    /* Nodes */
    for (el = 1; el <= num; el++) {
	if (!Vect_node_alive(Map, el))
	    continue;
	Vect_get_node_coor(Map, el, &xl, &yl, &zl);
	G_debug(3, "node = %d", el);

	sprintf(text, "%.2f", zl);
	show_label(&xl, &yl, lattr, text);
    }

    return 0;
}