Cpeak_interpBiLin.h

/**********************************************************************
 * 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_INTERPBILIN
#define CPEAK_INTERPBILIN

#include "Cpeak_subPix.h"


template<typename T> 
class Cpeak_interpBiLin : public Cpeak_subPix<T>
{
  
  public:
    
    std::string class_name;
    int m_neighb, m_interpFactor; 
    std::vector<T> m_posInt;
    
    Cpeak_interpBiLin() : Cpeak_subPix<T>(), m_neighb(3), m_interpFactor(3)
    {// constructor
      
      class_name = "Cpeak_interpBiLin";
      m_posInt.clear();
      
    }
    
    virtual void sub_pixel(CImg<T> &img)
    {
      
      //stores int peak location
      m_posInt.clear();
      for (int i=0; i<this->m_pos.size(); i++){m_posInt.push_back(this->m_pos[i]);}
      
      CImg<T> imgC = interp(img);

      Cpeak_max<T>::peak_detection(imgC);
      
      set_position();
      
    }

    CImg<T> interp(CImg<T> img)
    {// interp [n-neighbours][n][n+neighbours]
      
      if(!img.is_sameZ(1))
      {
        return img.crop(this->m_pos[0]-m_neighb, this->m_pos[1]-m_neighb, this->m_pos[2]-m_neighb,0, this->m_pos[0]+m_neighb, this->m_pos[1]+m_neighb, this->m_pos[2]+m_neighb,0, "false").resize(-m_interpFactor*100,-m_interpFactor*100,-m_interpFactor*100,-100,3);
      }
      else
      {
        return img.crop(this->m_pos[0]-m_neighb, this->m_pos[1]-m_neighb, 0,0,this->m_pos[0]+m_neighb, this->m_pos[1]+m_neighb, 0,0,"false").resize(-m_interpFactor*100,-m_interpFactor*100,-100,-100,3);
      }

    }

    void set_position()
    {//actualizes sub-pix position

      for (int i=0; i<m_posInt.size(); i++)
      {
        (*this).m_pos[i] = m_posInt[i]-m_neighb + (*this).m_pos[i]/(double)m_interpFactor;
      }
      
    }
   
    
};

#endif