如何使用boost :: geometry :: rtree和glm :: vec3作为自定义点类型?

时间:2017-12-05 12:13:44

标签: c++ templates boost glm-math boost-geometry

我正在使用

BOOST_GEOMETRY_REGISTER_POINT_3D(glm::vec3, float, boost::geometry::cs::cartesian, x, y, z);

将RTree定义为:

  using IndexedPoint = std::pair<glm::vec3, uint32_t>;
  using RTree = boost::geometry::index::rtree<IndexedPoint, boost::geometry::index::rstar<8>>;

当我尝试使用它运行最近邻居查询时,它无法编译:

auto it = rtree.qbegin(boost::geometry::index::nearest(glm::vec3(), 3))

错误是:

error C2664: 'int boost::mpl::assertion_failed<false>(boost::mpl::assert<false>::type)': cannot convert argument 1 from 'boost::mpl::failed ************(__cdecl boost::geometry::strategy::distance::services::default_strategy<boost::geometry::point_tag,boost::geometry::box_tag,glm::vec<3,float,0>,boost::geometry::model::point<float,3,boost::geometry::cs::cartesian>,boost::geometry::cartesian_tag,boost::geometry::cartesian_tag,void>::NOT_IMPLEMENTED_FOR_THIS_POINT_TYPE_COMBINATION::* ***********)(boost::mpl::assert_::types<Point1,Point2,CsTag1,CsTag2>)' to 'boost::mpl::assert<false>::type'
        with
        [
            Point1=glm::vec<3,float,0>,
            Point2=boost::geometry::model::point<float,3,boost::geometry::cs::cartesian>,
            CsTag1=boost::geometry::cartesian_tag,
            CsTag2=boost::geometry::cartesian_tag
        ]

似乎Matched_distance_result缺少vec3 vs boost :: geometry :: model :: point和boost :: geometry :: model :: box的特化。我曾尝试手动添加它们,但无法使其正常工作。如何添加所需的距离类型专精?

请注意,我可以使用相同的设置进行空间查询,因此看起来基本上是正确的。

2 个答案:

答案 0 :(得分:3)

我无法用GCC / Boost 1.65.1重现问题:

<强> Live¹ On Coliru

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

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

#include <glm/vec3.hpp>
BOOST_GEOMETRY_REGISTER_POINT_3D(glm::vec3, float, bg::cs::cartesian, x, y, z)

#include <iostream>
int main() {

    using IndexedPoint = std::pair<glm::vec3, uint32_t>;
    using RTree = boost::geometry::index::rtree<IndexedPoint, boost::geometry::index::rstar<8>>;

    RTree rtree;
    rtree.insert({glm::vec3(1,1,1), 1});
    rtree.insert({glm::vec3(2,2,2), 2});
    rtree.insert({glm::vec3(3,3,3), 3});
    rtree.insert({glm::vec3(4,4,4), 4});

    auto it = rtree.qbegin(bgi::nearest(glm::vec3(2.9, 2.9, 2.9), 99));

    auto p = it->first;
    std::cout << "Nearest: # " << it->second << " (" << p.x << ", " << p.y << " " << p.z << ")\n";
}

打印

Nearest: # 3 (3, 3 3)

¹Coliru没有libglm

答案 1 :(得分:0)

我只是想从对已接受答案的评论中(@BuschnicK)拉出实际上为我解决了这个问题的答案。

  

问题是我没有包括geometry.hpp,只有rtree,box和point。

在我的标题中添加以下包含项解决了上面显示的问题。

#include <boost/geometry.hpp>