MRSurfacePath.h

Go to the documentation of this file.

#pragma once
 
#include "MRExpected.h"
#include "MRMeshTriPoint.h"
#include "MREnums.h"
#include <vector>
#include <string>
 
namespace MR
{


struct GeodesicPath
{
    MeshTriPoint start; 
    SurfacePath mids;
    MeshTriPoint end;   
 
    [[nodiscard]] size_t numVertices() const { return start.valid() + mids.size() + end.valid(); }
 
    [[nodiscard]] MeshTriPoint operator[]( size_t i ) const
    {
        if ( start.valid() )
        {
            if ( i == 0 )
                return start;
            --i;
        }
        if ( i < mids.size() )
            return mids[i];
        assert( end.valid() );
        assert( i == mids.size() );
        return end;
    }
};
 
enum class PathError
{
    StartEndNotConnected, 
    InternalError         
};
 
inline std::string toString(PathError error)
{
    switch (error)
    {
        case PathError::StartEndNotConnected:   return "No path can be found from start to end, because they are not from the same connected component";
        case PathError::InternalError:          return "Report to developers for further investigations";
        default:                                return "Unknown error. Please, report to developers for further investigations";
    }
}

MRMESH_API Expected<SurfacePath, PathError> computeSurfacePath( const MeshPart & mp, 
    const MeshTriPoint & start, const MeshTriPoint & end, 
    int maxGeodesicIters = 5, 
    const VertBitSet* vertRegion = nullptr, VertScalars * outSurfaceDistances = nullptr );

MRMESH_API Expected<SurfacePath, PathError> computeGeodesicPath( const Mesh & mesh,
    const MeshTriPoint & start, const MeshTriPoint & end, GeodesicPathApprox atype = GeodesicPathApprox::FastMarching,
    int maxGeodesicIters = 100 ); 

MRMESH_API Expected<SurfacePath, PathError> computeGeodesicPathApprox( const Mesh & mesh,
    const MeshTriPoint & start, const MeshTriPoint & end, GeodesicPathApprox atype );

MRMESH_API Expected<SurfacePath, PathError> computeFastMarchingPath( const MeshPart & mp, 
    const MeshTriPoint & start, const MeshTriPoint & end, const VertBitSet* vertRegion = nullptr,
    VertScalars * outSurfaceDistances = nullptr );
 
struct ComputeSteepestDescentPathSettings
{
    MeshTriPoint end;
    VertId * outVertexReached = nullptr;
    EdgePoint * outBdReached = nullptr;
};

[[nodiscard]] MRMESH_API SurfacePath computeSteepestDescentPath( const MeshPart & mp, const VertScalars & field,
    const MeshTriPoint & start, const ComputeSteepestDescentPathSettings & settings = {} );

MRMESH_API void computeSteepestDescentPath( const MeshPart & mp, const VertScalars & field,
    const MeshTriPoint & start, SurfacePath * outPath, const ComputeSteepestDescentPathSettings & settings = {} );

[[nodiscard]] MRMESH_API MeshEdgePoint findSteepestDescentPoint( const MeshPart & mp, const VertScalars & field, VertId v );

[[nodiscard]] MRMESH_API MeshEdgePoint findSteepestDescentPoint( const MeshPart & mp, const VertScalars & field, const MeshEdgePoint & ep );

[[nodiscard]] MRMESH_API MeshEdgePoint findSteepestDescentPoint( const MeshPart & mp, const VertScalars & field, const MeshTriPoint & tp );
 
enum class ExtremeEdgeType
{
    Ridge, 
    Gorge  
};

[[nodiscard]] MRMESH_API UndirectedEdgeBitSet findExtremeEdges( const Mesh & mesh, const VertScalars & field, ExtremeEdgeType type );

[[nodiscard]] MRMESH_API HashMap<VertId, VertId> computeClosestSurfacePathTargets( const Mesh & mesh,
    const VertBitSet & starts, const VertBitSet & ends, const VertBitSet * vertRegion = nullptr,
    VertScalars * outSurfaceDistances = nullptr );

[[nodiscard]] MRMESH_API SurfacePaths getSurfacePathsViaVertices( const Mesh & mesh, const VertBitSet & vs );

[[nodiscard]] MRMESH_API float surfacePathLength( const Mesh& mesh, const SurfacePath& surfacePath );

[[nodiscard]] MRMESH_API float geodesicPathLength( const Mesh& mesh, const GeodesicPath& path );

[[nodiscard]] MRMESH_API Contour3f surfacePathToContour3f( const Mesh & mesh, const SurfacePath & line );
[[nodiscard]] MRMESH_API Contours3f surfacePathsToContours3f( const Mesh & mesh, const SurfacePaths & lines );

[[nodiscard]] MRMESH_API Contour3f geodesicPathToContour3f( const Mesh& mesh, const GeodesicPath& path );

 
}