57 #ifndef OPENMESH_POLYMESHT_HH
58 #define OPENMESH_POLYMESHT_HH
64 #include <OpenMesh/Core/System/config.h>
65 #include <OpenMesh/Core/Geometry/MathDefs.hh>
66 #include <OpenMesh/Core/Mesh/PolyConnectivity.hh>
67 #include <OpenMesh/Core/Mesh/FinalMeshItemsT.hh>
93 template <
class Kernel>
103 enum { IsPolyMesh = 1 };
105 enum { IsTriMesh = 0 };
112 typedef typename Kernel::Scalar
Scalar;
115 typedef typename Kernel::Point
Point;
131 typedef typename Kernel::Edge
Edge;
133 typedef typename Kernel::Face
Face;
140 typedef typename Kernel::HalfedgeHandle HalfedgeHandle;
141 typedef typename Kernel::EdgeHandle EdgeHandle;
142 typedef typename Kernel::FaceHandle FaceHandle;
146 typedef typename Kernel::VertexIter VertexIter;
147 typedef typename Kernel::HalfedgeIter HalfedgeIter;
148 typedef typename Kernel::EdgeIter EdgeIter;
149 typedef typename Kernel::FaceIter FaceIter;
151 typedef typename Kernel::ConstVertexIter ConstVertexIter;
152 typedef typename Kernel::ConstHalfedgeIter ConstHalfedgeIter;
153 typedef typename Kernel::ConstEdgeIter ConstEdgeIter;
154 typedef typename Kernel::ConstFaceIter ConstFaceIter;
203 {
return Kernel::new_vertex(); }
212 VertexHandle vh(Kernel::new_vertex());
213 this->set_point(vh, _p);
230 VertexHandle vh(Kernel::new_vertex_dirty());
231 this->set_point(vh, _p);
274 const Point& _p2)
const;
310 virtual Normal
calc_halfedge_normal(HalfedgeHandle _heh,
const double _feature_angle = 0.8)
const;
387 return this->point(this->to_vertex_handle(_heh)) -
388 this->point(this->from_vertex_handle(_heh));
392 Scalar calc_edge_length(EdgeHandle _eh)
const
393 {
return calc_edge_length(this->halfedge_handle(_eh,0)); }
398 {
return (Scalar)sqrt(calc_edge_sqr_length(_heh)); }
400 Scalar calc_edge_sqr_length(EdgeHandle _eh)
const
401 {
return calc_edge_sqr_length(this->halfedge_handle(_eh,0)); }
403 Scalar calc_edge_sqr_length(HalfedgeHandle _heh)
const
407 return edge_vec.sqrnorm();
429 Scalar denom = v0.norm()*v1.norm();
434 Scalar cos_a =
dot(v0 , v1) / denom;
435 if (this->is_boundary(_in_heh))
437 FaceHandle fh(this->face_handle(this->opposite_halfedge_handle(_in_heh)));
439 Scalar sign_a =
dot(
cross(v0, v1), f_n);
440 return angle(cos_a, sign_a);
474 _sector_normal =
cross(vec0, vec1);
482 Normal sector_normal;
484 return sector_normal.norm()/2;
492 assert(Kernel::has_face_normals());
494 if (this->is_boundary(this->edge_handle(_heh)))
498 const Normal& n0 = this->normal(this->face_handle(_heh));
499 const Normal& n1 = this->normal(this->face_handle(this->opposite_halfedge_handle(_heh)));
502 Scalar da_cos =
dot(n0, n1);
504 Scalar da_sin_sign =
dot(
cross(n0, n1), he);
505 return angle(da_cos, da_sin_sign);
514 Scalar calc_dihedral_angle(HalfedgeHandle _heh)
const
516 if (this->is_boundary(this->edge_handle(_heh)))
524 Scalar denom = n0.norm()*n1.norm();
529 Scalar da_cos =
dot(n0, n1)/denom;
531 Scalar da_sin_sign =
dot(
cross(n0, n1), he);
532 return angle(da_cos, da_sin_sign);
536 Scalar calc_dihedral_angle(EdgeHandle _eh)
const
537 {
return calc_dihedral_angle(this->halfedge_handle(_eh,0)); }
545 inline void split(FaceHandle _fh,
const Point& _p)
548 inline void split(FaceHandle _fh, VertexHandle _vh)
549 { Kernel::split(_fh, _vh); }
551 inline void split(EdgeHandle _eh,
const Point& _p)
554 inline void split(EdgeHandle _eh, VertexHandle _vh)
555 { Kernel::split_edge(_eh, _vh); }
558 struct PointIs3DTag {};
559 struct PointIsNot3DTag {};
560 Normal calc_face_normal_impl(FaceHandle, PointIs3DTag)
const;
561 Normal calc_face_normal_impl(FaceHandle, PointIsNot3DTag)
const;
562 Normal calc_face_normal_impl(
const Point&,
const Point&,
const Point&, PointIs3DTag)
const;
563 Normal calc_face_normal_impl(
const Point&,
const Point&,
const Point&, PointIsNot3DTag)
const;
591 template<
typename LHS,
typename KERNEL>
596 template<
typename LHS,
typename KERNEL>
598 return MeshCast<LHS, PolyMeshT<KERNEL>*>::cast(rhs);
601 template<
typename LHS,
typename KERNEL>
602 const LHS
mesh_cast(
const PolyMeshT<KERNEL> &rhs) {
603 return MeshCast<LHS, const PolyMeshT<KERNEL>&>::cast(rhs);
606 template<
typename LHS,
typename KERNEL>
607 const LHS
mesh_cast(
const PolyMeshT<KERNEL> *rhs) {
608 return MeshCast<LHS, const PolyMeshT<KERNEL>*>::cast(rhs);
615 #if defined(OM_INCLUDE_TEMPLATES) && !defined(OPENMESH_POLYMESH_C)
616 # define OPENMESH_POLYMESH_TEMPLATES
617 # include "PolyMeshT.cc"
620 #endif // OPENMESH_POLYMESHT_HH defined
Add normals to mesh item (vertices/faces)
Definition: Attributes.hh:87
void calc_vertex_normal_fast(VertexHandle _vh, Normal &_n) const
Different methods for calculation of the normal at _vh:
Definition: PolyMeshT.cc:387
Normal calc_edge_vector(HalfedgeHandle _heh) const
Calculates the edge vector as the difference of the the points defined by to_vertex_handle() and from...
Definition: PolyMeshT.hh:385
Kernel::Point Point
Coordinate type.
Definition: PolyMeshT.hh:115
Add colors to mesh item (vertices/faces/edges)
Definition: Attributes.hh:88
Kernel::VertexEdgeIter VertexEdgeIter
Circulator.
Definition: PolyMeshT.hh:168
LHS mesh_cast(PolyMeshT< KERNEL > &rhs)
Cast a mesh with different but identical traits into each other.
Definition: PolyMeshT.hh:592
Kernel::Normal Normal
Normal type.
Definition: PolyMeshT.hh:117
virtual Normal calc_face_normal(FaceHandle _fh) const
Calculate normal vector for face _fh.
Definition: PolyMeshT.cc:102
Kernel::ConstVertexEdgeIter ConstVertexEdgeIter
Circulator.
Definition: PolyMeshT.hh:178
Kernel::Halfedge Halfedge
Halfedge type.
Definition: PolyMeshT.hh:129
PolyMeshT< Kernel > This
Self type. Used to specify iterators/circulators.
Definition: PolyMeshT.hh:99
Kernel::Edge Edge
Edge type.
Definition: PolyMeshT.hh:131
static bool is_trimesh()
Determine whether this is a PolyMeshT or TriMeshT ( This function does not check the per face vertex ...
Definition: PolyMeshT.hh:107
Add 1D texture coordinates (vertices, halfedges)
Definition: Attributes.hh:91
Kernel::VertexVertexIter VertexVertexIter
Circulator.
Definition: PolyMeshT.hh:165
Kernel::ConstVertexOHalfedgeIter ConstVertexOHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:176
void calc_edge_vector(HalfedgeHandle _heh, Normal &_edge_vec) const
Calculates the edge vector as the difference of the the points defined by to_vertex_handle() and from...
Definition: PolyMeshT.hh:378
Kernel::VertexOHalfedgeIter VertexOHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:166
virtual Normal calc_halfedge_normal(HalfedgeHandle _heh, const double _feature_angle=0.8) const
Calculate halfedge normal for one specific halfedge.
Definition: PolyMeshT.cc:289
Kernel::TexCoord2D TexCoord2D
TexCoord2D type.
Definition: PolyMeshT.hh:123
Base type for a polygonal mesh.
Definition: PolyMeshT.hh:94
void calc_sector_vectors(HalfedgeHandle _in_heh, Normal &_vec0, Normal &_vec1) const
defines a consistent representation of a sector geometry: the halfedge _in_heh defines the sector ori...
Definition: PolyMeshT.hh:414
Kernel::Face Face
Face type.
Definition: PolyMeshT.hh:133
Kernel::ConstVertexVertexIter ConstVertexVertexIter
Circulator.
Definition: PolyMeshT.hh:175
VertexHandle new_vertex_dirty(const Point &_p)
Same as new_vertex(const Point&) but never shrinks, only enlarges the vertex property vectors...
Definition: PolyMeshT.hh:228
Kernel::ConstFaceVertexIter ConstFaceVertexIter
Circulator.
Definition: PolyMeshT.hh:180
Normal calc_vertex_normal(VertexHandle _vh) const
Calculate vertex normal for one specific vertex.
Definition: PolyMeshT.cc:373
Scalar calc_edge_length(HalfedgeHandle _heh) const
Calculates the length of the edge _heh.
Definition: PolyMeshT.hh:397
Kernel::Vertex Vertex
Vertex type.
Definition: PolyMeshT.hh:127
Kernel::FaceVertexIter FaceVertexIter
Circulator.
Definition: PolyMeshT.hh:170
VertexHandle add_vertex(const Point &_p)
Alias for new_vertex(const Point&).
Definition: PolyMeshT.hh:236
Normal calc_edge_vector(EdgeHandle _eh) const
Calculates the edge vector as the vector defined by the halfedge with id #0 (see below) ...
Definition: PolyMeshT.hh:371
Kernel::TexCoord1D TexCoord1D
TexCoord1D type.
Definition: PolyMeshT.hh:121
Kernel::ConstVertexIHalfedgeIter ConstVertexIHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:177
Scalar calc_sector_area(HalfedgeHandle _in_heh) const
calculates the area of the face sector defined by the angle <(_in_heh,next_halfedge(_in_heh)) NOTE: s...
Definition: PolyMeshT.hh:480
T sane_aarg(T _aarg)
Trigonometry/angles - related.
Definition: MathDefs.hh:127
Add 2D texture coordinates (vertices, halfedges)
Definition: Attributes.hh:92
Kernel::ConstFaceHalfedgeIter ConstFaceHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:181
Kernel::VertexHandle VertexHandle
Handle for referencing the corresponding item.
Definition: PolyMeshT.hh:139
void split(FaceHandle _fh, const Point &_p)
Face split (= 1-to-n split)
Definition: PolyMeshT.hh:545
VertexHandle new_vertex()
Uses default copy and assignment operator.
Definition: PolyMeshT.hh:202
void calc_vertex_normal_correct(VertexHandle _vh, Normal &_n) const
Compute normals for all primitives.
Definition: PolyMeshT.cc:397
Cast a mesh with different but identical traits into each other.
Definition: FinalMeshItemsT.hh:183
Contains all the mesh ingredients like the polygonal mesh, the triangle mesh, different mesh kernels ...
Definition: MeshItems.hh:64
void update_vertex_normals()
Update normal vectors for all vertices.
Definition: PolyMeshT.cc:448
Kernel::Scalar Scalar
Scalar type.
Definition: PolyMeshT.hh:113
bool is_estimated_feature_edge(HalfedgeHandle _heh, const double _feature_angle) const
identifies feature edges w.r.t.
Definition: PolyMeshT.cc:342
VertexHandle add_vertex_dirty(const Point &_p)
Alias for new_vertex_dirty().
Definition: PolyMeshT.hh:240
Scalar calc_dihedral_angle_fast(HalfedgeHandle _heh) const
calculates the dihedral angle on the halfedge _heh
Definition: PolyMeshT.hh:489
T angle(T _cos_angle, T _sin_angle)
returns the angle determined by its cos and the sign of its sin result is positive if the angle is in...
Definition: MathDefs.hh:145
void calc_face_centroid(FaceHandle _fh, Point &_pt) const
calculates the average of the vertices defining _fh
Definition: PolyMeshT.hh:276
bool is_zero(const T &_a, Real _eps)
comparison operators with user-selected precision control
Definition: MathDefs.hh:66
void update_halfedge_normals(const double _feature_angle=0.8)
Update normal vectors for all halfedges.
Definition: PolyMeshT.cc:274
void update_normal(FaceHandle _fh)
Update normal for face _fh.
Definition: PolyMeshT.hh:259
Kernel::FaceFaceIter FaceFaceIter
Circulator.
Definition: PolyMeshT.hh:173
Kernel::FaceHalfedgeIter FaceHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:171
Kernel::ConstFaceEdgeIter ConstFaceEdgeIter
Circulator.
Definition: PolyMeshT.hh:182
void update_face_normals()
Update normal vectors for all faces.
Definition: PolyMeshT.cc:259
osg::Vec3f cross(const osg::Vec3f &_v1, const osg::Vec3f &_v2)
Adapter for osg vector member computing a scalar product.
Definition: VectorAdapter.hh:196
Kernel::Color Color
Color type.
Definition: PolyMeshT.hh:119
Add 3D texture coordinates (vertices, halfedges)
Definition: Attributes.hh:93
Kernel::ConstVertexFaceIter ConstVertexFaceIter
Circulator.
Definition: PolyMeshT.hh:179
VertexHandle new_vertex(const Point &_p)
Adds a new vertex initialized to a custom position.
Definition: PolyMeshT.hh:210
void calc_vertex_normal_loop(VertexHandle _vh, Normal &_n) const
Compute normals for all primitives.
Definition: PolyMeshT.cc:425
Kernel::FaceEdgeIter FaceEdgeIter
Circulator.
Definition: PolyMeshT.hh:172
void calc_edge_vector(EdgeHandle _eh, Normal &_edge_vec) const
Calculates the edge vector as the vector defined by the halfedge with id #0 (see below) ...
Definition: PolyMeshT.hh:364
Scalar calc_dihedral_angle_fast(EdgeHandle _eh) const
calculates the dihedral angle on the edge _eh
Definition: PolyMeshT.hh:510
unsigned int find_feature_edges(Scalar _angle_tresh=OpenMesh::deg_to_rad(44.0))
tags an edge as a feature if its dihedral angle is larger than _angle_tresh returns the number of the...
Definition: PolyMeshT.cc:79
void update_normals()
Compute normals for all primitives.
Definition: PolyMeshT.cc:241
Kernel::VertexIHalfedgeIter VertexIHalfedgeIter
Circulator.
Definition: PolyMeshT.hh:167
void calc_sector_normal(HalfedgeHandle _in_heh, Normal &_sector_normal) const
calculates the normal (non-normalized) of the face sector defined by the angle <(_in_heh,next_halfedge(_in_heh))
Definition: PolyMeshT.hh:470
Kernel::ConstFaceFaceIter ConstFaceFaceIter
Circulator.
Definition: PolyMeshT.hh:183
void update_normal(VertexHandle _vh)
Update normal for vertex _vh.
Definition: PolyMeshT.hh:318
Scalar calc_sector_angle(HalfedgeHandle _in_heh) const
calculates the sector angle.
Definition: PolyMeshT.hh:425
Kernel::VertexFaceIter VertexFaceIter
Circulator.
Definition: PolyMeshT.hh:169
void update_normal(HalfedgeHandle _heh, const double _feature_angle=0.8)
Update normal for halfedge _heh.
Definition: PolyMeshT.hh:284
Kernel::TexCoord3D TexCoord3D
TexCoord3D type.
Definition: PolyMeshT.hh:125
static bool is_polymesh()
Determine whether this is a PolyMeshT or TriMeshT ( This function does not check the per face vertex ...
Definition: PolyMeshT.hh:106
osg::Vec3f::ValueType dot(const osg::Vec3f &_v1, const osg::Vec3f &_v2)
Adapter for osg vector member computing a scalar product.
Definition: VectorAdapter.hh:181