Ccorrelation_intercor_fftw_phase.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_FFTW_PHASE
#define CCORRELATION_INTERCOR_FFTW_PHASE



#include "Ccorrelation_intercor_fftw.h"

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

  int m_factor;
  
  public:
    
    Ccorrelation_intercor_fftw_phase() : Ccorrelation_intercor_fftw<T,Timg>()
    {// constructor
      
      this->class_name = "Correlation : Inter-correlation by FFTW-Phase";
      
    }

    virtual void init(CParameterNetCDF &fp)
    {// init from parameter file
    
      Ccorrelation_intercor<T,Timg>::init(fp);
      
      std::string attribute_name = "phase_threshold";
      int error = fp.loadAttribute(attribute_name,m_factor);
      printf("\t\t phase signal threshold -> %i\n", m_factor);

    }
    
    virtual void interCor(CImgList<T> &correlogram, CImgList<T> &F, CImgList<T> &G)
    {//inter-correlation
      
      cimg_forXYZ(F[0],x,y,z)
      {// CImgList correlogram is used to store ||F|| and ||G||, in order not to create other variables
        correlogram(0,x,y,z) = std::sqrt( std::pow(F(0,x,y,z),2) + std::pow(F(1,x,y,z),2) );
        correlogram(1,x,y,z) = std::sqrt( std::pow(G(0,x,y,z),2) + std::pow(G(1,x,y,z),2) );
      }
      
      normalized(correlogram, F, G);// compute F and G phases in place
      
      Ccorrelation_intercor<T,Timg>::interCor(correlogram, F, G);// classic inter-correlation
      
    }
    
    void normalized(CImgList<T> &correlogram, CImgList<T> &F, CImgList<T> &G)
    {// phase calcul from division by the norm
      
      float maxF(correlogram[0].max());
      
      cimglist_for(F,L)
      {//phase
        cimg_forXYZ(F[L],x,y,z)
        {

          if (correlogram(0,x,y,z)<maxF*m_factor/100 or correlogram(1,x,y,z)<maxF*m_factor/100)
          {//smoothing in Fourier space of low magnitude signal
            F(L,x,y,z) = 0;
            G(L,x,y,z) = 0;
          }
          else
          {
            F(L,x,y,z) = F(L,x,y,z)/correlogram(0,x,y,z);
            G(L,x,y,z) = G(L,x,y,z)/correlogram(1,x,y,z);
          }
          
        }
      }
      
    }

};

#endif