← Back to branch summary

~ubuntu-branches/ubuntu/saucy/merkaartor/saucy

« back to all changes in this revision

Viewing changes to include/ggl/projections/proj/mill.hpp

  • Committer: Bazaar Package Importer
  • Author(s): Bernd Zeimetz
  • Date: 2009-09-13 00:52:12 UTC
  • mto: (1.2.7 upstream) (0.1.3 upstream) (3.1.7 sid)
  • mto: This revision was merged to the branch mainline in revision 10.
  • Revision ID: james.westby@ubuntu.com-20090913005212-pjecal8zxm07x0fj
Import upstream version 0.14+svnfixes~20090912

Show diffs side-by-side

added added

removed removed

Lines of Context:
include/ggl/projections/proj/mill.hpp
 
1
#ifndef GGL_PROJECTIONS_MILL_HPP
 
2
#define GGL_PROJECTIONS_MILL_HPP
 
3
 
 
4
// Generic Geometry Library - projections (based on PROJ4)
 
5
// This file is automatically generated. DO NOT EDIT.
 
6
 
 
7
// Copyright Barend Gehrels (1995-2009), Geodan Holding B.V. Amsterdam, the Netherlands.
 
8
// Copyright Bruno Lalande (2008-2009)
 
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 Geometry Library by Barend Gehrels (Geodan, Amsterdam)
 
17
 
 
18
// Original copyright notice:
 
19
 
 
20
// Permission is hereby granted, free of charge, to any person obtaining a
 
21
// copy of this software and associated documentation files (the "Software"),
 
22
// to deal in the Software without restriction, including without limitation
 
23
// the rights to use, copy, modify, merge, publish, distribute, sublicense,
 
24
// and/or sell copies of the Software, and to permit persons to whom the
 
25
// Software is furnished to do so, subject to the following conditions:
 
26
 
 
27
// The above copyright notice and this permission notice shall be included
 
28
// in all copies or substantial portions of the Software.
 
29
 
 
30
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
 
31
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 
32
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
 
33
// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 
34
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
 
35
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
 
36
// DEALINGS IN THE SOFTWARE.
 
37
 
 
38
#include <boost/math/special_functions/hypot.hpp>
 
39
 
 
40
#include <ggl/projections/impl/base_static.hpp>
 
41
#include <ggl/projections/impl/base_dynamic.hpp>
 
42
#include <ggl/projections/impl/projects.hpp>
 
43
#include <ggl/projections/impl/factory_entry.hpp>
 
44
 
 
45
namespace ggl { namespace projection
 
46
{
 
47
    #ifndef DOXYGEN_NO_DETAIL
 
48
    namespace detail { namespace mill{ 
 
49
 
 
50
 
 
51
            // template class, using CRTP to implement forward/inverse
 
52
            template <typename Geographic, typename Cartesian, typename Parameters>
 
53
            struct base_mill_spheroid : public base_t_fi<base_mill_spheroid<Geographic, Cartesian, Parameters>,
 
54
                     Geographic, Cartesian, Parameters>
 
55
            {
 
56
 
 
57
                 typedef double geographic_type;
 
58
                 typedef double cartesian_type;
 
59
 
 
60
 
 
61
                inline base_mill_spheroid(const Parameters& par)
 
62
                    : base_t_fi<base_mill_spheroid<Geographic, Cartesian, Parameters>,
 
63
                     Geographic, Cartesian, Parameters>(*this, par) {}
 
64
 
 
65
                inline void fwd(geographic_type& lp_lon, geographic_type& lp_lat, cartesian_type& xy_x, cartesian_type& xy_y) const
 
66
                {
 
67
                        xy_x = lp_lon;
 
68
                        xy_y = log(tan(FORTPI + lp_lat * .4)) * 1.25;
 
69
                }
 
70
 
 
71
                inline void inv(cartesian_type& xy_x, cartesian_type& xy_y, geographic_type& lp_lon, geographic_type& lp_lat) const
 
72
                {
 
73
                        lp_lon = xy_x;
 
74
                        lp_lat = 2.5 * (atan(exp(.8 * xy_y)) - FORTPI);
 
75
                }
 
76
            };
 
77
 
 
78
            // Miller Cylindrical
 
79
            template <typename Parameters>
 
80
            void setup_mill(Parameters& par)
 
81
            {
 
82
                par.es = 0.;
 
83
                // par.inv = s_inverse;
 
84
                // par.fwd = s_forward;
 
85
            }
 
86
 
 
87
        }} // namespace detail::mill
 
88
    #endif // doxygen 
 
89
 
 
90
    /*!
 
91
        \brief Miller Cylindrical projection
 
92
        \ingroup projections
 
93
        \tparam Geographic latlong point type
 
94
        \tparam Cartesian xy point type
 
95
        \tparam Parameters parameter type
 
96
        \par Projection characteristics
 
97
         - Cylindrical
 
98
         - Spheroid
 
99
        \par Example
 
100
        \image html ex_mill.gif
 
101
    */
 
102
    template <typename Geographic, typename Cartesian, typename Parameters = parameters>
 
103
    struct mill_spheroid : public detail::mill::base_mill_spheroid<Geographic, Cartesian, Parameters>
 
104
    {
 
105
        inline mill_spheroid(const Parameters& par) : detail::mill::base_mill_spheroid<Geographic, Cartesian, Parameters>(par)
 
106
        {
 
107
            detail::mill::setup_mill(this->m_par);
 
108
        }
 
109
    };
 
110
 
 
111
    #ifndef DOXYGEN_NO_DETAIL
 
112
    namespace detail
 
113
    {
 
114
 
 
115
        // Factory entry(s)
 
116
        template <typename Geographic, typename Cartesian, typename Parameters>
 
117
        class mill_entry : public detail::factory_entry<Geographic, Cartesian, Parameters>
 
118
        {
 
119
            public :
 
120
                virtual projection<Geographic, Cartesian>* create_new(const Parameters& par) const
 
121
                {
 
122
                    return new base_v_fi<mill_spheroid<Geographic, Cartesian, Parameters>, Geographic, Cartesian, Parameters>(par);
 
123
                }
 
124
        };
 
125
 
 
126
        template <typename Geographic, typename Cartesian, typename Parameters>
 
127
        inline void mill_init(detail::base_factory<Geographic, Cartesian, Parameters>& factory)
 
128
        {
 
129
            factory.add_to_factory("mill", new mill_entry<Geographic, Cartesian, Parameters>);
 
130
        }
 
131
 
 
132
    } // namespace detail 
 
133
    #endif // doxygen
 
134
 
 
135
}} // namespace ggl::projection
 
136
 
 
137
#endif // GGL_PROJECTIONS_MILL_HPP
 
138