AssocGraph.h

Go to the documentation of this file.

/*************************************************************************************
**
** Copyright 2007, MeVis Medical Solutions AG
**
** The user may use this file in accordance with the license agreement provided with
** the Software or, alternatively, in accordance with the terms contained in a
** written agreement between the user and MeVis Medical Solutions AG.
**
** For further information use the contact form at https://www.mevislab.de/contact
**
**************************************************************************************/
 
#ifndef ASSOC_GRAPH_H
#define ASSOC_GRAPH_H
 

 
#include "mlVesselGraphSystem.h"
#include "mlMatrixTemplate.h"
 
#include "VesselGraphDisableWarnings.h"
#include <list>
#include <vector>
#include <valarray>
#include "VesselGraphRestoreWarnings.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