Cformat_char_double.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 CFORMAT_CHAR_FLOAT
#define CFORMAT_CHAR_FLOAT

#include "Cformat.h"

class Cformat_char_double
{

 
  Cformat_char_double()
  {// constructor
    
    class_name = "Format : char double";

  }    

  virtual void exec(int argc, char *argv[])
  {

    //output name
    std::string ofile = cimg_option("-o","./","ofile");
    
    //create image container object
    Cimage_factory<char> image_factory("image_list");
    Cimage<char>* pImage = image_factory.create(argc,argv);
  
    //create sequence correlation object
    Csequence<double,char>* pSequence = new Csequence<double,char>;
    pSequence->init(argc, argv);

    for (int s=0; s<pImage->num_samples(); s++)
    {

      frame_number = s;
    
      pImage->exec(argc, argv, s);//load images
    
      //Test if the image sequence is (A-B, B-C...) or (A-B, C-D...) as in PIV. In the last case, pSequence->m_dataField must be clear. Otherwise, the second correlation sequence will be treated with the previous field as init.
      if (pImage->m_mode != DEPENDENT){pSequence->m_modes.clear();}
      pSequence->exec(argc, argv, *pImage);//perform correlation sequences
    
      save(*pSequence, *pImage, ofile, s);
      
    }//end image list
  
    CImg<T> res(pSequence->m_meanRes.size(),1,1,1,0);
    cimg_forX(res,step){res[step] = pSequence->m_meanRes[step];}
    pImage->save_netcdf(res, ofile + "residue_mean", 0);
  
    std::cout<<"calcul done\n";

    delete pImage;
    delete pSequence;

  }
  
  
  
  
};

#endif