seaview::seaDBSCAN Class Reference

#include <seaDBSCAN.h>

Public Member Functions
	seaDBSCAN (double in_eps, int in_minpts)
	constructor More...
	~seaDBSCAN ()
	Default destructor. More...
std::vector< int >	Scan2D (std::vector< std::vector< double >> &pts)
std::vector< std::vector < double > >	GetNeighbours (size_t i, std::vector< std::vector< double >> &pts, bool)
int	UnionSets (std::vector< std::vector< double >> &seed, std::vector< std::vector< double >> &pts)
double	SimpleDist (double w1, double t1, double w2, double t2)

Public Attributes
double	m_eps
int	m_minpts

Detailed Description

Definition at line 23 of file seaDBSCAN.h.

Constructor & Destructor Documentation

seaview::seaDBSCAN::seaDBSCAN ( double  in_eps, int  in_minpts  )

inline

constructor

Definition at line 30 of file seaDBSCAN.h.

: m_eps(in_eps), m_minpts(in_minpts) {}

seaview::seaDBSCAN::~seaDBSCAN ( )

inline

Default destructor.

Definition at line 33 of file seaDBSCAN.h.

{}

Member Function Documentation

std::vector< std::vector< double > > seaview::seaDBSCAN::GetNeighbours	(	size_t	i,
		std::vector< std::vector< double >> &	pts,
		bool	include_self
	)

Definition at line 58 of file seaDBSCAN.cxx.

                                                                                                                        {
    std::vector<std::vector<double>> neighbours;
    std::vector<double> point = pts[ipoint];

    //VERY simple, will update soon to a DB 


    for(size_t ip=0; ip<pts.size(); ip++){
      std::vector<double> p = pts[ip];

      double dist = this->SimpleDist( p[0], p[1], point[0], point[1]);

      if(include_self){
        if(dist <= m_eps){
          std::vector<double> tp = {p[0],p[1],(double)ip};//Push back original index too
          neighbours.push_back(tp);
        }
      }else{
        if(dist <= m_eps && p != point ){
          std::vector<double> tp = {p[0],p[1],(double)ip};//Push back original index too
          neighbours.push_back(tp);
        }

      }
    }
    return neighbours;
  }

std::vector< int > seaview::seaDBSCAN::Scan2D

(

std::vector< std::vector< double >> &

pts

)

Definition at line 5 of file seaDBSCAN.cxx.

                                                                   {

    int cluster_count = 0;//we have no clusters
    size_t N = pts.size();
    int l_undef = -99;
    int l_noise = 0;
    std::vector<int> label(N,l_undef);

    for(size_t i=0; i<N; i++){
      if(label[i]!=l_undef) continue;
      std::vector<std::vector<double>> neighbours = this->GetNeighbours(i,pts,false);
      //std::cout<<i<<" has neightbours "<<neighbours.size()<<std::endl;

      if((int)neighbours.size()+1 < m_minpts){ // if there is less than minpts, its a noise point
        label[i]= l_noise;
        //  std::cout<<i<<" thats less than min, noise"<<std::endl;
        continue;
      }

      cluster_count+=1;
      label[i] = cluster_count;


      std::vector<std::vector<double>> seed_set = neighbours;
      for(size_t q=0; q<seed_set.size(); q++){
        size_t iq = (size_t)seed_set[q][2]; //This is original 

        if(label[iq]==l_noise){
          label[iq] = cluster_count;//Change noise to a border point;
        }

        if(label[iq]!=l_undef){
          continue; // previously processed, already identified to a cluster
        }

        // if label[iq] is l_undef
        label[iq]=cluster_count;// wasn't noise, new point, add to cluster

        std::vector<std::vector<double>> new_neighbours = this->GetNeighbours(iq,pts,true);//Get neighbours of this point, including itslef eh.

        if((int)new_neighbours.size() >= m_minpts ){ //expand the seed set
          //Guanqun: change vector while looping over its elements
          //new elements are pushed back to seed_set
          this->UnionSets(seed_set, new_neighbours); 
        }
      }
    }
    return label;

  }

double seaview::seaDBSCAN::SimpleDist	(	double	w1,
		double	t1,
		double	w2,
		double	t2
	)

Definition at line 110 of file seaDBSCAN.cxx.

                                                                        {
    // Guanqun: wire and tick conversion to distance??, if so should only be done on plane 2
    double wire_con = 0.3;
    double tick_con = 1.0/25.0;

    return  sqrt(pow(w1*wire_con-w2*wire_con,2)+pow(t1*tick_con-t2*tick_con,2));

  }

int seaview::seaDBSCAN::UnionSets	(	std::vector< std::vector< double >> &	seed,
		std::vector< std::vector< double >> &	pts
	)

Definition at line 86 of file seaDBSCAN.cxx.

                                                                                               {

    //VERY simple, will update soon if it works
    for(auto &p:pts){

      bool is_in=false;
      for(auto &s:seed){
        if(s==p){
          is_in = true;
          break;
        }
      }

      if(is_in == false){
        seed.push_back(p);
      }

    }



    return 0;
  }

Member Data Documentation

double seaview::seaDBSCAN::m_eps

Definition at line 26 of file seaDBSCAN.h.

int seaview::seaDBSCAN::m_minpts

Definition at line 27 of file seaDBSCAN.h.

---

The documentation for this class was generated from the following files:

• seaDBSCAN.h
• seaDBSCAN.cxx