example_funcarraylabeling.cpp

 
#include <Pds/Ra>
#include <Pds/Sp>
 
int main(void)
{
    Pds::Ra::Randomize();
    
    unsigned int L=512;
    Pds::Array<unsigned char> A(L,L);
    
    // Cargando datos de prueba
    Pds::Image::DataSet::LoadDataCircles(A,3,1);
    Pds::Matrix(A).Scale(0,255).ExportBmpFile(Pds::Colormap::Gray,"A-binary-circle.bmp");
 
    // Usando Pds::Image::IdMapUsingCC4() para obter la matriz ID
    Pds::Array<unsigned int> ID;
    
    Pds::Ra::Tic();
    ID=Pds::Image::IdMapUsingCC4(A);
    Pds::Ra::Toc();
 
    Pds::Matrix(ID).Scale(0,255).ExportBmpFile(Pds::Colormap::Steps,"A-id-circle.bmp");
    
    // Obtendo datos de los grupos en ID
    std::vector<std::vector<Pds::Position>> Data;
    Pds::Ra::Tic();
    Data=Pds::Image::VectorGroupFromIdMap(ID);
    Pds::Ra::Toc();
 
    for(unsigned int i=0;i<Data.size();i++)
    {
        double  Per=Pds::Image::PerimeterOfConvexGroup(Data[i]);
        unsigned int Area=Data[i].size();
        
        std::cout<<"Perimeter["<<i<<"]:\t"<<Per<<"\t";
        std::cout<<"     Area["<<i<<"]:\t"<<Area<<"\t";
        std::cout<<"    Gamma["<<i<<"]:\t"<<(Per/2.0)/sqrt(Area*Pds::Ra::Pi)<<"\n";
    }
    return 0;
}
 
/*
 
*/