dem.h

Go to the documentation of this file.

#ifndef GCOP_DEM_H
#define GCOP_DEM_H

namespace gcop {

  class Dem {
  public:
    Dem();

    Dem(const char *fname, double cs = 1.0, double ds = 1.0, const double *o = 0);

    Dem(double w, double h, double cs = 1.0, double ds = 1.0, const double* o = 0);
    
    Dem(const Dem &dem);
    
    virtual ~Dem();

    virtual bool IsValid(double x, double y) const;

    virtual double Get(double x, double y) const;

    double GetNormal(double n[3], double x, double y) const;

    const double* GetNormal(double x, double y) const;

    virtual void Get(double *p, int i, int j) const;    

    void Set(int i, int j, double z);

    void Set(double x, double y, double z, double s = 0);    

    void Clear();

    void Scale(double s);

    bool Inside(double x, double y, double z) const;

    void Load(const char *fname);

    const double* GetNormal(int i, int j) const;

    void ComputeNormals();

    static double bilinterp(const double* z, int w, int h, double xi, double yi, double eps = 1e-10);

    void Dilate(double r, bool cube = true);
    
    void AddBoundary(double h);

    void Convolve(double sigma, bool cn = true, double thresh = 0);
    
    void Point2Index(int &i, int &j, double x, double y) const;

    void Index2Point(double &x, double &y, int i, int j) const;


    double w;          
    double h;          
    double cs;         
    double ds;         
    double o[3];       
    int ni;            
    int nj;            
    double *data;      
    double *odata;     
        
    double *normals;   

    double eps;        
  };


  inline bool Dem::IsValid(double x, double y) const
    {
      return 
        x - o[0] >= -eps && x - o[0] <= w + eps && 
        y - o[1] >= -eps && y - o[1] <= h + eps;
    }
  
};



#endif