题
为什么查询box
是否在rtree
中的查询返回空结果?
说明
拿线串
LINESTRING(1 1, 2 2)
和两个多边形
POLYGON((0 0,0 1,1 1,1 0,0 0))
POLYGON((0 0,0 3,3 3,3 0,0 0))
我可以使用
boost::geometry::witin()
查询线串是否在任意一个多边形内(并且有效)但是,如果创建多边形的
rtree
,并在box
周围放置linestring
,则查询box
是否为within
时,rtree
返回空结果。例
这是同时显示有效
bg::within(line, polygon)
和无效rtree.query(bgi::within( line_box1 ), ...)
的示例void rtree_within() {
typedef bgm::point< double, 2, bg::cs::cartesian > point;
typedef bgm::box<point> box;
typedef bgm::linestring<point> line;
typedef bgm::polygon<point> polygon;
typedef std::pair<box, unsigned> value;
bgi::rtree<value, bgi::quadratic<16> > rtree;
std::vector<value> result_s;
polygon poly1;
polygon poly2;
line line1;
bg::read_wkt("POLYGON((0 0,0 1,1 1,1 0,0 0))", poly1);
bg::read_wkt("POLYGON((0 0,0 3,3 3,3 0,0 0))", poly2);
bg::read_wkt("LINESTRING(1 1, 2 2)", line1);
std::cout << "line2 in poly1: " << bg::within(line1, poly1) << std::endl;
std::cout << "line2 in poly3: " << bg::within(line1, poly2) << std::endl;
// boxes to insert into rtree
box poly_box1 = bg::return_envelope<box>( poly1 );
rtree.insert(std::make_pair(poly_box1, 0));
box poly_box2 = bg::return_envelope<box>( poly2 );
rtree.insert(std::make_pair(poly_box2, 2));
// box around the line
box line_box1 = bg::return_envelope<box>( line1 );
std::cout << "poly_box1: " << bg::wkt( poly_box1 ) << std::endl; // returns 0
std::cout << "poly_box2: " << bg::wkt( poly_box2 ) << std::endl; // returns 1
std::cout << "line_box1: " << bg::wkt( line_box1 ) << std::endl;
rtree.query(bgi::within( line_box1 ), std::back_inserter( result_s ));
std::cout << "line_box1 within rtree - size: " << result_s.size() << std::endl;
// result_s is empty (size == 0)
}
最佳答案
调用自由函数within(geom1,geom2)
时,如果geom1
在geom2
内部,则返回TRUE。但是,当您阅读有关使用within
作为谓词的参考时
生成定义值和几何关系的谓词。如果bg :: within(Indexable,Geometry)返回true,则查询将返回值。
因此,您尝试检查indexable
是否在geometry
内,并且结果为假[可索引为rect,几何为线]。您应该使用contains
谓词而不是使用rtree的within
。 contains
谓词与bg::within(Geometry, Indexable)
等效,那么您的代码将按预期工作。