Question

您是否知道如何在Matlab中获取由面和顶点集组成的两个多边形网格的交点（对应）体积/网格？

Answer 1

我不知道内置的MATLAB方法来做到这一点。但是通过mex函数，您可以使用ready-made C++ libraries。

以下示例使用 libigl's boolean mesh operation 。

它有一些依赖关系，我使用以下命令在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');

如何获得两个网格（面和顶点）的交点量？

1 个答案: