// Boost.Geometry (aka GGL, Generic Geometry Library) // Copyright (c) 2007-2016 Barend Gehrels, Amsterdam, the Netherlands. // This file was modified by Oracle on 2014, 2016. // Modifications copyright (c) 2014-2016 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_ANDOYER_HPP #define BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_ANDOYER_HPP #include #include #include #include #include #include #include #include #include #include namespace boost { namespace geometry { namespace strategy { namespace distance { /*! \brief Point-point distance approximation taking flattening into account \ingroup distance \tparam Spheroid The reference spheroid model \tparam CalculationType \tparam_calculation \author After Andoyer, 19xx, republished 1950, republished by Meeus, 1999 \note Although not so well-known, the approximation is very good: in all cases the results are about the same as Vincenty. In my (Barend's) testcases the results didn't differ more than 6 m \see http://nacc.upc.es/tierra/node16.html \see http://sci.tech-archive.net/Archive/sci.geo.satellite-nav/2004-12/2724.html \see http://home.att.net/~srschmitt/great_circle_route.html (implementation) \see http://www.codeguru.com/Cpp/Cpp/algorithms/article.php/c5115 (implementation) \see http://futureboy.homeip.net/frinksamp/navigation.frink (implementation) \see http://www.voidware.com/earthdist.htm (implementation) \see http://www.dtic.mil/docs/citations/AD0627893 \see http://www.dtic.mil/docs/citations/AD703541 */ template < typename Spheroid, typename CalculationType = void > class andoyer { public : template struct calculation_type : promote_floating_point < typename select_calculation_type < Point1, Point2, CalculationType >::type > {}; typedef Spheroid model_type; inline andoyer() : m_spheroid() {} explicit inline andoyer(Spheroid const& spheroid) : m_spheroid(spheroid) {} template inline typename calculation_type::type apply(Point1 const& point1, Point2 const& point2) const { return geometry::detail::andoyer_inverse < typename calculation_type::type, true, false >::apply(get_as_radian<0>(point1), get_as_radian<1>(point1), get_as_radian<0>(point2), get_as_radian<1>(point2), m_spheroid).distance; } inline Spheroid const& model() const { return m_spheroid; } private : Spheroid m_spheroid; }; #ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS namespace services { template struct tag > { typedef strategy_tag_distance_point_point type; }; template struct return_type, P1, P2> : andoyer::template calculation_type {}; template struct comparable_type > { typedef andoyer type; }; template struct get_comparable > { static inline andoyer apply(andoyer const& input) { return input; } }; template struct result_from_distance, P1, P2> { template static inline typename return_type, P1, P2>::type apply(andoyer const& , T const& value) { return value; } }; template struct default_strategy { typedef strategy::distance::andoyer < srs::spheroid < typename select_coordinate_type::type > > type; }; } // namespace services #endif // DOXYGEN_NO_STRATEGY_SPECIALIZATIONS }} // namespace strategy::distance }} // namespace boost::geometry #endif // BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_ANDOYER_HPP