// Boost.Geometry
|
|
// Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
|
|
// This file was modified by Oracle on 2014-2020.
|
// Modifications copyright (c) 2014-2020 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/radian_access.hpp>
|
#include <boost/geometry/core/radius.hpp>
|
|
#include <boost/geometry/formulas/spherical.hpp>
|
|
#include <boost/geometry/srs/spheroid.hpp>
|
|
#include <boost/geometry/util/math.hpp>
|
#include <boost/geometry/util/promote_floating_point.hpp>
|
#include <boost/geometry/util/select_calculation_type.hpp>
|
|
#include <boost/geometry/strategy/geographic/envelope.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/strategies/concepts/side_concept.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;
|
|
typedef strategy::envelope::geographic
|
<
|
FormulaPolicy,
|
Spheroid,
|
CalculationType
|
> envelope_strategy_type;
|
|
inline envelope_strategy_type get_envelope_strategy() const
|
{
|
return envelope_strategy_type(m_model);
|
}
|
|
typedef strategy::disjoint::segment_box_geographic
|
<
|
FormulaPolicy,
|
Spheroid,
|
CalculationType
|
> disjoint_strategy_type;
|
|
inline disjoint_strategy_type get_disjoint_strategy() const
|
{
|
return disjoint_strategy_type(m_model);
|
}
|
|
typedef strategy::within::spherical_point_point equals_point_point_strategy_type;
|
static inline equals_point_point_strategy_type get_equals_point_point_strategy()
|
{
|
return equals_point_point_strategy_type();
|
}
|
|
geographic()
|
{}
|
|
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 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
|
