VTK  9.2.6
vtkGreedyTerrainDecimation.h
Go to the documentation of this file.
1 /*=========================================================================
2 
3  Program: Visualization Toolkit
4  Module: vtkGreedyTerrainDecimation.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 =========================================================================*/
89 #ifndef vtkGreedyTerrainDecimation_h
90 #define vtkGreedyTerrainDecimation_h
91 
92 #include "vtkFiltersHybridModule.h" // For export macro
93 #include "vtkPolyDataAlgorithm.h"
94 
95 class vtkPriorityQueue;
96 class vtkDataArray;
97 class vtkPointData;
98 class vtkIdList;
99 class vtkDoubleArray;
100 class vtkFloatArray;
101 
102 // PIMPL Encapsulation for STL containers
103 class vtkGreedyTerrainDecimationTerrainInfoType;
104 class vtkGreedyTerrainDecimationPointInfoType;
105 
106 #define VTK_ERROR_NUMBER_OF_TRIANGLES 0
107 #define VTK_ERROR_SPECIFIED_REDUCTION 1
108 #define VTK_ERROR_ABSOLUTE 2
109 #define VTK_ERROR_RELATIVE 3
110 
111 class VTKFILTERSHYBRID_EXPORT vtkGreedyTerrainDecimation : public vtkPolyDataAlgorithm
112 {
113 public:
115  void PrintSelf(ostream& os, vtkIndent indent) override;
116 
121 
123 
129  vtkSetClampMacro(ErrorMeasure, int, VTK_ERROR_NUMBER_OF_TRIANGLES, VTK_ERROR_RELATIVE);
130  vtkGetMacro(ErrorMeasure, int);
132  {
133  this->SetErrorMeasure(VTK_ERROR_NUMBER_OF_TRIANGLES);
134  }
136  {
137  this->SetErrorMeasure(VTK_ERROR_SPECIFIED_REDUCTION);
138  }
139  void SetErrorMeasureToAbsoluteError() { this->SetErrorMeasure(VTK_ERROR_ABSOLUTE); }
140  void SetErrorMeasureToRelativeError() { this->SetErrorMeasure(VTK_ERROR_RELATIVE); }
142 
144 
150  vtkSetClampMacro(NumberOfTriangles, vtkIdType, 2, VTK_ID_MAX);
151  vtkGetMacro(NumberOfTriangles, vtkIdType);
153 
155 
160  vtkSetClampMacro(Reduction, double, 0.0, 1.0);
161  vtkGetMacro(Reduction, double);
163 
165 
170  vtkSetClampMacro(AbsoluteError, double, 0.0, VTK_DOUBLE_MAX);
171  vtkGetMacro(AbsoluteError, double);
173 
175 
181  vtkSetClampMacro(RelativeError, double, 0.0, VTK_DOUBLE_MAX);
182  vtkGetMacro(RelativeError, double);
184 
186 
190  vtkSetMacro(BoundaryVertexDeletion, vtkTypeBool);
191  vtkGetMacro(BoundaryVertexDeletion, vtkTypeBool);
192  vtkBooleanMacro(BoundaryVertexDeletion, vtkTypeBool);
194 
196 
199  vtkSetMacro(ComputeNormals, vtkTypeBool);
200  vtkGetMacro(ComputeNormals, vtkTypeBool);
201  vtkBooleanMacro(ComputeNormals, vtkTypeBool);
203 
204 protected:
207 
210 
213  void ComputePointNormal(int i, int j, float n[3]);
214 
215  // ivars that the API addresses
218  double Reduction;
221  vtkTypeBool BoundaryVertexDeletion; // Can we delete boundary vertices?
222 
223  // Used for convenience
230  double Tolerance;
232  int Dimensions[3];
233  double Origin[3];
234  double Spacing[3];
236  double Length;
237 
238  // Bookkeeping arrays
239  vtkPriorityQueue* TerrainError; // errors for each pt in height field
240  vtkGreedyTerrainDecimationTerrainInfoType* TerrainInfo; // owning triangle for each pt
241  vtkGreedyTerrainDecimationPointInfoType* PointInfo; // map mesh pt id to input pt id
242 
243  // Make a guess at initial allocation
244  void EstimateOutputSize(const vtkIdType numInputPts, vtkIdType& numPts, vtkIdType& numTris);
245 
246  // Returns non-zero if the error measure is satisfied.
247  virtual int SatisfiesErrorMeasure(double error);
248 
249  // Insert all the boundary vertices into the TIN
251 
252  // Insert a point into the triangulation; get a point from the triangulation
254  vtkIdType InsertNextPoint(vtkIdType inputPtId, double x[3]);
255  double* GetPoint(vtkIdType id);
256  void GetPoint(vtkIdType id, double x[3]);
257 
258  // Helper functions
259  void GetTerrainPoint(int i, int j, double x[3]);
260  void ComputeImageCoordinates(vtkIdType inputPtId, int ij[2]);
261  int InCircle(double x[3], double x1[3], double x2[3], double x3[3]);
262  vtkIdType FindTriangle(double x[3], vtkIdType ptIds[3], vtkIdType tri, double tol,
263  vtkIdType nei[3], vtkIdList* neighbors, int& status);
264  void CheckEdge(vtkIdType ptId, double x[3], vtkIdType p1, vtkIdType p2, vtkIdType tri, int depth);
265 
266  void UpdateTriangles(vtkIdType meshPtId); // update all points connected to this point
268  void UpdateTriangle(vtkIdType triId, int ij1[2], int ij2[2], int ij3[2], double h[4]);
269 
270  int CharacterizeTriangle(int ij1[2], int ij2[2], int ij[3], int*& min, int*& max, int*& midL,
271  int*& midR, int*& mid, int mid2[2], double h[3], double& hMin, double& hMax, double& hL,
272  double& hR);
273 
274 private:
276  void operator=(const vtkGreedyTerrainDecimation&) = delete;
277 };
278 
279 #endif
abstract superclass for arrays of numeric data
Definition: vtkDataArray.h:165
dynamic, self-adjusting array of double
dynamic, self-adjusting array of float
reduce height field (represented as image) to reduced TIN
int FillInputPortInformation(int port, vtkInformation *info) override
Fill the input port information objects for this algorithm.
int RequestData(vtkInformation *, vtkInformationVector **, vtkInformationVector *) override
This is called by the superclass.
void ComputeImageCoordinates(vtkIdType inputPtId, int ij[2])
void UpdateTriangle(vtkIdType triId, vtkIdType p1, vtkIdType p2, vtkIdType p3)
void GetTerrainPoint(int i, int j, double x[3])
int CharacterizeTriangle(int ij1[2], int ij2[2], int ij[3], int *&min, int *&max, int *&midL, int *&midR, int *&mid, int mid2[2], double h[3], double &hMin, double &hMax, double &hL, double &hR)
void SetErrorMeasureToNumberOfTriangles()
Specify how to terminate the algorithm: either as an absolute number of triangles,...
vtkGreedyTerrainDecimationTerrainInfoType * TerrainInfo
void ComputePointNormal(int i, int j, float n[3])
void SetErrorMeasureToAbsoluteError()
Specify how to terminate the algorithm: either as an absolute number of triangles,...
~vtkGreedyTerrainDecimation() override
vtkIdType InsertNextPoint(vtkIdType inputPtId, double x[3])
vtkIdType AddPointToTriangulation(vtkIdType inputPtId)
vtkIdType FindTriangle(double x[3], vtkIdType ptIds[3], vtkIdType tri, double tol, vtkIdType nei[3], vtkIdList *neighbors, int &status)
void PrintSelf(ostream &os, vtkIndent indent) override
Methods invoked by print to print information about the object including superclasses.
void EstimateOutputSize(const vtkIdType numInputPts, vtkIdType &numPts, vtkIdType &numTris)
vtkGreedyTerrainDecimationPointInfoType * PointInfo
void SetErrorMeasureToRelativeError()
Specify how to terminate the algorithm: either as an absolute number of triangles,...
void GetPoint(vtkIdType id, double x[3])
void UpdateTriangle(vtkIdType triId, int ij1[2], int ij2[2], int ij3[2], double h[4])
void SetErrorMeasureToSpecifiedReduction()
Specify how to terminate the algorithm: either as an absolute number of triangles,...
static vtkGreedyTerrainDecimation * New()
Instantiate the class.
virtual int SatisfiesErrorMeasure(double error)
int InCircle(double x[3], double x1[3], double x2[3], double x3[3])
void CheckEdge(vtkIdType ptId, double x[3], vtkIdType p1, vtkIdType p2, vtkIdType tri, int depth)
void UpdateTriangles(vtkIdType meshPtId)
double * GetPoint(vtkIdType id)
list of point or cell ids
Definition: vtkIdList.h:143
a simple class to control print indentation
Definition: vtkIndent.h:119
Store zero or more vtkInformation instances.
Store vtkAlgorithm input/output information.
represent and manipulate point attribute data
Definition: vtkPointData.h:151
Superclass for algorithms that produce only polydata as output.
concrete dataset represents vertices, lines, polygons, and triangle strips
Definition: vtkPolyData.h:200
a list of ids arranged in priority order
Computes the portion of a dataset which is inside a selection.
@ info
Definition: vtkX3D.h:382
@ port
Definition: vtkX3D.h:453
int vtkTypeBool
Definition: vtkABI.h:69
#define VTK_ERROR_SPECIFIED_REDUCTION
#define VTK_ERROR_ABSOLUTE
#define VTK_ERROR_RELATIVE
#define VTK_ERROR_NUMBER_OF_TRIANGLES
int vtkIdType
Definition: vtkType.h:332
#define VTK_ID_MAX
Definition: vtkType.h:336
#define VTK_DOUBLE_MAX
Definition: vtkType.h:165
#define max(a, b)