CNpeak.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 CNPEAK
#define CNPEAK



//-----------------------------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 "Cpeak_factory.h"
#include "CImg_math.h"
  
template<typename T> class CNpeak 
{
  
  public:
    
    std::string class_name;
    std::vector<Cpeak<T>*> m_Npeak;//attribute : vector of Cpeak object pointer
    CImgList<T> m_Nstat;// [flag,snr,x,y,z](0,0,0,peaks)
    float m_snr;// Signal Noise Ratio
    int m_nbPeak, m_dim; //number of peak required, spatial dimensions
    
    bool m_verbose;
    
    
    CNpeak(std::string &type, const int &number=3, const int &dim=2)
    {// constructor
      
      class_name = "CNpeak";
      
      m_nbPeak = number;//number of peak to detect
      m_dim = dim;// spatial dimensions
      m_snr = 0;// init of Signal Noise Ratio
      
      init(type);

    }
   
    ~CNpeak()
    {// destructor
      
      for (int i=0; i< m_Npeak.size(); i++)
      {
        delete m_Npeak[i];
      };
      
    }
   
    void init(std::string &type)
    {// init

      m_Npeak.clear();

      for (int i=0; i<=m_nbPeak; i++)
      {//m_nbPeak + 1 for SNR calcul
        Cpeak_factory<T> peak_factory("peak");
        Cpeak<T>* pPeak=peak_factory.create(type);
        pPeak->m_verbose = m_verbose;
        m_Npeak.push_back(pPeak);
      }

      m_Nstat.assign(m_dim+2,1,1,1,m_nbPeak+1,0); //[flag,snr,x,y,z](0,0,0,peaks)
      
    }
  
  
  
    void store(const int &i, const int &j)
    {// store peak stats
        m_Nstat(j,0,0,0,i) = (*m_Npeak[i]).content()[j];
    }

    void exec(CImg<T> &img)
    {// exec 
      
      for (int i=0; i<m_Npeak.size(); i++)
      {// ieme execution of Cpeak->exec() method
      
        m_Npeak[i]->exec(img);

        for (int j=0; j<m_Nstat.size(); j++)
        {// store results
          store(i,j);
        }
        
      }
      
      snr(m_Nstat);
      
    }
    
    void snr(CImgList<T> snRatio)
    {
     
      for (int i=0; i<m_nbPeak; i++)
      {// Signal Noise Ratio calcul -> 1 - peak_i/peak_i+1 (0 = poor snr, <<0 = good snr)
        snRatio(1,0,0,0,i) = snRatio(1,0,0,0,i+1)*100/snRatio(1,0,0,0,i);
      }
      
    }


};

#endif