1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
|
/*=========================================================================
Program: Visualization Toolkit
Module: vtkSelectEnclosedPoints.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
/**
* @class vtkSelectEnclosedPoints
* @brief mark points as to whether they are inside a closed surface
*
* vtkSelectEnclosedPoints is a filter that evaluates all the input points to
* determine whether they are in an enclosed surface. The filter produces a
* (0,1) mask (in the form of a vtkDataArray) that indicates whether points
* are outside (mask value=0) or inside (mask value=1) a provided surface.
* (The name of the output vtkDataArray is "SelectedPoints".)
*
* After running the filter, it is possible to query it as to whether a point
* is inside/outside by invoking the IsInside(ptId) method.
*
* @warning
* The filter assumes that the surface is closed and manifold. A boolean flag
* can be set to force the filter to first check whether this is true. If false,
* all points will be marked outside. Note that if this check is not performed
* and the surface is not closed, the results are undefined.
*
* @warning
* This filter produces and output data array, but does not modify the input
* dataset. If you wish to extract cells or points, various threshold filters
* are available (i.e., threshold the output array). Also, see the filter
* vtkExtractEnclosedPoints which operates on point clouds.
*
* @warning
* This class has been threaded with vtkSMPTools. Using TBB or other
* non-sequential type (set in the CMake variable
* VTK_SMP_IMPLEMENTATION_TYPE) may improve performance significantly.
*
* @sa
* vtkMaskPoints vtkExtractEnclosedPoints
*/
#ifndef vtkSelectEnclosedPoints_h
#define vtkSelectEnclosedPoints_h
#include "vtkDataSetAlgorithm.h"
#include "vtkFiltersModelingModule.h" // For export macro
#include "vtkIntersectionCounter.h" // to count intersections along ray
class vtkUnsignedCharArray;
class vtkAbstractCellLocator;
class vtkStaticCellLocator;
class vtkIdList;
class vtkGenericCell;
class vtkRandomPool;
class VTKFILTERSMODELING_EXPORT vtkSelectEnclosedPoints : public vtkDataSetAlgorithm
{
public:
//@{
/**
* Standard methods for type information and printing.
*/
vtkTypeMacro(vtkSelectEnclosedPoints, vtkDataSetAlgorithm);
void PrintSelf(ostream& os, vtkIndent indent) override;
//@}
/**
* Instantiate this class.
*/
static vtkSelectEnclosedPoints* New();
//@{
/**
* Set the surface to be used to test for containment. Two methods are
* provided: one directly for vtkPolyData, and one for the output of a
* filter.
*/
void SetSurfaceData(vtkPolyData* pd);
void SetSurfaceConnection(vtkAlgorithmOutput* algOutput);
//@}
//@{
/**
* Return a pointer to the enclosing surface.
*/
vtkPolyData* GetSurface();
vtkPolyData* GetSurface(vtkInformationVector* sourceInfo);
//@}
//@{
/**
* By default, points inside the surface are marked inside or sent to
* the output. If InsideOut is on, then the points outside the surface
* are marked inside.
*/
vtkSetMacro(InsideOut, vtkTypeBool);
vtkBooleanMacro(InsideOut, vtkTypeBool);
vtkGetMacro(InsideOut, vtkTypeBool);
//@}
//@{
/**
* Specify whether to check the surface for closure. If on, then the
* algorithm first checks to see if the surface is closed and manifold.
*/
vtkSetMacro(CheckSurface, vtkTypeBool);
vtkBooleanMacro(CheckSurface, vtkTypeBool);
vtkGetMacro(CheckSurface, vtkTypeBool);
//@}
/**
* Query an input point id as to whether it is inside or outside. Note that
* the result requires that the filter execute first.
*/
int IsInside(vtkIdType inputPtId);
//@{
/**
* Specify the tolerance on the intersection. The tolerance is expressed as
* a fraction of the diagonal of the bounding box of the enclosing surface.
*/
vtkSetClampMacro(Tolerance, double, 0.0, VTK_FLOAT_MAX);
vtkGetMacro(Tolerance, double);
//@}
//@{
/**
* This is a backdoor that can be used to test many points for containment.
* First initialize the instance, then repeated calls to IsInsideSurface()
* can be used without rebuilding the search structures. The Complete()
* method releases memory.
*/
void Initialize(vtkPolyData* surface);
int IsInsideSurface(double x[3]);
int IsInsideSurface(double x, double y, double z);
void Complete();
//@}
/**
* A static method for determining whether a point is inside a
* surface. This is the heart of the algorithm and is thread safe. The user
* must provide an input point x, the enclosing surface, the bounds of the
* enclosing surface, the diagonal length of the enclosing surface, an
* intersection tolerance, a cell locator for the surface, and two working
* objects (cellIds, genCell) to support computation. Finally, in threaded
* execution, generating random numbers is hard, so a precomputed random
* sequence can be provided with an index into the sequence.
*/
static int IsInsideSurface(double x[3], vtkPolyData* surface, double bds[6], double length,
double tol, vtkAbstractCellLocator* locator, vtkIdList* cellIds, vtkGenericCell* genCell,
vtkIntersectionCounter& counter, vtkRandomPool* poole = nullptr, vtkIdType seqIdx = 0);
/**
* A static method for determining whether a surface is closed. Provide as input
* a vtkPolyData. The method returns >0 is the surface is closed and manifold.
*/
static int IsSurfaceClosed(vtkPolyData* surface);
protected:
vtkSelectEnclosedPoints();
~vtkSelectEnclosedPoints() override;
vtkTypeBool CheckSurface;
vtkTypeBool InsideOut;
double Tolerance;
vtkUnsignedCharArray* InsideOutsideArray;
// Internal structures for accelerating the intersection test
vtkStaticCellLocator* CellLocator;
vtkIdList* CellIds;
vtkGenericCell* Cell;
vtkPolyData* Surface;
double Bounds[6];
double Length;
int RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*) override;
int FillInputPortInformation(int, vtkInformation*) override;
void ReportReferences(vtkGarbageCollector*) override;
private:
vtkSelectEnclosedPoints(const vtkSelectEnclosedPoints&) = delete;
void operator=(const vtkSelectEnclosedPoints&) = delete;
};
#endif
|
