---

gabor.h

Go to the documentation of this file.

#ifndef __GABOR__
#define __GABOR__

#include "base_2D.h"
#include "gaussian_2D.h"

namespace Function_2D{

  class Gabor:public Filter_interface_2D{

  public:
    Gabor::Gabor(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 orientation_angle = 0,
                 const real_type w_phase = 0,
                 const real_type w_length = 1);

    inline void set_phase(const real_type p);
    inline void set_wave_length(const real_type w);

  protected:

    real_type phase;

    real_type wave_length;

    
    virtual inline real_type base_value_proxy(const real_type x,
                                              const real_type y) const;


  };



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

  Gabor::Gabor(const real_type center_x,
               const real_type center_y,
               const real_type sigma_x,
               const real_type sigma_y,
               const real_type orientation_angle,
               const real_type w_phase,
               const real_type w_length):
    Filter_interface_2D(center_x,center_y,
                        sigma_x,sigma_y,
                        M_SQRT1_2,M_SQRT1_2,
                        orientation_angle,
                        0.5/M_PI),
    phase(w_phase),
    wave_length(w_length){}
  
  void Gabor::set_phase(const real_type p){
    phase = p;
  }

  void Gabor::set_wave_length(const real_type w){
    wave_length = w;
  }

  Gabor::real_type
  Gabor::base_value_proxy(const real_type x,
                          const real_type y) const{

    return cos(2*M_PI*x/wave_length + phase) * Gaussian::base_value(x,y);
  }
  
}

#endif

---