Cpeak_interpP.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_INTERPP
#define CPEAK_INTERPP

#include "Cpeak_interp.h"

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

      std::vector<T> cur, prev, next;
        cur.assign(3,0);
        prev.assign(3,0);
        next.assign(3,0);
      
      for (int i=0; i<this->m_pos.size(); i++)
      {// to keep only one formalisme for 2d and 3d
        cur[i] = this->m_pos[i];
        prev[i] = this->m_pos[i]; 
        next[i] = this->m_pos[i];
      }
      
      
      for (int i=0; i<this->m_pos.size(); i++)
      {
        
        if (cur[i]>0 and cur[i]<dim(img,i)-1)
        {//borders conditions

          prev[i] = prev[i]-1;
          next[i] = next[i]+1;
          
          this->interp(img(cur[0],cur[1],cur[2]), img(prev[0],prev[1],prev[2]), img(next[0],next[1],next[2]), this->m_pos[i]);
          
          prev[i] = cur[i];
          next[i] = cur[i];
          
        }
        //else one keeps integer value
        
      }

    }
    
    int dim(CImg<T> &img,const int &i)
    {
      
      switch(i)
      {
        case 0: return img.width(); break;
        case 1: return img.height(); break;
        case 2: return img.depth(); break;
        case 3: return img.spectrum(); break;
        default:
        {
          std::cerr<<class_name<<"::"<<__func__<<": error: id="<<i<<" is unknown as CImg dimensions\n"<<std::flush;
          return 0;
        }
      }
      
    }
    
};

#endif