splineresampler.h

Go to the documentation of this file.

/* This software and supporting documentation are distributed by
 *     Institut Federatif de Recherche 49
 *     CEA/NeuroSpin, Batiment 145,
 *     91191 Gif-sur-Yvette cedex
 *     France
 *
 * This software is governed by the CeCILL-B license under
 * French law and abiding by the rules of distribution of free software.
 * You can  use, modify and/or redistribute the software under the
 * terms of the CeCILL-B license as circulated by CEA, CNRS
 * and INRIA at the following URL "http://www.cecill.info".
 *
 * As a counterpart to the access to the source code and  rights to copy,
 * modify and redistribute granted by the license, users are provided only
 * with a limited warranty  and the software's author,  the holder of the
 * economic rights,  and the successive licensors  have only  limited
 * liability.
 *
 * In this respect, the user's attention is drawn to the risks associated
 * with loading,  using,  modifying and/or developing or reproducing the
 * software by the user in light of its specific status of free software,
 * that may mean  that it is complicated to manipulate,  and  that  also
 * therefore means  that it is reserved for developers  and  experienced
 * professionals having in-depth computer knowledge. Users are therefore
 * encouraged to load and test the software's suitability as regards their
 * requirements in conditions enabling the security of their systems and/or
 * data to be ensured and,  more generally, to use and operate it in the
 * same conditions as regards security.
 *
 * The fact that you are presently reading this means that you have had
 * knowledge of the CeCILL-B license and that you accept its terms.
 */


#ifndef AIMS_RESAMPLING_SPLINERESAMPLER_H
#define AIMS_RESAMPLING_SPLINERESAMPLER_H

#include <aims/resampling/resampler.h>
#include <cartobase/type/datatypetraits.h>

#include <cmath>
#include <vector>

namespace aims
{

// This implementation helper allows to change the implementation of methods
// depending if the data type is single-channel or multi-channel.
template<bool, typename>
struct MultiChannelResamplerSwitch;


template <class T>
class SplineResampler : public Resampler< T >
{
public:

  SplineResampler();
  ~SplineResampler() { }

  virtual int getOrder() const = 0;

  carto::VolumeRef<double> getSplineCoef( const carto::Volume< T >& inVolume,
                                          int t = 0,
                                          bool verbose = false );

  void reset();

  // Overridden for performance in case of multi-channel data
  void resample_inv_to_vox( const carto::Volume< T >& input_data,
                            const soma::Transformation3d&
                              inverse_transform_to_vox,
                            const T& background,
                            carto::Volume< T > & output_data,
                            bool verbose = false ) const CARTO_OVERRIDE;
  using Resampler<T>::resample_inv_to_vox;

  // Implements iteration over channels for multi-channel data types
  template<bool, typename> friend struct MultiChannelResamplerSwitch;

  typedef typename carto::DataTypeTraits<T>::ChannelType ChannelType;

  void resample_channel_inv_to_vox(
    const carto::Volume< ChannelType >& inVolume,
    const soma::Transformation3d& inverse_transform_to_vox,
    const ChannelType& outBackground,
    carto::Volume< ChannelType > & outVolume,
    bool verbose) const;
protected:

  void doResample( const carto::Volume< T > &inVolume,
                   const soma::Transformation3d &transform3d,
                   const T &outBackground,
                   const Point3df &outLocation,
                   T &outValue,
                   int t ) const CARTO_OVERRIDE;

  virtual void doResampleChannel( const carto::Volume< ChannelType > &inVolume,
                                  const soma::Transformation3d &transform3d,
                                  const ChannelType &outBackground,
                                  const Point3df &outLocation,
                                  ChannelType &outValue,
                                  int t ) const;

  void updateParameters( const carto::Volume< T >& inVolume, int t,
                         bool verbose ) const CARTO_OVERRIDE;
  virtual void updateParametersChannel( const carto::Volume< ChannelType >& inVolume, int t,
                         bool verbose ) const;
  void iirConvolveMirror( std::vector< double >& data ) const;

  int getFold( int i, int size ) const;

  virtual double getBSplineWeight( int i, double x ) const = 0;

  std::vector<double>  _poles;
  double               _gain;

  // These three mutable members handle the cache of spline coefficients
  mutable carto::VolumeRef<double>       _splineCoefficients;
  mutable const carto::Volume<ChannelType> *  _lastvolume;
  mutable int                            _lasttime;
};


  inline int mirrorCoeff( int i, int size )
  {
    i = std::abs( i );
    if ( i < size )
    {

      return i;

    }
    if ( size == 1 )
    {

      return 0;

    }
    int size2 = ( size * 2 ) - 2;
    int rem = i % size2;
    return ( rem < size ) ? rem : ( size2 - rem );
  }

} // namespace aims

#endif