123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165 |
- #ifndef BOOST_GEOMETRY_PROJECTIONS_RPOLY_HPP
- #define BOOST_GEOMETRY_PROJECTIONS_RPOLY_HPP
- #include <boost/geometry/srs/projections/impl/base_static.hpp>
- #include <boost/geometry/srs/projections/impl/base_dynamic.hpp>
- #include <boost/geometry/srs/projections/impl/factory_entry.hpp>
- #include <boost/geometry/srs/projections/impl/pj_param.hpp>
- #include <boost/geometry/srs/projections/impl/projects.hpp>
- namespace boost { namespace geometry
- {
- namespace projections
- {
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail { namespace rpoly
- {
- static const double epsilon = 1e-9;
- template <typename T>
- struct par_rpoly
- {
- T phi1;
- T fxa;
- T fxb;
- bool mode;
- };
- template <typename T, typename Parameters>
- struct base_rpoly_spheroid
- {
- par_rpoly<T> m_proj_parm;
-
-
- inline void fwd(Parameters const& par, T const& lp_lon, T const& lp_lat, T& xy_x, T& xy_y) const
- {
- T fa;
- if (this->m_proj_parm.mode)
- fa = tan(lp_lon * this->m_proj_parm.fxb) * this->m_proj_parm.fxa;
- else
- fa = 0.5 * lp_lon;
- if (fabs(lp_lat) < epsilon) {
- xy_x = fa + fa;
- xy_y = - par.phi0;
- } else {
- xy_y = 1. / tan(lp_lat);
- xy_x = sin(fa = 2. * atan(fa * sin(lp_lat))) * xy_y;
- xy_y = lp_lat - par.phi0 + (1. - cos(fa)) * xy_y;
- }
- }
- static inline std::string get_name()
- {
- return "rpoly_spheroid";
- }
- };
-
- template <typename Params, typename Parameters, typename T>
- inline void setup_rpoly(Params const& params, Parameters& par, par_rpoly<T>& proj_parm)
- {
- proj_parm.phi1 = fabs(pj_get_param_r<T, srs::spar::lat_ts>(params, "lat_ts", srs::dpar::lat_ts));
- if ((proj_parm.mode = (proj_parm.phi1 > epsilon)))
- {
- proj_parm.fxb = 0.5 * sin(proj_parm.phi1);
- proj_parm.fxa = 0.5 / proj_parm.fxb;
- }
- par.es = 0.;
- }
- }}
- #endif
-
- template <typename T, typename Parameters>
- struct rpoly_spheroid : public detail::rpoly::base_rpoly_spheroid<T, Parameters>
- {
- template <typename Params>
- inline rpoly_spheroid(Params const& params, Parameters & par)
- {
- detail::rpoly::setup_rpoly(params, par, this->m_proj_parm);
- }
- };
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail
- {
-
- BOOST_GEOMETRY_PROJECTIONS_DETAIL_STATIC_PROJECTION_F(srs::spar::proj_rpoly, rpoly_spheroid)
-
- BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_ENTRY_F(rpoly_entry, rpoly_spheroid)
- BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_BEGIN(rpoly_init)
- {
- BOOST_GEOMETRY_PROJECTIONS_DETAIL_FACTORY_INIT_ENTRY(rpoly, rpoly_entry)
- }
- }
- #endif
- }
- }}
- #endif
|