Geografic-Information-System/raster3d/r3.mkdspf/r3_data.c

#include "vizual.h"
#include <grass/raster3d.h>


#define NO_RESAMPLE


int r3read_level(void *g3map, RASTER3D_Region * g3reg, file_info * Headfax,
		 float *data, int n_lev)
{
#ifdef NO_RESAMPLE
    Rast3d_get_block(g3map, 0, 0, n_lev,
		 Headfax->xdim, Headfax->ydim, 1, (char *)data, FCELL_TYPE);
#else
    /* Rast3d_getBlockRegion might be handy */
    /*
       Rast3d_getAllignedVolume (map, originNorth, originWest, originBottom, 
       lengthNorth, lengthWest, lengthBottom, 
       nx, ny, nz, volumeBuf, type);
     */
    Rast3d_getAllignedVolume(g3map, g3reg->north, g3reg->west,
			  g3reg->top - n_lev * g3reg->tb_res,
			  g3reg->north - g3reg->south,
			  g3reg->east - g3reg->west,
			  g3reg->tb_res,
			  Headfax->xdim, Headfax->ydim, 1,
			  (char *)data, FCELL_TYPE);
#endif

    return 0;
}