test_funcclustering_kmeans2.cpp

 
#include <cmath>
#include <Pds/Ra>
#include <Pds/Ml>
    
int main(void)
{
    std::vector<Pds::Matrix> Block;
 
    Block=Pds::Matrix::ImportBmpFile("imagen1.bmp");
    unsigned int Nlin=Block[0].Nlin();
    unsigned int Ncol=Block[0].Ncol();
    Pds::Matrix X=Pds::Matrix::GetSamples(Block);
    Block[0].MakeEmpty();    Block[1].MakeEmpty();    Block[2].MakeEmpty();
 
    unsigned int K=256;//128;
    Pds::Matrix Colors;
    Pds::Array<unsigned int> MAPID;
    
    Pds::IterationConf Conf;
    Conf.SetMinError(1.0);
 
    Colors=Pds::Clustering::Kmeans(Conf,X,K,MAPID);
    MAPID.Reshape(Nlin,Ncol);
    MAPID.ExportBmpFile(Pds::Array<unsigned char>(Colors),"imagen1-kmeans.bmp");
       
    return 0;
}