vtkClipHyperOctree.h 7.54 KB
Newer Older
Francois Bertel's avatar
Francois Bertel committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
/*=========================================================================

  Program:   Visualization Toolkit
  Module:    vtkClipHyperOctree.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.

=========================================================================*/
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
/**
 * @class   vtkClipHyperOctree
 * @brief   clip an hyperoctree with user-specified implicit function or input scalar data
 *
 * vtkClipHyperOctree is a filter that clips an hyperoctree using either
 * any subclass of vtkImplicitFunction, or the input scalar
 * data. Clipping means that it actually "cuts" through the leaves (cells) of
 * the hyperoctree, returning everything inside of the specified implicit
 * function (or greater than the scalar value) including "pieces" of
 * a cell. (Compare this with vtkExtractGeometry, which pulls out
 * entire, uncut cells.) The output of this filter is an unstructured
 * grid.
 *
 * To use this filter, you must decide if you will be clipping with an
 * implicit function, or whether you will be using the input scalar
 * data.  If you want to clip with an implicit function, you must:
 * 1) define an implicit function
 * 2) set it with the SetClipFunction method
 * 3) apply the GenerateClipScalarsOn method
 * If a ClipFunction is not specified, or GenerateClipScalars is off
 * (the default), then the input's scalar data will be used to clip
 * the polydata.
 *
 * You can also specify a scalar value, which is used to decide what is
 * inside and outside of the implicit function. You can also reverse the
 * sense of what inside/outside is by setting the InsideOut instance
 * variable. (The clipping algorithm proceeds by computing an implicit
 * function value or using the input scalar data for each point in the
 * dataset. This is compared to the scalar value to determine
 * inside/outside.)
 *
 * This filter can be configured to compute a second output. The
 * second output is the part of the cell that is clipped away. Set the
 * GenerateClippedData boolean on if you wish to access this output data.
 *
 * @warning
 * vtkClipHyperOctree will triangulate all types of 3D cells (i.e., create
 * tetrahedra). This is true even if the cell is not actually cut. This
 * is necessary to preserve compatibility across face neighbors. 2D cells
 * will only be triangulated if the cutting function passes through them.
 *
 * @sa
 * vtkImplicitFunction vtkCutter vtkClipVolume vtkClipPolyData
*/
Francois Bertel's avatar
Francois Bertel committed
59

60
61
#ifndef vtkClipHyperOctree_h
#define vtkClipHyperOctree_h
Francois Bertel's avatar
Francois Bertel committed
62

63
#include "vtkFiltersHyperTreeModule.h" // For export macro
Francois Bertel's avatar
Francois Bertel committed
64
65
#include "vtkUnstructuredGridAlgorithm.h"

66
#if !defined(VTK_LEGACY_REMOVE)
Francois Bertel's avatar
Francois Bertel committed
67
class vtkImplicitFunction;
68
class vtkIncrementalPointLocator;
Francois Bertel's avatar
Francois Bertel committed
69
70
71
72
73
74
75
76
77
78
79
80
81
82
class vtkHyperOctreeCursor;
class vtkHyperOctree;
class vtkUnsignedCharArray;
class vtkIdTypeArray;
class vtkCellArray;
class vtkCellData;
class vtkPointData;
class vtkOrderedTriangulator;
class vtkDoubleArray;
class vtkTetra;
class vtkPoints;
class vtkPolygon;
class vtkHyperOctreeClipCutPointsGrabber;

83
class VTKFILTERSHYPERTREE_EXPORT vtkClipHyperOctree : public vtkUnstructuredGridAlgorithm
Francois Bertel's avatar
Francois Bertel committed
84
85
{
public:
86
  vtkTypeMacro(vtkClipHyperOctree,vtkUnstructuredGridAlgorithm);
87
  void PrintSelf(ostream& os, vtkIndent indent) override;
Francois Bertel's avatar
Francois Bertel committed
88

89
90
91
92
  /**
   * Construct with user-specified implicit function; InsideOut turned off;
   * value set to 0.0; and generate clip scalars turned off.
   */
Francois Bertel's avatar
Francois Bertel committed
93
94
  static vtkClipHyperOctree *New();

95
96
97
98
99
100
  //@{
  /**
   * Set the clipping value of the implicit function (if clipping with
   * implicit function) or scalar value (if clipping with
   * scalars). The default value is 0.0.
   */
Francois Bertel's avatar
Francois Bertel committed
101
102
  vtkSetMacro(Value,double);
  vtkGetMacro(Value,double);
103
104
105
106
107
108
109
110
111
112
113
  //@}

  //@{
  /**
   * Set/Get the InsideOut flag. When off, a vertex is considered
   * inside the implicit function if its value is greater than the
   * Value ivar. When InsideOutside is turned on, a vertex is
   * considered inside the implicit function if its implicit function
   * value is less than or equal to the Value ivar.  InsideOut is off
   * by default.
   */
Francois Bertel's avatar
Francois Bertel committed
114
115
116
  vtkSetMacro(InsideOut,int);
  vtkGetMacro(InsideOut,int);
  vtkBooleanMacro(InsideOut,int);
117
118
119
120
121
122
123
124
  //@}

  //@{
  /**
   * Specify the implicit function with which to perform the
   * clipping. If you do not define an implicit function,
   * then the selected input scalar data will be used for clipping.
   */
Francois Bertel's avatar
Francois Bertel committed
125
126
  virtual void SetClipFunction(vtkImplicitFunction*);
  vtkGetObjectMacro(ClipFunction,vtkImplicitFunction);
127
128
129
130
131
132
133
134
135
  //@}

  //@{
  /**
   * If this flag is enabled, then the output scalar values will be
   * interpolated from the implicit function values, and not the
   * input scalar data. If you enable this flag but do not provide an
   * implicit function an error will be reported.
   */
Francois Bertel's avatar
Francois Bertel committed
136
137
138
  vtkSetMacro(GenerateClipScalars,int);
  vtkGetMacro(GenerateClipScalars,int);
  vtkBooleanMacro(GenerateClipScalars,int);
139
  //@}
Francois Bertel's avatar
Francois Bertel committed
140

141
142
143
144
145
  //@{
  /**
   * Control whether a second output is generated. The second output
   * contains the polygonal data that's been clipped away.
   */
Francois Bertel's avatar
Francois Bertel committed
146
147
148
  vtkSetMacro(GenerateClippedOutput,int);
  vtkGetMacro(GenerateClippedOutput,int);
  vtkBooleanMacro(GenerateClippedOutput,int);
149
  //@}
150

151
152
153
  /**
   * Return the Clipped output.
   */
Francois Bertel's avatar
Francois Bertel committed
154
155
  vtkUnstructuredGrid *GetClippedOutput();

156
157
158
159
160
  //@{
  /**
   * Specify a spatial locator for merging points. By default, an
   * instance of vtkMergePoints is used.
   */
161
162
  void SetLocator(vtkIncrementalPointLocator *locator);
  vtkGetObjectMacro(Locator,vtkIncrementalPointLocator);
163
  //@}
Francois Bertel's avatar
Francois Bertel committed
164

165
166
167
168
  /**
   * Create default locator. Used to create one when none is specified. The
   * locator is used to merge coincident points.
   */
Francois Bertel's avatar
Francois Bertel committed
169
170
  void CreateDefaultLocator();

171
172
173
  /**
   * Return the mtime also considering the locator and clip function.
   */
174
  vtkMTimeType GetMTime() override;
Francois Bertel's avatar
Francois Bertel committed
175
176

protected:
177
  vtkClipHyperOctree(vtkImplicitFunction *cf=nullptr);
178
  ~vtkClipHyperOctree() override;
Francois Bertel's avatar
Francois Bertel committed
179

180
  int RequestData(vtkInformation *, vtkInformationVector **, vtkInformationVector *) override;
181

182
183
184
185
186
  /**
   * Clip the sub-hierarchy pointed by cursor.
   * \pre cursor_exists: cursor!=0
   * \pre positive_level: level>=0
   */
Francois Bertel's avatar
Francois Bertel committed
187
188
189
  void ClipNode(vtkHyperOctreeCursor *cursor,
                int level,
                double bounds[6]);
190

191
  int FillInputPortInformation(int port, vtkInformation *info) override;
Francois Bertel's avatar
Francois Bertel committed
192
  vtkImplicitFunction *ClipFunction;
193

194
195
  vtkIncrementalPointLocator *Locator;
  vtkIncrementalPointLocator *Locator2; // used for the clipped output
196

Francois Bertel's avatar
Francois Bertel committed
197
198
199
200
201
  int InsideOut;
  double Value;
  int GenerateClipScalars;

  int GenerateClippedOutput;
202

Francois Bertel's avatar
Francois Bertel committed
203
204
205
  vtkHyperOctree *Input;
  vtkUnstructuredGrid *Output;
  vtkUnstructuredGrid *ClippedOutput;
206

Francois Bertel's avatar
Francois Bertel committed
207
208
209
210
211
212
213
  vtkUnsignedCharArray *Types[2];
  vtkIdTypeArray *Locs[2];
  vtkCellArray *Conn[2];
  vtkCellData *InCD;
  vtkCellData *OutCD[2];
  vtkPointData *OutPD[2];
  vtkOrderedTriangulator *Triangulator;
214

Francois Bertel's avatar
Francois Bertel committed
215
  vtkHyperOctreeCursor *Sibling; // to avoid allocation in the loop
216

Francois Bertel's avatar
Francois Bertel committed
217
218
219
  vtkDoubleArray *CellScalars;
  vtkTetra *Tetra;
  vtkDoubleArray *TetScalars;
220

Francois Bertel's avatar
Francois Bertel committed
221
222
  vtkPoints *Pts;
  vtkPolygon *Polygon;
223

Francois Bertel's avatar
Francois Bertel committed
224
225
226
227
  vtkIdType CellTypeCounter[65536]; // up-to-65536 points per octant
  vtkIdType TotalCounter;
  vtkIdType TemplateCounter; // record the number of octants that succceed
  // to use the template triangulator
228

Francois Bertel's avatar
Francois Bertel committed
229
  vtkHyperOctreeClipCutPointsGrabber *Grabber;
230

Francois Bertel's avatar
Francois Bertel committed
231
private:
232
233
  vtkClipHyperOctree(const vtkClipHyperOctree&) = delete;
  void operator=(const vtkClipHyperOctree&) = delete;
Francois Bertel's avatar
Francois Bertel committed
234
};
235
#endif // LEGACY remove
Francois Bertel's avatar
Francois Bertel committed
236
237

#endif