TileSphere.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
**
**************************************************************************************/
 

 
#pragma once
 
#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(const float* point1, const float* point2)
  {
    const float d1 = * point1      - * point2;
    const float d2 = *(point1 + 1) - *(point2 + 1);
    const float d3 = *(point1 + 2) - *(point2 + 2);
    return sqrtf(d1*d1 + d2*d2 + d3*d3);
  }

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