C++ 在几何体中使用boost rtree查找结果点

我试图在我的一个项目中使用boost:：geometry的rtree DS，但我发现很难浏览文档。某些方法的文档记录很差，我找不到足够的例子。现在我正在尝试构建一个示例程序，以便进一步构建它

在下面的例子中，我有一个点和一个盒子的rtree，我需要找到盒子里的所有点。我想问的另一件事是，我找不到打包算法构造函数或方法，所以如何使用它。这是我到现在为止所做的-

#include <iostream>
#include <vector>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/register/point.hpp>
#include <boost/geometry/geometries/register/box.hpp>
#include <boost/geometry/index/rtree.hpp>


namespace bg = boost::geometry;
namespace bgi = boost::geometry::index;


struct my_point
{
    float x, y;
    my_point(float _x, float _y) : x(_x), y(_y) {}
};

struct my_box
{
    my_point ll, ur;
    my_box(float x1, float y1, float x2, float y2) : ll(x1,y1), ur(x2,y2) {}
};

// Register the point type
BOOST_GEOMETRY_REGISTER_POINT_2D(my_point, float, cs::cartesian, x, y)

// Register the box type, also notifying that it is based on "my_point"
BOOST_GEOMETRY_REGISTER_BOX(my_box, my_point, ll, ur)



int main()
{
    std::vector<std::pair<my_point, int>> pts;
    pts.emplace_back(std::make_pair(my_point(2,2), 5));
    pts.emplace_back(std::make_pair(my_point(3,3), 1));
    pts.emplace_back(std::make_pair(my_point(4,5), 3));
    pts.emplace_back(std::make_pair(my_point(4,4), 12));
    pts.emplace_back(std::make_pair(my_point(1,2), 50));
    // ....

    bgi::rtree<std::pair<my_point, int>, bgi::dynamic_rstar> rT(bgi::dynamic_rstar(pts.size()));
    rT.insert(pts.begin(), pts.end());

    my_box box1(1,1,4,4);
    // how to retrieve all points or their .second inside this box?

    return 0;
}

您需要在rtree上使用查询方法，并传递适当的查询策略，例如-

std::vector<std::pair<my_point, int>> result;
rT.query(bgi::covered_by(box1), std::back_inserter(result));

对于打包算法，虽然我不确定，但这应该有效-

bgi::rtree<std::pair<my_point, int>, bgi::dynamic_rstar> rT(pts, bgi::dynamic_rstar(pts.size()));

但是，我建议根据文档进行验证。这是我在我的机器上本地测试的完整程序-

#include <iostream>
#include <vector>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/register/point.hpp>
#include <boost/geometry/geometries/register/box.hpp>
#include <boost/geometry/index/rtree.hpp>


namespace bg = boost::geometry;
namespace bgi = boost::geometry::index;


struct my_point
{
    float x, y;
    my_point(float _x, float _y) : x(_x), y(_y) {}
};

struct my_box
{
    my_point ll, ur;
    my_box(float x1, float y1, float x2, float y2) : ll(x1,y1), ur(x2,y2) {}
};

// Register the point type
BOOST_GEOMETRY_REGISTER_POINT_2D(my_point, float, cs::cartesian, x, y)

// Register the box type, also notifying that it is based on "my_point"
BOOST_GEOMETRY_REGISTER_BOX(my_box, my_point, ll, ur)



int main()
{
    std::vector<std::pair<my_point, int>> pts;
    pts.emplace_back(std::make_pair(my_point(2,2), 5));
    pts.emplace_back(std::make_pair(my_point(3,3), 1));
    pts.emplace_back(std::make_pair(my_point(4,5), 3));
    pts.emplace_back(std::make_pair(my_point(4,4), 12));
    pts.emplace_back(std::make_pair(my_point(1,2), 50));
    // ....

    bgi::rtree<std::pair<my_point, int>, bgi::dynamic_rstar> rT(pts, bgi::dynamic_rstar(pts.size()));

    my_box box1(1,1,4,4);
    std::vector<std::pair<my_point, int>> result;
    rT.query(bgi::covered_by(box1), std::back_inserter(result));
    for(const auto &r: result){
        std::cout << r.second << ' ' << '\n';
    }
    return 0;
}

---

如果将范围传递给R树的构造函数，则使用打包算法：

可以使用查询方法或查询迭代器从R树检索值：

另见此

最后一句话。在代码中，将点数传递给bgi:：dynamic\u rstar。bgi:：dynamic_rstar的参数不是包含的值的数量。它是与持久/基于磁盘的R树的页面大小相对应的R树的单个节点中存储的最大值数。因此，如果传递bgi:：dynamic_rstarpts.size，R树只包含一个包含所有点的节点，这意味着没有节点层次结构，没有空间分区，并且rtree的行为实际上就像一个向量。如果值重叠较大，则使用例如bgi:：rstar或更大的数字，但在您的情况下不应该使用，因为您正在存储点。这就足够了：

bgi::rtree<std::pair<my_point, int>, bgi::rstar<4> > rT(pts);