Ccorrelation.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
#define CCORRELATION

//-----------------------------CIMG_NETCDF-------------------------------
#define cimg_plugin  "plugins/add_fileformat.h"
#define cimg_use_netcdf

#ifdef cimg_use_netcdf
 #include "../NetCDF.Tool/struct_parameter_NetCDF.h"
// #define     cimg_plugin1 "plugins/netcdf_file_format4CImg1.h"
// #define cimglist_plugin1 "plugins/netcdf_file_format4CImgList1.h"
 #define     cimg_plugin2 "plugins/netcdf_file_format4CImg2.h"
 #define cimglist_plugin2 "plugins/netcdf_file_format4CImgList2.h"
 #include "../CImg.Tool/CImg_NetCDF.h"
#endif
//-----------------------------CIMG_NETCDF-------------------------------


#include "Cimage.h"
#include "Cfield.h"
#include "Cmesh_factory.h"
#include "CNpeak.h"
#include "CshapeFunction_factory.h"

template<typename T, typename Timg> 
class Ccorrelation
{
  
  public:
  
    std::string class_name, m_correl_type, m_correl_name;
    T m_current_res;// mean residue between ref and corrected deformed image
    int m_threadNB, _3D;// storage mode activation, thread number & if 3D
    int m_medianFilter;

    bool m_verbose;

    
    Ccorrelation()
    {// constructor

      class_name = "Correlation : mother";
      
      #pragma omp parallel
      {//total number of thread required for reduction process
        m_threadNB = omp_get_num_threads();
      }
      
    }
 
 
 
    virtual void init(CParameterNetCDF &fp) = 0;
    
    virtual void exec(const Cimage<Timg> &oImage, const Cmesh<T,Timg> &oMesh, Cfield<T,Timg> &oField) = 0;//launch correlation method

    virtual void regularize(const Cmesh<T,Timg> &oMesh, const Cimage<Timg> &oImage, CImgList<T> &solution)
    {// apply median filter to intermediate fields
      
      cimglist_for(solution,var)
      {//median filter on fields
        solution[var].blur_median(m_medianFilter);
      }
      
    }
    
    void count(const int &increment, const int &dim, int &counter)
    {// counter

      if (100-increment*100/dim<=counter-10)
      {
        std::cout<<"remaining : "<<(int)(100.-increment*100./dim)<<"\%\n";
        counter = 100-increment*100/dim;
      }
      
    }
  

  
};

#endif