geoellipse.h

Go to the documentation of this file.

#pragma once

#include "cad/const.h"
#include "geocoordinate.h"
#include "cad/math/lcmath.h"
#include <memory>
#include "geobase.h"
#include <vector>
#include "cad/math/equation.h"

namespace lc {

    namespace geo {
        class Ellipse : public Base, virtual public Visitable {
            public:
                Ellipse(const Coordinate& center, const Coordinate& majorP, double minorRadius, double startAngle, double endAngle, bool reversed = false);
                const Coordinate center() const;

                const Coordinate majorP() const;

                double minorRadius() const;

                double startAngle() const;

                double endAngle() const;

                Coordinate getPoint(const double& angle) const;
                Coordinate startPoint() const;
                Coordinate endPoint() const;
                std::vector<Coordinate> findPotentialNearestPoints(const Coordinate& coord) const;
                Coordinate nearestPointOnPath(const Coordinate& coord) const;
                Coordinate nearestPointOnEntity(const Coordinate& coord) const;
                bool isArc() const;

                bool isReversed() const;

                double ratio() const;

                double majorRadius() const;

                double getAngle() const;

                double getEllipseAngle(const Coordinate& coord) const;

                bool isAngleBetween(double angle) const {
                    if (!isArc())
                        return true;
                    return Math::isAngleBetween(angle, _startAngle, _endAngle, !_isReversed);
                }

                Ellipse georotate(const Coordinate& center, const double rotation_angle) const;

                Ellipse geoscale(const Coordinate& center, const Coordinate &factor) const;

                const maths::Equation equation() const {
                    auto ce0 = _majorP.squared();
                    auto ce2 = this->ratio() * this->ratio() * ce0;

                    if (ce0 < LCTOLERANCE * LCTOLERANCE || ce2 < LCTOLERANCE * LCTOLERANCE) {
                        return maths::Equation();
                    }

                    return maths::Equation(1. / ce0, 0., 1. / ce2, 0., 0., -1.).rotate(getAngle()).move(_center);
                }

                virtual void accept(GeoEntityVisitor &v) const override { v.visit(*this); }

            private:
                friend std::ostream& operator<<(std::ostream& os, const Ellipse& e) {
                    os << "Ellipse(center=" << e._center << " majorP=" << e._majorP << " minorRadius=" << e._minorRadius << " startAngle=" << e._startAngle << " endAngle=" << e._endAngle << " isReversed=" << e._isReversed << ")";
                    return os;
                }
            private:
                const Coordinate _center;
                const Coordinate _majorP;
                const double _minorRadius;
                const double  _startAngle;
                const double _endAngle;
                const bool _isReversed;
        };
    }
}

// GEOELLIPSE_H