| OpenMesh

BaseDecimaterT_impl.hh
1 /* ========================================================================= *
2 * *
3 * OpenMesh *
4 * Copyright (c) 2001-2015, RWTH-Aachen University *
5 * Department of Computer Graphics and Multimedia *
6 * All rights reserved. *
7 * www.openmesh.org *
8 * *
9 *---------------------------------------------------------------------------*
10 * This file is part of OpenMesh. *
11 *---------------------------------------------------------------------------*
12 * *
13 * Redistribution and use in source and binary forms, with or without *
14 * modification, are permitted provided that the following conditions *
15 * are met: *
16 * *
17 * 1. Redistributions of source code must retain the above copyright notice, *
18 * this list of conditions and the following disclaimer. *
19 * *
20 * 2. Redistributions in binary form must reproduce the above copyright *
21 * notice, this list of conditions and the following disclaimer in the *
22 * documentation and/or other materials provided with the distribution. *
23 * *
24 * 3. Neither the name of the copyright holder nor the names of its *
25 * contributors may be used to endorse or promote products derived from *
26 * this software without specific prior written permission. *
27 * *
28 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS *
29 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED *
30 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A *
31 * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER *
32 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, *
33 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, *
34 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR *
35 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF *
36 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING *
37 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
38 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
39 * *
40 * ========================================================================= */
41
45 //=============================================================================
46 //
47 // CLASS DecimaterT - IMPLEMENTATION
48 //
49 //=============================================================================
50 #define OPENMESH_BASE_DECIMATER_DECIMATERT_CC
51
52 //== INCLUDES =================================================================
53
54 #include <vector>
55 #if defined(OM_CC_MIPS)
56 # include <float.h>
57 #else
58 # include <cfloat>
59 #endif
60
61 //== NAMESPACE ===============================================================
62
63 namespace OpenMesh {
64 namespace Decimater {
65
66 //== IMPLEMENTATION ==========================================================
67
68 template<class Mesh>
69 BaseDecimaterT<Mesh>::BaseDecimaterT(Mesh& _mesh) :
70 mesh_(_mesh), cmodule_(NULL), initialized_(false), observer_(NULL) {
71 // default properties
72 mesh_.request_vertex_status();
73 mesh_.request_edge_status();
74 mesh_.request_face_status();
75
76 }
77
78 //-----------------------------------------------------------------------------
79
80 template<class Mesh>
81 BaseDecimaterT<Mesh>::~BaseDecimaterT() {
82 // default properties
83 mesh_.release_vertex_status();
84 mesh_.release_edge_status();
85 mesh_.release_face_status();
86
87 // dispose of modules
88 {
89 set_uninitialized();
90 typename ModuleList::iterator m_it, m_end = all_modules_.end();
91 for (m_it = all_modules_.begin(); m_it != m_end; ++m_it)
92 delete *m_it;
93 all_modules_.clear();
94 }
95 }
96
97 //-----------------------------------------------------------------------------
98
99 template<class Mesh>
100 bool BaseDecimaterT<Mesh>::is_collapse_legal(const CollapseInfo& _ci) {
101 // std::clog << "McDecimaterT<>::is_collapse_legal()\n";
102
103 // locked ?
104 if (mesh_.status(_ci.v0).locked())
105 return false;
106
107 // this test checks:
108 // is v0v1 deleted?
109 // is v0 deleted?
110 // is v1 deleted?
111 // are both vlv0 and v1vl boundary edges?
112 // are both v0vr and vrv1 boundary edges?
113 // are vl and vr equal or both invalid?
114 // one ring intersection test
115 // edge between two boundary vertices should be a boundary edge
116 if (!mesh_.is_collapse_ok(_ci.v0v1))
117 return false;
118
119 if (_ci.vl.is_valid() && _ci.vr.is_valid()
120 && mesh_.find_halfedge(_ci.vl, _ci.vr).is_valid()
121 && mesh_.valence(_ci.vl) == 3 && mesh_.valence(_ci.vr) == 3) {
122 return false;
123 }
124 //--- feature test ---
125
126 if (mesh_.status(_ci.v0).feature()
127 && !mesh_.status(mesh_.edge_handle(_ci.v0v1)).feature())
128 return false;
129
130 //--- test boundary cases ---
131 if (mesh_.is_boundary(_ci.v0)) {
132
133 // don't collapse a boundary vertex to an inner one
134 if (!mesh_.is_boundary(_ci.v1))
135 return false;
136
137 // only one one ring intersection
138 if (_ci.vl.is_valid() && _ci.vr.is_valid())
139 return false;
140 }
141
142 // there have to be at least 2 incident faces at v0
143 if (mesh_.cw_rotated_halfedge_handle(
144 mesh_.cw_rotated_halfedge_handle(_ci.v0v1)) == _ci.v0v1)
145 return false;
146
147 // collapse passed all tests -> ok
148 return true;
149 }
150
151 //-----------------------------------------------------------------------------
152
153 template<class Mesh>
154 float BaseDecimaterT<Mesh>::collapse_priority(const CollapseInfo& _ci) {
155 typename ModuleList::iterator m_it, m_end = bmodules_.end();
156
157 for (m_it = bmodules_.begin(); m_it != m_end; ++m_it) {
158 if ((*m_it)->collapse_priority(_ci) < 0.0)
159 return ModBaseT< Mesh >::ILLEGAL_COLLAPSE;
160 }
161 return cmodule_->collapse_priority(_ci);
162 }
163
164 //-----------------------------------------------------------------------------
165
166 template<class Mesh>
167 void BaseDecimaterT<Mesh>::postprocess_collapse(CollapseInfo& _ci) {
168 typename ModuleList::iterator m_it, m_end = bmodules_.end();
169
170 for (m_it = bmodules_.begin(); m_it != m_end; ++m_it)
171 (*m_it)->postprocess_collapse(_ci);
172
173 cmodule_->postprocess_collapse(_ci);
174 }
175
176 //-----------------------------------------------------------------------------
177
178 template<class Mesh>
179 void BaseDecimaterT<Mesh>::preprocess_collapse(CollapseInfo& _ci) {
180 typename ModuleList::iterator m_it, m_end = bmodules_.end();
181
182 for (m_it = bmodules_.begin(); m_it != m_end; ++m_it)
183 (*m_it)->preprocess_collapse(_ci);
184
185 cmodule_->preprocess_collapse(_ci);
186 }
187
188 //-----------------------------------------------------------------------------
189
190 template<class Mesh>
191 void BaseDecimaterT<Mesh>::set_error_tolerance_factor(double _factor) {
192 if (_factor >= 0.0 && _factor <= 1.0) {
193 typename ModuleList::iterator m_it, m_end = bmodules_.end();
194
195 for (m_it = bmodules_.begin(); m_it != m_end; ++m_it)
196 (*m_it)->set_error_tolerance_factor(_factor);
197
198 cmodule_->set_error_tolerance_factor(_factor);
199 }
200 }
201
202 //-----------------------------------------------------------------------------
203
204 template<class Mesh>
205 void BaseDecimaterT<Mesh>::info(std::ostream& _os) {
206 if (initialized_) {
207 _os << "initialized : yes" << std::endl;
208 _os << "binary modules: " << bmodules_.size() << std::endl;
209 for (ModuleListIterator m_it = bmodules_.begin(); m_it != bmodules_.end();
210 ++m_it) {
211 _os << " " << (*m_it)->name() << std::endl;
212 }
213 _os << "priority module: " << cmodule_->name().c_str() << std::endl;
214 } else {
215 _os << "initialized : no" << std::endl;
216 _os << "available modules: " << all_modules_.size() << std::endl;
217 for (ModuleListIterator m_it = all_modules_.begin();
218 m_it != all_modules_.end(); ++m_it) {
219 _os << " " << (*m_it)->name() << " : ";
220 if ((*m_it)->is_binary()) {
221 _os << "binary";
222 if ((*m_it)->name() == "Quadric") {
223 _os << " and priority (special treatment)";
224 }
225 } else {
226 _os << "priority";
227 }
228 _os << std::endl;
229 }
230 }
231 }
232
233 //-----------------------------------------------------------------------------
234
235 template<class Mesh>
236 bool BaseDecimaterT<Mesh>::initialize() {
237 if (initialized_) {
238 return true;
239 }
240
241 // FIXME: quadric module shouldn't be treated specially.
242 // Q: Why?
243 // A: It isn't generic and breaks encapsulation. Also, using string
244 // name comparison is not reliable, since you can't guarantee that
245 // no one else will name their custom module "Quadric".
246 // Q: What should be done instead?
247 // A: ModBaseT API should support modules that can be both binary
248 // and priority, or BETTER YET, let the DecimaterT API specify the
249 // priority module explicitly.
250
251 // find the priority module: either the only non-binary module in the list, or "Quadric"
252 Module *quadric = NULL;
253 Module *pmodule = NULL;
254 for (ModuleListIterator m_it = all_modules_.begin(), m_end =
255 all_modules_.end(); m_it != m_end; ++m_it) {
256 if ((*m_it)->name() == "Quadric")
257 quadric = *m_it;
258
259 if (!(*m_it)->is_binary()) {
260 if (pmodule) {
261 // only one priority module allowed!
262 set_uninitialized();
263 return false;
264 }
265 pmodule = *m_it;
266 }
267 }
268
269 // Quadric is used as default priority module (even if it is set to be binary)
270 if (!pmodule && quadric) {
271 pmodule = quadric;
272 }
273
274 if (!pmodule) {
275 // At least one priority module required
276 set_uninitialized();
277 return false;
278 }
279
280 // set pmodule as the current priority module
281 cmodule_ = pmodule;
282
283 for (ModuleListIterator m_it = all_modules_.begin(), m_end =
284 all_modules_.end(); m_it != m_end; ++m_it) {
285 // every module gets initialized
286 (*m_it)->initialize();
287
288 if (*m_it != pmodule) {
289 // all other modules are binary, and go into bmodules_ list
290 bmodules_.push_back(*m_it);
291 }
292 }
293
294 return initialized_ = true;
295 }
296
297
298
299 //=============================================================================
300 }// END_NS_DECIMATER
301 } // END_NS_OPENMESH
302 //=============================================================================
303
OpenMesh::Decimater::ModBaseT
Base class for all decimation modules.
Definition: ModBaseT.hh:192
OpenMesh::Decimater::CollapseInfoT::vr
Mesh::VertexHandle vr
Right vertex.
Definition: CollapseInfoT.hh:98
OpenMesh::Decimater::CollapseInfoT::v0v1
Mesh::HalfedgeHandle v0v1
Halfedge to be collapsed.
Definition: CollapseInfoT.hh:89
OpenMesh::Decimater::BaseDecimaterT::set_uninitialized
void set_uninitialized()
Reset the initialized flag, and clear the bmodules_ and cmodule_.
Definition: BaseDecimaterT.hh:202
OpenMesh::Decimater::ModBaseT::name
virtual const std::string & name() const
Set module&#39;s name (using DECIMATER_MODNAME macro)
Definition: ModBaseT.hh:216
OpenMesh::Decimater::ModBaseT::collapse_priority
virtual float collapse_priority(const CollapseInfoT< MeshT > &)
Return collapse priority.
Definition: ModBaseT.hh:245
OpenMesh::Decimater::ModBaseT::postprocess_collapse
virtual void postprocess_collapse(const CollapseInfoT< MeshT > &)
After _from_vh has been collapsed into _to_vh, this method will be called.
Definition: ModBaseT.hh:257
OpenMesh::Decimater::BaseDecimaterT::set_error_tolerance_factor
void set_error_tolerance_factor(double _factor)
This provides a function that allows the setting of a percentage of the original constraint of the mo...
Definition: BaseDecimaterT_impl.hh:191
OpenMesh::Decimater::BaseDecimaterT::initialize
bool initialize()
Initialize decimater and decimating modules.
Definition: BaseDecimaterT_impl.hh:236
OpenMesh::Decimater::ModBaseT::preprocess_collapse
virtual void preprocess_collapse(const CollapseInfoT< MeshT > &)
Before _from_vh has been collapsed into _to_vh, this method will be called.
Definition: ModBaseT.hh:251
OpenMesh::Decimater::CollapseInfoT::vl
Mesh::VertexHandle vl
Left vertex.
Definition: CollapseInfoT.hh:97
OpenMesh::Decimater::BaseDecimaterT::collapse_priority
float collapse_priority(const CollapseInfo &_ci)
Calculate priority of an halfedge collapse (using the modules)
Definition: BaseDecimaterT_impl.hh:154
OpenMesh
Contains all the mesh ingredients like the polygonal mesh, the triangle mesh, different mesh kernels ...
Definition: MeshItems.hh:59
OpenMesh::Decimater::ModBaseT::set_error_tolerance_factor
virtual void set_error_tolerance_factor(double _factor)
This provides a function that allows the setting of a percentage of the original contraint.
Definition: ModBaseT.hh:268
OpenMesh::Decimater::CollapseInfoT::v1
Mesh::VertexHandle v1
Remaining vertex.
Definition: CollapseInfoT.hh:92
OpenMesh::Decimater::BaseDecimaterT::is_collapse_legal
bool is_collapse_legal(const CollapseInfo &_ci)
Is an edge collapse legal? Performs topological test only.
Definition: BaseDecimaterT_impl.hh:100
OpenMesh::Decimater::BaseDecimaterT::preprocess_collapse
void preprocess_collapse(CollapseInfo &_ci)
Pre-process a collapse.
Definition: BaseDecimaterT_impl.hh:179
OpenMesh::Decimater::ModBaseT::initialize
virtual void initialize()
Initialize module-internal stuff.
Definition: ModBaseT.hh:229
OpenMesh::Decimater::CollapseInfoT
Stores information about a halfedge collapse.
Definition: CollapseInfoT.hh:74
OpenMesh::Decimater::BaseDecimaterT::postprocess_collapse
void postprocess_collapse(CollapseInfo &_ci)
Post-process a collapse.
Definition: BaseDecimaterT_impl.hh:167
OpenMesh::Decimater::BaseDecimaterT::info
void info(std::ostream &_os)
Print information about modules to _os.
Definition: BaseDecimaterT_impl.hh:205
OpenMesh::Decimater::CollapseInfoT::v0
Mesh::VertexHandle v0
Vertex to be removed.
Definition: CollapseInfoT.hh:91