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"
 
namespace Ponca
{
    namespace internal
    {
        template < class DataPoint, class _NFilter, 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_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);
        };
    }
 
    template < class DataPoint, class _NFilter, typename T>
    class CurvatureEstimator : public internal::CurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
    public:
        PONCA_EXPLICIT_CAST_OPERATORS(CurvatureEstimator, curvatureEstimator)
    };
 
    template < class DataPoint, class _NFilter, int DiffType, typename T>
    class CurvatureEstimatorDer : public internal::CurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
    public:
        PONCA_EXPLICIT_CAST_OPERATORS_DER(CurvatureEstimatorDer, curvatureEstimator)
    };
 
} //namespace Ponca
 
#include "curvature.hpp"