projections/proj/tcea.hpp

0001 // Boost.Geometry - gis-projections (based on PROJ4)
0002
0003 // Copyright (c) 2008-2015 Barend Gehrels, Amsterdam, the Netherlands.
0004
0005 // This file was modified by Oracle on 2017, 2018, 2019.
0006 // Modifications copyright (c) 2017-2019, Oracle and/or its affiliates.
0007 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle.
0008
0009 // Use, modification and distribution is subject to the Boost Software License,
0010 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
0011 // http://www.boost.org/LICENSE_1_0.txt)
0012
0013 // This file is converted from PROJ4, http://trac.osgeo.org/proj
0014 // PROJ4 is originally written by Gerald Evenden (then of the USGS)
0015 // PROJ4 is maintained by Frank Warmerdam
0016 // PROJ4 is converted to Boost.Geometry by Barend Gehrels
0017
0018 // Last updated version of proj: 5.0.0
0019
0020 // Original copyright notice:
0021
0022 // Permission is hereby granted, free of charge, to any person obtaining a
0023 // copy of this software and associated documentation files (the "Software"),
0024 // to deal in the Software without restriction, including without limitation
0025 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
0026 // and/or sell copies of the Software, and to permit persons to whom the
0027 // Software is furnished to do so, subject to the following conditions:
0028
0029 // The above copyright notice and this permission notice shall be included
0030 // in all copies or substantial portions of the Software.
0031
0032 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
0033 // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
0034 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
0035 // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
0036 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
0037 // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
0038 // DEALINGS IN THE SOFTWARE.
0039
0040 #ifndef BOOST_GEOMETRY_PROJECTIONS_TCEA_HPP
0041 #define BOOST_GEOMETRY_PROJECTIONS_TCEA_HPP
0042
0043 #include <boost/geometry/srs/projections/impl/base_static.hpp>
0044 #include <boost/geometry/srs/projections/impl/base_dynamic.hpp>
0045 #include <boost/geometry/srs/projections/impl/projects.hpp>
0046 #include <boost/geometry/srs/projections/impl/factory_entry.hpp>
0047
0048 namespace boost { namespace geometry
0049 {
0050
0051 namespace projections
0052 {
0053     #ifndef DOXYGEN_NO_DETAIL
0054     namespace detail { namespace tcea
0055     {
0056             template <typename T, typename Parameters>
0057             struct base_tcea_spheroid
0058             {
0059                 // FORWARD(s_forward)  spheroid
0060                 // Project coordinates from geographic (lon, lat) to cartesian (x, y)
0061                 inline void fwd(Parameters const& par, T const& lp_lon, T const& lp_lat, T& xy_x, T& xy_y) const
0062                 {
0063                     xy_x = cos(lp_lat) * sin(lp_lon) / par.k0;
0064                     xy_y = par.k0 * (atan2(tan(lp_lat), cos(lp_lon)) - par.phi0);
0065                 }
0066
0067                 // INVERSE(s_inverse)  spheroid
0068                 // Project coordinates from cartesian (x, y) to geographic (lon, lat)
0069                 inline void inv(Parameters const& par, T xy_x, T xy_y, T& lp_lon, T& lp_lat) const
0070                 {
0071                     T t;
0072
0073                     xy_y = xy_y / par.k0 + par.phi0;
0074                     xy_x *= par.k0;
0075                     t = sqrt(1. - xy_x * xy_x);
0076                     lp_lat = asin(t * sin(xy_y));
0077                     lp_lon = atan2(xy_x, t * cos(xy_y));
0078                 }
0079
0080                 static inline std::string get_name()
0081                 {
0082                     return "tcea_spheroid";
0083                 }
0084
0085             };
0086
0087             // Transverse Cylindrical Equal Area
0088             template <typename Parameters>
0089             inline void setup_tcea(Parameters& par)
0090             {
0091                 par.es = 0.;
0092             }
0093
0094     }} // namespace detail::tcea
0095     #endif // doxygen
0096
0097     /*!
0098         \brief Transverse Cylindrical Equal Area projection
0099         \ingroup projections
0100         \tparam Geographic latlong point type
0101         \tparam Cartesian xy point type
0102         \tparam Parameters parameter type
0103         \par Projection characteristics
0104          - Cylindrical
0105          - Spheroid
0106         \par Example
0107         \image html ex_tcea.gif
0108     */
0109     template <typename T, typename Parameters>
0110     struct tcea_spheroid : public detail::tcea::base_tcea_spheroid<T, Parameters>
0111     {
0112         template <typename Params>
0113         inline tcea_spheroid(Params const& , Parameters & par)
0114         {
0115             detail::tcea::setup_tcea(par);
0116         }
0117     };
0118
0119     #ifndef DOXYGEN_NO_DETAIL
0120     namespace detail
0121     {
0122
0123         // Static projection
0124         BOOST_GEOMETRY_PROJECTIONS_DETAIL_STATIC_PROJECTION_FI(srs::spar::proj_tcea, tcea_spheroid)
0125
0126         // Factory entry(s)
0127         BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_ENTRY_FI(tcea_entry, tcea_spheroid)
0128
0129         BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_BEGIN(tcea_init)
0130         {
0131             BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_ENTRY(tcea, tcea_entry)
0132         }
0133
0134     } // namespace detail
0135     #endif // doxygen
0136
0137 } // namespace projections
0138
0139 }} // namespace boost::geometry
0140
0141 #endif // BOOST_GEOMETRY_PROJECTIONS_TCEA_HPP
0142