VTK
vtkOctreePointLocator.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkOctreePointLocator.h
5 
6  Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
7  All rights reserved.
8  See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
9 
10  This software is distributed WITHOUT ANY WARRANTY; without even
11  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12  PURPOSE. See the above copyright notice for more information.
13 
14 =========================================================================*/
15 /*----------------------------------------------------------------------------
16  Copyright (c) Sandia Corporation
17  See Copyright.txt or http://www.paraview.org/HTML/Copyright.html for details.
18 ----------------------------------------------------------------------------*/
19 
40 #ifndef vtkOctreePointLocator_h
41 #define vtkOctreePointLocator_h
42 
43 #include "vtkCommonDataModelModule.h" // For export macro
45 
46 class vtkCellArray;
47 class vtkIdTypeArray;
49 class vtkPoints;
50 class vtkPolyData;
51 
52 class VTKCOMMONDATAMODEL_EXPORT vtkOctreePointLocator : public vtkAbstractPointLocator
53 {
54 public:
56  void PrintSelf(ostream& os, vtkIndent indent) override;
57 
58  static vtkOctreePointLocator *New();
59 
61 
64  vtkSetMacro(MaximumPointsPerRegion, int);
65  vtkGetMacro(MaximumPointsPerRegion, int);
67 
69 
72  vtkSetMacro(CreateCubicOctants, int);
73  vtkGetMacro(CreateCubicOctants, int);
75 
77 
83  vtkGetMacro(FudgeFactor, double);
84  vtkSetMacro(FudgeFactor, double);
86 
88 
92  double *GetBounds() override;
93  void GetBounds(double *bounds) override;
95 
97 
100  vtkGetMacro(NumberOfLeafNodes, int);
102 
106  void GetRegionBounds(int regionID, double bounds[6]);
107 
111  void GetRegionDataBounds(int leafNodeID, double bounds[6]);
112 
116  int GetRegionContainingPoint(double x, double y, double z);
117 
123  void BuildLocator() override;
124 
126 
130  vtkIdType FindClosestPoint(const double x[3]) override;
131  vtkIdType FindClosestPoint(double x, double y, double z, double &dist2);
133 
140  double radius, const double x[3], double& dist2) override;
141 
143 
148  vtkIdType FindClosestPointInRegion(int regionId, double *x, double &dist2);
149  vtkIdType FindClosestPointInRegion(int regionId, double x, double y,
150  double z, double &dist2);
152 
158  double radius, const double x[3], vtkIdList *result) override;
159 
168  void FindClosestNPoints(int N, const double x[3],
169  vtkIdList *result) override;
170 
174  vtkIdTypeArray *GetPointsInRegion(int leafNodeId);
175 
179  void FreeSearchStructure() override;
180 
185  void GenerateRepresentation(int level, vtkPolyData *pd) override;
186 
193  void FindPointsInArea(double* area, vtkIdTypeArray* ids, bool clearArray = true);
194 
195 protected:
196 
198  ~vtkOctreePointLocator() override;
199 
201  vtkOctreePointLocatorNode **LeafNodeList; // indexed by region/node ID
202 
203  void BuildLeafNodeList(vtkOctreePointLocatorNode* node, int & index);
204 
206 
210  int FindRegion(vtkOctreePointLocatorNode* node, float x, float y, float z);
211  int FindRegion(vtkOctreePointLocatorNode* node, double x, double y, double z);
213 
214  static void SetDataBoundsToSpatialBounds(vtkOctreePointLocatorNode *node);
215 
216  static void DeleteAllDescendants(vtkOctreePointLocatorNode* octant);
217 
223  void FindPointsWithinRadius(vtkOctreePointLocatorNode* node, double radiusSquared,
224  const double x[3], vtkIdList* ids);
225 
226  // Recursive helper for public FindPointsWithinRadius
227  void AddAllPointsInRegion(vtkOctreePointLocatorNode* node, vtkIdList* ids);
228 
229  // Recursive helper for public FindPointsInArea
230  void FindPointsInArea(vtkOctreePointLocatorNode* node, double* area, vtkIdTypeArray* ids);
231 
232  // Recursive helper for public FindPointsInArea
233  void AddAllPointsInRegion(vtkOctreePointLocatorNode* node, vtkIdTypeArray* ids);
234 
235  void DivideRegion(vtkOctreePointLocatorNode *node, int* ordering, int level);
236 
237  int DivideTest(int size, int level);
238 
239  void AddPolys(vtkOctreePointLocatorNode *node, vtkPoints *pts, vtkCellArray *polys);
240 
245  int _FindClosestPointInRegion(int leafNodeId, double x, double y,
246  double z, double &dist2);
247 
255  int FindClosestPointInSphere(double x, double y, double z, double radius,
256  int skipRegion, double &dist2);
257 
259 
265 
266  double FudgeFactor; // a very small distance, relative to the dataset's size
270 
271  float MaxWidth;
272 
281 
283  void operator=(const vtkOctreePointLocator&) = delete;
284 };
285 #endif
virtual void BuildLocator()=0
Build the locator from the input dataset.
an octree spatial decomposition of a set of points
virtual vtkIdType FindClosestPointWithinRadius(double radius, const double x[3], double &dist2)=0
Given a position x and a radius r, return the id of the point closest to the point in that radius...
int CreateCubicOctants
If CreateCubicOctants is non-zero, the bounding box of the points will be expanded such that all octa...
virtual double * GetBounds()
Provide an accessor to the bounds.
int MaximumPointsPerRegion
The maximum number of points in a region/octant before it is subdivided.
dynamic, self-adjusting array of vtkIdType
int vtkIdType
Definition: vtkType.h:347
concrete dataset represents vertices, lines, polygons, and triangle strips
Definition: vtkPolyData.h:85
virtual void FreeSearchStructure()=0
Free the memory required for the spatial data structure.
int NumberOfLeafNodes
The maximum number of points in a region/octant before it is subdivided.
vtkOctreePointLocatorNode ** LeafNodeList
a simple class to control print indentation
Definition: vtkIndent.h:39
void PrintSelf(ostream &os, vtkIndent indent) override
Standard type and print methods.
abstract class to quickly locate points in 3-space
list of point or cell ids
Definition: vtkIdList.h:36
virtual void FindPointsWithinRadius(double R, const double x[3], vtkIdList *result)=0
Find all points within a specified radius R of position x.
virtual void FindClosestNPoints(int N, const double x[3], vtkIdList *result)=0
Find the closest N points to a position.
object to represent cell connectivity
Definition: vtkCellArray.h:50
Octree node that has 8 children each of equal size.
vtkOctreePointLocatorNode * Top
virtual vtkIdType FindClosestPoint(const double x[3])=0
Given a position x, return the id of the point closest to it.
static vtkObject * New()
Create an object with Debug turned off, modified time initialized to zero, and reference counting on...
virtual void GenerateRepresentation(int level, vtkPolyData *pd)=0
Method to build a representation at a particular level.
represent and manipulate 3D points
Definition: vtkPoints.h:39