example_line2d_create.cpp

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    double a;    double b;    double c;
    Pds::Octave::DAspect=true;
    Pds::Matrix X;
    
    double angle=Pds::Ra::Pi/4;
    Pds::Line2D Line1(angle,1);
    Line1.Print("Line1: ");
    Line1.GetLinearForm(a,b,c);
    std::cout<<"a: "<<a<<"\tb: "<<b<<"\tc: "<<c<<"\n\n";
    
    
    Pds::Line2D Line2(a,b,c);
    Line2.Print("Line2: ");
    
    X=Line1.RandSamplesAngle(32,angle-Pds::Ra::Pi/4,angle+Pds::Ra::Pi/4);
    Pds::Octave::Plot::PointsX2D(X,"testando.m","example_line2d_create_randsamples.png");
    
    std::cout<<std::endl;
    
    Pds::Line2D Line3=Line1.GetLineCovEigMethod(X);
    Pds::Matrix P=Line3.RandSamplesAngle(32,angle-Pds::Ra::Pi/4,angle+Pds::Ra::Pi/4);
    Line3.Print("Line3: ");
    
    Pds::Octave::Plot::PointsX2DY(  Pds::MergeVer({P,X}),
                                    Pds::MergeVer({Pds::Ones(P.Nlin(),1),Pds::Zeros(X.Nlin(),1)}),
                                    "testando.m","example_line2d_create_getline.png");
}