Geografic-Information-System/vector/v.lidar.growing/growing.c

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include "growing.h"

void P_Aux_to_Coor(struct Map_info *In, struct Map_info *Out,
		   dbDriver * driver, FILE * fsite)
{

    int cont = 0;
    int more, ltype, line_num, ID_type, Interp_type;
    double quotaZ;

    struct line_pnts *point;
    struct line_cats *cat;
    dbTable *table;
    dbColumn *ID_column, *Interp_column;
    dbValue *ID_value, *Interp_value;
    dbCursor cursor;
    dbString sql;

    point = Vect_new_line_struct();
    cat = Vect_new_cats_struct();
    db_init_string(&sql);
    db_zero_string(&sql);

    db_append_string(&sql,
		     "select ID, sum(Interp) from Auxiliar_table group by ID");
    db_open_select_cursor(driver, &sql, &cursor, DB_SEQUENTIAL);

    while (db_fetch(&cursor, DB_NEXT, &more) == DB_OK && more) {
	cont++;

	table = db_get_cursor_table(&cursor);

	ID_column = db_get_table_column(table, 0);
	Interp_column = db_get_table_column(table, 1);

	ID_type = db_sqltype_to_Ctype(db_get_column_sqltype(ID_column));
	Interp_type =
	    db_sqltype_to_Ctype(db_get_column_sqltype(Interp_column));

	if (ID_type == DB_C_TYPE_INT)
	    ID_value = db_get_column_value(ID_column);
	else
	    continue;

	if (Interp_type == DB_C_TYPE_DOUBLE)
	    Interp_value = db_get_column_value(Interp_column);
	else
	    continue;

	line_num = db_get_value_int(ID_value);
	quotaZ = db_get_value_double(Interp_value);

	ltype = Vect_read_line(In, point, cat, line_num);

	if (!(ltype & GV_POINT))
	    continue;

	point->z[0] = quotaZ;
	Vect_write_line(Out, ltype, point, cat);
    }
    return;
}