Function getPointsForRefineFeature

source/MRMesh/MRFeatureRefine.cpp:96–122

}

VertBitSet getPointsForRefineFeature( const PointCloud& pointCloud, const FeatureObject& feature, float distanceEps, float normalEps )
{
    VertBitSet detectedVerts;
    detectedVerts.resize( pointCloud.validPoints.size() + 1, false );

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

    BitSetParallelFor( pointCloud.validPoints, [&] ( VertId v )
    {
        const auto proj = feature.projectPoint( pointCloud.points[v] );
        const auto lenSq = ( proj.point - pointCloud.points[v] ).lengthSq();
        if ( lenSq >= distanceEpsSq )
            return;

        if ( pointCloud.hasNormals() && proj.normal )
        {
            const auto& worldNormal = *proj.normal;
            if ( std::abs( dot( worldNormal, pointCloud.normals[v] ) ) < cosNormalEps )
                return;
        }

        detectedVerts.set( v, true );
    } );

    return detectedVerts;
}

VertBitSet filterDisjointPoints( const Mesh& mesh, const VertBitSet& selectedPoints, const FaceBitSet& referenceFaces )
{