1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
#include <vector>
using namespace std;
// VCG headers for triangular mesh processing
#include<vcg/simplex/edge/base.h>
#include<vcg/simplex/edge/component.h>
#include<vcg/simplex/vertex/base.h>
#include<vcg/simplex/face/base.h>
#include <vcg/complex/complex.h>
#include <vcg/complex/algorithms/update/topology.h>
#include <vcg/complex/algorithms/update/edges.h>
#include <vcg/complex/algorithms/update/bounding.h>
#include <vcg/complex/algorithms/update/flag.h>
#include <vcg/complex/algorithms/clean.h>
#include <vcg/complex/algorithms/intersection.h>
#include <vcg/space/index/grid_static_ptr.h>
//#include <vcg/simplex/edge/with/ae.h>
//#include <vcg/complex/edgemesh/base.h>
//#include <vcg/complex/edgemesh/allocate.h>
//#include <vcg/complex/edgemesh/update/bounding.h>
// VCG File Format Importer/Exporter
#include <wrap/io_trimesh/import.h>
#include <wrap/io_edgemesh/export_svg.h>
#include <wrap/io_edgemesh/export_dxf.h>
// VCG Vertex
// VCG Faces
using namespace vcg;
class MyFace;
class MyEdge;
class MyVertex;
struct MyUsedTypes : public UsedTypes< Use<MyVertex> ::AsVertexType,
Use<MyEdge> ::AsEdgeType,
Use<MyFace> ::AsFaceType>{};
class MyVertex : public Vertex< MyUsedTypes, vertex::Coord3f, vertex::BitFlags, vertex::Normal3f, vertex::Mark>{};
class MyEdge : public Edge< MyUsedTypes, edge::VertexRef, edge::EVAdj> {};
class MyFace : public Face <MyUsedTypes, face::VertexRef,face::FFAdj, face::BitFlags, face::Normal3f> {};
class MyEdgeMesh: public tri::TriMesh< vector<MyVertex>, vector<MyEdge> > {};
class MyMesh : public tri::TriMesh< vector<MyVertex>, vector<MyFace > >{};
typedef vcg::GridStaticPtr<MyMesh::FaceType, MyMesh::ScalarType> TriMeshGrid;
int main(int argc,char ** argv)
{
if (argc<6)
{
printf("\n");
printf(" Usage: trimesh_intersection <filename> <a> <b> <c> <d>\n\n");
printf(" <filename> Mesh model to intersect (PLY format).\n");
printf(" <a> <b> <c> <d> The coefficients that specifying a plane in the form:\n\n");
printf(" a*x + b*y + c*z + d = 0\n\n\n");
printf(" Example: trimesh_intersection bunny.ply 0.0 1.0 0.0 0.0\n\n");
return 0;
}
MyMesh m;
// open a mesh
int err = tri::io::Importer<MyMesh>::Open(m,argv[1]);
if(err)
{
printf("Error in reading %s: '%s'\n",argv[1],tri::io::Importer<MyMesh>::ErrorMsg(err));
exit(-1);
}
// some cleaning to get rid of bad file formats like stl that duplicate vertexes..
int dup = tri::Clean<MyMesh>::RemoveDuplicateVertex(m);
int unref = tri::Clean<MyMesh>::RemoveUnreferencedVertex(m);
if (dup > 0 || unref > 0)
printf("Removed %i duplicate and %i unreferenced vertices from mesh %s\n",dup,unref,argv[1]);
// Compute cross-intersection with the given plane
/////////////////////////////////////////////////////////
MyMesh::ScalarType a = static_cast<MyMesh::ScalarType>(atof(argv[2]));
MyMesh::ScalarType b = static_cast<MyMesh::ScalarType>(atof(argv[3]));
MyMesh::ScalarType c = static_cast<MyMesh::ScalarType>(atof(argv[4]));
MyMesh::ScalarType d = static_cast<MyMesh::ScalarType>(atof(argv[5]));
vcg::Point3<MyMesh::ScalarType> direction(a, b, c);
MyMesh::ScalarType distance = -d / direction.Norm();
direction.Normalize();
vcg::Plane3<MyMesh::ScalarType> plane(distance, direction);
MyEdgeMesh edge_mesh; // returned EdgeMesh (i.e. the cross-section)
vcg::IntersectionPlaneMesh<MyMesh, MyEdgeMesh, MyMesh::ScalarType>(m, plane, edge_mesh);
// Compute bounding box
vcg::tri::UpdateBounding<MyEdgeMesh>::Box(edge_mesh);
// export the cross-section
tri::io::SVGProperties pro;
if (tri::io::ExporterSVG<MyEdgeMesh>::Save(edge_mesh, "out.svg",pro))
printf(" The cross-intersection has been successfully saved (OUT.SVG).\n");
else
printf(" The cross-intersection cannot be saved.\n");
return 0;
}
|
