decimation_modules.h

Go to the documentation of this file.

#ifndef COB_DECIMATION_MODULES_H
#define COB_DECIMATION_MODULES_H

#include <OpenMesh/Tools/Decimater/ModBaseT.hh>

namespace OpenMesh
{
  namespace Decimater
  {
    template<typename MeshT>
      class ModNormalQuadricT : public OpenMesh::Decimater::ModBaseT<MeshT>
    {
    public:
      // Defines the types Self, Handle, Base, Mesh, and CollapseInfo
      // and the memberfunction name()
      DECIMATING_MODULE( ModNormalQuadricT, MeshT, NormalQuadric );

    public:
    ModNormalQuadricT( MeshT &_mesh )
      : Base(_mesh, false)
      {
        unset_max_err();
        Base::mesh().add_property( quadrics_ );
      }

      virtual ~ModNormalQuadricT()
      {
        Base::mesh().remove_property(quadrics_);
      }

    public: // inherited

      virtual void initialize(void);

      virtual float collapse_priority(const CollapseInfo& _ci)
      {
        typedef Geometry::QuadricT<double> Q;

        Q q = Base::mesh().property(quadrics_, _ci.v0);
        q += Base::mesh().property(quadrics_, _ci.v1);
        /*Vec3f v = (_ci.p1 + _ci.p0);
          v = v * 0.5;
          double err = q(v);*/
        double err = q(_ci.p1);
        //std::cout << v << " Err: " << err << std::endl;

        return float( (err < max_err_) ? err : Base::ILLEGAL_COLLAPSE );
      }

      virtual void postprocess_collapse(const CollapseInfo& _ci)
      {
        Base::mesh().property(quadrics_, _ci.v1) +=
          Base::mesh().property(quadrics_, _ci.v0);

        //Base::mesh().point(_ci.v1) += Base::mesh().point(_ci.v0);
        //Base::mesh().point(_ci.v1) *= 0.5;
      }

      void set_error_tolerance_factor(double _factor);

    public: // specific methods

      void set_max_err(double _err, bool _binary=true)
      {
        max_err_ = _err;
        Base::set_binary(_binary);
      }

      void set_normal_property(const VPropHandleT<Vec4f>& normals)
      {
        normals_ = normals;
      }

      void set_label_property(const VPropHandleT<int>& labels)
      {
        labels_ = labels;
      }

      void unset_max_err(void) { max_err_ = DBL_MAX; Base::set_binary(false); }
      double max_err() const { return max_err_; }

    private:
      double max_err_;

      VPropHandleT< Geometry::QuadricT<double> >  quadrics_;
      VPropHandleT< Vec4f > normals_;
      VPropHandleT< int > labels_;
    };


//=============================================================================
  } // END_NS_DECIMATER
} // END_NS_OPENMESH
//=============================================================================

#endif