我正在创建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?
我该如何做到这一点?
答案 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