Function getComplexFillMetric

source/MRMesh/MRMeshMetrics.cpp:207–252

}

FillHoleMetric getComplexFillMetric( const Mesh& mesh, EdgeId e0 )
{
    float maxEdgeLengthSq = 0.0f;
    for ( auto e : leftRing( mesh.topology, e0 ) )
        maxEdgeLengthSq = std::max( maxEdgeLengthSq, mesh.edgeLengthSq( e ) );

    assert( maxEdgeLengthSq > 0.0f );

    float reverseCharacteristicTriArea{ 0.0f };
    if ( maxEdgeLengthSq <= 0.0f )
        reverseCharacteristicTriArea = 1.0f;
    else
        reverseCharacteristicTriArea = 1.0f / maxEdgeLengthSq;

    FillHoleMetric metric;
    metric.triangleMetric = [&mesh, reverseCharacteristicTriArea] ( VertId a, VertId b, VertId c )
    {
        double aspectRatio = triangleAspectRatio( mesh.points[a], mesh.points[b], mesh.points[c] );
        if ( aspectRatio > BadTriangulationMetric )
            return BadTriangulationMetric;

        double normedArea = TriangleAreaModifier * cross( mesh.points[b] - mesh.points[a], mesh.points[c] - mesh.points[a] ).length() * reverseCharacteristicTriArea;

        return aspectRatio + normedArea;
    };
    metric.edgeMetric = [&mesh] ( VertId a, VertId b, VertId l, VertId r )
    {
        auto abVec = mesh.points[b] - mesh.points[a];
        auto bcVec = -abVec;
        auto normA = cross( mesh.points[r] - mesh.points[b], bcVec );
        auto normC = cross( mesh.points[l] - mesh.points[a], abVec );

        auto s_Abc_double = normA.length();
        auto s_abC_double = normC.length();
        auto denom = s_Abc_double * s_abC_double;
        if ( denom == 0.0f )
            return BadTriangulationMetric;
        auto cosAC = dot( normA, normC ) / ( s_Abc_double * s_abC_double );

        if ( cosAC <= -1.0f )
            return BadTriangulationMetric;

        return double( sqr( sqr( ( 1.0f - cosAC ) / ( 1.0f + cosAC ) ) ) );
    };
    return metric;
}

FillHoleMetric getParallelPlaneFillMetric( const Mesh& mesh, EdgeId e0, const Plane3f* plane /*= nullptr */ )
{