Geografic-Information-System/lib/raster3d/open2.c

#include <grass/gis.h>
#include <grass/G3d.h>

/*----------------------------------------------------------------------------*/

/*!
 * \brief 
 *
 * Opens new g3d-file with <em>name</em> in the current mapset. Tiles
 * are stored in memory with <em>typeIntern</em> which must be one of FCELL_TYPE,
 * DCELL_TYPE, or G3D_TILE_SAME_AS_FILE. <em>cache</em> specifies the
 * cache-mode used and must be either G3D_NO_CACHE, G3D_USE_CACHE_DEFAULT,
 * G3D_USE_CACHE_X, G3D_USE_CACHE_Y, G3D_USE_CACHE_Z,
 * G3D_USE_CACHE_XY, G3D_USE_CACHE_XZ, G3D_USE_CACHE_YZ,
 * G3D_USE_CACHE_XYZ, the result of <tt>G3d_cacheSizeEncode ()</tt>
 * (cf.{g3d:G3d.cacheSizeEncode}), or any positive integer which
 * specifies the number of tiles buffered in the cache.  <em>region</em> specifies
 * the 3d region.  
 * The map is created using the <em>type</em> which must be of FCELL_TYPE or DCELL_TYPE.
 * The methods for compression can be specified LZW or RLE. The digits of the floating point mantissa
 * can be specified. In case of FCELL_TYPE 0-23 and 0-52 in case of DCELL_TYPE. 
 * The number cells in X, Y and Z direction defines the size of each tile.
 * Returns a pointer to the cell structure ... if successful,
 * NULL ... otherwise.
 *
 *  \param name The name of the map
 *  \param typeIntern The internal storage type for in memory chached tiles
 *  \param cache The type of the caching
 *  \param region The region of the map
 *  \param type The type of the map (FCELL_TYPE, or DCELL_TYPE)
 *  \param doLzw Use the LZW compression algorithm  
 *  \param doRle Use the Run-Length-Encoding algroithm for compression
 *  \param precision The precision used for the mantissa (0 - 52) or G3D_MAX_PRECISION
 *  \param tileX The number of cells in X direction of a tile
 *  \param tileY The number of cells in Y direction of a tile
 *  \param tileZ The number of cells in Z direction of a tile
 *  \return void * 
 */


void *G3d_openNewParam(const char *name, int typeIntern, int cache,
		       G3D_Region * region, int type, int doLzw, int doRle,
		       int precision, int tileX, int tileY, int tileZ)
{
    void *map;
    int oldCompress, oldLzw, oldRle, oldPrecision, oldTileX, oldTileY,
	oldTileZ;
    int oldType;

    G3d_initDefaults();

    G3d_getCompressionMode(&oldCompress, &oldLzw, &oldRle, &oldPrecision);
    G3d_setCompressionMode(oldCompress, doLzw, doRle, precision);

    G3d_getTileDimension(&oldTileX, &oldTileY, &oldTileZ);
    G3d_setTileDimension(tileX, tileY, tileZ);

    oldType = G3d_getFileType();
    G3d_setFileType(type);

    map = G3d_openCellNew(name, typeIntern, cache, region);

    G3d_setCompressionMode(oldCompress, oldLzw, oldRle, oldPrecision);
    G3d_setTileDimension(oldTileX, oldTileY, oldTileZ);
    G3d_setFileType(oldType);

    return map;
}

/*----------------------------------------------------------------------------*/

/*!
 * \brief 
 *
 * Opens new g3d-file with <em>name</em> in the current mapset. This method tries to compute 
 * optimal tile size based on the number of rows, cols and depths and the maximum allowed tile size in KB. 
 * Tiles are stored in memory using G3D_TILE_SAME_AS_FILE method. <em>cache</em> specifies the
 * cache-mode used and must be either G3D_NO_CACHE, G3D_USE_CACHE_DEFAULT,
 * G3D_USE_CACHE_X, G3D_USE_CACHE_Y, G3D_USE_CACHE_Z,
 * G3D_USE_CACHE_XY, G3D_USE_CACHE_XZ, G3D_USE_CACHE_YZ,
 * G3D_USE_CACHE_XYZ, the result of <tt>G3d_cacheSizeEncode ()</tt>
 * (cf.{g3d:G3d.cacheSizeEncode}), or any positive integer which
 * specifies the number of tiles buffered in the cache.  <em>region</em> specifies
 * the 3d region.  
 * The map is created using the <em>type</em> which must be of FCELL_TYPE or DCELL_TYPE.
 * Returns a pointer to the cell structure ... if successful,
 * NULL ... otherwise.
 *
 *  \param name The name of the map
 *  \param cache The type of the caching
 *  \param region The region of the map
 *  \param type The type of the map (FCELL_TYPE, or DCELL_TYPE)
 *  \param maxSize The maximum size of a tile in kilo bytes
 *  \return void * 
 */


void *G3d_openNewOptTileSize(const char *name, int cache, G3D_Region * region, int type, int maxSize)
{
    void *map;
    int oldTileX, oldTileY, oldTileZ, oldType;
    int tileX, tileY, tileZ;

    G3d_initDefaults();


    G3d_getTileDimension(&oldTileX, &oldTileY, &oldTileZ);
    
    G3d_computeOptimalTileDimension(region, type, &tileX, &tileY, &tileZ, maxSize);

    G3d_setTileDimension(tileX, tileY, tileZ);

    oldType = G3d_getFileType();
    G3d_setFileType(type);

    map = G3d_openCellNew(name, G3D_TILE_SAME_AS_FILE, cache, region);

    G3d_setTileDimension(oldTileX, oldTileY, oldTileZ);
    G3d_setFileType(oldType);

    return map;
}