1 // Boost.Geometry (aka GGL, Generic Geometry Library)
3 // Copyright (c) 2012-2014 Barend Gehrels, Amsterdam, the Netherlands.
5 // Use, modification and distribution is subject to the Boost Software License,
6 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
7 // http://www.boost.org/LICENSE_1_0.txt)
9 #ifndef BOOST_GEOMETRY_STRATEGIES_CARTESIAN_BUFFER_JOIN_MITER_HPP
10 #define BOOST_GEOMETRY_STRATEGIES_CARTESIAN_BUFFER_JOIN_MITER_HPP
12 #include <boost/geometry/core/assert.hpp>
13 #include <boost/geometry/core/cs.hpp>
14 #include <boost/geometry/policies/compare.hpp>
15 #include <boost/geometry/util/math.hpp>
16 #include <boost/geometry/util/select_most_precise.hpp>
18 #include <boost/geometry/strategies/buffer.hpp>
21 namespace boost { namespace geometry
24 namespace strategy { namespace buffer
28 \brief Let the buffer create sharp corners
30 \details This strategy can be used as JoinStrategy for the buffer algorithm.
31 It creates a sharp corners around each convex vertex. It can be applied
32 for (multi)linestrings and (multi)polygons.
33 If corners are sharp by themselves, the miters might become very long. Therefore
34 there is a limit (miter_limit), in terms of the used distance, which limits
35 their length. The miter is not changed to a bevel form (as done in some
36 other software), it is just adapted to the specified miter_limit but keeps
38 If the buffer distance is 5.0, and the miter limit is 2.0, generated points
39 will be located at a distance of at most 10.0 (2*5) units.
40 This strategy is only applicable for Cartesian coordinate systems.
46 [$img/strategies/buffer_join_miter.png]
48 \* [link geometry.reference.algorithms.buffer.buffer_7_with_strategies buffer (with strategies)]
49 \* [link geometry.reference.strategies.strategy_buffer_join_round join_round]
56 //! \brief Constructs the strategy
57 //! \param miter_limit The miter limit, to avoid excessively long miters around sharp corners
58 explicit inline join_miter(double miter_limit = 5.0)
59 : m_miter_limit(valid_limit(miter_limit))
62 #ifndef DOXYGEN_SHOULD_SKIP_THIS
63 //! Fills output_range with a sharp shape around a vertex
64 template <typename Point, typename DistanceType, typename RangeOut>
65 inline bool apply(Point const& ip, Point const& vertex,
66 Point const& perp1, Point const& perp2,
67 DistanceType const& buffer_distance,
68 RangeOut& range_out) const
70 geometry::equal_to<Point> equals;
71 if (equals(ip, vertex))
75 if (equals(perp1, perp2))
80 typedef typename coordinate_type<Point>::type coordinate_type;
81 typedef typename geometry::select_most_precise
85 >::type promoted_type;
89 // Check the distance ip-vertex (= miter distance)
90 // (We calculate it manually (not using Pythagoras strategy) to reuse
92 coordinate_type const dx = get<0>(p) - get<0>(vertex);
93 coordinate_type const dy = get<1>(p) - get<1>(vertex);
95 promoted_type const distance = geometry::math::sqrt(dx * dx + dy * dy);
97 promoted_type const max_distance
98 = m_miter_limit * geometry::math::abs(buffer_distance);
100 if (distance > max_distance)
102 BOOST_GEOMETRY_ASSERT(distance != 0.0);
104 promoted_type const proportion = max_distance / distance;
105 set<0>(p, get<0>(vertex) + dx * proportion);
106 set<1>(p, get<1>(vertex) + dy * proportion);
109 range_out.push_back(perp1);
110 range_out.push_back(p);
111 range_out.push_back(perp2);
115 template <typename NumericType>
116 inline NumericType max_distance(NumericType const& distance) const
118 return distance * m_miter_limit;
121 #endif // DOXYGEN_SHOULD_SKIP_THIS
124 double valid_limit(double miter_limit) const
126 if (miter_limit < 1.0)
128 // It should always exceed the buffer distance
134 double m_miter_limit;
137 }} // namespace strategy::buffer
140 }} // namespace boost::geometry
142 #endif // BOOST_GEOMETRY_STRATEGIES_CARTESIAN_BUFFER_JOIN_MITER_HPP