#define NOMINMAX #include "OpenVDBUtils.hpp" #ifdef _MSC_VER // Suppress warning C4146 in OpenVDB: unary minus operator applied to unsigned type, result still unsigned #pragma warning(push) #pragma warning(disable : 4146) #endif // _MSC_VER #include #ifdef _MSC_VER #pragma warning(pop) #endif // _MSC_VER #include #include #include #include //#include "MTUtils.hpp" namespace Slic3r { class TriangleMeshDataAdapter { public: const indexed_triangle_set &its; float voxel_scale; size_t polygonCount() const { return its.indices.size(); } size_t pointCount() const { return its.vertices.size(); } size_t vertexCount(size_t) const { return 3; } // Return position pos in local grid index space for polygon n and vertex v // The actual mesh will appear to openvdb as scaled uniformly by voxel_size // And the voxel count per unit volume can be affected this way. void getIndexSpacePoint(size_t n, size_t v, openvdb::Vec3d& pos) const { auto vidx = size_t(its.indices[n](Eigen::Index(v))); Slic3r::Vec3d p = its.vertices[vidx].cast() * voxel_scale; pos = {p.x(), p.y(), p.z()}; } TriangleMeshDataAdapter(const indexed_triangle_set &m, float voxel_sc = 1.f) : its{m}, voxel_scale{voxel_sc} {}; }; openvdb::FloatGrid::Ptr mesh_to_grid(const indexed_triangle_set & mesh, const openvdb::math::Transform &tr, float voxel_scale, float exteriorBandWidth, float interiorBandWidth) { // Might not be needed but this is now proven to be working openvdb::initialize(); std::vector meshparts = its_split(mesh); auto it = std::remove_if(meshparts.begin(), meshparts.end(), [](auto &m) { return its_volume(m) < EPSILON; }); meshparts.erase(it, meshparts.end()); openvdb::FloatGrid::Ptr grid; for (auto &m : meshparts) { auto subgrid = openvdb::tools::meshToVolume( TriangleMeshDataAdapter{m, voxel_scale}, tr, 1.f, 1.f); if (grid && subgrid) openvdb::tools::csgUnion(*grid, *subgrid); else if (subgrid) grid = std::move(subgrid); } if(meshparts.empty()) { // Splitting failed, fall back to hollow the original mesh grid = openvdb::tools::meshToVolume( TriangleMeshDataAdapter{mesh}, tr, 1.f, 1.f); } constexpr int DilateIterations = 1; grid = openvdb::tools::dilateSdf( *grid, interiorBandWidth, openvdb::tools::NN_FACE_EDGE, DilateIterations, openvdb::tools::FastSweepingDomain::SWEEP_LESS_THAN_ISOVALUE); grid = openvdb::tools::dilateSdf( *grid, exteriorBandWidth, openvdb::tools::NN_FACE_EDGE, DilateIterations, openvdb::tools::FastSweepingDomain::SWEEP_GREATER_THAN_ISOVALUE); grid->insertMeta("voxel_scale", openvdb::FloatMetadata(voxel_scale)); return grid; } indexed_triangle_set grid_to_mesh(const openvdb::FloatGrid &grid, double isovalue, double adaptivity, bool relaxDisorientedTriangles) { openvdb::initialize(); std::vector points; std::vector triangles; std::vector quads; openvdb::tools::volumeToMesh(grid, points, triangles, quads, isovalue, adaptivity, relaxDisorientedTriangles); float scale = 1.; try { scale = grid.template metaValue("voxel_scale"); } catch (...) { } indexed_triangle_set ret; ret.vertices.reserve(points.size()); ret.indices.reserve(triangles.size() + quads.size() * 2); for (auto &v : points) ret.vertices.emplace_back(to_vec3f(v) / scale); for (auto &v : triangles) ret.indices.emplace_back(to_vec3i(v)); for (auto &quad : quads) { ret.indices.emplace_back(quad(0), quad(1), quad(2)); ret.indices.emplace_back(quad(2), quad(3), quad(0)); } return ret; } openvdb::FloatGrid::Ptr redistance_grid(const openvdb::FloatGrid &grid, double iso, double er, double ir) { auto new_grid = openvdb::tools::levelSetRebuild(grid, float(iso), float(er), float(ir)); // Copies voxel_scale metadata, if it exists. new_grid->insertMeta(*grid.deepCopyMeta()); return new_grid; } openvdb::FloatGrid::Ptr redistance_grid(const openvdb::FloatGrid &grid, double iso) { auto new_grid = openvdb::tools::levelSetRebuild(grid, float(iso)); // Copies voxel_scale metadata, if it exists. new_grid->insertMeta(*grid.deepCopyMeta()); return new_grid; } } // namespace Slic3r