connectivity.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_CONNECTIVITY_CONNECTCLASS_H
#define AIMS_CONNECTIVITY_CONNECTCLASS_H

#include <string>

#include <aims/config/aimsdata_config.h>
#include <aims/vector/vector.h>
#include <aims/bucket/bucketMap.h>  // for CoordinatesLess operator
#include <set>


namespace aims
{

  class Connectivity
  {
  public:
    enum Type
      {
        CONNECTIVITY_4_XY,
        CONNECTIVITY_4_XZ,
        CONNECTIVITY_4_YZ,
        CONNECTIVITY_6_XYZ,
        CONNECTIVITY_8_XY,
        CONNECTIVITY_8_XZ,
        CONNECTIVITY_8_YZ,
        CONNECTIVITY_18_XYZ,
        CONNECTIVITY_26_XYZ,
        CONNECTIVITY_9_XY_Zminus,
        CONNECTIVITY_9_XY_Zplus,
        CONNECTIVITY_9_XZ_Yminus,
        CONNECTIVITY_9_XZ_Yplus,
        CONNECTIVITY_9_YZ_Xminus,
        CONNECTIVITY_9_YZ_Xplus,
        CONNECTIVITY_5_XminusY,
        CONNECTIVITY_5_XplusY,
        CONNECTIVITY_5_XYminus,
        CONNECTIVITY_5_XYplus,
        CONNECTIVITY_5_XminusZ,
        CONNECTIVITY_5_XplusZ,
        CONNECTIVITY_5_XZminus,
        CONNECTIVITY_5_XZplus,
        CONNECTIVITY_5_YminusZ,
        CONNECTIVITY_5_YplusZ,
        CONNECTIVITY_5_YZminus,
        CONNECTIVITY_5_YZplus,
        CONNECTIVITY_8_XYZ,
        CONNECTIVITY_4_XYdiag,
        CONNECTIVITY_4_XZdiag,
        CONNECTIVITY_4_YZdiag,
      };

    static Type type_from_string(const std::string &);
    static std::string type_to_string(Type);

    Connectivity( int oline, int oslice, Type type);
    virtual ~Connectivity();
    
    Type type() const { return _type; }
    std::string type_string() const { return type_to_string(_type); }
    int nbNeighbors() const { return _nbNeighbors; }
    int offset( int n ) const { return _offset[n]; }
    const Point3d& xyzOffset( int n ) const { return _xyzOffset[n]; }
    const Point3df& dir( int n ) const { return _dir[n]; }
    bool isNeighbor( const Point3d & offset )
    { return _xyzOffsets.find( offset ) != _xyzOffsets.end(); }

  Connectivity& 
  operator=(const Connectivity& other)
  {
    if( &other == this )
      return *this ;
    
    _xSize = other._xSize ;
    _sliceSize = other._sliceSize ;
    _type = other._type ;
    
    if( _xyzOffset ){
      delete[] _xyzOffset;
      _xyzOffset = 0 ;
    }
    
    if(_offset){
      delete[] _offset;
      _offset = 0 ;
    }
    
    if( _dir ){
      delete[] _dir;
      _dir = 0 ;
    }
    setConnectivity() ;
    return *this ;
  }

  protected: 
    int _xSize;
    int _sliceSize;
    Type _type;
    int _nbNeighbors;
    int *_offset;
    Point3d* _xyzOffset;
    std::set<Point3d, BucketMapLess> _xyzOffsets;
    Point3df* _dir;

    void setConnectivity();
    void setConnectivity_4_XY();
    void setConnectivity_4_XZ();
    void setConnectivity_4_YZ();
    void setConnectivity_6_XYZ();
    void setConnectivity_8_XY();
    void setConnectivity_8_XZ();
    void setConnectivity_8_YZ();
    void setConnectivity_18_XYZ();
    void setConnectivity_26_XYZ();
    void setConnectivity_9_XY_Zminus();
    void setConnectivity_9_XY_Zplus();
    void setConnectivity_9_XZ_Yminus();
    void setConnectivity_9_XZ_Yplus();
    void setConnectivity_9_YZ_Xminus();
    void setConnectivity_9_YZ_Xplus();
    void setConnectivity_5_XminusY();
    void setConnectivity_5_XplusY();
    void setConnectivity_5_XYminus();
    void setConnectivity_5_XYplus();
    void setConnectivity_5_XminusZ();
    void setConnectivity_5_XplusZ();
    void setConnectivity_5_XZminus();
    void setConnectivity_5_XZplus();
    void setConnectivity_5_YminusZ();
    void setConnectivity_5_YplusZ();
    void setConnectivity_5_YZminus();
    void setConnectivity_5_YZplus();
  };

}


#endif