Chris@16
|
1 // Boost.Geometry (aka GGL, Generic Geometry Library)
|
Chris@16
|
2
|
Chris@16
|
3 // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
|
Chris@16
|
4 // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
|
Chris@16
|
5 // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
|
Chris@16
|
6
|
Chris@16
|
7 // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
|
Chris@16
|
8 // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
|
Chris@16
|
9
|
Chris@16
|
10 // Use, modification and distribution is subject to the Boost Software License,
|
Chris@16
|
11 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
|
Chris@16
|
12 // http://www.boost.org/LICENSE_1_0.txt)
|
Chris@16
|
13
|
Chris@16
|
14 #ifndef BOOST_GEOMETRY_MULTI_STRATEGIES_CARTESIAN_CENTROID_AVERAGE_HPP
|
Chris@16
|
15 #define BOOST_GEOMETRY_MULTI_STRATEGIES_CARTESIAN_CENTROID_AVERAGE_HPP
|
Chris@16
|
16
|
Chris@16
|
17
|
Chris@16
|
18 #include <boost/numeric/conversion/cast.hpp>
|
Chris@16
|
19 #include <boost/type_traits.hpp>
|
Chris@16
|
20
|
Chris@16
|
21 #include <boost/geometry/core/coordinate_type.hpp>
|
Chris@16
|
22 #include <boost/geometry/core/point_type.hpp>
|
Chris@16
|
23 #include <boost/geometry/arithmetic/arithmetic.hpp>
|
Chris@16
|
24 #include <boost/geometry/strategies/centroid.hpp>
|
Chris@16
|
25
|
Chris@16
|
26
|
Chris@16
|
27 namespace boost { namespace geometry
|
Chris@16
|
28 {
|
Chris@16
|
29
|
Chris@16
|
30 namespace strategy { namespace centroid
|
Chris@16
|
31 {
|
Chris@16
|
32
|
Chris@16
|
33
|
Chris@16
|
34 /*!
|
Chris@16
|
35 \brief Centroid calculation taking average of points
|
Chris@16
|
36 \ingroup strategies
|
Chris@16
|
37 */
|
Chris@16
|
38 template
|
Chris@16
|
39 <
|
Chris@16
|
40 typename PointCentroid,
|
Chris@16
|
41 typename Point = PointCentroid
|
Chris@16
|
42 >
|
Chris@16
|
43 class average
|
Chris@16
|
44 {
|
Chris@16
|
45 private :
|
Chris@16
|
46
|
Chris@16
|
47 /*! subclass to keep state */
|
Chris@16
|
48 class sum
|
Chris@16
|
49 {
|
Chris@16
|
50 friend class average;
|
Chris@16
|
51 int count;
|
Chris@16
|
52 PointCentroid centroid;
|
Chris@16
|
53
|
Chris@16
|
54 public :
|
Chris@16
|
55 inline sum()
|
Chris@16
|
56 : count(0)
|
Chris@16
|
57 {
|
Chris@16
|
58 assign_zero(centroid);
|
Chris@16
|
59 }
|
Chris@16
|
60 };
|
Chris@16
|
61
|
Chris@16
|
62 public :
|
Chris@16
|
63 typedef sum state_type;
|
Chris@16
|
64 typedef PointCentroid centroid_point_type;
|
Chris@16
|
65 typedef Point point_type;
|
Chris@16
|
66
|
Chris@16
|
67 static inline void apply(Point const& p, sum& state)
|
Chris@16
|
68 {
|
Chris@16
|
69 add_point(state.centroid, p);
|
Chris@16
|
70 state.count++;
|
Chris@16
|
71 }
|
Chris@16
|
72
|
Chris@16
|
73 static inline void result(sum const& state, PointCentroid& centroid)
|
Chris@16
|
74 {
|
Chris@16
|
75 centroid = state.centroid;
|
Chris@16
|
76 divide_value(centroid, state.count);
|
Chris@16
|
77 }
|
Chris@16
|
78
|
Chris@16
|
79 };
|
Chris@16
|
80
|
Chris@16
|
81
|
Chris@16
|
82 #ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
|
Chris@16
|
83
|
Chris@16
|
84
|
Chris@16
|
85 namespace services
|
Chris@16
|
86 {
|
Chris@16
|
87
|
Chris@16
|
88 template <typename Point, std::size_t DimensionCount, typename Geometry>
|
Chris@16
|
89 struct default_strategy
|
Chris@16
|
90 <
|
Chris@16
|
91 cartesian_tag,
|
Chris@16
|
92 pointlike_tag,
|
Chris@16
|
93 DimensionCount,
|
Chris@16
|
94 Point,
|
Chris@16
|
95 Geometry
|
Chris@16
|
96 >
|
Chris@16
|
97 {
|
Chris@16
|
98 typedef average
|
Chris@16
|
99 <
|
Chris@16
|
100 Point,
|
Chris@16
|
101 typename point_type<Geometry>::type
|
Chris@16
|
102 > type;
|
Chris@16
|
103 };
|
Chris@16
|
104
|
Chris@16
|
105 } // namespace services
|
Chris@16
|
106
|
Chris@16
|
107 #endif
|
Chris@16
|
108
|
Chris@16
|
109
|
Chris@16
|
110 }} // namespace strategy::centroid
|
Chris@16
|
111
|
Chris@16
|
112
|
Chris@16
|
113 }} // namespace boost::geometry
|
Chris@16
|
114
|
Chris@16
|
115
|
Chris@16
|
116 #endif // BOOST_GEOMETRY_MULTI_STRATEGIES_CARTESIAN_CENTROID_AVERAGE_HPP
|