ceph/src/boost/boost/geometry/srs/projections/proj/vandg4.hpp

   1 // Boost.Geometry - gis-projections (based on PROJ4)
   2
   3 // Copyright (c) 2008-2015 Barend Gehrels, Amsterdam, the Netherlands.
   4
   5 // This file was modified by Oracle on 2017, 2018, 2019.
   6 // Modifications copyright (c) 2017-2019, Oracle and/or its affiliates.
   7 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle.
   8
   9 // Use, modification and distribution is subject to the Boost Software License,
  10 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  11 // http://www.boost.org/LICENSE_1_0.txt)
  12
  13 // This file is converted from PROJ4, http://trac.osgeo.org/proj
  14 // PROJ4 is originally written by Gerald Evenden (then of the USGS)
  15 // PROJ4 is maintained by Frank Warmerdam
  16 // PROJ4 is converted to Boost.Geometry by Barend Gehrels
  17
  18 // Last updated version of proj: 5.0.0
  19
  20 // Original copyright notice:
  21
  22 // Permission is hereby granted, free of charge, to any person obtaining a
  23 // copy of this software and associated documentation files (the "Software"),
  24 // to deal in the Software without restriction, including without limitation
  25 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
  26 // and/or sell copies of the Software, and to permit persons to whom the
  27 // Software is furnished to do so, subject to the following conditions:
  28
  29 // The above copyright notice and this permission notice shall be included
  30 // in all copies or substantial portions of the Software.
  31
  32 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
  33 // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  34 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
  35 // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  36 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  37 // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
  38 // DEALINGS IN THE SOFTWARE.
  39
  40 #ifndef BOOST_GEOMETRY_PROJECTIONS_VANDG4_HPP
  41 #define BOOST_GEOMETRY_PROJECTIONS_VANDG4_HPP
  42
  43 #include <boost/geometry/util/math.hpp>
  44
  45 #include <boost/geometry/srs/projections/impl/base_static.hpp>
  46 #include <boost/geometry/srs/projections/impl/base_dynamic.hpp>
  47 #include <boost/geometry/srs/projections/impl/projects.hpp>
  48 #include <boost/geometry/srs/projections/impl/factory_entry.hpp>
  49
  50 namespace boost { namespace geometry
  51 {
  52
  53 namespace projections
  54 {
  55     #ifndef DOXYGEN_NO_DETAIL
  56     namespace detail { namespace vandg4
  57     {
  58
  59             static const double tolerance = 1e-10;
  60
  61             template <typename T, typename Parameters>
  62             struct base_vandg4_spheroid
  63             {
  64                 // FORWARD(s_forward)  spheroid
  65                 // Project coordinates from geographic (lon, lat) to cartesian (x, y)
  66                 inline void fwd(Parameters const& , T const& lp_lon, T const& lp_lat, T& xy_x, T& xy_y) const
  67                 {
  68                     static const T half_pi = detail::half_pi<T>();
  69                     static const T two_div_pi = detail::two_div_pi<T>();
  70
  71                     T x1, t, bt, ct, ft, bt2, ct2, dt, dt2;
  72
  73                     if (fabs(lp_lat) < tolerance) {
  74                         xy_x = lp_lon;
  75                         xy_y = 0.;
  76                     } else if (fabs(lp_lon) < tolerance || fabs(fabs(lp_lat) - half_pi) < tolerance) {
  77                         xy_x = 0.;
  78                         xy_y = lp_lat;
  79                     } else {
  80                         bt = fabs(two_div_pi * lp_lat);
  81                         bt2 = bt * bt;
  82                         ct = 0.5 * (bt * (8. - bt * (2. + bt2)) - 5.)
  83                             / (bt2 * (bt - 1.));
  84                         ct2 = ct * ct;
  85                         dt = two_div_pi * lp_lon;
  86                         dt = dt + 1. / dt;
  87                         dt = sqrt(dt * dt - 4.);
  88                         if ((fabs(lp_lon) - half_pi) < 0.) dt = -dt;
  89                         dt2 = dt * dt;
  90                         x1 = bt + ct; x1 *= x1;
  91                         t = bt + 3.*ct;
  92                         ft = x1 * (bt2 + ct2 * dt2 - 1.) + (1.-bt2) * (
  93                             bt2 * (t * t + 4. * ct2) +
  94                             ct2 * (12. * bt * ct + 4. * ct2) );
  95                         x1 = (dt*(x1 + ct2 - 1.) + 2.*sqrt(ft)) /
  96                             (4.* x1 + dt2);
  97                         xy_x = half_pi * x1;
  98                         xy_y = half_pi * sqrt(1. + dt * fabs(x1) - x1 * x1);
  99                         if (lp_lon < 0.) xy_x = -xy_x;
 100                         if (lp_lat < 0.) xy_y = -xy_y;
 101                     }
 102                 }
 103
 104                 static inline std::string get_name()
 105                 {
 106                     return "vandg4_spheroid";
 107                 }
 108
 109             };
 110
 111             // van der Grinten IV
 112             template <typename Parameters>
 113             inline void setup_vandg4(Parameters& par)
 114             {
 115                 par.es = 0.;
 116             }
 117
 118     }} // namespace detail::vandg4
 119     #endif // doxygen
 120
 121     /*!
 122         \brief van der Grinten IV projection
 123         \ingroup projections
 124         \tparam Geographic latlong point type
 125         \tparam Cartesian xy point type
 126         \tparam Parameters parameter type
 127         \par Projection characteristics
 128          - Miscellaneous
 129          - Spheroid
 130          - no inverse
 131         \par Example
 132         \image html ex_vandg4.gif
 133     */
 134     template <typename T, typename Parameters>
 135     struct vandg4_spheroid : public detail::vandg4::base_vandg4_spheroid<T, Parameters>
 136     {
 137         template <typename Params>
 138         inline vandg4_spheroid(Params const& , Parameters & par)
 139         {
 140             detail::vandg4::setup_vandg4(par);
 141         }
 142     };
 143
 144     #ifndef DOXYGEN_NO_DETAIL
 145     namespace detail
 146     {
 147
 148         // Static projection
 149         BOOST_GEOMETRY_PROJECTIONS_DETAIL_STATIC_PROJECTION_F(srs::spar::proj_vandg4, vandg4_spheroid)
 150
 151         // Factory entry(s)
 152         BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_ENTRY_F(vandg4_entry, vandg4_spheroid)
 153
 154         BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_BEGIN(vandg4_init)
 155         {
 156             BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_ENTRY(vandg4, vandg4_entry)
 157         }
 158
 159     } // namespace detail
 160     #endif // doxygen
 161
 162 } // namespace projections
 163
 164 }} // namespace boost::geometry
 165
 166 #endif // BOOST_GEOMETRY_PROJECTIONS_VANDG4_HPP
 167