Geografic-Information-System/raster/r.profile/read_rast.c

/*
 * Copyright (C) 2000 by the GRASS Development Team
 * Author: Bob Covill <bcovill@tekmap.ns.ca>
 * 
 * This Program is free software under the GPL (>=v2)
 * Read the file COPYING coming with GRASS for details
 *
 */

#include <grass/gis.h>
#include <grass/raster.h>
#include <grass/glocale.h>
#include "local_proto.h"

int read_rast(double east, double north, double dist, int fd, int coords,
	      RASTER_MAP_TYPE data_type, FILE * fp, char *null_string)
{
    static DCELL *dcell;
    static int cur_row = -1;
    static CELL nullcell;
    static int nrows, ncols;
    static struct Cell_head window;
    int row, col;
    int outofbounds = FALSE;

    if (!dcell) {
	Rast_set_c_null_value(&nullcell, 1);
	dcell = Rast_allocate_d_buf();
	G_get_window(&window);
	nrows = window.rows;
	ncols = window.cols;
    }

    row = (window.north - north) / window.ns_res;
    col = (east - window.west) / window.ew_res;
    G_debug(4, "row=%d:%d  col=%d:%d", row, nrows, col, ncols);

    if ((row < 0) || (row >= nrows) || (col < 0) || (col >= ncols))
	outofbounds = TRUE;

    if (!outofbounds) {
	if (row != cur_row)
	    Rast_get_d_row(fd, dcell, row);
	cur_row = row;
    }

    if (coords)
	fprintf(fp, "%f %f", east, north);

    fprintf(fp, " %f", dist);

    if (outofbounds || Rast_is_d_null_value(&dcell[col]))
	fprintf(fp, " %s", null_string);
    else {
	if (data_type == CELL_TYPE)
	    fprintf(fp, " %d", (int) dcell[col]);
	else
	    fprintf(fp, " %f", dcell[col]);
    }

    if (clr) {
	int red, green, blue;

	if (outofbounds)
	    Rast_get_c_color(&nullcell, &red, &green, &blue, &colors);
	else
	    Rast_get_d_color(&dcell[col], &red, &green, &blue,
				 &colors);

	fprintf(fp, " %03d:%03d:%03d", red, green, blue);
    }

    fprintf(fp, "\n");

    return 0;
}