Geografic-Information-System/raster/r.recode/recode.c

#include <stdio.h>
#include "global.h"

#define F2I(map_type) \
	(map_type == CELL_TYPE ? 0 : (map_type == FCELL_TYPE ? 1 : 2))

static void process_row_ii(int), process_row_if(int), process_row_id(int);
static void process_row_fi(int), process_row_ff(int), process_row_fd(int);
static void process_row_di(int), process_row_df(int), process_row_dd(int);

static void (*process_row_FtypeOtype[3][3]) () = { {
process_row_ii, process_row_if, process_row_id}, {
process_row_fi, process_row_ff, process_row_fd}, {
process_row_di, process_row_df, process_row_dd}};

#define PROCESS_ROW \
     (process_row_FtypeOtype [F2I (in_type)] [F2I (out_type)])

static void *in_rast, *out_rast;
static int nrows, ncols;

int do_recode(void)
{
    struct Cell_head window, cellhd;
    int row, i;
    struct History hist;

    /* set the window from the header for the input file */
    if (align_wind) {
	G_get_window(&window);
	if (G_get_cellhd(name, "", &cellhd) >= 0) {
	    G_align_window(&window, &cellhd);
	    G_set_window(&window);
	}
    }

    G_get_set_window(&window);

    nrows = G_window_rows();
    ncols = G_window_cols();

    /* open the input file for reading */
    in_fd = G_open_cell_old(name, "");
    if (in_fd < 0)
	G_fatal_error("Can't open input map");

    out_fd = G_open_raster_new(result, out_type);

    out_rast = G_allocate_raster_buf(out_type);
    in_rast = G_allocate_raster_buf(in_type);

    for (row = 0; row < nrows; row++) {
	G_percent(row, nrows, 2);
	PROCESS_ROW(row);
    }
    G_percent(row, nrows, 2);
    G_close_cell(in_fd);
    G_close_cell(out_fd);

    /* writing history file */
    G_short_history(result, "raster", &hist);
    sprintf(hist.edhist[0], "recode of raster map %s", name);
    /* if there are more rules than history lines allocated, write only 
       MAXEDLINES-1 rules , and "...." as a last rule */
    for (i = 0; (i < nrules) && (i < MAXEDLINES - 2); i++)
	sprintf(hist.edhist[i + 1], "%s", rules[i]);
    if (nrules > MAXEDLINES - 1) {
	sprintf(hist.edhist[MAXEDLINES - 1], "...");
	hist.edlinecnt = MAXEDLINES;
    }
    else
	hist.edlinecnt = nrules + 1;
    sprintf(hist.datsrc_1, "raster map %s", name);
    G_write_history(result, &hist);

    return 0;
}

static void process_row_ii(int row)
{
    if (no_mask)
	G_get_c_raster_row_nomask(in_fd, (CELL *) in_rast, row);
    else
	G_get_c_raster_row(in_fd, (CELL *) in_rast, row);
    G_fpreclass_perform_ii(&rcl_struct, (CELL *) in_rast, (CELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (CELL *) out_rast, CELL_TYPE);
}

static void process_row_if(int row)
{
    if (no_mask)
	G_get_c_raster_row_nomask(in_fd, (CELL *) in_rast, row);
    else
	G_get_c_raster_row(in_fd, (CELL *) in_rast, row);
    G_fpreclass_perform_if(&rcl_struct, (CELL *) in_rast, (FCELL *) out_rast,
			   ncols);
    G_put_f_raster_row(out_fd, (FCELL *) out_rast);
}

static void process_row_id(int row)
{
    if (no_mask)
	G_get_c_raster_row_nomask(in_fd, (CELL *) in_rast, row);
    else
	G_get_c_raster_row(in_fd, (CELL *) in_rast, row);
    G_fpreclass_perform_id(&rcl_struct, (CELL *) in_rast, (DCELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (DCELL *) out_rast, DCELL_TYPE);
}

static void process_row_fi(int row)
{
    if (no_mask)
	G_get_f_raster_row_nomask(in_fd, (FCELL *) in_rast, row);
    else
	G_get_f_raster_row(in_fd, (FCELL *) in_rast, row);
    G_fpreclass_perform_fi(&rcl_struct, (FCELL *) in_rast, (CELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (CELL *) out_rast, CELL_TYPE);
}

static void process_row_ff(int row)
{
    if (no_mask)
	G_get_f_raster_row_nomask(in_fd, (FCELL *) in_rast, row);
    else
	G_get_f_raster_row(in_fd, (FCELL *) in_rast, row);
    G_fpreclass_perform_ff(&rcl_struct, (FCELL *) in_rast, (FCELL *) out_rast,
			   ncols);
    G_put_f_raster_row(out_fd, (FCELL *) out_rast);
}

static void process_row_fd(int row)
{
    if (no_mask)
	G_get_f_raster_row_nomask(in_fd, (FCELL *) in_rast, row);
    else
	G_get_f_raster_row(in_fd, (FCELL *) in_rast, row);
    G_fpreclass_perform_fd(&rcl_struct, (FCELL *) in_rast, (DCELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (DCELL *) out_rast, DCELL_TYPE);
}

static void process_row_di(int row)
{
    if (no_mask)
	G_get_d_raster_row_nomask(in_fd, (DCELL *) in_rast, row);
    else
	G_get_d_raster_row(in_fd, (DCELL *) in_rast, row);
    G_fpreclass_perform_di(&rcl_struct, (DCELL *) in_rast, (CELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (CELL *) out_rast, CELL_TYPE);
}

static void process_row_df(int row)
{
    if (no_mask)
	G_get_d_raster_row_nomask(in_fd, (DCELL *) in_rast, row);
    else
	G_get_d_raster_row(in_fd, (DCELL *) in_rast, row);
    G_fpreclass_perform_df(&rcl_struct, (DCELL *) in_rast, (FCELL *) out_rast,
			   ncols);
    G_put_f_raster_row(out_fd, (FCELL *) out_rast);
}

static void process_row_dd(int row)
{
    if (no_mask)
	G_get_d_raster_row_nomask(in_fd, (DCELL *) in_rast, row);
    else
	G_get_d_raster_row(in_fd, (DCELL *) in_rast, row);
    G_fpreclass_perform_dd(&rcl_struct, (DCELL *) in_rast, (DCELL *) out_rast,
			   ncols);
    G_put_raster_row(out_fd, (DCELL *) out_rast, DCELL_TYPE);
}