GdalRaster.hh

Go to the documentation of this file.

#ifndef _GDAL_RASTERHH_
#define _GDAL_RASTERHH_

#include <openmodeller/om_defs.hh>
#include <openmodeller/env_io/Header.hh>
#include <openmodeller/env_io/Raster.hh>

#include <string>
#include <vector>

class GDALDataset;
class MapFormat;

/****************************************************************/
/************************** Raster Gdal *************************/

class dllexp GdalRaster : public Raster
{
public:

  // Empty constructor
  GdalRaster(): f_data(0), f_size(0), f_format(-1), f_currentRow(-1), f_changed(0), f_warped_ds(0) {};

  void createRaster( const std::string& file, int categ = 0 );

  #ifdef MPI_FOUND

  void createRaster( const std::string& output_file_source, const std::string& file, const MapFormat& format );
  #else

  void createRaster( const std::string& file, const MapFormat& format );
  #endif

  ~GdalRaster();

  int iget( int x, int y, Scalar *val );
  int iput( int x, int y, Scalar val );

  int get( Coord x, Coord y, Scalar *val );

  int put( Coord x, Coord y, Scalar val );

  int put( Coord x, Coord y );

  int getMinMax( Scalar *min, Scalar *max );

  int calcMinMax( int band=0 );

  void finish();

  int deleteRaster();

  bool hasCustomGeotransform() { return (f_warped_ds == 0) ? false: true; }

  int getExtentInStandardCs( Coord *xmin, Coord *ymin, Coord *xmax, Coord *ymax );

private:

  void open( char mode );

  void create( int format );

  void initBuffer();

  static void initGdal();

  void read ( Scalar *buf, int first_row, int num_rows );

  void write( Scalar *buf, int first_row, int num_rows );

  void loadRow( int row, bool writeOperation=false );

  void saveRow(); // Save the current row.

  GDALDataset *f_ds;

  Scalar *f_data; // One line data for all bands.
  int     f_size; // Size of one line.

  int     f_format; // File format used to create the raster (MapFormat::getFormat())

  int f_currentRow;
  int f_changed;

  // Disable copying.
  GdalRaster( const GdalRaster& );
  GdalRaster& operator=( const GdalRaster& );

  // Warped data set converting to standard openModeller coordinate system.
  // Only used when the raster cs is different from openModeller's default.
  GDALDataset *f_warped_ds;
};

#endif