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<AimsData<T> > : public BorderIterator
  {
  public:
    Border3dIteratorOf(const AimsData<T> & data, const Point3d & width);
    Border3dIteratorOf(const AimsData<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:
    AimsData<T> _data;
    Point3d _width;
    Point3d _current;
    Point3d _lastPoint;
  };

  // Border extraction
  template< class T >
  inline AimsData< T > border( AimsData< T > &data, Point3d border )
  {
    std::vector< Point3d > v;

    Border3dIteratorOf<AimsData<T> > b(data, border);
    while(b.isValid())
    {
      v.push_back(b.value());
      b.next();
    }

    AimsData< 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 AimsData< T > border( AimsData< T > &data, uint32_t width )
  {
    return border(data, Point3d(width, width, width));
  }
}

#endif