MeVisLab/Standard/Sources/ML/MLVesselGraph/AssocGraph.h

Go to the documentation of this file.

// **InsertLicense** code

// $Id: AssocGraph.h 50846 2007-02-13 17:59:51Z okonrad $
// $Source$
//--------------------------------------------------------------------------------------
#ifndef __AssocGraph_h__
#define __AssocGraph_h__


#include "mlVesselGraphSystem.h"
#include "mlMatrixTemplate.h"
#include "mlSystemWarningsDisable.h"
#include <list>

#include <valarray>
#include "mlSystemWarningsRestore.h"

ML_START_NAMESPACE



typedef std::vector<double>       CDouble1Vec;
typedef std::vector<CDouble1Vec>  CDouble2Vec;
typedef std::vector<CDouble2Vec>  CDoubleArray;  

#define N_XPPOS   0
#define N_YPPOS   1
#define N_ZPPOS   2
#define N_RL      3

#define N_EDGEA   0
#define N_EDGEB   1
#define N_EDGENO  2
#define N_DIST    3

#define DELTA_POS 0.1

//------------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------------
class PointMap
{
public:
  int _nP1;                         
  int _nP2;                         


  inline bool operator<(const PointMap& pm) const { if(_nP1 != pm._nP1){return _nP1 < pm._nP1;} else {return _nP2 < pm._nP2;} }
  inline bool operator==(const PointMap& pm) const { return ((_nP1 == pm._nP1) && (_nP2 == pm._nP2)); }
};


//------------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------------
class PathInfo
{
public:
  std::valarray<bool> _vPath;       
  double              _dLength;     
  int                 _nIdx;        

  void           *    _pNode;       

  PathInfo(void);                   
  PathInfo(size_t size);            

};




//------------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------------
class IsomorphismMap
{
  size_t _N1;                       
  size_t _N2;                       

public:
  std::valarray<bool>  _w1;         
  std::valarray<bool>  _w2;         
  std::list<PointMap>  _LMap;       


  IsomorphismMap(void);
  IsomorphismMap(size_t N1, size_t N2);
  IsomorphismMap(size_t N1, size_t N2, const PointMap&);
  IsomorphismMap(const IsomorphismMap&);


  IsomorphismMap& operator=(const IsomorphismMap&);       
  bool operator==(const IsomorphismMap&) const;                 
  bool operator<(const IsomorphismMap&) const;                  



  bool addPointMap          (const PointMap& pm);         
  bool addPointMapAscending (const PointMap& pm, int n);  
  void removePointMap       (const PointMap& pm);         
  void clear ();                                          
  void resize (size_t N1, size_t N2);                     

};



//------------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------------
class VESSELGRAPH_EXPORT AssocGraph
{
  MatrixSizedTemplate<double, 0, 0> _mTAG;        
  std::list<PointMap>       _LMapTable;           

  std::list<IsomorphismMap> _LIsoMap;             

  size_t _N1;                                     
  size_t _N2;                                     


public:
  //--------------------------------------------------------------------------------------------

  typedef std::list<PointMap>::iterator MapIterator;
  typedef std::list<PointMap>::const_iterator constMapIterator;
  typedef std::list<IsomorphismMap>::iterator IsoIterator;
  typedef std::list<IsomorphismMap>::const_iterator constIsoIterator;
  //--------------------------------------------------------------------------------------------


  //--------------------------------------------------------------------------------------------






  bool calcAssocMatrix(const MatrixTemplate<double>& mDist1,  
                       const MatrixTemplate<double>& mDist2,  
                       double dSigma,                         
                       double dScale                          
                       );

  bool calcAssocMatrix(const CDoubleArray & AClass1,  
                       const CDoubleArray & AClass2,  
                       double dSigma,                 
                       double dScale                  
                      );

  bool calcMaxClique (double dAssocLevel);

  bool calcAutoMaxClique(double *pdAssocLevel, double dStep=0.1);

  bool expandMap(int nIdx);

  bool expandMap(int nIdx, double dAssocLevel);
  //--------------------------------------------------------------------------------------------


  //--------------------------------------------------------------------------------------------

  std::string printMapInformation ();                                   

  std::list<PointMap>       * getPointMap () { return &_LMapTable; }    
  std::list<IsomorphismMap> * getIsoMap   () { return &_LIsoMap; }      

  //--------------------------------------------------------------------------------------------

};



ML_END_NAMESPACE

#endif