MRFreeFormDeformer.h

Go to the documentation of this file.

#pragma once
#include "MRMeshFwd.h"
#include "MRVector3.h"
#include "MRBox.h"
#include "MRVector.h"
#include "MRMeshEigen.h"
 
namespace MR
{

class FreeFormDeformer
{
public:
    MRMESH_API FreeFormDeformer( VertCoords& coords, const VertBitSet& valid );
    MRMESH_API FreeFormDeformer( Mesh& mesh, const VertBitSet* region = nullptr );
    MRMESH_API void init( const Vector3i& resolution = Vector3i::diagonal( 2 ), const Box3f& initialBox = Box3f() );
    MRMESH_API void setRefGridPointPosition( const Vector3i& coordOfPointInGrid, const Vector3f& newPos );
    MRMESH_API const Vector3f& getRefGridPointPosition( const Vector3i& coordOfPointInGrid ) const;
    MRMESH_API void apply() const;
    MRMESH_API Vector3f applySinglePoint( const Vector3f& point ) const;
    MRMESH_API int getIndex( const Vector3i& coordOfPointInGrid ) const;
    MRMESH_API Vector3i getCoord( int index ) const;
 
    const std::vector<Vector3f>& getAllRefGridPositions() const { return refPointsGrid_; }
    void setAllRefGridPositions( const std::vector<Vector3f>& refPoints ) { refPointsGrid_ = refPoints; }
    const Vector3i& getResolution() const { return resolution_; }
private:
    VertCoords& coords_;
    const VertBitSet& validPoints_;
    std::vector<Vector3f> refPointsGrid_;
    Box3f initialBox_;
    VertCoords normedCoords_;
    Vector3i resolution_;
 
    Vector3f applyToNormedPoint_( const Vector3f& normedPoint, std::vector<Vector3f>& xPlaneCache, std::vector<Vector3f>& yLineCache, std::vector<Vector3f>& tempPoints ) const;
};

class MRMESH_CLASS FreeFormBestFit
{
public:
    MRMESH_API FreeFormBestFit( const Box3d& box, const Vector3i& resolution = Vector3i::diagonal( 2 ) );

    MRMESH_API void addPair( const Vector3d& src, const Vector3d& tgt, double w = 1.0 );
    void addPair( const Vector3f& src, const Vector3f& tgt, float w = 1.0f ) { addPair( Vector3d( src ), Vector3d( tgt ), double( w ) ); }

    MRMESH_API void addOther( const FreeFormBestFit& other );

    void setStabilizer( double stabilizer ) { stabilizer_ = stabilizer; }
    double getStabilizer() const { return stabilizer_; }

    [[nodiscard]] MRMESH_API std::vector<Vector3f> findBestDeformationReferenceGrid();
private:
    Box3d box_;
    Vector3i resolution_;
    size_t resXY_{ 0 };
    size_t size_{ 0 };
    double sumWeight_{ 0.0 };
 
    Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> accumA_;
    Eigen::Matrix<double, Eigen::Dynamic, 3> accumB_;
 
    std::vector<int> pascalLineX_;
    std::vector<int> pascalLineY_;
    std::vector<int> pascalLineZ_;
 
    Vector3d reverseDiagonal_;
 
    double stabilizer_{ 0.1 };
    void stabilize_();
};

MRMESH_API std::vector<Vector3f> makeFreeFormOriginGrid( const Box3f& box, const Vector3i& resolution );

MRMESH_API std::vector<Vector3f> findBestFreeformDeformation( const Box3f& box, const std::vector<Vector3f>& source, const std::vector<Vector3f>& target,
                                                              const Vector3i& resolution = Vector3i::diagonal( 2 ), const AffineXf3f* samplesToBox = nullptr );
 
}