/*****************************************************************************
 * MeshLab                                                           o o     *
 * An extendible mesh processor                                    o     o   *
 *                                                                _   O  _   *
 * Copyright(C) 2005-2021                                           \/)\/    *
 * 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 "cleanfilter.h"

#include <QCoreApplication>
#include <vcg/complex/algorithms/clean.h>
#include <vcg/complex/algorithms/create/ball_pivoting.h>
#include <vcg/complex/algorithms/create/platonic.h>
#include <vcg/complex/algorithms/stat.h>
#include <vcg/complex/algorithms/update/texture.h>

using namespace std;
using namespace vcg;

int SnapVertexBorder(CMeshO& m, Scalarm threshold, vcg::CallBackPos* cb);

CleanFilter::CleanFilter()
{
	typeList = {
		FP_BALL_PIVOTING,
		FP_REMOVE_WRT_Q,
		FP_REMOVE_ISOLATED_COMPLEXITY,
		FP_REMOVE_ISOLATED_DIAMETER,
		FP_REMOVE_TVERTEX,
		FP_SNAP_MISMATCHED_BORDER,
		FP_REMOVE_DUPLICATE_FACE,
		FP_REMOVE_FOLD_FACE,
		FP_REPAIR_NON_MANIF_EDGE,
		FP_REMOVE_NON_MANIF_VERT,
		FP_REMOVE_UNREFERENCED_VERTEX,
		FP_REMOVE_DUPLICATED_VERTEX,
		FP_REMOVE_FACE_ZERO_AREA,
		FP_MERGE_CLOSE_VERTEX,
		FP_MERGE_WEDGE_TEX};

	for (ActionIDType tt : types())
		actionList.push_back(new QAction(filterName(tt), this));

	QCoreApplication* app = QCoreApplication::instance();
	if (app != nullptr)
		getFilterAction(FP_SNAP_MISMATCHED_BORDER)->setShortcut(QKeySequence("ALT+`"));
}

CleanFilter::~CleanFilter()
{
}

QString CleanFilter::pluginName() const
{
	return "FilterClean";
}

QString CleanFilter::filterName(ActionIDType filter) const
{
	switch (filter) {
	case FP_BALL_PIVOTING: return QString("Surface Reconstruction: Ball Pivoting");
	case FP_REMOVE_WRT_Q: return QString("Remove Vertices wrt Quality");
	case FP_REMOVE_ISOLATED_DIAMETER: return QString("Remove Isolated pieces (wrt Diameter)");
	case FP_REMOVE_ISOLATED_COMPLEXITY: return QString("Remove Isolated pieces (wrt Face Num.)");
	case FP_REMOVE_TVERTEX: return QString("Remove T-Vertices");
	case FP_SNAP_MISMATCHED_BORDER: return QString("Snap Mismatched Borders");
	case FP_MERGE_CLOSE_VERTEX: return QString("Merge Close Vertices");
	case FP_MERGE_WEDGE_TEX: return QString("Merge Wedge Texture Coord");
	case FP_REMOVE_DUPLICATE_FACE: return QString("Remove Duplicate Faces");
	case FP_REMOVE_FOLD_FACE: return QString("Remove Isolated Folded Faces by Edge Flip");
	case FP_REPAIR_NON_MANIF_EDGE: return QString("Repair non Manifold Edges");
	case FP_REMOVE_NON_MANIF_VERT: return QString("Repair non Manifold Vertices by splitting");
	case FP_REMOVE_UNREFERENCED_VERTEX: return QString("Remove Unreferenced Vertices");
	case FP_REMOVE_DUPLICATED_VERTEX: return QString("Remove Duplicate Vertices");
	case FP_REMOVE_FACE_ZERO_AREA: return QString("Remove Zero Area Faces");
	default: assert(0);
	}
	return QString("error!");
}

QString CleanFilter::filterInfo(ActionIDType filterId) const
{
	switch (filterId) {
	case FP_BALL_PIVOTING:
		return QString(
			"Given a point cloud with normals it reconstructs a surface using the <b>Ball Pivoting "
			"Algorithm</b>."
			"Starting with a seed triangle, the BPA algorithm  pivots a ball of the given radius "
			"around the already formed edges until it touches another point, forming another "
			"triangle. The process continues until all reachable edges have been tried. This "
			"surface reconstruction algorithm uses the existing points without creating new ones. "
			"Works better with uniformly sampled point clouds. If needed first perform a poisson "
			"disk subsampling of the point cloud. <br>"
			"Bernardini F., Mittleman J., Rushmeier H., Silva C., Taubin G.<br>"
			"<b>The ball-pivoting algorithm for surface reconstruction.</b><br>"
			"IEEE TVCG 1999");
	case FP_REMOVE_ISOLATED_COMPLEXITY:
		return QString(
			"Delete isolated connected components composed by a limited number of triangles");
	case FP_REMOVE_ISOLATED_DIAMETER:
		return QString(
			"Delete isolated connected components whose diameter is smaller than the specified "
			"constant");
	case FP_REMOVE_WRT_Q:
		return QString(
			"Delete all the vertices with a quality lower smaller than the specified constant");
	case FP_REMOVE_TVERTEX:
		return QString(
			"Delete t-vertices from the mesh by edge collapse (collapsing the shortest of the "
			"incident edges) or edge flip (flipping the opposite edge on the degenerate face if "
			"the triangulation quality improves).");
	case FP_SNAP_MISMATCHED_BORDER:
		return QString(
			"Try to snap together adjacent borders that are slightly mismatched.<br>"
			"This situation can happen on badly triangulated adjacent patches defined by high "
			"order surfaces.<br>For each border vertex the filter snap it onto the closest "
			"boundary edge only if it is closest of <i>edge_length*threshold</i>. When vertex is "
			"snapped the corresponding face is split and a new vertex is created.");
	case FP_MERGE_CLOSE_VERTEX:
		return QString(
			"Merge together all the vertices that are nearer than the specified threshold. Like a "
			"unify duplicated vertices but with some tolerance.");
	case FP_MERGE_WEDGE_TEX:
		return QString(
			"Merge together per-wedge texture coords that are very close. Used to correct apparent "
			"texture seams that can arise from numerical approximations when saving in ascii "
			"formats.");
	case FP_REMOVE_DUPLICATE_FACE:
		return QString(
			"Delete all the duplicate faces. Two faces are considered equal if they are composed "
			"by the same set of vertices, regardless of the order of the vertices.");
	case FP_REMOVE_FOLD_FACE:
		return QString(
			"Delete all the single folded faces. A face is considered folded if its normal is "
			"opposite to all the adjacent faces. It is removed by flipping it against the face f "
			"adjacent along the edge e such that the vertex opposite to e fall inside f");
	case FP_REPAIR_NON_MANIF_EDGE:
		return QString(
			"Remove non-manifold edges by removing faces (for each non Manifold edge it "
			"iteratively deletes the smallest area face until it becomes 2-Manifold) or by "
			"splitting vertices (each non manifold edges chain will become a border).");
	case FP_REMOVE_NON_MANIF_VERT:
		return QString("Split non Manifold vertices until it becomes 2-Manifold.");
	case FP_REMOVE_UNREFERENCED_VERTEX:
		return QString(
			"Check for every vertex on the mesh: if it is NOT referenced by a face, removes it");
	case FP_REMOVE_DUPLICATED_VERTEX:
		return QString(
			"Check for every vertex on the mesh: if there are two vertices with same coordinates "
			"they are merged into a single one.");
	case FP_REMOVE_FACE_ZERO_AREA:
		return QString("Remove null faces (the one with area equal to zero)");
	default: assert(0);
	}
	return QString("error!");
}

CleanFilter::FilterClass CleanFilter::getClass(const QAction* a) const
{
	switch (ID(a)) {
	case FP_REMOVE_WRT_Q:
	case FP_REMOVE_ISOLATED_DIAMETER:
	case FP_REMOVE_ISOLATED_COMPLEXITY:
	case FP_REMOVE_TVERTEX:
	case FP_REMOVE_FOLD_FACE:
	case FP_MERGE_CLOSE_VERTEX:
	case FP_REMOVE_DUPLICATE_FACE:
	case FP_SNAP_MISMATCHED_BORDER:
	case FP_REPAIR_NON_MANIF_EDGE:
	case FP_REMOVE_NON_MANIF_VERT:
	case FP_REMOVE_FACE_ZERO_AREA:
	case FP_REMOVE_UNREFERENCED_VERTEX:
	case FP_REMOVE_DUPLICATED_VERTEX:
	case FP_BALL_PIVOTING: return FilterPlugin::Remeshing;
	case FP_MERGE_WEDGE_TEX:
		return FilterPlugin::FilterClass(FilterPlugin::Cleaning + FilterPlugin::Texture);
	default: assert(0);
	}
	return FilterPlugin::Generic;
}

int CleanFilter::getRequirements(const QAction* action)
{
	switch (ID(action)) {
	case FP_REMOVE_WRT_Q:
	case FP_BALL_PIVOTING: return MeshModel::MM_VERTMARK;
	case FP_REMOVE_ISOLATED_COMPLEXITY:
	case FP_REMOVE_ISOLATED_DIAMETER: return MeshModel::MM_FACEFACETOPO | MeshModel::MM_FACEMARK;
	case FP_REMOVE_TVERTEX: return MeshModel::MM_FACEFACETOPO | MeshModel::MM_VERTMARK;
	case FP_REPAIR_NON_MANIF_EDGE: return MeshModel::MM_FACEFACETOPO | MeshModel::MM_VERTMARK;
	case FP_REMOVE_NON_MANIF_VERT: return MeshModel::MM_FACEFACETOPO | MeshModel::MM_VERTMARK;
	case FP_SNAP_MISMATCHED_BORDER:
		return MeshModel::MM_FACEFACETOPO | MeshModel::MM_VERTMARK | MeshModel::MM_FACEMARK;
	case FP_REMOVE_FOLD_FACE: return MeshModel::MM_FACEFACETOPO | MeshModel::MM_VERTMARK;
	case FP_MERGE_CLOSE_VERTEX:
	case FP_REMOVE_DUPLICATE_FACE: return MeshModel::MM_NONE;
	case FP_MERGE_WEDGE_TEX: return MeshModel::MM_VERTFACETOPO | MeshModel::MM_WEDGTEXCOORD;
	case FP_REMOVE_UNREFERENCED_VERTEX: return MeshModel::MM_NONE;
	case FP_REMOVE_DUPLICATED_VERTEX: return MeshModel::MM_NONE;
	case FP_REMOVE_FACE_ZERO_AREA: return MeshModel::MM_NONE;
	default: assert(0);
	}
	return 0;
}

int CleanFilter::postCondition(const QAction* action) const
{
	switch (ID(action)) {
	case FP_BALL_PIVOTING:
	case FP_REMOVE_WRT_Q:
	case FP_REMOVE_ISOLATED_DIAMETER:
	case FP_REMOVE_ISOLATED_COMPLEXITY:
	case FP_REMOVE_TVERTEX:
	case FP_SNAP_MISMATCHED_BORDER:
	case FP_MERGE_CLOSE_VERTEX:
	case FP_MERGE_WEDGE_TEX:
	case FP_REMOVE_DUPLICATE_FACE:
	case FP_REMOVE_FOLD_FACE:
	case FP_REPAIR_NON_MANIF_EDGE:
	case FP_REMOVE_NON_MANIF_VERT:
	case FP_REMOVE_UNREFERENCED_VERTEX:
	case FP_REMOVE_DUPLICATED_VERTEX:
	case FP_REMOVE_FACE_ZERO_AREA: return MeshModel::MM_GEOMETRY_AND_TOPOLOGY_CHANGE;
	}
	return MeshModel::MM_ALL;
}

RichParameterList CleanFilter::initParameterList(const QAction* action, const MeshDocument& md)
{
	RichParameterList  parlst;
	pair<float, float> qualityRange;
	switch (ID(action)) {
	case FP_BALL_PIVOTING:
		parlst.addParam(RichAbsPerc(
			"BallRadius",
			0.0f,
			0.0f,
			md.mm()->cm.bbox.Diag(),
			"Pivoting Ball radius (0 autoguess)",
			"The radius of the ball pivoting (rolling) over the set of points. Gaps that are "
			"larger than the ball radius will not be filled; similarly the small pits that are "
			"smaller than the ball radius will be filled."));
		parlst.addParam(RichFloat(
			"Clustering",
			20.0f,
			"Clustering radius (% of ball radius)",
			"To avoid the creation of too small triangles, if a vertex is found too close to a "
			"previous one, it is clustered/merged with it."));
		parlst.addParam(RichFloat(
			"CreaseThr",
			90.0f,
			"Angle Threshold (degrees)",
			"If we encounter a crease angle that is too large we should stop the ball rolling"));
		parlst.addParam(RichBool(
			"DeleteFaces",
			false,
			"Delete initial set of faces",
			"if true all the initial faces of the mesh are deleted and the whole surface is "
			"rebuilt from scratch. Otherwise the current faces are used as a starting point. "
			"Useful if you run the algorithm multiple times with an increasing ball radius."));
		break;
	case FP_REMOVE_ISOLATED_DIAMETER:
		parlst.addParam(RichAbsPerc(
			"MinComponentDiag",
			md.mm()->cm.bbox.Diag() / 10.0f,
			0.0f,
			md.mm()->cm.bbox.Diag(),
			"Enter max diameter of isolated pieces",
			"Delete all the connected components (floating pieces) with a diameter smaller than "
			"the specified one"));
		parlst.addParam(RichBool(
			"removeUnref",
			true,
			"Remove unfreferenced vertices",
			"if true, the unreferenced vertices remaining after the face deletion are removed."));
		break;
	case FP_REMOVE_ISOLATED_COMPLEXITY:
		parlst.addParam(RichInt(
			"MinComponentSize",
			25,
			"Enter minimum conn. comp size:",
			"Delete all the connected components (floating pieces) composed by a number of "
			"triangles smaller than the specified one"));
		parlst.addParam(RichBool(
			"removeUnref",
			true,
			"Remove unfreferenced vertices",
			"if true, the unreferenced vertices remaining after the face deletion are removed."));
		break;
	case FP_REMOVE_WRT_Q:
		qualityRange = tri::Stat<CMeshO>::ComputePerVertexQualityMinMax(md.mm()->cm);
		parlst.addParam(RichAbsPerc(
			"MaxQualityThr",
			(float) 1.0,
			qualityRange.first,
			qualityRange.second,
			"Delete all vertices with quality under:"));
		break;
	case FP_MERGE_CLOSE_VERTEX:
		parlst.addParam(RichAbsPerc(
			"Threshold",
			md.mm()->cm.bbox.Diag() / 10000.0f,
			0.0f,
			md.mm()->cm.bbox.Diag() / 100.0f,
			"Merging distance",
			"All the vertices that closer than this threshold are merged together. Use very small "
			"values, default values is 1/10000 of bounding box diagonal. "));
		break;
	case FP_MERGE_WEDGE_TEX:
		parlst.addParam(RichFloat(
			"MergeThr",
			1.0f / 10000.0f,
			"Merging Threshold",
			"All the per-wedge texture coords that are on the same vertex and are distant less "
			"then the given threshold are merged together. It can be used to remove the fake "
			"texture seams that arise from error. Distance is in texture space (the default, 1e-4, "
			"corresponds to one texel on a 10kx10x texture) "));
		break;
	case FP_SNAP_MISMATCHED_BORDER:
		parlst.addParam(RichFloat(
			"EdgeDistRatio",
			1 / 100.0f,
			"Edge Distance Ratio",
			"Collapse edge when the edge / distance ratio is greater than this value. E.g. for "
			"default value 1000 two straight border edges are collapsed if the central vertex dist "
			"from the straight line composed by the two edges less than a 1/1000 of the sum of the "
			"edges length. Larger values enforce that only vertices very close to the line are "
			"removed."));
		parlst.addParam(RichBool(
			"UnifyVertices",
			true,
			"UnifyVertices",
			"if true the snap vertices are weld together."));
		break;
	case FP_REMOVE_TVERTEX:
		parlst.addParam(RichEnum(
			"method",
			0,
			{"Edge Collapse", "Edge Flip"},
			"Method",
			"Selects wether to remove t-vertices by edge collapse or edge flip."));
		parlst.addParam(RichFloat(
			"Threshold",
			40,
			"Ratio",
			"Detects faces where the base/height ratio is lower than this value"));
		parlst.addParam(RichBool(
			"Repeat",
			true,
			"Iterate until convergence",
			"Iterates the algorithm until it reaches convergence"));
		break;
	case FP_REMOVE_NON_MANIF_VERT:
		parlst.addParam(RichFloat(
			"VertDispRatio",
			0,
			"Vertex Displacement Ratio",
			"This parameter denote the ratio ⍺ of displacement of a vertex. When a vertex <i>v</i> "
			"is split, it is moved towards the barycenter <i>b</i> of the FF connected faces "
			"sharing it of a (<i>v</i>-<i>b</i>)*⍺. When ⍺ is zero vertex is not displaced. When ⍺ "
			"is 0.5 the new vertex is half away toward the barycenter of the face. Reasonable "
			"values are in the [0 .. 0.1] range."));
		break;
	case FP_REPAIR_NON_MANIF_EDGE:
		parlst.addParam(RichEnum(
			"method",
			0,
			{"Remove Faces", "Split Vertices"},
			"Method",
			"Selects wether to remove non manifold edges by removing faces or by splitting "
			"vertices."));
	default: break; // do not add any parameter for the other filters
	}
	return parlst;
}

std::map<std::string, QVariant> CleanFilter::applyFilter(
	const QAction*           filter,
	const RichParameterList& par,
	MeshDocument&            md,
	unsigned int&            postConditionMask,
	vcg::CallBackPos*        cb)
{
	MeshModel& m = *(md.mm());
	switch (ID(filter)) {
	case FP_BALL_PIVOTING: {
		vcg::tri::Allocator<CMeshO>::CompactEveryVector(m.cm);
		Scalarm Radius      = par.getAbsPerc("BallRadius");
		Scalarm Clustering  = par.getFloat("Clustering") / 100.0f;
		Scalarm CreaseThr   = math::ToRad(par.getFloat("CreaseThr"));
		bool    DeleteFaces = par.getBool("DeleteFaces");
		if (DeleteFaces) {
			m.cm.fn = 0;
			m.cm.face.resize(0);
		}
		m.updateDataMask(MeshModel::MM_VERTFACETOPO);
		int                       startingFn = m.cm.fn;
		tri::BallPivoting<CMeshO> pivot(m.cm, Radius, Clustering, CreaseThr);
		// the main processing
		pivot.BuildMesh(cb);
		m.clearDataMask(MeshModel::MM_FACEFACETOPO);
		log("Reconstructed surface. Added %i faces", m.cm.fn - startingFn);
	} break;

	case FP_REMOVE_ISOLATED_DIAMETER: {
		Scalarm             minCC = par.getAbsPerc("MinComponentDiag");
		std::pair<int, int> delInfo =
			tri::Clean<CMeshO>::RemoveSmallConnectedComponentsDiameter(m.cm, minCC);
		log("Removed %i connected components out of %i", delInfo.second, delInfo.first);
		if (par.getBool("removeUnref")) {
			int delvert = tri::Clean<CMeshO>::RemoveUnreferencedVertex(m.cm);
			log("Removed %d unreferenced vertices", delvert);
		}
		m.updateBoxAndNormals();
	} break;
	case FP_REMOVE_ISOLATED_COMPLEXITY: {
		int                 minCC = par.getInt("MinComponentSize");
		std::pair<int, int> delInfo =
			tri::Clean<CMeshO>::RemoveSmallConnectedComponentsSize(m.cm, minCC);
		log("Removed %i connected components out of %i", delInfo.second, delInfo.first);
		if (par.getBool("removeUnref")) {
			int delvert = tri::Clean<CMeshO>::RemoveUnreferencedVertex(m.cm);
			log("Removed %d unreferenced vertices", delvert);
		}
		m.updateBoxAndNormals();
	} break;

	case FP_REMOVE_WRT_Q: {
		int                    deletedFN = 0;
		int                    deletedVN = 0;
		Scalarm                val       = par.getAbsPerc("MaxQualityThr");
		CMeshO::VertexIterator vi;
		for (vi = m.cm.vert.begin(); vi != m.cm.vert.end(); ++vi)
			if (!(*vi).IsD() && (*vi).Q() < val) {
				tri::Allocator<CMeshO>::DeleteVertex(m.cm, *vi);
				deletedVN++;
			}

		CMeshO::FaceIterator fi;
		for (fi = m.cm.face.begin(); fi != m.cm.face.end(); ++fi)
			if (!(*fi).IsD())
				if ((*fi).V(0)->IsD() || (*fi).V(1)->IsD() || (*fi).V(2)->IsD()) {
					tri::Allocator<CMeshO>::DeleteFace(m.cm, *fi);
					deletedFN++;
				}

		m.clearDataMask(MeshModel::MM_FACEFACETOPO);
		log("Deleted %i vertices and %i faces with a quality lower than %f",
			deletedVN,
			deletedFN,
			val);
		m.updateBoxAndNormals();
	} break;

	case FP_REMOVE_TVERTEX: {
		int     total     = 0;
		Scalarm threshold = par.getFloat("Threshold");
		bool    repeat    = par.getBool("Repeat");
		if (par.getEnum("method") == 0) { // edge collapse
			bool hadFF = false;
			bool hadVF = false;
			if (m.hasDataMask(MeshModel::MM_FACEFACETOPO)) {
				m.clearDataMask(MeshModel::MM_FACEFACETOPO);
				hadFF = true;
			}
			if (m.hasDataMask(MeshModel::MM_VERTFACETOPO)) {
				m.clearDataMask(MeshModel::MM_VERTFACETOPO);
				hadVF = true;
			}
			total = tri::Clean<CMeshO>::RemoveTVertexByCollapse(m.cm, threshold, repeat);

			if (hadFF)
				m.updateDataMask(MeshModel::MM_FACEFACETOPO);
			if (hadVF)
				m.updateDataMask(MeshModel::MM_VERTFACETOPO);
		}
		else { // edge flip
			if (vcg::tri::Clean<CMeshO>::CountNonManifoldEdgeFF(m.cm) > 0 ||
				vcg::tri::Clean<CMeshO>::CountNonManifoldVertexFF(m.cm) > 0) {
				throw MLException("Non manifold mesh. Please clean the mesh first.");
			}
			total = tri::Clean<CMeshO>::RemoveTVertexByFlip(m.cm, threshold, repeat);
		}
		log("Successfully removed %d t-vertices", total);
	} break;

	case FP_MERGE_WEDGE_TEX: {
		Scalarm threshold = par.getFloat("MergeThr");
		tri::UpdateTopology<CMeshO>::VertexFace(m.cm);
		int total = tri::UpdateTexture<CMeshO>::WedgeTexMergeClose(m.cm, threshold);
		log("Successfully merged %d wedge tex coord distant less than %f", total, threshold);
	} break;

	case FP_MERGE_CLOSE_VERTEX: {
		Scalarm threshold = par.getAbsPerc("Threshold");
		int     total     = tri::Clean<CMeshO>::MergeCloseVertex(m.cm, threshold);
		log("Successfully merged %d vertices", total);
	} break;

	case FP_REMOVE_DUPLICATE_FACE: {
		int total = tri::Clean<CMeshO>::RemoveDuplicateFace(m.cm);
		log("Successfully deleted %d duplicated faces", total);
	} break;

	case FP_REMOVE_FOLD_FACE: {
		m.updateDataMask(MeshModel::MM_FACECOLOR);
		int total = tri::Clean<CMeshO>::RemoveFaceFoldByFlip(m.cm);
		log("Successfully flipped %d folded faces", total);
		m.updateBoxAndNormals();
	} break;

	case FP_REPAIR_NON_MANIF_EDGE: {
		int total = 0;
		if (par.getEnum("method") == 0) { // removing faces
			total = tri::Clean<CMeshO>::RemoveNonManifoldFace(m.cm);
			log("Successfully removed %d non-manifold faces", total);
			m.updateBoxAndNormals();
		}
		else { // splitting vertices
			total = tri::Clean<CMeshO>::SplitManifoldComponents(m.cm);
			log("Successfully split the mesh into %d edge manifold components", total);
			m.updateBoxAndNormals();
			if (m.hasDataMask(MeshModel::MM_WEDGTEXCOORD))
				postConditionMask =
					MeshModel::MM_GEOMETRY_AND_TOPOLOGY_CHANGE | MeshModel::MM_WEDGTEXCOORD;
		}
	} break;

	case FP_REMOVE_NON_MANIF_VERT: {
		Scalarm threshold = par.getFloat("VertDispRatio");
		int     total     = tri::Clean<CMeshO>::SplitNonManifoldVertex(m.cm, threshold);
		log("Successfully split %d non manifold vertices faces", total);
		m.updateBoxAndNormals();
	} break;

	case FP_REMOVE_FACE_ZERO_AREA: {
		int nullFaces = tri::Clean<CMeshO>::RemoveFaceOutOfRangeArea(m.cm, 0);
		log("Removed %d null faces", nullFaces);
		m.clearDataMask(MeshModel::MM_FACEFACETOPO);
	} break;

	case FP_REMOVE_UNREFERENCED_VERTEX: {
		int delvert = tri::Clean<CMeshO>::RemoveUnreferencedVertex(m.cm);
		log("Removed %d unreferenced vertices", delvert);
		if (delvert != 0)
			m.updateBoxAndNormals();
	} break;

	case FP_REMOVE_DUPLICATED_VERTEX: {
		int delvert = tri::Clean<CMeshO>::RemoveDuplicateVertex(m.cm);
		log("Removed %d duplicated vertices", delvert);
		if (delvert != 0)
			m.updateBoxAndNormals();
		m.clearDataMask(MeshModel::MM_FACEFACETOPO);
		m.clearDataMask(MeshModel::MM_VERTFACETOPO);
	} break;

	case FP_SNAP_MISMATCHED_BORDER: {
		Scalarm threshold = par.getFloat("EdgeDistRatio");
		int     total     = SnapVertexBorder(m.cm, threshold, cb);
		log("Successfully Split %d faces to snap", total);
		m.clearDataMask(MeshModel::MM_FACEFACETOPO);
		m.clearDataMask(MeshModel::MM_VERTFACETOPO);
	} break;

	default: wrongActionCalled(filter); // unknown filter;
	}
	return std::map<std::string, QVariant>();
}

int SnapVertexBorder(CMeshO& m, Scalarm threshold, vcg::CallBackPos* cb)
{
	tri::Allocator<CMeshO>::CompactEveryVector(m);

	tri::UpdateTopology<CMeshO>::FaceFace(m);
	tri::UpdateFlags<CMeshO>::FaceBorderFromFF(m);
	tri::UpdateFlags<CMeshO>::VertexBorderFromFaceBorder(m);
	tri::UpdateNormal<CMeshO>::PerVertexNormalizedPerFaceNormalized(m);
	typedef GridStaticPtr<CMeshO::FaceType, CMeshO::ScalarType> MetroMeshFaceGrid;
	MetroMeshFaceGrid                                           unifGridFace;
	typedef tri::FaceTmark<CMeshO>                              MarkerFace;
	MarkerFace                                                  markerFunctor(&m);
	vcg::face::PointDistanceBaseFunctor<CMeshO::ScalarType>     PDistFunct;
	tri::UpdateFlags<CMeshO>::FaceClearV(m);
	unifGridFace.Set(m.face.begin(), m.face.end());

	int                         faceFound;
	int                         K = 20;
	Point3m                     startPt;
	float                       maxDist = m.bbox.Diag() / 20;
	vector<Point3m>             splitVertVec;
	vector<CMeshO::FacePointer> splitFaceVec;
	vector<int>                 splitEdgeVec;
	for (CMeshO::VertexIterator vi = m.vert.begin(); vi != m.vert.end(); ++vi)
		if ((*vi).IsB()) {
			cb((int(tri::Index(m, *vi)) * 100) / m.vn, "Snapping vertices");
			vector<CMeshO::FacePointer> faceVec;
			vector<float>               distVec;
			vector<Point3m>             pointVec;
			Point3m                     u;
			startPt   = (*vi).P();
			faceFound = unifGridFace.GetKClosest(
				PDistFunct, markerFunctor, K, startPt, maxDist, faceVec, distVec, pointVec);

			CMeshO::FacePointer bestFace = 0;
			float               localThr, bestDist = std::numeric_limits<float>::max();
			Point3m             bestPoint;
			int                 bestEdge;
			//        qDebug("Found %i face for vertex %i",faceFound,vi-m.vert.begin());
			for (int i = 0; i < faceFound; ++i) {
				const float         epsilonSmall = float(1e-5);
				const float         epsilonBig   = float(1e-2);
				CMeshO::FacePointer fp           = faceVec[i];
				InterpolationParameters(*fp, fp->cN(), pointVec[i], u);
				//          qDebug(" face %i face for vertex %5.3f %5.3f %5.3f  dist %5.3f  (%c %c
				//          %c)",fp-&*m.face.begin(),u[0],u[1],u[2],distVec[i],IsBorder(*fp,0)?'b':'
				//          ',IsBorder(*fp,1)?'b':' ',IsBorder(*fp,2)?'b':' ');
				for (int j = 0; j < 3; ++j) {
					if (IsBorder(*fp, j) && !fp->IsV()) {
						if (u[(j + 0) % 3] > epsilonBig && u[(j + 1) % 3] > epsilonBig &&
							u[(j + 2) % 3] < epsilonSmall) {
							if (distVec[i] < bestDist) {
								bestDist = distVec[i];
								// bestPoint=pointVec[i];
								bestPoint = (*vi).cP();
								bestFace  = fp;
								bestEdge  = j;
							}
						}
					}
				}
			} // end for each faceFound

			if (bestFace) {
				localThr = threshold * Distance(bestFace->P0(bestEdge), bestFace->P1(bestEdge));
				if (bestDist < localThr && !bestFace->IsV()) {
					bestFace->SetV();
					(*vi).C() = Color4b::Blue;
					// bestFace->C()=Color4b::LightBlue;
					(*vi).SetS();
					splitVertVec.push_back(bestPoint);
					splitEdgeVec.push_back(bestEdge);
					splitFaceVec.push_back(bestFace);
				}
			}
		} // end for all border vertices
	tri::Allocator<CMeshO>::PointerUpdater<CMeshO::FacePointer> pu;
	CMeshO::VertexIterator firstVert = tri::Allocator<CMeshO>::AddVertices(m, splitVertVec.size());
	CMeshO::FaceIterator   firstface = tri::Allocator<CMeshO>::AddFaces(m, splitVertVec.size(), pu);
	//
	//             ^                           ^
	//           /   \                       / | \          .
	//         /       \                   /   |   \        .
	//       /           \               /     |     \      .
	//     /       fp      \           /       |       \    .
	//   /                   \       /    fp   |   ff    \  .
	//  V0 ------------------V2     V0 -------fv---------V2
	//       i

	for (size_t i = 0; i < splitVertVec.size(); ++i) {
		firstVert->P()           = splitVertVec[i];
		int                 eInd = splitEdgeVec[i];
		CMeshO::FacePointer fp   = splitFaceVec[i];
		pu.Update(fp);
		firstface->V(0) = &*firstVert;
		firstface->V(1) = fp->V2(eInd);
		firstface->V(2) = fp->V0(eInd);
		//      firstface->C()=Color4b::LightBlue;

		fp->V0(eInd) = &*firstVert;

		++firstface;
		++firstVert;
	}
	tri::UpdateNormal<CMeshO>::PerVertexNormalizedPerFaceNormalized(m);
	return int(splitVertVec.size());
}

//
//             ^-----------------
//           / | \             /
//         /   |   \     fj   /
//       /  i+1|     \ij     /
//     /       |       \    /
//   /    fi   |   fadj  \ /
//  V0 --------V1--------V2
//       i
//

int DeleteCollinearBorder(CMeshO& m, float threshold)
{
	int                  total = 0;
	CMeshO::FaceIterator fi;
	for (fi = m.face.begin(); fi != m.face.end(); ++fi) {
		if (!(*fi).IsD()) {
			for (int i = 0; i < 3; ++i) {
				if (face::IsBorder(*fi, i) && !face::IsBorder(*fi, (i + 1) % 3)) {
					CMeshO::VertexPointer V0         = (*fi).V0(i);
					CMeshO::VertexPointer V1         = (*fi).V1(i);
					CMeshO::VertexPointer V2         = 0;
					CMeshO::FacePointer   fadj       = (*fi).FFp((i + 1) % 3);
					int                   adjBordInd = (*fi).FFi((i + 1) % 3);
					if (fadj->V1(adjBordInd) == V1)
						V2 = fadj->V2(adjBordInd);
					else
						continue; // non coerent face ordering.
					if (face::IsBorder(*fadj, (adjBordInd + 1) % 3)) {
						// the colinearity test;
						Point3m            pp;
						CMeshO::ScalarType dist;
						SegmentPointDistance(Segment3m(V0->cP(), V2->cP()), V1->cP(), pp, dist);
						if (dist * threshold < Distance(V0->cP(), V2->cP())) {
							(*fi).V1(i) = V2;
							if (face::IsBorder(*fadj, (adjBordInd + 2) % 3)) {
								(*fi).FFp((i + 1) % 3) = &*fi;
								(*fi).FFi((i + 1) % 3) = (i + 1) % 3;
							}
							else {
								CMeshO::FacePointer fj = fadj->FFp((adjBordInd + 2) % 3);
								int                 ij = fadj->FFi((adjBordInd + 2) % 3);
								(*fi).FFp((i + 1) % 3) = fj;
								(*fi).FFi((i + 1) % 3) = ij;
								fj->FFp(ij)            = &*fi;
								fj->FFi(ij)            = (i + 1) % 3;
							}
							tri::Allocator<CMeshO>::DeleteFace(m, *fadj);
							total++;
						}
					}
				}
			}
		}
	}
	return total;
}

MESHLAB_PLUGIN_NAME_EXPORTER(CleanFilter)
