example_ellipse_getellipse.cpp

 
#include <cmath>
#include <Pds/Ra>
    
    
    
int main(void)
{
    Pds::Ra::Randomize();
    Pds::Octave::DAspect=true;
    Pds::Octave::MarkerSize=1;
    unsigned int L=1000;
    
    Pds::Ellipse E(3,1.5,5,3,-Pds::Ra::Pi/6);
    std::cout<<"  a: "<<E.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<E.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<E.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<" Xc: "<<E.GetCentroid().T();
    
    Pds::Matrix Xeq=E.RandSamples(L);
    
    Pds::Octave::Plot::PointsX2D(Xeq,"testando.m","example_ellipse_getellipse_elipse1.png");
    
    std::cout<<std::endl;
 
    Pds::Ellipse E2=Pds::Ellipse::GetEllipse(Xeq);
    
    std::cout<<"  a: "<<E2.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<E2.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<E2.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<" Xc: "<<E2.GetCentroid().T();
    
    Pds::Matrix P=E2.RandPerimeter(L);
    Pds::Octave::Plot::PointsX2DY(  Pds::MergeVer({P,Xeq}),
                                    Pds::MergeVer({Pds::Ones(P.Nlin(),1),Pds::Zeros(Xeq.Nlin(),1)}),
                                    "testando.m","example_ellipse_getellipse_elipse2.png");
 
    return 0;
}