Fade2D Documentation pages v2.03
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:
81  explicit Fade_2D(unsigned numExpectedVertices=3)
82  {
83  initFade(numExpectedVertices);
84  }
87 
88 
89 
90 
108  bool saveTriangulation(const char* filename,std::vector<Zone2*>& vSaveZones);
109 
127  bool saveTriangulation(std::ostream& stream,std::vector<Zone2*>& vSaveZones);
128 
129 
152  bool saveZones(const char* filename,std::vector<Zone2*>& vSaveZones);
153 
154 
177  bool saveZones(std::ostream& stream,std::vector<Zone2*>& vSaveZones);
178 
179 
191  bool load(const char* filename,std::vector<Zone2*>& vZones);
192 
193 
205  bool load(std::istream& stream,std::vector<Zone2*>& vZones);
206 
207 
217  void exportTriangulation(FadeExport& fadeExport,bool bWithCustomIndices,bool bClear);
218 
219 
232  bool checkValidity(bool bCheckEmptyCircleProperty,const char* msg) const;
233 
240  int setNumCPU(int numCPU);
241 
248 
263  void setFastMode(bool bFast);
264 
265 
270  void statistics(const char* s) const;
271 
282  void show(const char* postscriptFilename,bool bWithConstraints=true) const;
283 
284 
297  void show(Visualizer2* pVis,bool bWithConstraints=true) const;
298 
299 #if GEOM_PSEUDO3D==GEOM_TRUE
307  void showGeomview(const char* filename, const char* color="1 1 1 0.5") const;
308 
316  void showGeomview(Visualizer3* pVis, const char* color="1 1 1 0.5") const;
317 
318 
319 #endif
320 
330  void remove(Point2* pVertex);
331 
332 
338  void reset();
339 
340 
353  void getConvexHull(bool bAllVertices,std::vector<Point2*>& vConvexHullPointsOut);
354 
355 #if GEOM_PSEUDO3D==GEOM_TRUE
373  void insert(CloudPrepare* pCloudPrepare,bool bClear=true);
374 #endif
375 
389  Point2* insert(const Point2& p);
390 
399  void insert(const std::vector<Point2>& vInputPoints);
400 
415  void insert(const std::vector<Point2>& vInputPoints,std::vector<Point2*>& vHandles);
416 
417 #if GEOM_PSEUDO3D==GEOM_TRUE
426 #else
435 #endif
436  void insert(int numPoints,double * aCoordinates,Point2 ** aHandles);
437 
438 
455  double measureTriangulationTime(std::vector<Point2>& vPoints);
456 
457 
458 
473  Triangle2* locate(const Point2& p);
474 
475 
485 
486 #if GEOM_PSEUDO3D==GEOM_TRUE
512  bool getHeight(double x,double y,double& heightOut,Triangle2* pApproxT=NULL,double tolerance=0) const;
513 #endif
514 
530  void refine( Zone2* pZone,
531  double minAngleDegree,
532  double minEdgeLength,
533  double maxEdgeLength,
534  bool bAllowConstraintSplitting
535  );
536 
572  void refineAdvanced(MeshGenParams* pParameters);
573 
574 
575 
576 
577 
578 
588  size_t numberOfPoints() const;
589 
596  size_t numberOfTriangles() const;
597 
598 
610  void getTrianglePointers(std::vector<Triangle2*>& vAllTriangles) const;
611 
625  void getVertexPointers(std::vector<Point2*>& vAllPoints) const;
626 
627 
628 
639 
656  bool hasArea() const;
657 
663  bool is2D() const;
664 
665 
721 #if GEOM_PSEUDO3D==GEOM_TRUE
722  ConstraintGraph2* createConstraint(
723  std::vector<Segment2>& vSegments,
725  bool bOrientedSegments=false,
726  bool bUseHeightOfLatest=false
727  );
728 #else
730  std::vector<Segment2>& vSegments,
732  bool bOrientedSegments=false
733  );
734 #endif
735 
736 
737 
738 
739 
762  Zone2* createZone(ConstraintGraph2* pConstraintGraph,ZoneLocation zoneLoc,bool bVerbose=true);
763 
764 
776  Zone2* createZone(const std::vector<ConstraintGraph2*>& vConstraintGraphs,ZoneLocation zoneLoc,const Point2& startPoint,bool bVerbose=true);
777 
790  Zone2* createZone(ConstraintGraph2* pConstraintGraph,ZoneLocation zoneLoc,const Point2& startPoint,bool bVerbose=true);
791 
802  Zone2* createZone( std::vector<Triangle2*>& vTriangles,bool bVerbose=true );
803 
804 
805 
818  void deleteZone(Zone2* pZone);
819 
820 
841 
842 
850 
858  bool isConstraint(Triangle2* pT,int ith) const;
859 
862  void getAliveConstraintSegments(std::vector<ConstraintSegment2*>& vAliveConstraintSegments) const;
863 
866  void getAliveAndDeadConstraintSegments(std::vector<ConstraintSegment2*>& vAllConstraintSegments) const;
867 
876 
877 
882  void getIncidentTriangles(Point2* pVtx,std::vector<Triangle2*>& vIncidentT) const;
883 
888  void getIncidentVertices(Point2* pVtx,std::vector<Point2*>& vIncidentVertices) const;
889 
890 
898  void writeObj(const char* filename) const;
899 
906  void writeObj(const char* filename,Zone2* pZone) const;
907 
908 
914  void writeWebScene(const char* path) const;
915 
921  void writeWebScene(const char* path,Zone2* pZone) const;
922 
928 void subscribe(MsgType msgType,MsgBase* pMsg);
929 
935 void unsubscribe(MsgType msgType,MsgBase* pMsg);
936 
937 
944  bool isConstraint(Point2* p0,Point2* p1) const;
945 
951  bool isConstraint(Point2* pVtx) const;
952 
954  void printLicense() const;
955 
957  bool checkZoneQuality(Zone2* pZone,double minAngle,const char* name,const AcceptExperimentalFeature& accept);
959  void setName(const char* s);
961  std::string getName() const;
962 
963 
985  Zone2* importTriangles( std::vector<Point2>& vPoints,
986  bool bReorientIfNeeded,
987  bool bCreateExtendedBoundingBox
988  );
989 
994  Orientation2 getOrientation(const Point2& p0,const Point2& p1,const Point2& p2);
995 
1026  void cutTriangles( const Point2& knifeStart,
1027  const Point2& knifeEnd,
1028  bool bTurnEdgesIntoConstraints);
1029 
1042  void cutTriangles( std::vector<Segment2>& vSegments,
1043  bool bTurnEdgesIntoConstraints);
1044 
1045 
1065  Zone2* createZone_cookieCutter(std::vector<Segment2>& vSegments,bool bProtectEdges);
1066 
1067 
1113 #if GEOM_PSEUDO3D==GEOM_TRUE
1114  bool drape( std::vector<Segment2>& vSegmentsIn,
1115  std::vector<Segment2>& vSegmentsOut,
1116  double zTolerance) const;
1117 #else
1118  bool drape( std::vector<Segment2>& vSegmentsIn,
1119  std::vector<Segment2>& vSegmentsOut) const;
1120 #endif
1121 
1122 
1123 
1124 
1130 void getDirectedEdges(std::vector<Edge2>& vDirectedEdgesOut) const;
1131 
1137 void getUndirectedEdges(std::vector<Edge2>& vUndirectedEdgesOut) const;
1138 
1139 
1140 
1142  /* Deprecated: Use setNumCPU() instead. This method is kept for
1143  * compatibility with existing applications. Internally it calls
1144  * setNumCPU(0) to automatically determine and use the number of
1145  * available CPU cores.
1146  */
1147  void enableMultithreading();
1148  // Development functions, not for public use
1150  void internal(int au,int fu,const char* s="");
1152  Dt2* getImpl();
1154  void setDev(const char* s,int ival,double dval);
1155 
1156 protected:
1158  void initFade(unsigned numExpectedVertices);
1160  Fade_2D(const Fade_2D&); // No copy constructor
1162  Fade_2D& operator=(const Fade_2D&); // No assignment allowed
1164  Dt2* pImpl;
1165 }; // end of class
1166 
1167 
1168 
1169 
1170 
1171 } // (namespace)
1172 
ConstraintInsertionStrategy
Constraint Insertion Strategy determines how a constraint edge shall be inserted:
Definition: ConstraintSegment2.h:52
Bbox2 is an axis aligned 2D bounding box.
Definition: Bbox2.h:37
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 main class.
Definition: Fade_2D.h:73
Triangle2 * getAdjacentTriangle(Point2 *p0, Point2 *p1) const
Get adjacent triangle.
void insert(const std::vector< Point2 > &vInputPoints, std::vector< Point2 * > &vHandles)
Insert points from a std::vector and store pointers in vHandles.
void unsubscribe(MsgType msgType, MsgBase *pMsg)
Unregister a message receiver.
int setNumCPU(int numCPU)
Set the number CPU cores for multithreading (deprecated)
void setFastMode(bool bFast)
Set fast mode.
bool drape(std::vector< Segment2 > &vSegmentsIn, std::vector< Segment2 > &vSegmentsOut) const
Drape segments along a surface.
void getConvexHull(bool bAllVertices, std::vector< Point2 * > &vConvexHullPointsOut)
Compute the convex hull.
void subscribe(MsgType msgType, MsgBase *pMsg)
Register a message receiver.
bool hasArea() const
Check if the triangulation contains triangles (which is the case if at least 3 non-collinear points e...
Orientation2 getOrientation(const Point2 &p0, const Point2 &p1, const Point2 &p2)
Compute the orientation of 3 points.
void show(Visualizer2 *pVis, bool bWithConstraints=true) const
Draws the triangulation as postscript file using an existing Visualizer2 object.
void getUndirectedEdges(std::vector< Edge2 > &vUndirectedEdgesOut) const
Get undirected edges.
ConstraintGraph2 * createConstraint(std::vector< Segment2 > &vSegments, ConstraintInsertionStrategy cis, bool bOrientedSegments=false)
Add constraint edges (edges, polyline, polygon) .
void insert(const std::vector< Point2 > &vInputPoints)
Insert a vector of points.
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.
void writeWebScene(const char *path, Zone2 *pZone) const
Write a zone to an *.obj file.
void writeObj(const char *filename, Zone2 *pZone) const
Write a zone to an *.obj file.
bool load(std::istream &stream, std::vector< Zone2 * > &vZones)
Load a triangulation.
void getVertexPointers(std::vector< Point2 * > &vAllPoints) const
Get pointers to all vertices.
void writeObj(const char *filename) const
Write the current triangulation to an *.obj file.
Bbox2 computeBoundingBox() const
Compute the axis-aligned bounding box of the points.
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 getTrianglePointers(std::vector< Triangle2 * > &vAllTriangles) const
Get pointers to all triangles.
Point2 * insert(const Point2 &p)
Insert a single point.
void printLicense() const
Prints license information.
Fade_2D(unsigned numExpectedVertices=3)
Constructor of the main triangulation class.
Definition: Fade_2D.h:81
ConstraintSegment2 * getConstraintSegment(Point2 *p0, Point2 *p1) const
Retrieve a ConstraintSegment2.
void remove(Point2 *pVertex)
Remove a single vertex.
void getIncidentTriangles(Point2 *pVtx, std::vector< Triangle2 * > &vIncidentT) const
Get incident triangles.
double measureTriangulationTime(std::vector< Point2 > &vPoints)
Measure the Delaunay triangulation time.
bool saveTriangulation(const char *filename, std::vector< Zone2 * > &vSaveZones)
Save a triangulation.
Zone2 * createZone(ConstraintGraph2 *pConstraintGraph, ZoneLocation zoneLoc, bool bVerbose=true)
Create a zone.
Voronoi2 * getVoronoiDiagram()
Get the Voronoi diagram.
Point2 * getNearestNeighbor(const Point2 &p)
Get nearest neighbor.
bool saveZones(std::ostream &stream, std::vector< Zone2 * > &vSaveZones)
Save zones.
bool saveTriangulation(std::ostream &stream, std::vector< Zone2 * > &vSaveZones)
Save a triangulation.
void applyConstraintsAndZones()
Apply conforming constraints and zones (deprecated!)
void reset()
Reset.
size_t numberOfPoints() const
Number of points.
~Fade_2D()
Destructor.
void getAliveAndDeadConstraintSegments(std::vector< ConstraintSegment2 * > &vAllConstraintSegments) const
Get all (alive and dead) constraint segments.
void getDirectedEdges(std::vector< Edge2 > &vDirectedEdgesOut) const
Get directed edges Edges are counterclockwise oriented around their triangle. The present method retu...
void refineAdvanced(MeshGenParams *pParameters)
Delaunay refinement and grid meshing.
bool isConstraint(Point2 *pVtx) const
Check if a vertex is a constraint vertex.
bool load(const char *filename, std::vector< Zone2 * > &vZones)
Load a triangulation.
void cutTriangles(const Point2 &knifeStart, const Point2 &knifeEnd, bool bTurnEdgesIntoConstraints)
Cut through a triangulation.
void exportTriangulation(FadeExport &fadeExport, bool bWithCustomIndices, bool bClear)
Export triangulation data from Fade.
void statistics(const char *s) const
Statistics.
void getAliveConstraintSegments(std::vector< ConstraintSegment2 * > &vAliveConstraintSegments) const
Get active (alive) constraint segments.
void refine(Zone2 *pZone, double minAngleDegree, double minEdgeLength, double maxEdgeLength, bool bAllowConstraintSplitting)
Delaunay refinement.
void cutTriangles(std::vector< Segment2 > &vSegments, bool bTurnEdgesIntoConstraints)
Cut through a triangulation.
bool checkValidity(bool bCheckEmptyCircleProperty, const char *msg) const
Checks if a triangulation is valid.
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.
bool saveZones(const char *filename, std::vector< Zone2 * > &vSaveZones)
Save zones.
bool isConstraint(Triangle2 *pT, int ith) const
Check if an edge is a constraint edge.
bool isConstraint(Point2 *p0, Point2 *p1) const
Check if an edge is a constraint edge.
void writeWebScene(const char *path) const
Write the current triangulation to an *.obj file.
void getIncidentVertices(Point2 *pVtx, std::vector< Point2 * > &vIncidentVertices) const
Get incident vertices.
void insert(int numPoints, double *aCoordinates, Point2 **aHandles)
Insert points from an array.
void show(const char *postscriptFilename, bool bWithConstraints=true) const
Draws the triangulation as postscript file.
Zone2 * importTriangles(std::vector< Point2 > &vPoints, bool bReorientIfNeeded, bool bCreateExtendedBoundingBox)
Import triangles.
size_t numberOfTriangles() const
Number of triangles.
void deleteZone(Zone2 *pZone)
Delete a Zone2 object.
Zone2 * createZone(std::vector< Triangle2 * > &vTriangles, bool bVerbose=true)
Create a zone defined by a vector of triangles.
Triangle2 * locate(const Point2 &p)
Locate a triangle which contains p.
Parameters for the mesh generator.
Definition: MeshGenParams.h:58
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
Voronoi diagram.
Definition: Voronoi2.h:56
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