Fade2.5D Documentation pages v2.07
Delaunay Features
Fade_2D.h
Go to the documentation of this file.
1 // Copyright (C) Geom Software e.U, Bernhard Kornberger, Graz/Austria
2 //
3 // This file is part of the Fade2D library. The student license is free
4 // of charge and covers personal non-commercial research. Licensees
5 // holding a commercial license may use this file in accordance with
6 // the Commercial License Agreement.
7 //
8 // This software is provided AS IS with NO WARRANTY OF ANY KIND,
9 // INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS
10 // FOR A PARTICULAR PURPOSE.
11 //
12 // Please contact the author if any conditions of this licensing are
13 // not clear to you.
14 //
15 // Author: Bernhard Kornberger, bkorn (at) geom.at
16 // http://www.geom.at
17 
18 
24 #pragma once
25 
26 #include "common.h"
27 #include "freeFunctions.h"
28 #ifndef FADE2D_EXPORT
29  #include "License.h"
30 #endif
31 #include "Point2.h"
32 #include "Triangle2.h"
34 #include "Visualizer2.h"
35 #include "Zone2.h"
36 #include "ConstraintGraph2.h"
37 #include "Performance.h"
38 #include "MeshGenParams.h"
39 #include "MsgBase.h"
40 #include "SegmentChecker.h"
41 #include "testDataGenerators.h"
42 #include "FadeExport.h"
43 #include "Voronoi2.h"
44 
45 #if GEOM_PSEUDO3D==GEOM_TRUE
46  #include "IsoContours.h"
47  #include "EfficientModel.h"
48  #include "CutAndFill.h"
49  #include "CloudPrepare.h"
50 #endif
51 
52 
53 
54 #if GEOM_PSEUDO3D==GEOM_TRUE
55  namespace GEOM_FADE25D {
56 #elif GEOM_PSEUDO3D==GEOM_FALSE
57  namespace GEOM_FADE2D {
58 #else
59  #error GEOM_PSEUDO3D is not defined
60 #endif
61 
62 
63 
64 class Dt2; // Forward
65 class Visualizer3; // Forward
66 
67 
72 class CLASS_DECLSPEC Fade_2D
73 {
74 public:
84  explicit Fade_2D(unsigned numExpectedVertices=3)
85  {
86  initFade(numExpectedVertices);
87  }
90 
91 
92 
93 
111  bool saveTriangulation(const char* filename,std::vector<Zone2*>& vSaveZones);
112 
130  bool saveTriangulation(std::ostream& stream,std::vector<Zone2*>& vSaveZones);
131 
132 
155  bool saveZones(const char* filename,std::vector<Zone2*>& vSaveZones);
156 
157 
180  bool saveZones(std::ostream& stream,std::vector<Zone2*>& vSaveZones);
181 
182 
194  bool load(const char* filename,std::vector<Zone2*>& vZones);
195 
196 
208  bool load(std::istream& stream,std::vector<Zone2*>& vZones);
209 
210 
220  void exportTriangulation(FadeExport& fadeExport,bool bWithCustomIndices,bool bClear);
221 
222 
235  bool checkValidity(bool bCheckEmptyCircleProperty,const char* msg) const;
236 
243  int setNumCPU(int numCPU);
244 
250  Voronoi2* getVoronoiDiagram();
251 
266  void setFastMode(bool bFast);
267 
268 
273  void statistics(const char* s) const;
274 
285  void show(const char* postscriptFilename,bool bWithConstraints=true) const;
286 
287 
300  void show(Visualizer2* pVis,bool bWithConstraints=true) const;
301 
302 #if GEOM_PSEUDO3D==GEOM_TRUE
310  void showGeomview(const char* filename, const char* color="1 1 1 0.5") const;
311 
319  void showGeomview(Visualizer3* pVis, const char* color="1 1 1 0.5") const;
320 
321 
322 #endif
323 
333  void remove(Point2* pVertex);
334 
342  void remove(std::vector<Point2*>& vPoints);
343 
344 
350  void reset();
351 
352 
365  void getConvexHull(bool bAllVertices,std::vector<Point2*>& vConvexHullPointsOut);
366 
367 #if GEOM_PSEUDO3D==GEOM_TRUE
385  void insert(CloudPrepare* pCloudPrepare,bool bClear=true);
386 #endif
387 
401  Point2* insert(const Point2& p);
402 
411  void insert(const std::vector<Point2>& vInputPoints);
412 
427  void insert(const std::vector<Point2>& vInputPoints,std::vector<Point2*>& vHandles);
428 
429 #if GEOM_PSEUDO3D==GEOM_TRUE
438 #else
447 #endif
448  void insert(int numPoints,double * aCoordinates,Point2 ** aHandles);
449 
450 
467  double measureTriangulationTime(std::vector<Point2>& vPoints);
468 
469 
470 
485  Triangle2* locate(const Point2& p);
486 
487 
497 
498 #if GEOM_PSEUDO3D==GEOM_TRUE
524  bool getHeight(double x,double y,double& heightOut,Triangle2* pApproxT=NULL,double tolerance=0) const;
525 #endif
526 
542  void refine( Zone2* pZone,
543  double minAngleDegree,
544  double minEdgeLength,
545  double maxEdgeLength,
546  bool bAllowConstraintSplitting
547  );
548 
584  void refineAdvanced(MeshGenParams* pParameters);
585 
586 
587 
588 
589 
590 
600  size_t numberOfPoints() const;
601 
608  size_t numberOfTriangles() const;
609 
610 
622  void getTrianglePointers(std::vector<Triangle2*>& vAllTriangles) const;
623 
637  void getVertexPointers(std::vector<Point2*>& vAllPoints) const;
638 
639 
640 
651 
668  bool hasArea() const;
669 
675  bool is2D() const;
676 
677 
733 #if GEOM_PSEUDO3D==GEOM_TRUE
735  std::vector<Segment2>& vSegments,
737  bool bOrientedSegments=false,
738  bool bUseHeightOfLatest=false
739  );
740 #else
741  ConstraintGraph2* createConstraint(
742  std::vector<Segment2>& vSegments,
744  bool bOrientedSegments=false
745  );
746 #endif
747 
748 
749 
750 
751 
774  Zone2* createZone(ConstraintGraph2* pConstraintGraph,ZoneLocation zoneLoc,bool bVerbose=true);
775 
776 
788  Zone2* createZone(const std::vector<ConstraintGraph2*>& vConstraintGraphs,ZoneLocation zoneLoc,const Point2& startPoint,bool bVerbose=true);
789 
802  Zone2* createZone(ConstraintGraph2* pConstraintGraph,ZoneLocation zoneLoc,const Point2& startPoint,bool bVerbose=true);
803 
814  Zone2* createZone( std::vector<Triangle2*>& vTriangles,bool bVerbose=true );
815 
816 
817 
830  void deleteZone(Zone2* pZone);
831 
832 
853 
854 
862 
870  bool isConstraint(Triangle2* pT,int ith) const;
871 
874  void getAliveConstraintSegments(std::vector<ConstraintSegment2*>& vAliveConstraintSegments) const;
875 
878  void getAliveAndDeadConstraintSegments(std::vector<ConstraintSegment2*>& vAllConstraintSegments) const;
879 
888 
889 
894  void getIncidentTriangles(Point2* pVtx,std::vector<Triangle2*>& vIncidentT) const;
895 
900  void getIncidentVertices(Point2* pVtx,std::vector<Point2*>& vIncidentVertices) const;
901 
902 
910  void writeObj(const char* filename) const;
911 
918  void writeObj(const char* filename,Zone2* pZone) const;
919 
920 
926  void writeWebScene(const char* path) const;
927 
933  void writeWebScene(const char* path,Zone2* pZone) const;
934 
940 void subscribe(MsgType msgType,MsgBase* pMsg);
941 
947 void unsubscribe(MsgType msgType,MsgBase* pMsg);
948 
949 
956  bool isConstraint(Point2* p0,Point2* p1) const;
957 
963  bool isConstraint(Point2* pVtx) const;
964 
966  void printLicense() const;
967 
969  bool checkZoneQuality(Zone2* pZone,double minAngle,const char* name,const AcceptExperimentalFeature& accept);
971  void setName(const char* s);
973  std::string getName() const;
974 
975 
997  Zone2* importTriangles( std::vector<Point2>& vPoints,
998  bool bReorientIfNeeded,
999  bool bCreateExtendedBoundingBox
1000  );
1001 
1006  Orientation2 getOrientation(const Point2& p0,const Point2& p1,const Point2& p2);
1007 
1038  void cutTriangles( const Point2& knifeStart,
1039  const Point2& knifeEnd,
1040  bool bTurnEdgesIntoConstraints);
1041 
1054  void cutTriangles( std::vector<Segment2>& vSegments,
1055  bool bTurnEdgesIntoConstraints);
1056 
1057 
1077  Zone2* createZone_cookieCutter(std::vector<Segment2>& vSegments,bool bProtectEdges);
1078 
1079 
1125 #if GEOM_PSEUDO3D==GEOM_TRUE
1126  bool drape( std::vector<Segment2>& vSegmentsIn,
1127  std::vector<Segment2>& vSegmentsOut,
1128  double zTolerance) const;
1129 #else
1130  bool drape( std::vector<Segment2>& vSegmentsIn,
1131  std::vector<Segment2>& vSegmentsOut) const;
1132 #endif
1133 
1134 
1135 
1136 
1142 void getDirectedEdges(std::vector<Edge2>& vDirectedEdgesOut) const;
1143 
1149 void getUndirectedEdges(std::vector<Edge2>& vUndirectedEdgesOut) const;
1150 
1152  /* Deprecated: Use setNumCPU() instead. This method is kept for
1153  * compatibility with existing applications. Internally it calls
1154  * setNumCPU(0) to automatically determine and use the number of
1155  * available CPU cores.
1156  */
1157  void enableMultithreading();
1158  // Development functions, not for public use
1160  void internal(int au,int fu,const char* s="");
1162  Dt2* getImpl();
1164  void setDev(const char* s,int ival,double dval);
1165 
1166 protected:
1168  void initFade(unsigned numExpectedVertices);
1170  Fade_2D(const Fade_2D&); // No copy constructor
1172  Fade_2D& operator=(const Fade_2D&); // No assignment allowed
1174  Dt2* pImpl;
1175 }; // end of class
1176 
1177 
1185 inline void copyFade(Fade_2D& dt0,std::vector<Zone2*>& vZones0,Fade_2D& dt1,std::vector<Zone2*>& vZones1)
1186 {
1187  std::stringstream ss;
1188  dt0.saveTriangulation(ss,vZones0);
1189  dt1.load(ss,vZones1);
1190 }
1191 
1192 
1193 
1194 
1195 } // (namespace)
1196 
ConstraintInsertionStrategy
Constraint Insertion Strategy determines how a constraint edge shall be inserted:
Definition: ConstraintSegment2.h:52
void copyFade(Fade_2D &dt0, std::vector< Zone2 * > &vZones0, Fade_2D &dt1, std::vector< Zone2 * > &vZones1)
Copy a Fade_2D object and selected Zone2 objects.
Definition: Fade_2D.h:1185
Bbox2 is an axis aligned 2D bounding box.
Definition: Bbox2.h:37
CloudPrepare simplifies overdense point clouds and helps to avoid memory-usage-peaks during data tran...
Definition: CloudPrepare.h:87
ConstraintGraph2 is a set of Constraint Edges (ConstraintSegment2)
Definition: ConstraintGraph2.h:52
A ConstraintSegment2 represents a Constraint Edge.
Definition: ConstraintSegment2.h:69
Fade_2D is the Delaunay triangulation class.
Definition: Fade_2D.h:73
void getTrianglePointers(std::vector< Triangle2 * > &vAllTriangles) const
Get pointers to all triangles.
void writeWebScene(const char *path, Zone2 *pZone) const
Write a zone to an *.obj file.
bool getHeight(double x, double y, double &heightOut, Triangle2 *pApproxT=NULL, double tolerance=0) const
Compute the height of a certain point.
bool load(const char *filename, std::vector< Zone2 * > &vZones)
Load a triangulation.
void getVertexPointers(std::vector< Point2 * > &vAllPoints) const
Get pointers to all vertices.
void statistics(const char *s) const
Statistics.
void cutTriangles(const Point2 &knifeStart, const Point2 &knifeEnd, bool bTurnEdgesIntoConstraints)
Cut through a triangulation.
void getConvexHull(bool bAllVertices, std::vector< Point2 * > &vConvexHullPointsOut)
Compute the convex hull.
void getAliveConstraintSegments(std::vector< ConstraintSegment2 * > &vAliveConstraintSegments) const
Get active (alive) constraint segments.
bool drape(std::vector< Segment2 > &vSegmentsIn, std::vector< Segment2 > &vSegmentsOut, double zTolerance) const
Drape segments along a surface.
void insert(const std::vector< Point2 > &vInputPoints)
Insert a vector of points.
void showGeomview(Visualizer3 *pVis, const char *color="1 1 1 0.5") const
Draws the triangulation in 3D.
void writeWebScene(const char *path) const
Write the current triangulation to an *.obj file.
Voronoi2 * getVoronoiDiagram()
Get the Voronoi diagram.
void unsubscribe(MsgType msgType, MsgBase *pMsg)
Unregister a message receiver.
void setFastMode(bool bFast)
Set fast mode.
void remove(std::vector< Point2 * > &vPoints)
Remove vertices.
void insert(CloudPrepare *pCloudPrepare, bool bClear=true)
Insert points from a CloudPrepare object.
void insert(int numPoints, double *aCoordinates, Point2 **aHandles)
Insert points from an array.
void refineAdvanced(MeshGenParams *pParameters)
Delaunay refinement and grid meshing.
size_t numberOfTriangles() const
Number of triangles.
void writeObj(const char *filename, Zone2 *pZone) const
Write a zone to an *.obj file.
void insert(const std::vector< Point2 > &vInputPoints, std::vector< Point2 * > &vHandles)
Insert points from a std::vector and store pointers in vHandles.
ConstraintGraph2 * createConstraint(std::vector< Segment2 > &vSegments, ConstraintInsertionStrategy cis, bool bOrientedSegments=false, bool bUseHeightOfLatest=false)
Add constraint edges (edges, polyline, polygon) .
ConstraintSegment2 * getConstraintSegment(Point2 *p0, Point2 *p1) const
Retrieve a ConstraintSegment2.
bool saveZones(std::ostream &stream, std::vector< Zone2 * > &vSaveZones)
Save zones.
bool isConstraint(Point2 *p0, Point2 *p1) const
Check if an edge is a constraint edge.
Point2 * insert(const Point2 &p)
Insert a single point.
void applyConstraintsAndZones()
Apply conforming constraints and zones (deprecated!)
bool hasArea() const
Check if the triangulation contains triangles (which is the case if at least 3 non-collinear points e...
void refine(Zone2 *pZone, double minAngleDegree, double minEdgeLength, double maxEdgeLength, bool bAllowConstraintSplitting)
Delaunay refinement.
void remove(Point2 *pVertex)
Remove a single vertex.
Zone2 * createZone(ConstraintGraph2 *pConstraintGraph, ZoneLocation zoneLoc, bool bVerbose=true)
Create a zone.
void show(const char *postscriptFilename, bool bWithConstraints=true) const
Draws the triangulation as postscript file.
bool saveTriangulation(const char *filename, std::vector< Zone2 * > &vSaveZones)
Save a triangulation.
void showGeomview(const char *filename, const char *color="1 1 1 0.5") const
Draws the triangulation in 3D.
Point2 * getNearestNeighbor(const Point2 &p)
Get nearest neighbor.
void deleteZone(Zone2 *pZone)
Delete a Zone2 object.
size_t numberOfPoints() const
Number of points.
bool saveZones(const char *filename, std::vector< Zone2 * > &vSaveZones)
Save zones.
void getIncidentTriangles(Point2 *pVtx, std::vector< Triangle2 * > &vIncidentT) const
Get incident triangles.
Zone2 * createZone(const std::vector< ConstraintGraph2 * > &vConstraintGraphs, ZoneLocation zoneLoc, const Point2 &startPoint, bool bVerbose=true)
Create a zone limited by multiple ConstraintGraph2 objects by growing from a start point.
void show(Visualizer2 *pVis, bool bWithConstraints=true) const
Draws the triangulation as postscript file using an existing Visualizer2 object.
Zone2 * importTriangles(std::vector< Point2 > &vPoints, bool bReorientIfNeeded, bool bCreateExtendedBoundingBox)
Import triangles.
Triangle2 * getAdjacentTriangle(Point2 *p0, Point2 *p1) const
Get adjacent triangle.
Triangle2 * locate(const Point2 &p)
Locate a triangle which contains p.
bool isConstraint(Point2 *pVtx) const
Check if a vertex is a constraint vertex.
void exportTriangulation(FadeExport &fadeExport, bool bWithCustomIndices, bool bClear)
Export triangulation data from Fade.
Fade_2D(unsigned numExpectedVertices=3)
Constructor of the triangulation class.
Definition: Fade_2D.h:84
void getUndirectedEdges(std::vector< Edge2 > &vUndirectedEdgesOut) const
Get undirected edges.
Zone2 * createZone(std::vector< Triangle2 * > &vTriangles, bool bVerbose=true)
Create a zone defined by a vector of triangles.
void getDirectedEdges(std::vector< Edge2 > &vDirectedEdgesOut) const
Get directed edges Edges are counterclockwise oriented around their triangle. The present method retu...
bool checkValidity(bool bCheckEmptyCircleProperty, const char *msg) const
Checks if a triangulation is valid.
bool isConstraint(Triangle2 *pT, int ith) const
Check if an edge is a constraint edge.
void printLicense() const
Prints license information.
int setNumCPU(int numCPU)
Set the number CPU cores for multithreading (deprecated)
void subscribe(MsgType msgType, MsgBase *pMsg)
Register a message receiver.
void writeObj(const char *filename) const
Write the current triangulation to an *.obj file.
bool saveTriangulation(std::ostream &stream, std::vector< Zone2 * > &vSaveZones)
Save a triangulation.
void cutTriangles(std::vector< Segment2 > &vSegments, bool bTurnEdgesIntoConstraints)
Cut through a triangulation.
Zone2 * createZone(ConstraintGraph2 *pConstraintGraph, ZoneLocation zoneLoc, const Point2 &startPoint, bool bVerbose=true)
Create a zone limited by a ConstraintGraph by growing from a start point.
bool load(std::istream &stream, std::vector< Zone2 * > &vZones)
Load a triangulation.
Zone2 * createZone_cookieCutter(std::vector< Segment2 > &vSegments, bool bProtectEdges)
Cookie Cutter The Cookie Cutter cuts out a part of a triangulation and returns it as a Zone2 object.
Orientation2 getOrientation(const Point2 &p0, const Point2 &p1, const Point2 &p2)
Compute the orientation of 3 points.
void getIncidentVertices(Point2 *pVtx, std::vector< Point2 * > &vIncidentVertices) const
Get incident vertices.
~Fade_2D()
Destructor.
void getAliveAndDeadConstraintSegments(std::vector< ConstraintSegment2 * > &vAllConstraintSegments) const
Get all (alive and dead) constraint segments.
double measureTriangulationTime(std::vector< Point2 > &vPoints)
Measure the Delaunay triangulation time.
Bbox2 computeBoundingBox() const
Compute the axis-aligned bounding box of the points.
Parameters for the mesh generator.
Definition: MeshGenParams.h:59
MsgBase, a base class for message subscriber classes.
Definition: MsgBase.h:47
Point.
Definition: Point2.h:43
Triangle.
Definition: Triangle2.h:60
Visualizer2 is a general Postscript writer. It draws the objects Point2, Segment2,...
Definition: Visualizer2.h:56
Visualizer3 is a 3D scene writer for the Geomview viewer.
Definition: Visualizer3.h:37
Zone2 is a certain defined area of a triangulation.
Definition: Zone2.h:92
FadeExport is a simple struct to export triangulation data.
Definition: FadeExport.h:43