Cpeak_interpG.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_INTERPG
#define CPEAK_INTERPG

#include "Cpeak_interpP.h"

template<typename T>
class Cpeak_interpG : public Cpeak_interpP<T>
{

  public:
    
    std::string class_name;
    
    Cpeak_interpG() : Cpeak_interpP<T>()
    {// constructor
      
      class_name = "Cpeak_interpG";
      
    }
    
    virtual void sub_pixel(CImg<T> &img)
    {
      
      std::vector<T> cur;
        cur.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];
      }
      
      
      if (img.get_crop(cur[0]-1,cur[1]-1,cur[2]-1,0,cur[0]+1,cur[1]+1,cur[2]+1,0).get_stats()[0]<=0)// min(img)<=0
      {
        Cpeak_interpP<T>::sub_pixel(img);
        printf("\t\t\t\t\t\t\t\tWarning, image value incompatible with gaussian fit, polynomial one is used");
      }
      else
      {
        
        std::vector<T> prev, next;
          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
          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]<this->dim(img,i)-1)
          {

            prev[i] = prev[i]-1;
            next[i] = next[i]+1;
          
            this->interp(log(img(cur[0],cur[1],cur[2])), log(img(prev[0],prev[1],prev[2])), log(img(next[0],next[1],next[2])), this->m_pos[i]);
          
            prev[i] = cur[i];
            next[i] = cur[i];
          }
          
        }
        
      }

    }// sub_pixel (gaussian)
    
};

#endif