ModNormalDeviationT.hh 9.99 KB
Newer Older
Jan Möbius's avatar
Jan Möbius committed
1
/* ========================================================================= *
Jan Möbius's avatar
Jan Möbius committed
2 3
 *                                                                           *
 *                               OpenMesh                                    *
Jan Möbius's avatar
Jan Möbius committed
4
 *           Copyright (c) 2001-2015, RWTH-Aachen University                 *
Jan Möbius's avatar
Typo  
Jan Möbius committed
5
 *           Department of Computer Graphics and Multimedia                  *
Jan Möbius's avatar
Jan Möbius committed
6 7
 *                          All rights reserved.                             *
 *                            www.openmesh.org                               *
Jan Möbius's avatar
Jan Möbius committed
8 9
 *                                                                           *
 *---------------------------------------------------------------------------*
Jan Möbius's avatar
Jan Möbius committed
10 11
 * This file is part of OpenMesh.                                            *
 *---------------------------------------------------------------------------*
Jan Möbius's avatar
Jan Möbius committed
12
 *                                                                           *
Jan Möbius's avatar
Jan Möbius committed
13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
 * Redistribution and use in source and binary forms, with or without        *
 * modification, are permitted provided that the following conditions        *
 * are met:                                                                  *
 *                                                                           *
 * 1. Redistributions of source code must retain the above copyright notice, *
 *    this list of conditions and the following disclaimer.                  *
 *                                                                           *
 * 2. Redistributions in binary form must reproduce the above copyright      *
 *    notice, this list of conditions and the following disclaimer in the    *
 *    documentation and/or other materials provided with the distribution.   *
 *                                                                           *
 * 3. Neither the name of the copyright holder nor the names of its          *
 *    contributors may be used to endorse or promote products derived from   *
 *    this software without specific prior written permission.               *
 *                                                                           *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS       *
 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED *
 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A           *
 * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER *
 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,  *
 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,       *
 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR        *
 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF    *
 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING      *
 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS        *
 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.              *
Jan Möbius's avatar
Jan Möbius committed
39 40
 *                                                                           *
 * ========================================================================= */
Jan Möbius's avatar
Jan Möbius committed
41 42 43

/*===========================================================================*\
 *                                                                           *
Jan Möbius's avatar
Jan Möbius committed
44 45
 *   $Revision$                                                        *
 *   $Date$                   *
Jan Möbius's avatar
Jan Möbius committed
46 47 48 49 50 51 52 53 54 55 56 57 58
 *                                                                           *
\*===========================================================================*/

/** \file ModNormalDeviationT.hh
 */

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


59 60
#ifndef OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH
#define OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH
Jan Möbius's avatar
Jan Möbius committed
61 62 63 64 65 66


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

#include <OpenMesh/Tools/Decimater/ModBaseT.hh>
#include <OpenMesh/Core/Utils/Property.hh>
67
#include <OpenMesh/Core/Geometry/NormalConeT.hh>
Jan Möbius's avatar
Jan Möbius committed
68 69 70 71 72 73 74 75 76 77 78


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

namespace OpenMesh  {
namespace Decimater {


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


79 80 81 82 83 84 85 86 87 88 89 90 91 92
/** \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
  */
93 94
template <class MeshT>
class ModNormalDeviationT : public ModBaseT< MeshT >
Jan Möbius's avatar
Jan Möbius committed
95 96
{
public:
97

98
  DECIMATING_MODULE( ModNormalDeviationT, MeshT, NormalDeviation );
Jan Möbius's avatar
Jan Möbius committed
99 100 101 102 103 104 105

  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;
106
  typedef NormalConeT<Scalar>                       NormalCone;
Jan Möbius's avatar
Jan Möbius committed
107 108 109 110 111 112



public:

  /// Constructor
113 114
  ModNormalDeviationT(MeshT& _mesh, float _max_dev = 180.0)
  : Base(_mesh, true), mesh_(Base::mesh())
Jan Möbius's avatar
Jan Möbius committed
115 116 117
  {
    set_normal_deviation(_max_dev);
    mesh_.add_property(normal_cones_);
118 119 120 121 122 123

    const bool mesh_has_normals = _mesh.has_face_normals();
    _mesh.request_face_normals();

    if (!mesh_has_normals)
    {
Jan Möbius's avatar
Jan Möbius committed
124
      omerr() << "Mesh has no face normals. Compute them automatically." << std::endl;
125 126
      _mesh.update_face_normals();
    }
Jan Möbius's avatar
Jan Möbius committed
127
  }
128

Jan Möbius's avatar
Jan Möbius committed
129 130 131

  /// Destructor
  ~ModNormalDeviationT() {
Matthias Möller's avatar
Matthias Möller committed
132
    mesh_.remove_property(normal_cones_);
133
    mesh_.release_face_normals();
Jan Möbius's avatar
Jan Möbius committed
134 135 136
  }


137
  /// Get normal deviation ( 0 .. 360 )
Jan Möbius's avatar
Jan Möbius committed
138 139 140 141
  Scalar normal_deviation() const {
    return normal_deviation_ / M_PI * 180.0;
  }

142
  /// Set normal deviation ( 0 .. 360 )
Jan Möbius's avatar
Jan Möbius committed
143
  void set_normal_deviation(Scalar _s) {
Jan Möbius's avatar
Jan Möbius committed
144
    normal_deviation_ = _s / static_cast<Scalar>(180.0) * static_cast<Scalar>(M_PI);
Jan Möbius's avatar
Jan Möbius committed
145 146 147 148 149 150 151 152 153
  }


  /// 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(),
154 155
        f_end = mesh_.faces_end();

Jan Möbius's avatar
Jan Möbius committed
156
    for (; f_it != f_end; ++f_it)
Jan Möbius's avatar
Jan Möbius committed
157
      mesh_.property(normal_cones_, *f_it) = NormalCone(mesh_.normal(*f_it));
Jan Möbius's avatar
Jan Möbius committed
158 159
  }

160 161 162 163 164 165 166 167 168 169 170 171 172 173
  /** \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
174 175 176 177 178 179 180 181 182 183 184 185
  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
186
    for (; vf_it.is_valid(); ++vf_it) {
Jan Möbius's avatar
Jan Möbius committed
187
      fh = *vf_it;
Jan Möbius's avatar
Jan Möbius committed
188 189 190 191 192 193 194 195 196
      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();
197
          if (max_angle > 0.5 * normal_deviation_)
Jan Möbius's avatar
Jan Möbius committed
198 199 200 201 202 203 204 205 206 207
            break;
        }
      }
    }


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


208
    return (max_angle < 0.5 * normal_deviation_ ? max_angle : float( Base::ILLEGAL_COLLAPSE ));
Jan Möbius's avatar
Jan Möbius committed
209 210
  }

211 212 213 214 215 216
  /// 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
Jan Möbius's avatar
Jan Möbius committed
217 218
      Scalar normal_deviation = normal_deviation_ *  static_cast<Scalar>(  180.0 / M_PI * _factor / this->error_tolerance_factor_);

219 220 221 222 223
      set_normal_deviation(normal_deviation);
      this->error_tolerance_factor_ = _factor;
    }
  }

Jan Möbius's avatar
Jan Möbius committed
224 225 226 227

  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
228
    for (; vf_it.is_valid(); ++vf_it)
Jan Möbius's avatar
Jan Möbius committed
229 230
      mesh_.property(normal_cones_, *vf_it).
      merge(NormalCone(mesh_.normal(*vf_it)));
Jan Möbius's avatar
Jan Möbius committed
231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249


    // 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));
    }
  }
250

Jan Möbius's avatar
Jan Möbius committed
251 252 253 254 255 256 257 258 259 260 261 262 263 264


private:

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


//=============================================================================
} // END_NS_DECIMATER
} // END_NS_OPENMESH
//=============================================================================
265
#endif // OPENMESH_DECIMATER_MODNORMALDEVIATIONT_HH defined
Jan Möbius's avatar
Jan Möbius committed
266 267
//=============================================================================