curvature.h

/*
 This Source Code Form is subject to the terms of the Mozilla Public
 License, v. 2.0. If a copy of the MPL was not distributed with this
 file, You can obtain one at http://mozilla.org/MPL/2.0/.
*/
 
#pragma once
 
#include "./defines.h"
 
#include <Eigen/Eigenvalues>
#include <Eigen/Core>
 
namespace Ponca
{
    template < class DataPoint, class _NFilter, int DiffType, typename T>
    class CurvatureEstimatorBase : public T
    {
    PONCA_FITTING_DECLARE_DEFAULT_TYPES
    PONCA_FITTING_DECLARE_MATRIX_TYPE
 
    protected:
        enum
        {
            PROVIDES_PRINCIPAL_CURVATURES
        };
 
    private:
        Scalar m_kmin {0},
        m_kmax {0};
        VectorType m_vmin {VectorType::Zero()},
        m_vmax {VectorType::Zero()};
 
        bool m_isValid {false};
 
        static_assert ( DataPoint::Dim == 3, "CurvatureEstimatorBase is only valid in 3D");
 
    public:
        PONCA_EXPLICIT_CAST_OPERATORS_DER(CurvatureEstimatorBase, curvatureEstimatorBase)
        PONCA_FITTING_DECLARE_INIT
 
        PONCA_MULTIARCH [[nodiscard]] inline bool isValid() const { return m_isValid; }
 
        /**************************************************************************/
        /* Use results                                                            */
        /**************************************************************************/
        PONCA_MULTIARCH [[nodiscard]] inline Scalar kmin() const { return m_kmin; }
 
        PONCA_MULTIARCH [[nodiscard]] inline Scalar kmax() const { return m_kmax; }
 
        PONCA_MULTIARCH [[nodiscard]] inline VectorType kminDirection() const { return m_vmin; }
 
        PONCA_MULTIARCH [[nodiscard]] inline VectorType kmaxDirection() const { return m_vmax; }
 
        PONCA_MULTIARCH [[nodiscard]] inline Scalar kMean() const { return (m_kmin + m_kmax)/Scalar(2);}
 
        PONCA_MULTIARCH [[nodiscard]] inline Scalar GaussianCurvature() const { return m_kmin * m_kmax;}
 
    protected:
        PONCA_MULTIARCH inline void setCurvatureValues(Scalar kmin, Scalar kmax, const VectorType& vmin, const VectorType& vmax);
    };
 
} //namespace Ponca
 
#include "curvature.hpp"