提示框的交叉點給出與段的錯誤交集

Question

Boost rtree給出了與段查詢的某些交集的錯誤交集結果。 在這種情況下，邊界框在y = 0時是一個y平面的10×10平方。我用從（2，1，0）到（2，1，10）的z對齊線查詢。有趣的是，如果我用查詢框而不是段，那麼它按預期工作。當盒子不是平面的時候，這個行爲也是存在的，只是將最小角落移動到（0，-5，0）並且它仍然發生。

我使用這個錯誤還是它在提升錯誤？

編輯：在Boost 1.56和1.59上試過了。

#include <vector> 
#include "gtest/gtest.h" 
#include "gmock/gmock.h" 
#include <boost/geometry.hpp> 
#include <boost/geometry/geometries/point.hpp> 
#include <boost/geometry/geometries/box.hpp> 
#include <boost/geometry/geometries/segment.hpp> 
#include <boost/geometry/index/rtree.hpp> 
#include <vector> 
#include <iterator> 
#include <memory> 

namespace bg = boost::geometry; 
namespace bgi = boost::geometry::index; 
typedef bg::model::point<double, 3, bg::cs::cartesian> point_def; 
typedef bg::model::box<point_def> box; 
typedef bg::model::segment<point_def> segment; 
typedef std::pair<box, size_t> tri_box; 
typedef bgi::rtree< tri_box, bgi::linear<8>> tree_type; 

using namespace std; 

TEST(boost_rtree, cant_intersect_box_with_segment) { 
    vector<tri_box> buff(1); 
    buff[0].first = box{point_def{0, 0, 0}, point_def{10, 0, 10}}; 
    buff[0].second = 1; 
    tree_type tree(buff); 

    segment query{point_def{2, 1, 0}, point_def{2, 1, 10}}; 
// box query{point_def{2, 1, 0}, point_def{2, 1, 10}}; 
    vector<tri_box> out; 

    size_t count = tree.query(bgi::intersects(query), back_inserter(out)); 

    ASSERT_EQ(0, count); // fails here 
    ASSERT_EQ(0, out.size()); 
} 

編輯：問題正在採取行動，擴大郵件列表：lists.boost.org/geometry/2015/09/3472.php

Answer 1

不太可能，因爲它似乎，這在我看來是一個錯誤。

甚至編譯這個的第一個版本是Boost 1.56。所有以前的版本中失敗，並

BOOST_MPL_ASSERT_MSG 
    (
     false, NOT_OR_NOT_YET_IMPLEMENTED_FOR_THIS_GEOMETRY_TYPE 
     , (types<Geometry>) 
    ); 

但是，即使代碼被編譯，它似乎並不正確......：在intersects呼叫構成了這一查詢謂詞本身返回「假陽性」看來。

大大簡化：Live On Coliru

#include <boost/geometry.hpp> 
#include <boost/geometry/geometries/box.hpp> 
#include <boost/geometry/geometries/segment.hpp> 

namespace bg = boost::geometry; 

typedef bg::model::point<int, 3, bg::cs::cartesian> point; 
typedef bg::model::box<point>  box; 
typedef bg::model::segment<point> segment; 

int main() { 
    box y0rect = box{point{0, 0, 0}, point{10, 0, 10}}; 
    segment seg{point{2, 1, 0}, point{2, 1, 10}}; 

    bg::correct(y0rect); 
    bg::correct(seg); 
    assert(!bg::intersects(seg, y0rect)); 
} 

UPDATE

有趣的是，它似乎正常工作有時爲2D。我不知道結果是不是簡單的不確定......

Live On Coliru

#include <boost/geometry.hpp> 
#include <boost/geometry/geometries/box.hpp> 
#include <boost/geometry/geometries/segment.hpp> 

namespace bg = boost::geometry; 

typedef bg::model::point<int, 4, bg::cs::cartesian> point; 
typedef bg::model::box<point>  box; 
typedef bg::model::segment<point> segment; 

int main() { 
    box y0rect = box{point{0, 0}, point{10, 10}}; 
    bg::correct(y0rect); 

    { 
     segment seg{point{12, 0}, point{20, 10}}; 
     bg::correct(seg); 
     assert(!bg::intersects(seg, y0rect)); 
    } 
    { 
     segment seg{point{2, 0}, point{8, 6}}; 
     bg::correct(seg); 
     assert(bg::intersects(seg, y0rect)); 
    } 
    { 
     segment seg{point{2, 0}, point{18, 6}}; 
     bg::correct(seg); 
     assert(bg::intersects(seg, y0rect)); // OOPS BREAKS? 
    } 
}