Ccorrelation_intercor_fft.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 CCORRELATION_INTERCOR_FFT
#define CCORRELATION_INTERCOR_FFT



#include "Ccorrelation_intercor_fftw_phase.h"



template<typename T,typename Timg> 
class Ccorrelation_intercor_fft : public Ccorrelation_intercor_fftw_phase<T,Timg>
{

  public:
    
    
    Ccorrelation_intercor_fft() : Ccorrelation_intercor_fftw_phase<T,Timg>()
    {// constructor
      
      this->class_name = "Correlation : Inter-correlation by FFT";
      
    }

    void correl(CImgList<T> &correlogram, CImgList<T> &F, CImgList<T> &G)
    {// intercorrelation + zeroPadd
      
      zeroPadd(correlogram, F, G);
      Ccorrelation_intercor_fftw<T,Timg>::correl(correlogram, F, G);
      
    }
    
    void zeroPadd(CImgList<T> &correlogram, CImgList<T> &F, CImgList<T> &G)
    {
      
      // nearest power of 2
      int dimX(-100), dimY(-100), dimZ(-100);
      
      dimX = cimg::nearest_pow2(correlogram[0].width());
      if(!correlogram[0].is_sameY(1)){dimY = cimg::nearest_pow2(F[0].height()); };
      if(!correlogram[0].is_sameZ(1)){dimZ = cimg::nearest_pow2(F[0].depth()); };
      
      cimglist_for(F,i)
      {
        correlogram[i].resize(dimX,dimY,dimZ,-100, 0, 0, 0.5, 0.5, 0.5, 0);
        F[i].assign(dimX,dimY,dimZ);
        G[i].assign(dimX,dimY,dimZ);
      }
      
    }// zeroPadd    

    virtual void interCor(CImgList<T> &correlogram, CImgList<T> &F, CImgList<T> &G)
    {//inter-correlation
      
      Ccorrelation_intercor<T,Timg>::interCor(correlogram, F, G);
      
    }
    
};

#endif