]>
Commit | Line | Data |
---|---|---|
92f5a8d4 | 1 | // Boost.Geometry - gis-projections (based on PROJ4) |
11fdf7f2 TL |
2 | |
3 | // Copyright (c) 2008-2015 Barend Gehrels, Amsterdam, the Netherlands. | |
4 | ||
92f5a8d4 TL |
5 | // This file was modified by Oracle on 2017, 2018, 2019. |
6 | // Modifications copyright (c) 2017-2019, Oracle and/or its affiliates. | |
11fdf7f2 TL |
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 | ||
92f5a8d4 | 18 | // Last updated version of proj: 5.0.0 |
11fdf7f2 TL |
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 | ||
92f5a8d4 TL |
40 | #ifndef BOOST_GEOMETRY_PROJECTIONS_PUTP4P_HPP |
41 | #define BOOST_GEOMETRY_PROJECTIONS_PUTP4P_HPP | |
11fdf7f2 TL |
42 | |
43 | #include <boost/geometry/srs/projections/impl/base_static.hpp> | |
44 | #include <boost/geometry/srs/projections/impl/base_dynamic.hpp> | |
45 | #include <boost/geometry/srs/projections/impl/projects.hpp> | |
46 | #include <boost/geometry/srs/projections/impl/factory_entry.hpp> | |
47 | #include <boost/geometry/srs/projections/impl/aasincos.hpp> | |
48 | ||
49 | namespace boost { namespace geometry | |
50 | { | |
51 | ||
11fdf7f2 TL |
52 | namespace projections |
53 | { | |
54 | #ifndef DOXYGEN_NO_DETAIL | |
55 | namespace detail { namespace putp4p | |
56 | { | |
57 | template <typename T> | |
58 | struct par_putp4p | |
59 | { | |
60 | T C_x, C_y; | |
61 | }; | |
62 | ||
92f5a8d4 TL |
63 | template <typename T, typename Parameters> |
64 | struct base_putp4p_spheroid | |
11fdf7f2 | 65 | { |
92f5a8d4 | 66 | par_putp4p<T> m_proj_parm; |
11fdf7f2 TL |
67 | |
68 | // FORWARD(s_forward) spheroid | |
69 | // Project coordinates from geographic (lon, lat) to cartesian (x, y) | |
92f5a8d4 | 70 | inline void fwd(Parameters const& , T const& lp_lon, T lp_lat, T& xy_x, T& xy_y) const |
11fdf7f2 | 71 | { |
92f5a8d4 TL |
72 | static T const third = detail::third<T>(); |
73 | ||
11fdf7f2 TL |
74 | lp_lat = aasin(0.883883476 * sin(lp_lat)); |
75 | xy_x = this->m_proj_parm.C_x * lp_lon * cos(lp_lat); | |
92f5a8d4 | 76 | xy_x /= cos(lp_lat *= third); |
11fdf7f2 TL |
77 | xy_y = this->m_proj_parm.C_y * sin(lp_lat); |
78 | } | |
79 | ||
80 | // INVERSE(s_inverse) spheroid | |
81 | // Project coordinates from cartesian (x, y) to geographic (lon, lat) | |
92f5a8d4 | 82 | inline void inv(Parameters const& , T const& xy_x, T const& xy_y, T& lp_lon, T& lp_lat) const |
11fdf7f2 TL |
83 | { |
84 | lp_lat = aasin(xy_y / this->m_proj_parm.C_y); | |
85 | lp_lon = xy_x * cos(lp_lat) / this->m_proj_parm.C_x; | |
86 | lp_lat *= 3.; | |
87 | lp_lon /= cos(lp_lat); | |
88 | lp_lat = aasin(1.13137085 * sin(lp_lat)); | |
89 | } | |
90 | ||
91 | static inline std::string get_name() | |
92 | { | |
93 | return "putp4p_spheroid"; | |
94 | } | |
95 | ||
96 | }; | |
97 | ||
11fdf7f2 TL |
98 | |
99 | // Putnins P4' | |
100 | template <typename Parameters, typename T> | |
101 | inline void setup_putp4p(Parameters& par, par_putp4p<T>& proj_parm) | |
102 | { | |
103 | proj_parm.C_x = 0.874038744; | |
104 | proj_parm.C_y = 3.883251825; | |
92f5a8d4 TL |
105 | |
106 | par.es = 0.; | |
11fdf7f2 TL |
107 | } |
108 | ||
109 | // Werenskiold I | |
110 | template <typename Parameters, typename T> | |
111 | inline void setup_weren(Parameters& par, par_putp4p<T>& proj_parm) | |
112 | { | |
113 | proj_parm.C_x = 1.; | |
114 | proj_parm.C_y = 4.442882938; | |
92f5a8d4 TL |
115 | |
116 | par.es = 0.; | |
11fdf7f2 TL |
117 | } |
118 | ||
119 | }} // namespace detail::putp4p | |
120 | #endif // doxygen | |
121 | ||
122 | /*! | |
123 | \brief Putnins P4' projection | |
124 | \ingroup projections | |
125 | \tparam Geographic latlong point type | |
126 | \tparam Cartesian xy point type | |
127 | \tparam Parameters parameter type | |
128 | \par Projection characteristics | |
129 | - Pseudocylindrical | |
130 | - Spheroid | |
131 | \par Example | |
132 | \image html ex_putp4p.gif | |
133 | */ | |
92f5a8d4 TL |
134 | template <typename T, typename Parameters> |
135 | struct putp4p_spheroid : public detail::putp4p::base_putp4p_spheroid<T, Parameters> | |
11fdf7f2 | 136 | { |
92f5a8d4 TL |
137 | template <typename Params> |
138 | inline putp4p_spheroid(Params const& , Parameters & par) | |
11fdf7f2 | 139 | { |
92f5a8d4 | 140 | detail::putp4p::setup_putp4p(par, this->m_proj_parm); |
11fdf7f2 TL |
141 | } |
142 | }; | |
143 | ||
144 | /*! | |
145 | \brief Werenskiold I projection | |
146 | \ingroup projections | |
147 | \tparam Geographic latlong point type | |
148 | \tparam Cartesian xy point type | |
149 | \tparam Parameters parameter type | |
150 | \par Projection characteristics | |
151 | - Pseudocylindrical | |
152 | - Spheroid | |
153 | \par Example | |
154 | \image html ex_weren.gif | |
155 | */ | |
92f5a8d4 TL |
156 | template <typename T, typename Parameters> |
157 | struct weren_spheroid : public detail::putp4p::base_putp4p_spheroid<T, Parameters> | |
11fdf7f2 | 158 | { |
92f5a8d4 TL |
159 | template <typename Params> |
160 | inline weren_spheroid(Params const& , Parameters & par) | |
11fdf7f2 | 161 | { |
92f5a8d4 | 162 | detail::putp4p::setup_weren(par, this->m_proj_parm); |
11fdf7f2 TL |
163 | } |
164 | }; | |
165 | ||
166 | #ifndef DOXYGEN_NO_DETAIL | |
167 | namespace detail | |
168 | { | |
169 | ||
170 | // Static projection | |
92f5a8d4 TL |
171 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_STATIC_PROJECTION_FI(srs::spar::proj_putp4p, putp4p_spheroid) |
172 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_STATIC_PROJECTION_FI(srs::spar::proj_weren, weren_spheroid) | |
11fdf7f2 TL |
173 | |
174 | // Factory entry(s) | |
92f5a8d4 TL |
175 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_ENTRY_FI(putp4p_entry, putp4p_spheroid) |
176 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_ENTRY_FI(weren_entry, weren_spheroid) | |
11fdf7f2 | 177 | |
92f5a8d4 | 178 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_BEGIN(putp4p_init) |
11fdf7f2 | 179 | { |
92f5a8d4 TL |
180 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_ENTRY(putp4p, putp4p_entry) |
181 | BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_ENTRY(weren, weren_entry) | |
11fdf7f2 TL |
182 | } |
183 | ||
184 | } // namespace detail | |
185 | #endif // doxygen | |
186 | ||
187 | } // namespace projections | |
188 | ||
189 | }} // namespace boost::geometry | |
190 | ||
191 | #endif // BOOST_GEOMETRY_PROJECTIONS_PUTP4P_HPP | |
192 |