如何获得两个网格(面&顶点)

Question

你知道如何在Matlab中获得由面和顶点集合组成的两个多边形网格的相交(对应的)体积/网格吗?

Answer 1

我不知道一个内置的MATLAB方法来做到这一点。但是通过mex函数，你可以使用现成的c++库。

下面的例子使用了libigl的布尔网格操作。

它有一些依赖，我安装在Ubuntu Linux上使用:

$ sudo apt-get install libcgal-dev git mercurial
$ cd /opt/
$ hg clone https://bitbucket.org/eigen/eigen/
$ git clone https://github.com/libigl/libigl.git

---

get_intersection_mesh.cpp

#include "mex.h"
#include <Eigen/Core>
#include <igl/matlab/prepare_lhs.h>
#include <igl/matlab/parse_rhs.h>
#include <igl/boolean/mesh_boolean.h>
void mexFunction(
     int          nlhs,
     mxArray      *plhs[],
     int          nrhs,
     const mxArray *prhs[]
     )
{
  using namespace Eigen;
  if (nrhs != 4) 
  {
    mexErrMsgIdAndTxt("MATLAB:mexcpp:nargin",
        "get_intersection_mesh requires 4 input parameters: F1, V1, F2, V2");
  }
  if (nlhs != 2) 
  {
    mexErrMsgIdAndTxt("MATLAB:mexcpp:nargout",
        "get_intersection_mesh requires 2 output parameters: F, V");
  }
  MatrixXd V1, V2, V_i;
  MatrixXi F1, F2, F_i;
  // read input meshes
  igl::matlab::parse_rhs_index(prhs+0, F1);
  igl::matlab::parse_rhs_double(prhs+1,V1);
  igl::matlab::parse_rhs_index(prhs+2, F2);
  igl::matlab::parse_rhs_double(prhs+3,V2);
  // calculate intersection
  igl::boolean::mesh_boolean(V1,F1,V2,F2,igl::boolean::MESH_BOOLEAN_TYPE_INTERSECT,V_i,F_i);
  // write output
  igl::matlab::prepare_lhs_index(F_i,plhs+1);
  igl::matlab::prepare_lhs_double(V_i,plhs);
}

---

在MATLAB中使用:

编译它

mex get_intersection_mesh.cpp -I/opt/libigl/include -I/opt/eigen -lboost_thread -lCGAL -lmpfr -lgmp -lboost_system

---

你可以这样使用它:

[V_i, F_i] = get_intersection_mesh(fv1.faces, fv1.vertices, fv2.faces, fv2.vertices);
% plot intersection mesh
patch('Faces',F_i,'Vertices',V_i, 'FaceColor','w');