137 lines
4.0 KiB
C++
137 lines
4.0 KiB
C++
// Boost.Geometry
|
|
|
|
// Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
|
|
|
|
// This file was modified by Oracle on 2014-2021.
|
|
// Modifications copyright (c) 2014-2021 Oracle and/or its affiliates.
|
|
|
|
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
|
|
|
|
// Use, modification and distribution is subject to the Boost Software License,
|
|
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
|
|
// http://www.boost.org/LICENSE_1_0.txt)
|
|
|
|
#ifndef BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP
|
|
#define BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP
|
|
|
|
#include <boost/geometry/core/cs.hpp>
|
|
#include <boost/geometry/core/access.hpp>
|
|
#include <boost/geometry/core/coordinate_promotion.hpp>
|
|
#include <boost/geometry/core/radian_access.hpp>
|
|
#include <boost/geometry/core/radius.hpp>
|
|
|
|
#include <boost/geometry/formulas/spherical.hpp>
|
|
|
|
#include <boost/geometry/srs/spheroid.hpp>
|
|
|
|
//#include <boost/geometry/strategies/concepts/side_concept.hpp>
|
|
#include <boost/geometry/strategies/geographic/disjoint_segment_box.hpp>
|
|
#include <boost/geometry/strategies/geographic/parameters.hpp>
|
|
#include <boost/geometry/strategies/side.hpp>
|
|
#include <boost/geometry/strategies/spherical/point_in_point.hpp>
|
|
|
|
#include <boost/geometry/strategy/geographic/envelope.hpp>
|
|
|
|
#include <boost/geometry/util/math.hpp>
|
|
#include <boost/geometry/util/select_calculation_type.hpp>
|
|
|
|
namespace boost { namespace geometry
|
|
{
|
|
|
|
|
|
namespace strategy { namespace side
|
|
{
|
|
|
|
|
|
/*!
|
|
\brief Check at which side of a segment a point lies
|
|
left of segment (> 0), right of segment (< 0), on segment (0)
|
|
\ingroup strategies
|
|
\tparam FormulaPolicy Geodesic solution formula policy.
|
|
\tparam Spheroid Reference model of coordinate system.
|
|
\tparam CalculationType \tparam_calculation
|
|
|
|
\qbk{
|
|
[heading See also]
|
|
[link geometry.reference.srs.srs_spheroid srs::spheroid]
|
|
}
|
|
*/
|
|
template
|
|
<
|
|
typename FormulaPolicy = strategy::andoyer,
|
|
typename Spheroid = srs::spheroid<double>,
|
|
typename CalculationType = void
|
|
>
|
|
class geographic
|
|
{
|
|
public:
|
|
typedef geographic_tag cs_tag;
|
|
|
|
geographic() = default;
|
|
|
|
explicit geographic(Spheroid const& model)
|
|
: m_model(model)
|
|
{}
|
|
|
|
template <typename P1, typename P2, typename P>
|
|
inline int apply(P1 const& p1, P2 const& p2, P const& p) const
|
|
{
|
|
typedef strategy::within::spherical_point_point equals_point_point_strategy_type;
|
|
if (equals_point_point_strategy_type::apply(p, p1)
|
|
|| equals_point_point_strategy_type::apply(p, p2)
|
|
|| equals_point_point_strategy_type::apply(p1, p2))
|
|
{
|
|
return 0;
|
|
}
|
|
|
|
typedef typename promote_floating_point
|
|
<
|
|
typename select_calculation_type_alt
|
|
<
|
|
CalculationType,
|
|
P1, P2, P
|
|
>::type
|
|
>::type calc_t;
|
|
|
|
typedef typename FormulaPolicy::template inverse
|
|
<calc_t, false, true, false, false, false> inverse_formula;
|
|
|
|
calc_t a1p = azimuth<calc_t, inverse_formula>(p1, p, m_model);
|
|
calc_t a12 = azimuth<calc_t, inverse_formula>(p1, p2, m_model);
|
|
|
|
return formula::azimuth_side_value(a1p, a12);
|
|
}
|
|
|
|
Spheroid const& model() const
|
|
{
|
|
return m_model;
|
|
}
|
|
|
|
private:
|
|
template <typename ResultType,
|
|
typename InverseFormulaType,
|
|
typename Point1,
|
|
typename Point2,
|
|
typename ModelT>
|
|
static inline ResultType azimuth(Point1 const& point1, Point2 const& point2,
|
|
ModelT const& model)
|
|
{
|
|
return InverseFormulaType::apply(get_as_radian<0>(point1),
|
|
get_as_radian<1>(point1),
|
|
get_as_radian<0>(point2),
|
|
get_as_radian<1>(point2),
|
|
model).azimuth;
|
|
}
|
|
|
|
Spheroid m_model;
|
|
};
|
|
|
|
|
|
}} // namespace strategy::side
|
|
|
|
|
|
}} // namespace boost::geometry
|
|
|
|
|
|
#endif // BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP
|