Function getVertsForRefineFeature

source/MRMesh/MRFeatureRefine.cpp:45–94

}

VertBitSet getVertsForRefineFeature( const Mesh& mesh, const FeatureObject& feature, float distanceEps, float normalEps )
{
    VertBitSet detectedVerts;
    detectedVerts.resize( mesh.topology.lastValidVert() + 1, false );

    const auto distanceEpsSq = distanceEps * distanceEps;
    const auto cosNormalEps = std::cos( normalEps / 180.0f * PI_F );

    BitSetParallelForAll( mesh.topology.getValidVerts(), [&] ( VertId v )
    {
        const auto proj = feature.projectPoint( mesh.points[v] );
        const auto lenSq = ( proj.point - mesh.points[v] ).lengthSq();
        if ( lenSq >= distanceEpsSq )
            return;

        if ( !proj.normal )
        {
            detectedVerts.set( v, true );  // normals are not support for current feature type;
            return;
        }

        const auto& worldNormal = *proj.normal;

        // check vert normal.
        if ( std::abs( dot( worldNormal, mesh.normal( v ) ) ) >= cosNormalEps )
        {
            detectedVerts.set( v, true );
            return;
        }

        // Also check the face normals surrounding the current vertex.
        // If at least one face has a correct normal, add the vertex to the list.
        // It is required for working with meshes created by CAD (like STEP files or create primitives),
        // where almost all verts located at the sharp edged and have normals which is not correlated with surrounded faces.
        for ( const auto e : orgRing( mesh.topology, v ) )
        {
            const auto f = mesh.topology.left( e );
            if ( !f.valid() )
                continue;

            if ( std::abs( dot( worldNormal, mesh.normal( f ) ) ) >= cosNormalEps )
            {
                detectedVerts.set( v, true );
                return;
            }
        }
    } );

    return detectedVerts;
}

VertBitSet getPointsForRefineFeature( const PointCloud& pointCloud, const FeatureObject& feature, float distanceEps, float normalEps )
{