From a53ac8cd4933f2fbc2c03e1937d414d0b0f7af47 Mon Sep 17 00:00:00 2001
From: alemuntoni <muntoni.alessandro@gmail.com>
Date: Mon, 13 Dec 2021 13:07:32 +0100
Subject: [PATCH] lscm parametrization filter

---
 src/common/utilities/eigen_mesh_conversions.h |  9 ++--
 .../filter_parametrization.cpp                | 52 +++++++++++++++++--
 .../filter_parametrization.h                  |  4 +-
 3 files changed, 57 insertions(+), 8 deletions(-)

diff --git a/src/common/utilities/eigen_mesh_conversions.h b/src/common/utilities/eigen_mesh_conversions.h
index 7bf55913e..4ad5a0e06 100644
--- a/src/common/utilities/eigen_mesh_conversions.h
+++ b/src/common/utilities/eigen_mesh_conversions.h
@@ -31,9 +31,12 @@
 typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 1>      EigenVectorXm;
 typedef Eigen::Matrix<bool, Eigen::Dynamic, 1>         EigenVectorXb;
 typedef Eigen::Matrix<unsigned int, Eigen::Dynamic, 1> EigenVectorXui;
-typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 2>      EigenMatrixX2m;
-typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 3>      EigenMatrixX3m;
-typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 4>      EigenMatrixX4m;
+
+typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 2> EigenMatrixX2m;
+typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 3> EigenMatrixX3m;
+typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, 4> EigenMatrixX4m;
+
+typedef Eigen::Matrix<Scalarm, Eigen::Dynamic, Eigen::Dynamic> EigenMatrixXm;
 
 namespace meshlab {
 
diff --git a/src/meshlabplugins/filter_parametrization/filter_parametrization.cpp b/src/meshlabplugins/filter_parametrization/filter_parametrization.cpp
index e45af7f56..e70540ed2 100644
--- a/src/meshlabplugins/filter_parametrization/filter_parametrization.cpp
+++ b/src/meshlabplugins/filter_parametrization/filter_parametrization.cpp
@@ -27,11 +27,12 @@
 
 #include <igl/boundary_loop.h>
 #include <igl/harmonic.h>
+#include <igl/lscm.h>
 #include <igl/map_vertices_to_circle.h>
 
 FilterParametrizationPlugin::FilterParametrizationPlugin()
-{ 
-	typeList = { FP_HARMONIC_PARAM};
+{
+	typeList = { FP_HARMONIC_PARAM, FP_LEAST_SQUARES_PARAM};
 
 	for(const ActionIDType& tt : typeList)
 		actionList.push_back(new QAction(filterName(tt), this));
@@ -52,6 +53,8 @@ QString FilterParametrizationPlugin::filterName(ActionIDType filterId) const
 	switch(filterId) {
 	case FP_HARMONIC_PARAM :
 		return "Harmonic Parametrization";
+	case FP_LEAST_SQUARES_PARAM:
+		return "Least Squares Conformal Maps Parametrization";
 	default :
 		assert(0);
 		return "";
@@ -63,6 +66,8 @@ QString FilterParametrizationPlugin::filterInfo(ActionIDType filterId) const
 	switch(filterId) {
 	case FP_HARMONIC_PARAM :
 		return "";
+	case FP_LEAST_SQUARES_PARAM:
+		return "";
 	default :
 		assert(0);
 		return "Unknown Filter";
@@ -73,6 +78,7 @@ FilterParametrizationPlugin::FilterClass FilterParametrizationPlugin::getClass(c
 {
 	switch(ID(a)) {
 	case FP_HARMONIC_PARAM :
+	case FP_LEAST_SQUARES_PARAM:
 		return FilterPlugin::Texture;
 	default :
 		assert(0);
@@ -100,12 +106,15 @@ int FilterParametrizationPlugin::postCondition(const QAction*) const
 	return MeshModel::MM_VERTTEXCOORD;
 }
 
-RichParameterList FilterParametrizationPlugin::initParameterList(const QAction *action, const MeshModel &m)
+RichParameterList FilterParametrizationPlugin::initParameterList(const QAction *action, const MeshModel &)
 {
 	RichParameterList parlst;
 	switch(ID(action)) {
 	case FP_HARMONIC_PARAM :
 		parlst.addParam(RichInt("harm_function", 1,"N-Harmonic Function", "1 denotes harmonic function, 2 biharmonic, 3 triharmonic, etc."));
+		break;
+	case FP_LEAST_SQUARES_PARAM:
+
 		break;
 	default :
 		assert(0);
@@ -134,7 +143,10 @@ std::map<std::string, QVariant> FilterParametrizationPlugin::applyFilter(
 		Eigen::VectorXi bnd;
 
 		igl::boundary_loop(faces,bnd);
-		igl::map_vertices_to_circle(verts,bnd,bnd_uv);
+		if (bnd.size() == 0)
+			throw MLException(
+				"Harmonic Parametrization can be applied only on meshes that have a boundary.");
+		igl::map_vertices_to_circle(verts, bnd, bnd_uv);
 		igl::harmonic(verts,faces,bnd,bnd_uv,1,V_uv);
 
 		unsigned int i = 0;
@@ -146,6 +158,38 @@ std::map<std::string, QVariant> FilterParametrizationPlugin::applyFilter(
 
 		break;
 	}
+	case FP_LEAST_SQUARES_PARAM: {
+		EigenMatrixX3m v = meshlab::vertexMatrix(md.mm()->cm);
+		Eigen::MatrixX3d verts = v.cast<double>();
+		Eigen::MatrixX3i faces = meshlab::faceMatrix(md.mm()->cm);
+		Eigen::VectorXi bnd, boundaryPoints(2, 1);
+
+		Eigen::MatrixXd V_uv;
+
+		igl::boundary_loop(faces,bnd);
+		if (bnd.size() == 0)
+			throw MLException(
+				"Least Squares Conformal Maps Parametrization can be applied only on meshes that "
+				"have a boundary.");
+
+		boundaryPoints(0) = bnd(0);
+		boundaryPoints(1) = bnd(bnd.size()/2);
+
+		Eigen::MatrixXd bc(2,2);
+		bc<<0,0,1,0;
+
+		// LSCM parametrization
+		igl::lscm(verts,faces,boundaryPoints,bc,V_uv);
+
+		unsigned int i = 0;
+		for (auto& v : md.mm()->cm.vert){
+			v.T().u() =V_uv(i, 0);
+			v.T().v() =V_uv(i, 1);
+			i++;
+		}
+
+		break;
+	}
 	default :
 		wrongActionCalled(action);
 	}
diff --git a/src/meshlabplugins/filter_parametrization/filter_parametrization.h b/src/meshlabplugins/filter_parametrization/filter_parametrization.h
index b9fc1476c..753d77403 100644
--- a/src/meshlabplugins/filter_parametrization/filter_parametrization.h
+++ b/src/meshlabplugins/filter_parametrization/filter_parametrization.h
@@ -35,7 +35,9 @@ class FilterParametrizationPlugin : public QObject, public FilterPlugin
 
 public:
 	//enum used to give an ID to every filter implemented in the plugin
-	enum FileterIds {FP_HARMONIC_PARAM};
+	enum FileterIds {
+		FP_HARMONIC_PARAM,
+		FP_LEAST_SQUARES_PARAM};
 
 	FilterParametrizationPlugin();