example_matrix_kernel_gaussian.cpp

Programa para el testeo de las funciones.

Programa para el testeo de las funciones.Para compilar o código example_matrix_kernel_gaussian.cpp:

g++ -static -o example_matrix_kernel_gaussian example_matrix_kernel_gaussian.cpp -lpdsramm 

Para executar o programa:

./example_matrix_kernel_gaussian

Retornando por consola:

G.Size():       7       7
Gd.Size():      7       7
Sigma:4.2

imagesc_gaussian.png

imagesc_gaussiandiff.png

imagesc_gaussian_sigma_fit.png

Código example_matrix_kernel_gaussian.cpp:

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Octave::XLabel=" ";
    Pds::Octave::YLabel=" ";
    
    Pds::Matrix G;
    Pds::Matrix Gd;
    
    unsigned int K=3;
    double Sigma=1.0;
    
    G=Pds::Kernel2D::Gaussian( K,Sigma);
    G.Size().Print("G.Size():\t");
    Pds::Octave::Plot::ImageSc( G,"testandog.m","imagesc_gaussian.png");
    
    Gd=Pds::Kernel2D::GaussianDiff( K,Sigma,1.2*Sigma);
    Gd.Size().Print("Gd.Size():\t");
    Pds::Octave::Plot::ImageSc( Gd,"testandogq.m","imagesc_gaussiandiff.png");
    
    K=10;
    Sigma=Pds::Kernel2D::GaussianSigmaFit(7,K);
    std::cout<<"Sigma:"<<Sigma<<"\n";
    G=Pds::Kernel2D::Gaussian(K,Sigma);
    Gd=G-G.Mean();
    Pds::Octave::Plot::ImageSc(Gd.Geq(0),"testandogsf.m","imagesc_gaussian_sigma_fit.png");
    
    
    return 0;
}