Function meshDenoiseViaNormals

source/MRMesh/MRNormalDenoising.cpp:201–260

}

Expected<void> meshDenoiseViaNormals( Mesh & mesh, const DenoiseViaNormalsSettings & settings )
{
    MR_TIMER;
    if ( settings.normalIters <= 0 || settings.pointIters <= 0 )
    {
        assert( false );
        return unexpected( "Bad parameters" );
    }

    if ( !reportProgress( settings.cb, 0.0f ) )
        return unexpectedOperationCanceled();

    auto fnormals0 = computePerFaceNormals( mesh );
    Vector<float, UndirectedEdgeId> v( mesh.topology.undirectedEdgeSize(), 1 );

    if ( !reportProgress( settings.cb, 0.05f ) )
        return unexpectedOperationCanceled();

    auto sp = subprogress( settings.cb, 0.05f, 0.95f );
    FaceNormals fnormals;
    for ( int i = 0; i < settings.normalIters; ++i )
    {
        fnormals = fnormals0;
        denoiseNormals( mesh, fnormals, v, settings.gamma );
        if ( !reportProgress( sp, float( 2 * i ) / ( 2 * settings.normalIters ) ) )
            return unexpectedOperationCanceled();

        if ( settings.fastIndicatorComputation )
            updateIndicatorFast( mesh.topology, v, fnormals, settings.beta, settings.gamma );
        else
            updateIndicator( mesh, v, fnormals, settings.beta, settings.gamma );
        if ( !reportProgress( sp, float( 2 * i + 1 ) / ( 2 * settings.normalIters ) ) )
            return unexpectedOperationCanceled();
    }

    if ( settings.outCreases )
    {
        settings.outCreases->clear();
        settings.outCreases->resize( mesh.topology.undirectedEdgeSize() );
        BitSetParallelForAll( *settings.outCreases, [&]( UndirectedEdgeId ue )
        {
            if ( v[ue] < 0.5f )
                settings.outCreases->set( ue );
        } );
    }

    if ( !reportProgress( settings.cb, 0.95f ) )
        return unexpectedOperationCanceled();

    const auto guide = mesh.points;
    NormalsToPoints n2p;
    n2p.prepare( mesh.topology, settings.guideWeight );
    auto maxInitialDistSq = settings.limitNearInitial ? sqr( settings.maxInitialDist )
        : std::numeric_limits<float>::infinity();
    for ( int i = 0; i < settings.pointIters; ++i )
        n2p.run( guide, fnormals, mesh.points, maxInitialDistSq );

    reportProgress( settings.cb, 1.0f );