MeVisLab/Standard/Sources/Shared/MLPointCloudUtils/MLMinimalDistancePointClouds/TileSphere.h

Go to the documentation of this file.

// **InsertLicense** code
//----------------------------------------------------------------------------------

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


#ifndef __TileSphere_H
#define __TileSphere_H

#include "../MLPointCloudUtilsSystem.h"

#include <mlVector3.h>


ML_START_NAMESPACE

class MinimalDistancePointClouds;

class TileSphere 
{
public:

  TileSphere();  
  TileSphere(MinimalDistancePointClouds* minimalDistance, int maxPartition = 2, int numElements = 1103, float error = 0);
  virtual ~TileSphere();

  void setParameter(MinimalDistancePointClouds* minimalDistance, int maxPartition, int numElements, float error);

  void setBB(float  minX, float  maxX, float  minY, float  maxY, float  minZ, float  maxZ);
  void addPoint(float* position);
  void getDeFactoSize();
  float computeDistance(TileSphere* tileSphere, float*& point1, float*& point2);

  float computeDistance(TileSphere* tileSphere, unsigned int& node1, unsigned int& nodet2);


  unsigned int getSphereIndex();

  void getStatistics(int &_treeDepth, int &usedSpheres, int &leafSpheres, int &_maxEnt);


private:

  void _initSphere(MinimalDistancePointClouds* minimalDistance, int partition, int numElements, float error);
  MinimalDistancePointClouds* _minimalDistance;
  TileSphere* _tileSpheres;

  unsigned int _sphereIndex;

  float _error;

  bool _hasSubSpheresFlag;
  unsigned int _partition;
  unsigned int _cubicPartition;


  float** _subset;
  unsigned int _numEntries;
  int _maxEntries;

  float _minX, _maxX, _minY, _maxY, _minZ, _maxZ;
  float _position[3];
  float _radius;

  void _getBB(float& minX, float& maxX, float& minY, float& maxY, float& minZ, float& maxZ);
  float* _getPosition();
  float _getRadius();
  void _getSubset(float** &subset, int &size);
  bool _hasSubSpheres();

  TileSphere* _getPreciserSphere(TileSphere* referenceSphere, float distance);

  float _distance(TileSphere* tileSphere1, TileSphere* tileSphere2);

  inline float _distance(register float* point1, register float* point2) {
    register float d1 = * point1      - * point2;
    register float d2 = *(point1 + 1) - *(point2 + 1);
    register float d3 = *(point1 + 2) - *(point2 + 2);
    return sqrtf(d1*d1 + d2*d2 + d3*d3);
  }

  inline float _fastDistance(register float* point1, register float* point2) {
    register float d1 = * point1      - * point2;
    register float d2 = *(point1 + 1) - *(point2 + 1);
    register float d3 = *(point1 + 2) - *(point2 + 2);
    return (d1*d1 + d2*d2 + d3*d3);
  }
};

ML_END_NAMESPACE

#endif // __TileSphere_H