---

base_2D.h

Go to the documentation of this file.

\endverbatim

#ifndef __BASE_2D__
#define __BASE_2D__


namespace Function_2D{



  class Base_2D{

  public:
    typedef double real_type;

    inline Base_2D(const real_type in_offset_x = 0,
                   const real_type in_offset_y = 0,
                   const real_type in_scale_x = 1,
                   const real_type in_scale_y = 1,
                   const real_type in_rotation_angle = 0,
                   const real_type out_offset = 0,
                   const real_type out_scale = 1);
  
    virtual ~Base_2D();

    inline real_type operator()(const real_type x,
                                const real_type y) const;
    
  
  protected:
    virtual inline real_type base_value_proxy(const real_type x,
                                              const real_type y) const=0;

    inline void update_input_rotation(const real_type a);
    
  


  
    real_type input_offset_x,input_offset_y;



  
    real_type input_scale_x,input_scale_y;

    real_type input_rotation_angle;


  
    real_type cos_a,sin_a;

    real_type output_scale;

    real_type output_offset;
  

  };

  

  class Direct_interface_2D:public Base_2D{

  public:

    inline Direct_interface_2D(const real_type in_offset_x = 0,
                               const real_type in_offset_y = 0,
                               const real_type in_scale_x = 1,
                               const real_type in_scale_y = 1,
                               const real_type in_rotation_angle = 0,
                               const real_type out_offset = 0,
                               const real_type out_scale = 1);
    
    inline void set_input_offset(const real_type x,const real_type y);
    inline void set_input_scale(const real_type x,const real_type y);
    inline void set_input_rotation(const real_type a);

    inline void set_output_offset(const real_type o);
    inline void set_output_scale(const real_type s);

  };





  class Filter_interface_2D:public Base_2D{

  public:

    inline Filter_interface_2D(const real_type center_x = 0,
                               const real_type center_y = 0,
                               const real_type sigma_x = 1,
                               const real_type sigma_y = 1,
                               const real_type sig_scale_x = 1,
                               const real_type sig_scale_y = 1,
                               const real_type orientation_angle = 0,
                               const real_type out_sigma_scale = 1);


    inline void set_center(const real_type x,const real_type y);
    inline void set_orientation(const real_type o);
    inline void set_sigma(const real_type x,const real_type y);

  protected:

    
    real_type sigma_scale_x,sigma_scale_y;

    real_type output_sigma_scale;
};




/*
  ###########################################################
  ###########################################################
  ###########################################################
  ##############                               ##############
  ##############  I M P L E M E N T A T I O N  ##############
  ##############                               ##############
  ###########################################################
  ###########################################################
  ###########################################################
*/


  /*

  ###########
  # Base_2D #
  ###########

  */

  Base_2D::Base_2D(const real_type in_offset_x,
                   const real_type in_offset_y,
                   const real_type in_scale_x,
                   const real_type in_scale_y,
                   const real_type in_rotation_angle,
                   const real_type out_offset,
                   const real_type out_scale):
  
    input_offset_x(in_offset_x),
    input_offset_y(in_offset_y),
    input_scale_x(in_scale_x),
    input_scale_y(in_scale_y),
    input_rotation_angle(in_rotation_angle),
    cos_a(cos(in_rotation_angle)),
    sin_a(sin(in_rotation_angle)),
    output_scale(out_scale),
    output_offset(out_offset){
  }

  Base_2D::~Base_2D(){
  }

  Base_2D::real_type Base_2D::operator()(const real_type x,
                                         const real_type y) const{


    const real_type off_x = x + input_offset_x;
    const real_type off_y = y + input_offset_y;
    
    const real_type input_x = (off_x*cos_a + off_y*sin_a)*input_scale_x;
    const real_type input_y =(-off_x*sin_a + off_y*cos_a)*input_scale_y;
    
    return base_value_proxy(input_x,input_y)*output_scale + output_offset;
  }

  void Base_2D::update_input_rotation(const real_type a){
    input_rotation_angle = a;
    cos_a = cos(input_rotation_angle);
    sin_a = sin(input_rotation_angle);
  }


  
/*

  #######################
  # Direct_interface_2D #
  #######################

*/

  Direct_interface_2D::Direct_interface_2D(const real_type in_offset_x,
                                           const real_type in_offset_y,
                                           const real_type in_scale_x,
                                           const real_type in_scale_y,
                                           const real_type in_rotation_angle,
                                           const real_type out_offset,
                                           const real_type out_scale):
    Base_2D(in_offset_x,in_offset_y,
            in_scale_x,in_scale_y,
            in_rotation_angle,
            out_offset,out_scale){}


  
  void Direct_interface_2D::set_input_offset(const real_type x,
                                             const real_type y){

    input_offset_x = x;
    input_offset_y = y;
  }

  void Direct_interface_2D::set_input_scale(const real_type x,
                                            const real_type y){
  
    input_scale_x = x;
    input_scale_y = y;
  }

  void Direct_interface_2D::set_input_rotation(const real_type a){
  
    update_input_rotation(a);
  }

  void Direct_interface_2D::set_output_scale(const real_type s){

    output_scale = s;
  }
  
  void Direct_interface_2D::set_output_offset(const real_type o){

    output_offset = o;
  }


  
/*

  #######################
  # Filter_interface_2D #
  #######################

*/


  Filter_interface_2D::Filter_interface_2D(const real_type center_x,
                                           const real_type center_y,
                                           const real_type sigma_x,
                                           const real_type sigma_y,
                                           const real_type sig_scale_x,
                                           const real_type sig_scale_y,
                                           const real_type orientation_angle,
                                           const real_type out_sigma_scale):
    Base_2D(-center_x,-center_y,
            sig_scale_x/sigma_x,
            sig_scale_y/sigma_y,
            -orientation_angle,
            0,
            out_sigma_scale/(sigma_x*sigma_y)),
  sigma_scale_x(sig_scale_x),
  sigma_scale_y(sig_scale_y),
  output_sigma_scale(out_sigma_scale){
  }


  void Filter_interface_2D::set_center(const real_type x,const real_type y){
    
    input_offset_x = -x;
    input_offset_y = -y;
  }
  
  void Filter_interface_2D::set_orientation(const real_type o){
    
    update_input_rotation(-o);
  }

  void Filter_interface_2D::set_sigma(const real_type x,const real_type y){

    input_scale_x = sigma_scale_x / x;
    input_scale_y = sigma_scale_y / y;
    output_scale  = output_sigma_scale / (x*y);
  }


}
  
#endif

---