Cpeak_interp.h

Go to the documentation of this file.

/**********************************************************************
 * Copyright (C) 2012, The YaDICs Project Developers. 
 * See the COPYRIGHT file at the top-level directory of this distribution ./COPYRIGHT.
 * See ./COPYING file for copying and redistribution conditions.
 * 
 * This file is part of YaDICs.
 * YaDICs is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * YaDICs is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with YaDICs.  If not, see <http://www.gnu.org/licenses/>.
 * 
 * Information about how to use the software are provided at http://yadic.univ-lille1.fr/
 **********************************************************************/


#ifndef CPEAK_INTERP
#define CPEAK_INTERP

#include "Cpeak_subPix.h"


template<typename T> 
class Cpeak_interp : public Cpeak_subPix<T>
{
  
  public:
    
    std::string class_name;
    
    Cpeak_interp() : Cpeak_subPix<T>()
    {// constructor
      
      class_name = "Cpeak_interp";
      
    }
    
    virtual void sub_pixel(CImg<T> &img) = 0;

    void interp(const T &current, const T &previous, const T &next, T &intPos)
    {// interp [n-1][n][n+1]
    
      //sub_pix position = integer one + sub_position (interpolated at pos +- 1)
      if (2*previous+2*next-4*current!=0)
      {
        intPos += (previous-next)/(2*previous+2*next-4*current);
      }

    }
    
};

#endif