无法通过其索引

时间:2016-01-22 16:27:35

标签: c++ boost r-tree boost-geometry

我正在创建boost geometry rtree

在我创建的示例页面中:

typedef bg::model::point<float, 2, bg::cs::cartesian> point;
typedef bg::model::box<point> box;
typedef std::pair<box, unsigned> value;

我想在树中添加value元素:

bgi::rtree< value, bgi::quadratic<16> > rtree;
box b(point(2.0f, 2.0f), point(2.5f, 2.5f));
// insert new value
rtree.insert(std::make_pair(b, 2));

现在我想知道是否可以通过知道box标识符(或者如果可以设置更多具有相同id的元素)来删除元素。

我想要做的就是删除我在上面添加的元素,比如调用类似的东西:

rtree.remove(2); // why I can't do this?

我该如何做到这一点?

1 个答案:

答案 0 :(得分:2)

您必须使用迭代器。可能的重载是

remove(value_type const &)
remove(Iterator, Iterator)
remove(ConvertibleOrRange const &)

首先,您应该通过搜索找到要删除的值(可能会运行查询,但我认为您不需要这样做,因为您的问题不存在。)

因此,您可能希望使用std::find_if。我不确定如何改进性能(不使用反向映射)。我想你的探查器应该用于在那里提供保证。

更新:概念证明

天真的方法失败了:

// UNDEFINED BEHAVIOUR:
template <typename Rtree, typename Id>
size_t remove_ids_loop(Rtree& rtree, Id const& id) {
    using V = typename Rtree::value_type;
    static_assert(sizeof(V) == 0, "don't use; UNDEFINED BEHAVIOUR!");

    size_t removed = 0;
    std::for_each(rtree.begin(), rtree.end(), [&](V const& v) { 
            if (id == v.second)
                removed += rtree.remove(v);
        });

    return removed;
}

这是因为as documented

  

⚠警告

     

修改rtree可能会使迭代器失效。

简单方法:

template <typename Rtree, typename Id>
size_t remove_ids_bulk(Rtree& rtree, Id const& id) {
    using V = typename Rtree::value_type;
    std::vector<V> v;
    std::copy_if(rtree.begin(), rtree.end(), back_inserter(v), [id](V const& v) { return v.second == id; });

    return rtree.remove(v.begin(), v.end());
}

通过与删除分开执行查询来避开竞争。请注意,这可能更有效,因为我们会立即传递整个范围。

  

注意:对于重复的条目,这应该没问题,因为rtree::remove一次只删除一个值。

现场演示

<强> Live On Coliru

#include <boost/geometry.hpp>
#include <boost/geometry/index/rtree.hpp>

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

typedef bg::model::point<float, 2, bg::cs::cartesian> point;
typedef bg::model::box<point> box;

typedef std::pair<box, unsigned> value;

template <typename Rtree, typename Id>
size_t remove_ids_bulk(Rtree& rtree, Id const& id) {
    using V = typename Rtree::value_type;
    std::vector<V> v;
    std::copy_if(rtree.begin(), rtree.end(), back_inserter(v), [id](V const& v) { return v.second == id; });

    return rtree.remove(v.begin(), v.end());
}

// UNDEFINED BEHAVIOUR:
template <typename Rtree, typename Id>
size_t remove_ids_loop(Rtree& rtree, Id const& id) {
    using V = typename Rtree::value_type;
    static_assert(sizeof(V) == 0, "don't use; UNDEFINED BEHAVIOUR!");

    size_t removed = 0;
    std::for_each(rtree.begin(), rtree.end(), [&](V const& v) { 
            if (id == v.second)
                removed += rtree.remove(v);
        });

    return removed;
}

int main() {
    //I want to add value elements to the tree:

    bgi::rtree< value, bgi::quadratic<16> > rtree;
    box b(point(2.0f, 2.0f), point(2.5f, 2.5f));
    // insert new value
    rtree.insert(std::make_pair(b, 2));
    std::cout << "Elements: " << rtree.size() << "\n";

    // remove id 2
    remove_ids_bulk(rtree, 2u);
    std::cout << "Elements: " << rtree.size() << "\n";
}

打印

Elements: 1
Elements: 0