2d aabbtree в CGAL с пользовательской картой свойств ⇐ C++

[Anonymous]

Я пытаюсь создать aabbtree для проверки, в каком треугольнике 2D -сетки некоторые точки. Для этого я использую CGAL, в частности, пакет Polygon_mesh_processing и функции местоположения. Тем не менее, это кажется невозможной задачей для меня. Для этого мне нужна функция pmp :: build_aabb_tree , которая получает сетку и aabbtree .
В этом конкретном случае сетка имеет тип cgal :: surface_mesh , с Simple_carters kernel. Точки запроса также 2D. Согласно документации, мне нужно предоставить читаемый PropertyMap, который преобразует мои 2D -точки в 3D -точки. < /P>
Я делаю это со следующим: < /p>
struct Point3VPM {
using key_type = Mesh::Vertex_index;
using value_type = VNCS::Sim2D::Kernel::K::Point_3;
using reference = value_type;
using category = boost::readable_property_map_tag;

Point3VPM(const Mesh &m)
: mesh(std::cref(m))
{
}

friend Point3VPM::value_type get(const Point3VPM &map, Mesh::Vertex_index idx)
{
Point p = map.mesh.get().point(idx);
return {p[0], p[1], 0};
}

std::reference_wrapper mesh;
};
< /code>
Здесь начинает мою первую проблему. Мой тип aabbtree - следующее: < /p>
using AABBTreeTraits =
CGAL::AABB_traits;
using AABBTree = CGAL::AABB_tree;

Как я могу построить свой aabbtree ? Моя точка3VPM нуждается в Surfacemesh , чтобы иметь возможность преобразовать vertex_index в point_3 , но я не могу найти способ создания aabbtreetraits с помощью surfacemesh . CGAL ожидает конструктора по умолчанию для карты свойств, но, поскольку моя карта свойств не может работать без сетки, мне нужно как -то предоставить сетку. Как я могу это исправить?#include
#include
#include
#include
#include
#include

namespace
{
using K = CGAL::Simple_cartesian;
using Point2 = K::Point_2;
using Point3 = K::Point_3;
using Mesh = CGAL::Surface_mesh;
namespace PMP = CGAL::Polygon_mesh_processing;

struct Point3VPM {
using key_type = Mesh::Vertex_index;
using value_type = Point3;
using reference = value_type;
using category = boost::readable_property_map_tag;

Point3VPM(const Mesh &m)
: mesh(std::cref(m))
{
}

friend Point3VPM::value_type get(const Point3VPM &map, Mesh::Vertex_index idx)
{
Point2 p = map.mesh.get().point(idx);
return {p[0], p[1], 0};
}

std::reference_wrapper mesh;
};

using AABBTreeTraits = CGAL::AABB_traits;
using AABBTree = CGAL::AABB_tree;
} // namespace

int main()
{
Mesh mesh;

AABBTree tree;
PMP::build_AABB_tree(mesh, tree);
auto location = PMP::locate_with_AABB_tree({0, 0}, tree, mesh);
}


Подробнее здесь: https://stackoverflow.com/questions/663 ... operty-map