border.h

Go to the documentation of this file.

/* Copyright (C) 2000-2013 CEA
 *
 * This software and supporting documentation were developed by
 *     bioPICSEL
 *     CEA/DSV/I²BM/MIRCen/LMN, Batiment 61,
 *     18, route du Panorama
 *     92265 Fontenay-aux-Roses
 *     France
 */
 
#ifndef BRAINRAT_UTILITY_BORDER_H
#define BRAINRAT_UTILITY_BORDER_H
 
#include <vector>
#include <cartobase/smart/rcptr.h>
#include <aims/data/data_g.h>
#include <aims/data/data_g.h>
 
#define DEFAULT_BORDER     10        //   Around
 
namespace brainrat
{
 
  class BorderIterator : public carto::RCObject
  {
  public:
    virtual inline ~BorderIterator() {}
 
    // Returns the current point in border
    virtual const Point3d & value() const = 0;
    virtual void next() = 0;
    virtual bool isValid() const = 0;
    virtual void restart() = 0;
 
    // Check if a point (in voxel referential) belongs to the border
    virtual bool contains( const Point3d & ) const = 0;
    virtual Point3d & width() = 0;
  };
 
  //---------------------------------------------------------------------------
  template <class T>
  class Border3dIteratorOf : public BorderIterator
  {
  };
 
  //---------------------------------------------------------------------------
  template< class T >
  class Border3dIteratorOf<carto::VolumeRef<T> > : public BorderIterator
  {
  public:
    Border3dIteratorOf(const carto::VolumeRef<T> & data, const Point3d & width);
    Border3dIteratorOf(const carto::VolumeRef<T> & data, const int32_t width);
    virtual inline ~Border3dIteratorOf() {}
 
    // Returns the current point in border
    virtual const Point3d & value() const;
    virtual void next();
    virtual bool isValid() const;
    virtual void restart();
 
    // Check if a point (in voxel referential) belongs to the border
    virtual bool contains( const Point3d & ) const;
    virtual Point3d & width();
 
  private:
    carto::VolumeRef<T> _data;
    Point3d _width;
    Point3d _current;
    Point3d _lastPoint;
  };
 
  // Border extraction
  template< class T >
  inline carto::VolumeRef< T > border( carto::VolumeRef< T > &data, Point3d border )
  {
    std::vector< Point3d > v;
 
    Border3dIteratorOf<carto::VolumeRef<T> > b(data, border);
    while(b.isValid())
    {
      v.push_back(b.value());
      b.next();
    }
 
    carto::VolumeRef<T > r(v.size());
 
    //std::cout << "border size : " << carto::toString(v.size()) << std::endl << std::flush;
    for ( int32_t i = 0; i < v.size(); i++ )
      r(i) = data(v[i]);
 
    return r;
  }
 
  template< class T >
  inline carto::VolumeRef< T > border( carto::VolumeRef< T > &data, uint32_t width )
  {
    return border(data, Point3d(width, width, width));
  }
}
 
#endif