Boost几何体内用户定义的空间谓词

Question

我需要确定两个几何体的内部是否相交。一个InteriorsIntersect谓语不被OGC或升压几何规定，而是由DE-9IM矩阵（see also）定义：

T * * 
* * * 
* * * 

我已经使用升压几何的relate功能创建了自己的谓语。

namespace bgr = boost::geometry::detail::relate; 
using InteriorsIntersectMask = bgr::static_mask<'T','*','*','*','*','*','*','*','*'>; 

template<typename Geom1, typename Geom2> 
inline bool interiorsIntersect(const Geom1& geom1, const Geom2& geom2) 
{ 
    return bgr::relate<InteriorsIntersectMask>(geom1, geom2); 
} 

This works great。我唯一担心的是relate函数和static_mask类型没有作为Boost Geometry API的一部分记录下来，并且就我所知道的实现细节而言。以这种方式使用relate安全吗？使用Boost几何实现相同的目标是否有其他选择？理想情况下，我希望relate是boost/geometry/algorithms内的算法。

Answer 1

relate()和relation()函数计划在Boost 1.59中发布。接口比在问题中提到的一个silghtly不同：

namespace bg = boost::geometry; 
using II = bg::de9im::static_mask<'T','*','*','*','*','*','*','*','*'>; 
bool check1 = bg::relate(geom1, geom2, II()); 

bg::de9im::mask ii("T********"); 
bool check2 = bg::relate(geom1, geom2, ii); 

bg::de9im::matrix m = bg::relation(geom1, geom2); 
std::cout << m.str(); 

也可以通过更加复杂的面具：

bg::de9im::mask ii1("1********"); 
bg::de9im::mask ii2("2********"); 
// check if the intersection of interiors is linear or areal 
bool check2 = bg::relate(geom1, geom2, ii1 || ii2);