Function pointGridSampling

source/MRMesh/MRGridSampling.cpp:166–201

}

std::optional<VertBitSet> pointGridSampling( const PointCloudPart& pcp, float voxelSize, const ProgressCallback & cb )
{
    MR_TIMER;
    const auto & ps = pcp.cloud.getVertIds( pcp.region );
    if ( voxelSize <= 0 )
        return ps;

    const auto bbox = pcp.cloud.computeBoundingBox( pcp.region );
    const auto bboxSz = bbox.max - bbox.min;
    constexpr float maxVoxelsInOneDim = 1 << 10;
    const Vector3i dims
    {
        ( int )std::clamp( std::ceil( bboxSz.x / voxelSize ),1.0f, maxVoxelsInOneDim ),
        ( int )std::clamp( std::ceil( bboxSz.y / voxelSize ),1.0f, maxVoxelsInOneDim ),
        ( int )std::clamp( std::ceil( bboxSz.z / voxelSize ),1.0f, maxVoxelsInOneDim )
    };

    Grid grid( bbox, dims );
    if ( cb && !cb( 0.1f ) )
        return {};

    int counter = 0;
    int size = int( ps.count() );
    for ( auto v : ps )
    {
        grid.addVertex( pcp.cloud.points[v], v );
        if ( !reportProgress( cb, [&]{ return 0.1f + 0.8f * float( counter ) / float( size ); }, counter++, 1024 ) )
            return {};
    }

    const auto res = grid.getSamples();
    if ( cb && !cb( 1.0f ) )
        return {};

    return res;
}

std::optional<std::vector<ObjVertId>> multiModelGridSampling( const Vector<ModelPointsData, ObjId>& models, float voxelSize, const ProgressCallback& cb )
{