trimesh_intersection.cpp

Go to the documentation of this file.

/****************************************************************************
* VCGLib                                                            o o     *
* Visual and Computer Graphics Library                            o     o   *
*                                                                _   O  _   *
* Copyright(C) 2004-2012                                           \/)\/    *
* Visual Computing Lab                                            /\/|      *
* ISTI - Italian National Research Council                           |      *
*                                                                    \      *
* All rights reserved.                                                      *
*                                                                           *
* This program is free software; you can redistribute it and/or modify      *
* it under the terms of the GNU General Public License as published by      *
* the Free Software Foundation; either version 2 of the License, or         *
* (at your option) any later version.                                       *
*                                                                           *
* This program is distributed in the hope that it will be useful,           *
* but WITHOUT ANY WARRANTY; without even the implied warranty of            *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the             *
* GNU General Public License (http://www.gnu.org/licenses/gpl.txt)          *
* for more details.                                                         *
*                                                                           *
****************************************************************************/

#include <vcg/complex/complex.h>
#include <vcg/complex/algorithms/update/topology.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>

// 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>

using namespace std;
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;
}