CircularArcs.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2024 ISciences, LLC
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Public Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 **********************************************************************/
 
#pragma once
 
#include <geos/export.h>
#include <geos/geom/Coordinate.h>
 
#include <array>
#include <optional>
 
namespace geos {
namespace geom {
class CoordinateSequence;
class Envelope;
}
}
 
namespace geos {
namespace algorithm {
 
class GEOS_DLL CircularArcs {
public:

    static geom::CoordinateXY getCenter(const geom::CoordinateXY& p0, const geom::CoordinateXY& p1,
                                        const geom::CoordinateXY& p2);
 
    static double getAngle(const geom::CoordinateXY& pt, const geom::CoordinateXY& center);
 
    static geom::CoordinateXY getDirectionPoint(const geom::CoordinateXY& center, double radius, double theta, bool isCCW);
 
    static double getMidpointAngle(double theta0, double theta2, bool isCCW);
 
    static geom::CoordinateXY getMidpoint(const geom::CoordinateXY& p0, const geom::CoordinateXY& p2, const geom::CoordinateXY& center, double radius, bool isCCW);

    static void expandEnvelope(geom::Envelope& e, const geom::CoordinateXY& p0, const geom::CoordinateXY& p1,
                               const geom::CoordinateXY& p2);

    static geom::CoordinateXY createPoint(const geom::CoordinateXY& center, double radius, double theta);

    static void interpolateZM(const geom::CoordinateSequence &seq, size_t i0, const geom::CoordinateXY &center, bool isCCW, geom::CoordinateXY &pt, double
                              &z, double &m);

    static int
    circleIntersectsLine(const geom::CoordinateXY& center, double r,
                         const geom::CoordinateXY& p0, const geom::CoordinateXY& p1,
                         geom::CoordinateXY& isect0, geom::CoordinateXY& isect1);

    static int
    circleIntersectsSegment(const geom::CoordinateXY& center, double r,
                            const geom::CoordinateXY& p0, const geom::CoordinateXY& p1,
                            geom::CoordinateXY& isect0, geom::CoordinateXY& isect1);

    static bool
    containsAngle(double t0, double t2, bool isCCW, double theta);

    static bool
    containsPointOnCircle(const geom::CoordinateXY& center,
                          const geom::CoordinateXY& p0, const geom::CoordinateXY& p2, bool isCCW,
                          const geom::CoordinateXY& q);

    static double
    distanceArcPoint(const geom::CoordinateXY& c, double r,
                     const geom::CoordinateXY& p0, const geom::CoordinateXY& p2, bool isCCW,
                     const geom::CoordinateXY& q);

    static double
    distanceArcSegment(const geom::CoordinateXY& c, double r,
                       const geom::CoordinateXY& p0, const geom::CoordinateXY& p2, bool isCCW,
                       const geom::CoordinateXY& q0, const geom::CoordinateXY& q1);

    static double
    distanceArcArc(const geom::CoordinateXY& ca, double ra,
                   const geom::CoordinateXY& a0, const geom::CoordinateXY& a2, bool aCCW,
                   const geom::CoordinateXY& cb, double rb,
                   const geom::CoordinateXY& b0, const geom::CoordinateXY& b2, bool bCCW);

    static geom::CoordinateXY
    closestPointArcPoint(const geom::CoordinateXY& c, double r,
                         const geom::CoordinateXY& p0, const geom::CoordinateXY& p2,
                         bool isCCW, const geom::CoordinateXY &q);

    static std::array<geom::CoordinateXY, 2>
    closestPointsArcSegment(const geom::CoordinateXY& c, double r,
                            const geom::CoordinateXY& p0, const geom::CoordinateXY& p2, bool isCCW, const geom::CoordinateXY& q0, const geom::CoordinateXY& q1);
 

    static std::array<geom::CoordinateXY, 2> closestPointsArcArc(const geom::CoordinateXY& ca,
        double ra,
        const geom::CoordinateXY& a0,
        const geom::CoordinateXY& a2,
        bool aCCW,
        const geom::CoordinateXY& cb,
        double rb,
        const geom::CoordinateXY& b0,
        const geom::CoordinateXY& b2,
        bool bCCW);

    static std::optional<geom::CoordinateXY>
    cocircularArcsIntersectionPoint(const geom::CoordinateXY& center, double radius,
                                    const geom::CoordinateXY& a0, const geom::CoordinateXY& a2, bool aCCW,
                                    const geom::CoordinateXY& b0, const geom::CoordinateXY& b2, bool bCCW);

    static std::optional<geom::CoordinateXY>
    arcIntersectionPoint(const geom::CoordinateXY& ca, double ra,
                         const geom::CoordinateXY& a0, const geom::CoordinateXY& a2, bool aCCW,
                         const geom::CoordinateXY& cb, double rb,
                         const geom::CoordinateXY& b0, const geom::CoordinateXY& b2, bool bCCW);
 
};
}
}