Question

我试图评估几个RTree实现并遇到了boost::geometry::rtree::contains（升级版本1.55）的问题。我正在做的是尝试获取包含一个点的所有盒子的列表。它返回正确的方框，但它会多次返回相同的方框。我不确定为什么会这样，libspatialindex不会这样做。

这是我的代码：

void testBoostRTree(const Polygons& polygons)
{
    using namespace boost::geometry;

    typedef model::point<double, 2, cs::cartesian> BoostPoint;
    typedef model::box<BoostPoint> BoostBox;
    typedef model::polygon<BoostPoint, true, true> BoostPolygon; // clockwise, closed.
    typedef std::pair<BoostBox, unsigned int> RTreeValue;

    index::rtree<RTreeValue, index::rstar<16, 4>> rtree;
    std::vector<BoostPoint> centrePoints;

    for (const auto& p : polygons)
    {
        BoostPolygon bp;

        for (const auto& point : p.m_points)
        {
             bp.outer().push_back(BoostPoint(point.first, point.second));
        } 

        BoostBox box = return_envelope<BoostBox>(bp);
        rtree.insert(std::make_pair(box, p.m_id));

        centrePoints.push_back(return_centroid<BoostPoint>(box));
    }

    std::vector<RTreeValue> hits;
    for (const auto& cp : centrePoints)
    {
         std::cout << "* Query point: "  << get<0>(cp) << ", " << get<1>(cp) << "\n";
         hits.clear();

         rtree.query(index::contains(cp), std::back_inserter(hits));

         for (const auto& r : hits)
         {
             std::cout << r.second << "\n";
         }
    }
}

我已经确认要检查的点是正确的。为了使用std::set而不是std::vector，我也不想真正做迭代或过时。

以下是一些示例结果：

* Query point: 51.4181, 0.20462
278566
278566
278566
278566
278566
278566
261819
261821
261819
261820
261820
261821
13741
278566
278566
...

将这些与libspatialindex：

的结果进行比较

* Query 51.4181 0.20462 
261819
261820
261821
13741
278566

我已经搜索了文档和源代码，但我找不到任何明显的错误。

Answer 1

由于你没有写出准确插入rtree的内容，我无法重现这个问题。

在当前实现的rtree（Boost 1.55）版本中，您插入的每个值都被编入索引。我想你可以认为它更像std :: multiset而不是std :: set。因此，如果由于某种原因生成相同的值并将其插入到rtree中，您可能会注意到输出中的重复值。

boost :: geometry :: rtree :: contains返回多个相同的结果

1 个答案: