packages/vtk/src/bbvtkPiecewiseFunction.cxx

   1 /*
   2  # ---------------------------------------------------------------------
   3  #
   4  # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
   5  #                        pour la SantÈ)
   6  # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
   7  # Previous Authors : Laurent Guigues, Jean-Pierre Roux
   8  # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
   9  #
  10  #  This software is governed by the CeCILL-B license under French law and
  11  #  abiding by the rules of distribution of free software. You can  use,
  12  #  modify and/ or redistribute the software under the terms of the CeCILL-B
  13  #  license as circulated by CEA, CNRS and INRIA at the following URL
  14  #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  15  #  or in the file LICENSE.txt.
  16  #
  17  #  As a counterpart to the access to the source code and  rights to copy,
  18  #  modify and redistribute granted by the license, users are provided only
  19  #  with a limited warranty  and the software's author,  the holder of the
  20  #  economic rights,  and the successive licensors  have only  limited
  21  #  liability.
  22  #
  23  #  The fact that you are presently reading this means that you have had
  24  #  knowledge of the CeCILL-B license and that you accept its terms.
  25  # ------------------------------------------------------------------------ */
  26
  27
  28 /*=========================================================================
  29   Program:   bbtk
  30   Module:    $RCSfile: bbvtkPiecewiseFunction.cxx,v $
  31   Language:  C++
  32   Date:      $Date: 2012/11/16 08:51:58 $
  33   Version:   $Revision: 1.3 $
  34 =========================================================================*/
  35
  36 /**
  37  *  \file
  38  *  \brief
  39  */
  40
  41
  42
  43 #ifdef _USE_VTK_
  44
  45
  46 #include "bbvtkPiecewiseFunction.h"
  47 #include "bbvtkPackage.h"
  48
  49 namespace bbvtk
  50 {
  51
  52
  53    BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,PiecewiseFunction)
  54    BBTK_BLACK_BOX_IMPLEMENTATION(PiecewiseFunction,bbtk::AtomicBlackBox);
  55
  56
  57         //---------------------------------------------------------------------
  58    void PiecewiseFunction::bbUserSetDefaultValues()
  59    {
  60      // Identity function
  61      std::vector<float> x,y;
  62      x.push_back(0);
  63      x.push_back(1);
  64      y.push_back(0);
  65      y.push_back(1);
  66      bbSetInputX(x);
  67      bbSetInputY(y);
  68          bbSetOutputOut(0);
  69    }
  70
  71
  72         //---------------------------------------------------------------------
  73    void PiecewiseFunction::bbUserInitializeProcessing()
  74    {
  75      bbSetOutputOut(vtkPiecewiseFunction::New());
  76    }
  77
  78         //---------------------------------------------------------------------
  79    void PiecewiseFunction::bbUserFinalizeProcessing()
  80    {
  81            if (bbGetOutputOut()!=NULL)
  82            {
  83                    bbGetOutputOut()->Delete();
  84                    bbSetOutputOut(0);
  85            }
  86    }
  87 //---------------------------------------------------------------------
  88
  89    void PiecewiseFunction::Process()
  90    {
  91      if ( ( bbGetInputStatus("X") != bbtk::UPTODATE ) ||
  92           ( bbGetInputStatus("Y") != bbtk::UPTODATE ) )
  93        {
  94          if  ( bbGetInputX().size() != bbGetInputY().size() )
  95            bbtkError(bbGetFullName()
  96                      <<" : input vectors X and Y do not have the same size");
  97          if  ( bbGetInputX().size() < 2 )
  98            bbtkError(bbGetFullName()
  99                      <<" : input vectors X and Y have a size < 2");
 100          bbGetOutputOut()->RemoveAllPoints ();
 101          //      std::cout << bbGetInputX().size()<< ","<< bbGetInputY().size()<< std::endl;
 102          std::vector<float>::const_iterator x,y;
 103          std::vector<float> X = bbGetInputX();
 104          std::vector<float> Y = bbGetInputY();
 105          for (x = X.begin(), y = Y.begin();
 106               x != X.end(),  y != Y.end();
 107               ++x, ++y )
 108            {
 109              //              std::cout << *x << " -> " << *y << std::endl;
 110              bbGetOutputOut()->AddPoint(*x,*y);
 111            }
 112        }
 113    }
 114
 115 }//namespace bbtk
 116
 117 #endif // _USE_VTK_
 118