ModNormalDeviationT.hh 9.1 KB
Newer Older
Jan Möbius's avatar
Jan Möbius committed
1
2
3
/*===========================================================================*\
 *                                                                           *
 *                               OpenMesh                                    *
Jan Möbius's avatar
Jan Möbius committed
4
 *      Copyright (C) 2001-2012 by Computer Graphics Group, RWTH Aachen      *
Jan Möbius's avatar
Jan Möbius committed
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
 *                           www.openmesh.org                                *
 *                                                                           *
 *---------------------------------------------------------------------------*
 *  This file is part of OpenMesh.                                           *
 *                                                                           *
 *  OpenMesh is free software: you can redistribute it and/or modify         *
 *  it under the terms of the GNU Lesser General Public License as           *
 *  published by the Free Software Foundation, either version 3 of           *
 *  the License, or (at your option) any later version with the              *
 *  following exceptions:                                                    *
 *                                                                           *
 *  If other files instantiate templates or use macros                       *
 *  or inline functions from this file, or you compile this file and         *
 *  link it with other files to produce an executable, this file does        *
 *  not by itself cause the resulting executable to be covered by the        *
 *  GNU Lesser General Public License. This exception does not however       *
 *  invalidate any other reasons why the executable file might be            *
 *  covered by the GNU Lesser General Public License.                        *
 *                                                                           *
 *  OpenMesh is distributed in the hope that it will be useful,              *
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of           *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the            *
 *  GNU Lesser General Public License for more details.                      *
 *                                                                           *
 *  You should have received a copy of the GNU LesserGeneral Public          *
 *  License along with OpenMesh.  If not,                                    *
 *  see <http://www.gnu.org/licenses/>.                                      *
 *                                                                           *
\*===========================================================================*/

/*===========================================================================*\
 *                                                                           *
Jan Möbius's avatar
Jan Möbius committed
37
38
 *   $Revision$                                                        *
 *   $Date$                   *
Jan Möbius's avatar
Jan Möbius committed
39
40
41
42
43
44
45
46
47
48
49
50
51
 *                                                                           *
\*===========================================================================*/

/** \file ModNormalDeviationT.hh
 */

//=============================================================================
//
//  CLASS ModNormalDeviationT
//
//=============================================================================


52
53
#ifndef OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH
#define OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH
Jan Möbius's avatar
Jan Möbius committed
54
55
56
57
58
59


//== INCLUDES =================================================================

#include <OpenMesh/Tools/Decimater/ModBaseT.hh>
#include <OpenMesh/Core/Utils/Property.hh>
60
#include <OpenMesh/Core/Geometry/NormalConeT.hh>
Jan Möbius's avatar
Jan Möbius committed
61
62
63
64
65
66
67
68
69
70
71


//== NAMESPACES ===============================================================

namespace OpenMesh  {
namespace Decimater {


//== CLASS DEFINITION =========================================================


72
73
74
75
76
77
78
79
80
81
82
83
84
85
/** \brief Use Normal deviation to control decimation
  *
  * The module tracks the normals while decimating
  * a normal cone consisting of all normals of the
  * faces collapsed together is computed and if
  * a collapse would increase the size of
  * the cone to a value greater than the given value
  * the collapse will be illegal.
  *
  * In binary and mode, the collapse is legal if:
  *  - The normal deviation after the collapse is lower than the given value
  *
  * In continuous mode the maximal deviation is returned
  */
86
87
template <class MeshT>
class ModNormalDeviationT : public ModBaseT< MeshT >
Jan Möbius's avatar
Jan Möbius committed
88
89
{
public:
90

91
  DECIMATING_MODULE( ModNormalDeviationT, MeshT, NormalDeviation );
Jan Möbius's avatar
Jan Möbius committed
92
93
94
95
96
97
98

  typedef typename Mesh::Scalar                     Scalar;
  typedef typename Mesh::Point                      Point;
  typedef typename Mesh::Normal                     Normal;
  typedef typename Mesh::VertexHandle               VertexHandle;
  typedef typename Mesh::FaceHandle                 FaceHandle;
  typedef typename Mesh::EdgeHandle                 EdgeHandle;
99
  typedef NormalConeT<Scalar>                       NormalCone;
Jan Möbius's avatar
Jan Möbius committed
100
101
102
103
104
105



public:

  /// Constructor
106
107
  ModNormalDeviationT(MeshT& _mesh, float _max_dev = 180.0)
  : Base(_mesh, true), mesh_(Base::mesh())
Jan Möbius's avatar
Jan Möbius committed
108
109
110
111
  {
    set_normal_deviation(_max_dev);
    mesh_.add_property(normal_cones_);
  }
112

Jan Möbius's avatar
Jan Möbius committed
113
114
115
116
117
118
119

  /// Destructor
  ~ModNormalDeviationT() {
    mesh_.add_property(normal_cones_);
  }


120
  /// Get normal deviation ( 0 .. 360 )
Jan Möbius's avatar
Jan Möbius committed
121
122
123
124
  Scalar normal_deviation() const {
    return normal_deviation_ / M_PI * 180.0;
  }

125
  /// Set normal deviation ( 0 .. 360 )
Jan Möbius's avatar
Jan Möbius committed
126
127
128
129
130
131
132
133
134
135
136
  void set_normal_deviation(Scalar _s) {
    normal_deviation_ = _s / 180.0 * M_PI;
  }


  /// Allocate and init normal cones
  void  initialize() {
    if (!normal_cones_.is_valid())
      mesh_.add_property(normal_cones_);

    typename Mesh::FaceIter f_it  = mesh_.faces_begin(),
137
138
        f_end = mesh_.faces_end();

Jan Möbius's avatar
Jan Möbius committed
139
    for (; f_it != f_end; ++f_it)
Jan Möbius's avatar
Jan Möbius committed
140
      mesh_.property(normal_cones_, *f_it) = NormalCone(mesh_.normal(*f_it));
Jan Möbius's avatar
Jan Möbius committed
141
142
  }

143
144
145
146
147
148
149
150
151
152
153
154
155
156
  /** \brief Control normals when Decimating
   *
   * Binary and Cont. mode.
   *
   * The module tracks the normals while decimating
   * a normal cone consisting of all normals of the
   * faces collapsed together is computed and if
   * a collapse would increase the size of
   * the cone to a value greater than the given value
   * the collapse will be illegal.
   *
   * @param _ci Collapse info data
   * @return Half of the normal cones size (radius in radians)
   */
Jan Möbius's avatar
Jan Möbius committed
157
158
159
160
161
162
163
164
165
166
167
168
  float collapse_priority(const CollapseInfo& _ci) {
    // simulate collapse
    mesh_.set_point(_ci.v0, _ci.p1);


    typename Mesh::Scalar               max_angle(0.0);
    typename Mesh::ConstVertexFaceIter  vf_it(mesh_, _ci.v0);
    typename Mesh::FaceHandle           fh, fhl, fhr;

    if (_ci.v0vl.is_valid())  fhl = mesh_.face_handle(_ci.v0vl);
    if (_ci.vrv0.is_valid())  fhr = mesh_.face_handle(_ci.vrv0);

Jan Möbius's avatar
Jan Möbius committed
169
    for (; vf_it.is_valid(); ++vf_it) {
Jan Möbius's avatar
Jan Möbius committed
170
      fh = *vf_it;
Jan Möbius's avatar
Jan Möbius committed
171
172
173
174
175
176
177
178
179
      if (fh != _ci.fl && fh != _ci.fr) {
        NormalCone nc = mesh_.property(normal_cones_, fh);

        nc.merge(NormalCone(mesh_.calc_face_normal(fh)));
        if (fh == fhl) nc.merge(mesh_.property(normal_cones_, _ci.fl));
        if (fh == fhr) nc.merge(mesh_.property(normal_cones_, _ci.fr));

        if (nc.angle() > max_angle) {
          max_angle = nc.angle();
180
          if (max_angle > 0.5 * normal_deviation_)
Jan Möbius's avatar
Jan Möbius committed
181
182
183
184
185
186
187
188
189
190
            break;
        }
      }
    }


    // undo simulation changes
    mesh_.set_point(_ci.v0, _ci.p0);


191
    return (max_angle < 0.5 * normal_deviation_ ? max_angle : float( Base::ILLEGAL_COLLAPSE ));
Jan Möbius's avatar
Jan Möbius committed
192
193
  }

194
195
196
197
198
199
200
201
202
203
204
205
  /// set the percentage of normal deviation
  void set_error_tolerance_factor(double _factor) {
    if (_factor >= 0.0 && _factor <= 1.0) {
      // the smaller the factor, the smaller normal_deviation_ gets
      // thus creating a stricter constraint
      // division by error_tolerance_factor_ is for normalization
      Scalar normal_deviation = (normal_deviation_ * 180.0/M_PI) * _factor / this->error_tolerance_factor_;
      set_normal_deviation(normal_deviation);
      this->error_tolerance_factor_ = _factor;
    }
  }

Jan Möbius's avatar
Jan Möbius committed
206
207
208
209

  void  postprocess_collapse(const CollapseInfo& _ci) {
    // account for changed normals
    typename Mesh::VertexFaceIter vf_it(mesh_, _ci.v1);
Jan Möbius's avatar
Jan Möbius committed
210
    for (; vf_it.is_valid(); ++vf_it)
Jan Möbius's avatar
Jan Möbius committed
211
212
      mesh_.property(normal_cones_, *vf_it).
      merge(NormalCone(mesh_.normal(*vf_it)));
Jan Möbius's avatar
Jan Möbius committed
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231


    // normal cones of deleted triangles
    typename Mesh::FaceHandle fh;

    if (_ci.vlv1.is_valid()) {
      fh = mesh_.face_handle(mesh_.opposite_halfedge_handle(_ci.vlv1));
      if (fh.is_valid())
        mesh_.property(normal_cones_, fh).
        merge(mesh_.property(normal_cones_, _ci.fl));
    }

    if (_ci.v1vr.is_valid()) {
      fh = mesh_.face_handle(mesh_.opposite_halfedge_handle(_ci.v1vr));
      if (fh.is_valid())
        mesh_.property(normal_cones_, fh).
        merge(mesh_.property(normal_cones_, _ci.fr));
    }
  }
232

Jan Möbius's avatar
Jan Möbius committed
233
234
235
236
237
238
239
240
241
242
243
244
245
246


private:

  Mesh&                               mesh_;
  Scalar                              normal_deviation_;
  OpenMesh::FPropHandleT<NormalCone>  normal_cones_;
};


//=============================================================================
} // END_NS_DECIMATER
} // END_NS_OPENMESH
//=============================================================================
247
#endif // OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH defined
Jan Möbius's avatar
Jan Möbius committed
248
249
//=============================================================================