boost/geometry
2
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2026-02-11 11:52:11 +00:00
Code Issues Projects Releases Wiki Activity

relate(L,L) recently introduced turn position info used to optimize boundary points detection

Browse Source
This commit is contained in:
Adam Wulkiewicz
2014-02-25 17:37:00 +01:00
parent 736f26db4a
commit 37db735b38
1 changed files with 69 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

86
include/boost/geometry/algorithms/detail/relate/linear_linear.hpp
View File

@@ -571,13 +571,21 @@ struct linear_linear
// interiors overlaps
update<interior, interior, '1', transpose_result>(res);
bool this_b = is_ip_on_boundary<boundary_front>(it->point,
it->operations[op_id],
boundary_checker,
seg_id);
// going inside on boundary point
if ( boundary_checker.template is_boundary<boundary_front>(it->point, seg_id) )
// may be front only
if ( this_b )
{
bool other_b =
it->operations[other_op_id].operation == overlay::operation_blocked ?
other_boundary_checker.template is_endpoint_boundary<boundary_back>(it->point) :
other_boundary_checker.template is_boundary<boundary_any>(it->point, other_id);
// may be front and back
bool other_b = is_ip_on_boundary<boundary_any>(it->point,
it->operations[other_op_id],
other_boundary_checker,
other_id);
// it's also the boundary of the other geometry
if ( other_b )
{
@@ -623,10 +631,11 @@ struct linear_linear
// check if this is indeed the boundary point
if ( boundary_checker.template is_endpoint_boundary<boundary_back>(it->point) )
{
bool other_b =
it->operations[other_op_id].operation == overlay::operation_blocked ?
other_boundary_checker.template is_endpoint_boundary<boundary_back>(it->point) :
other_boundary_checker.template is_boundary<boundary_any>(it->point, other_id);
// may be front and back
bool other_b = is_ip_on_boundary<boundary_any>(it->point,
it->operations[other_op_id],
other_boundary_checker,
other_id);
// it's also the boundary of the other geometry
if ( other_b )
{
@@ -663,15 +672,15 @@ struct linear_linear
// method other than crosses, check more conditions
else
{
bool this_b =
op_blocked ?
boundary_checker.template is_endpoint_boundary<boundary_back>(it->point) :
boundary_checker.template is_boundary<boundary_front>(it->point, seg_id);
bool this_b = is_ip_on_boundary<boundary_any>(it->point,
it->operations[op_id],
boundary_checker,
seg_id);
bool other_b =
it->operations[other_op_id].operation == overlay::operation_blocked ?
other_boundary_checker.template is_endpoint_boundary<boundary_back>(it->point) :
other_boundary_checker.template is_boundary<boundary_any>(it->point, other_id);
bool other_b = is_ip_on_boundary<boundary_any>(it->point,
it->operations[other_op_id],
other_boundary_checker,
other_id);
// if current IP is on boundary of the geometry
if ( this_b )
@@ -738,6 +747,49 @@ struct linear_linear
}
}
template <boundary_query BoundaryQuery,
typename IntersectionPoint,
typename OperationInfo,
typename BoundaryChecker>
static inline
bool is_ip_on_boundary(IntersectionPoint const& ip,
OperationInfo const& operation_info,
BoundaryChecker & boundary_checker,
segment_identifier const& seg_id)
{
boost::ignore_unused_variable_warning(seg_id);
bool res = false;
if ( (BoundaryQuery == boundary_back || BoundaryQuery == boundary_any)
&& operation_info.operation == overlay::operation_blocked )
{
BOOST_ASSERT(operation_info.position == overlay::position_back);
res = boundary_checker.template is_endpoint_boundary<boundary_back>(ip);
#ifdef BOOST_GEOMETRY_DEBUG_RELATE_LINEAR_LINEAR
BOOST_ASSERT(res == boundary_checker.template is_boundary<boundary_back>(ip, seg_id));
#endif
}
else if ( (BoundaryQuery == boundary_front || BoundaryQuery == boundary_any)
&& operation_info.position == overlay::position_front )
{
res = boundary_checker.template is_endpoint_boundary<boundary_front>(ip);
#ifdef BOOST_GEOMETRY_DEBUG_RELATE_LINEAR_LINEAR
BOOST_ASSERT(res == boundary_checker.template is_boundary<boundary_front>(ip, seg_id));
#endif
}
else
{
#ifdef BOOST_GEOMETRY_DEBUG_RELATE_LINEAR_LINEAR
BOOST_ASSERT(res == boundary_checker.template is_boundary<boundary_any>(ip, seg_id));
#endif
}
return res;
}
private:
exit_watcher<turn_point_type> m_exit_watcher;
segment_watcher m_seg_watcher;