1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2
3 // Copyright (c) 2007-2014 Barend Gehrels, Amsterdam, the Netherlands.
4 // Copyright (c) 2008-2014 Bruno Lalande, Paris, France.
5 // Copyright (c) 2009-2014 Mateusz Loskot, London, UK.
6
7 // This file was modified by Oracle on 2014-2021.
8 // Modifications copyright (c) 2014-2021, Oracle and/or its affiliates.
9 // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
10 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
11
12 // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
13 // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
14
15 // Use, modification and distribution is subject to the Boost Software License,
16 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
17 // http://www.boost.org/LICENSE_1_0.txt)
18
19 #ifndef BOOST_GEOMETRY_ALGORITHMS_APPEND_HPP
20 #define BOOST_GEOMETRY_ALGORITHMS_APPEND_HPP
21
22
23 #include <boost/range/begin.hpp>
24 #include <boost/range/end.hpp>
25 #include <boost/range/value_type.hpp>
26
27 #include <boost/geometry/algorithms/num_interior_rings.hpp>
28 #include <boost/geometry/algorithms/detail/convert_point_to_point.hpp>
29 #include <boost/geometry/algorithms/detail/signed_size_type.hpp>
30 #include <boost/geometry/core/access.hpp>
31 #include <boost/geometry/core/mutable_range.hpp>
32 #include <boost/geometry/core/point_type.hpp>
33 #include <boost/geometry/core/tags.hpp>
34 #include <boost/geometry/core/visit.hpp>
35 #include <boost/geometry/geometries/adapted/boost_variant.hpp> // for backward compatibility
36 #include <boost/geometry/geometries/concepts/check.hpp>
37 #include <boost/geometry/util/range.hpp>
38
39
40 namespace boost { namespace geometry
41 {
42
43
44 #ifndef DOXYGEN_NO_DETAIL
45 namespace detail { namespace append
46 {
47
48 struct append_no_action
49 {
50 template <typename Geometry, typename Point>
applyboost::geometry::detail::append::append_no_action51 static inline void apply(Geometry& , Point const& ,
52 signed_size_type = -1, signed_size_type = 0)
53 {
54 }
55 };
56
57 struct to_range_point
58 {
59 template <typename Geometry, typename Point>
applyboost::geometry::detail::append::to_range_point60 static inline void apply(Geometry& geometry, Point const& point,
61 signed_size_type = -1, signed_size_type = 0)
62 {
63 typename geometry::point_type<Geometry>::type copy;
64 geometry::detail::conversion::convert_point_to_point(point, copy);
65 traits::push_back<Geometry>::apply(geometry, copy);
66 }
67 };
68
69
70 struct to_range_range
71 {
72 template <typename Geometry, typename Range>
applyboost::geometry::detail::append::to_range_range73 static inline void apply(Geometry& geometry, Range const& range,
74 signed_size_type = -1, signed_size_type = 0)
75 {
76 using point_type = typename boost::range_value<Range>::type;
77
78 auto const end = boost::end(range);
79 for (auto it = boost::begin(range); it != end; ++it)
80 {
81 to_range_point::apply<Geometry, point_type>(geometry, *it);
82 }
83 }
84 };
85
86
87 struct to_polygon_point
88 {
89 template <typename Polygon, typename Point>
applyboost::geometry::detail::append::to_polygon_point90 static inline void apply(Polygon& polygon, Point const& point,
91 signed_size_type ring_index, signed_size_type = 0)
92 {
93 using ring_type = typename ring_type<Polygon>::type;
94 using exterior_ring_type = typename ring_return_type<Polygon>::type;
95 using interior_ring_range_type = typename interior_return_type<Polygon>::type;
96
97 if (ring_index == -1)
98 {
99 exterior_ring_type ext_ring = exterior_ring(polygon);
100 to_range_point::apply<ring_type, Point>(ext_ring, point);
101 }
102 else if (ring_index < signed_size_type(num_interior_rings(polygon)))
103 {
104 interior_ring_range_type int_rings = interior_rings(polygon);
105 to_range_point::apply<ring_type, Point>(range::at(int_rings, ring_index), point);
106 }
107 }
108 };
109
110
111 struct to_polygon_range
112 {
113 template <typename Polygon, typename Range>
applyboost::geometry::detail::append::to_polygon_range114 static inline void apply(Polygon& polygon, Range const& range,
115 signed_size_type ring_index, signed_size_type = 0)
116 {
117 using ring_type = typename ring_type<Polygon>::type;
118 using exterior_ring_type = typename ring_return_type<Polygon>::type;
119 using interior_ring_range_type = typename interior_return_type<Polygon>::type;
120
121 if (ring_index == -1)
122 {
123 exterior_ring_type ext_ring = exterior_ring(polygon);
124 to_range_range::apply<ring_type, Range>(ext_ring, range);
125 }
126 else if (ring_index < signed_size_type(num_interior_rings(polygon)))
127 {
128 interior_ring_range_type int_rings = interior_rings(polygon);
129 to_range_range::apply<ring_type, Range>(range::at(int_rings, ring_index), range);
130 }
131 }
132 };
133
134
135 template <typename Policy>
136 struct to_multigeometry
137 {
138 template <typename MultiGeometry, typename RangeOrPoint>
applyboost::geometry::detail::append::to_multigeometry139 static inline void apply(MultiGeometry& multigeometry,
140 RangeOrPoint const& range_or_point,
141 signed_size_type ring_index, signed_size_type multi_index)
142 {
143 Policy::template apply
144 <
145 typename boost::range_value<MultiGeometry>::type,
146 RangeOrPoint
147 >(range::at(multigeometry, multi_index), range_or_point, ring_index);
148 }
149 };
150
151
152 }} // namespace detail::append
153 #endif // DOXYGEN_NO_DETAIL
154
155
156 #ifndef DOXYGEN_NO_DISPATCH
157 namespace dispatch
158 {
159
160 template
161 <
162 typename Geometry,
163 typename RangeOrPoint,
164 typename Tag = typename geometry::tag<Geometry>::type,
165 typename OtherTag = typename geometry::tag<RangeOrPoint>::type
166 >
167 struct append
168 : detail::append::append_no_action
169 {};
170
171 template <typename Geometry, typename Point>
172 struct append<Geometry, Point, linestring_tag, point_tag>
173 : detail::append::to_range_point
174 {};
175
176 template <typename Geometry, typename Point>
177 struct append<Geometry, Point, ring_tag, point_tag>
178 : detail::append::to_range_point
179 {};
180
181 template <typename Polygon, typename Point>
182 struct append<Polygon, Point, polygon_tag, point_tag>
183 : detail::append::to_polygon_point
184 {};
185
186 template <typename Geometry, typename Range, typename RangeTag>
187 struct append<Geometry, Range, linestring_tag, RangeTag>
188 : detail::append::to_range_range
189 {};
190
191 template <typename Geometry, typename Range, typename RangeTag>
192 struct append<Geometry, Range, ring_tag, RangeTag>
193 : detail::append::to_range_range
194 {};
195
196 template <typename Polygon, typename Range, typename RangeTag>
197 struct append<Polygon, Range, polygon_tag, RangeTag>
198 : detail::append::to_polygon_range
199 {};
200
201
202 template <typename Geometry, typename Point>
203 struct append<Geometry, Point, multi_point_tag, point_tag>
204 : detail::append::to_range_point
205 {};
206
207 template <typename Geometry, typename Range, typename RangeTag>
208 struct append<Geometry, Range, multi_point_tag, RangeTag>
209 : detail::append::to_range_range
210 {};
211
212 template <typename MultiGeometry, typename Point>
213 struct append<MultiGeometry, Point, multi_linestring_tag, point_tag>
214 : detail::append::to_multigeometry<detail::append::to_range_point>
215 {};
216
217 template <typename MultiGeometry, typename Range, typename RangeTag>
218 struct append<MultiGeometry, Range, multi_linestring_tag, RangeTag>
219 : detail::append::to_multigeometry<detail::append::to_range_range>
220 {};
221
222 template <typename MultiGeometry, typename Point>
223 struct append<MultiGeometry, Point, multi_polygon_tag, point_tag>
224 : detail::append::to_multigeometry<detail::append::to_polygon_point>
225 {};
226
227 template <typename MultiGeometry, typename Range, typename RangeTag>
228 struct append<MultiGeometry, Range, multi_polygon_tag, RangeTag>
229 : detail::append::to_multigeometry<detail::append::to_polygon_range>
230 {};
231
232
233 template <typename Geometry, typename RangeOrPoint, typename OtherTag>
234 struct append<Geometry, RangeOrPoint, dynamic_geometry_tag, OtherTag>
235 {
applyboost::geometry::dispatch::append236 static inline void apply(Geometry& geometry,
237 RangeOrPoint const& range_or_point,
238 signed_size_type ring_index, signed_size_type multi_index)
239 {
240 traits::visit<Geometry>::apply([&](auto & g)
241 {
242 append
243 <
244 std::remove_reference_t<decltype(g)>, RangeOrPoint
245 >::apply(g, range_or_point, ring_index, multi_index);
246 }, geometry);
247 }
248 };
249
250 // TODO: It's unclear how append should work for GeometryCollection because
251 // it can hold multiple different geometries.
252
253 } // namespace dispatch
254 #endif // DOXYGEN_NO_DISPATCH
255
256
257 /*!
258 \brief Appends one or more points to a linestring, ring, polygon, multi-geometry
259 \ingroup append
260 \tparam Geometry \tparam_geometry
261 \tparam RangeOrPoint Either a range or a point, fullfilling Boost.Range concept or Boost.Geometry Point Concept
262 \param geometry \param_geometry
263 \param range_or_point The point or range to add
264 \param ring_index The index of the ring in case of a polygon:
265 exterior ring (-1, the default) or interior ring index
266 \param multi_index The index of the geometry to which the points are appended
267
268 \qbk{[include reference/algorithms/append.qbk]}
269 }
270 */
271 template <typename Geometry, typename RangeOrPoint>
append(Geometry & geometry,RangeOrPoint const & range_or_point,signed_size_type ring_index=-1,signed_size_type multi_index=0)272 inline void append(Geometry& geometry, RangeOrPoint const& range_or_point,
273 signed_size_type ring_index = -1, signed_size_type multi_index = 0)
274 {
275 concepts::check<Geometry>();
276
277 dispatch::append
278 <
279 Geometry, RangeOrPoint
280 >::apply(geometry, range_or_point, ring_index, multi_index);
281 }
282
283
284 }} // namespace boost::geometry
285
286
287 #endif // BOOST_GEOMETRY_ALGORITHMS_APPEND_HPP
288