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;
}