An Application of Manifold Learning in Global Shape Descriptors

With the rapid expansion of applied 3D computational vision, shape descriptors have become increasingly important for a wide variety of applications and objects from molecules to planets. Appropriate shape descriptors are critical for accurate (and efficient) shape retrieval and 3D model classification. Several spectral-based shape descriptors have been introduced by solving various physical equations over a 3D surface model. In this paper, for the first time, we incorporate a specific group of techniques in statistics and machine learning, known as manifold learning, to develop a global shape descriptor in the computer graphics domain. The proposed descriptor utilizes the Laplacian Eigenmap technique in which the Laplacian eigenvalue problem is discretized using an exponential weighting scheme. As a result, our descriptor eliminates the limitations tied to the existing spectral descriptors, namely dependency on triangular mesh representation and high intra-class quality of 3D models. We also present a straightforward normalization method to obtain a scale-invariant descriptor. The extensive experiments performed in this study show that the present contribution provides a highly discriminative and robust shape descriptor under the presence of a high level of noise, random scale variations, and low sampling rate, in addition to the known isometric-invariance property of the Laplace-Beltrami operator. The proposed method significantly outperforms state-of-the-art algorithms on several non-rigid shape retrieval benchmarks.


Introduction
Three-dimensional models are ubiquitous data in the form of 3D surface meshes, point clouds, volumetric data, etc. in a wide variety of domains such as material and mechanical engineering [1], genetics [2], molecular biology [3], entomology [4], and dentistry [5,6], to name a few. Processing such large datasets (e.g., shape retrieval, matching, or recognition) is computationally expensive and memory intensive. For example, to query against a large database of 3D models to find the closest match for a 3D model of interest, one needs to develop an appropriate similarity measure as well as an efficient algorithm for search and retrieval. Shape descriptors assist with the example problem by providing discriminating feature vectors for shape retrieval [7,8] and play a fundamental role when dealing with shape analysis problems such as shape matching [9,10] and classification [11].
In general, there are two types of shape descriptors: local descriptors, also called point signatures, and global descriptors, referred to as shape fingerprints. A local shape descriptor computes a feature vector for every point of a 3D model. On the other hand, a global shape descriptor represents the whole 3D shape model in the form of a low-dimension vector. A descriptor that is informative and concise captures as much information as possible from the 3D shape including the geometric and topological features. Such a vector drastically lowers the shape analysis burdens in terms of both computational intensity and memory.
While a large number of successful non-spectral shape descriptors have been proposed in the literature, spectral descriptors have proved to be beneficial in many applications [12,13]. The spectral methods take advantage of eigendecomposition of the Laplace-Beltrami (LB) operator applied on the shapes and construct their informative descriptors using the eigenvalues and eigenvectors. These methods have found successful applications in graph processing [14], computational biology [15], and point-to-point correspondence [16].
One of the first spectral descriptors introduced to the computer graphics community is Shape-DNA, developed by Reuter et al. in 2006 [17]. Shape-DNA attracted a great deal of attention for its unique isometric and rotation invariant features [17]. Since then, several local as well as global shape descriptors have been introduced in accordance with Shape-DNA such as Heat Kernel Signature (HKS) [18], Wave Kernel Signature (WKS) [19], and Global Point Signature (GPS) [20]. The common ground between these methods is the discretization approach used to solve the Laplacian eigenvalue problem, which uses a cotangent weighting scheme along with area normalization.
Although there are many advantages of using variations of the cotangent scheme, there are several limitations. First, by their nature, they are limited to triangulated meshes. Second, they do not perform well when dealing with degenerate and non-uniform sampled meshes [21,22]. Also, their convergence error depends on factors such as the linearity of the function on the surface [22]. One possible approach to address these limitations is through the use of manifold learning.
Nonlinear dimensionality reduction techniques, known as manifold learning, assume the existence of a low-dimensional space, which represents a highdimensional manifold without much loss of information [23]. Similar to global descriptors, manifold learning methods attempt to learn the geometry of a manifold in order to extract a low dimensional vector of features that is informative and discriminative. However, unlike shape descriptors, the number of dimensions of a space does not confine manifold learning methods. To the best of our knowledge, the application of manifold learning, an active research topic in statistics and machine learning, has not been investigated in the computer graphics community for extracting global shape descriptors. This motivates the primary aim of this research, which is to explore the effectiveness of a manifold learning method, more specifically Laplacian Eigenmap [24], in representing a 3D model with a low-dimensional vector. Our work introduces a novel Laplacian Eigenmap-based global shape descriptor and provides a straightforward normalization method that significantly outperforms existing state-of-the-art approaches.
In our first contribution, inspired by the idea of Laplacian Eigenmaps [24], we learn the manifold of a 3D model and then, analogous to the approach taken by Shape-DNA, we make use of the spectrum of the embedded manifold to build the global shape descriptor. This approach has two main advantages. First, it relies on the adjacency of the nodes, disregarding the fine details of the mesh structure. Therefore, it can be used for degenerate or non-uniform sampled meshes. Second, as manifold learning does not rely on the mesh structure and is not limited to a specific type of meshes, e.g., triangulated meshes, it can be applied easily to any other mesh types such as quadrilateral meshes.
In our second contribution, we present a simple and straightforward normalization technique (motivated by [17,25,26]) to obtain a scale-invariant global shape descriptor that is more robust to noise. To this end, we propose to subtract the first non-zero eigenvalue from the shape descriptor after taking the logarithm of the spectrum. One advantage of our approach over the idea of Bronstein et al. [25] is that we avoid taking the direct derivative; this advantage is significant since the differential operator amplifies the noise. Taking the logarithm additionally helps to suppress the effect of the noise that is present in higher order elements of the spectrum.
The remainder of this paper is organized as follows. In Section 2, we briefly overview spectral shape analysis and manifold learning. Then in Section 3, we introduce the proposed shape descriptor along with some technical background. In Section 4, the performance of the proposed method, as well as the robustness of the algorithm are examined and compared with multiple well-known shape descriptors by performing several qualitative and quantitative experiments using widely used 3D model datasets. Section 5 discusses the results in more detail and draws conclusions.

Background
In this section, we first review spectral shape analysis, more specifically global shape fingerprints, and different discretization methods of the LB operator. Then, we briefly review manifold learning, more specifically Laplacian Eigenmap, to provide the necessary foundation for developing our proposed Laplacian Eigenmap based scale-invariant shape descriptor, which from now on we call LESI.

Spectral Shape Analysis
The Laplace-Beltrami operator ∆ is a linear differential operator defined on the differentiable manifold M as the divergence of the gradient of a function f as the following form [17,27]: ∆f = div(grad(f )). (1) Lévy [28] noted that the eigenfunctions φ i of the continuous LB operator, which are the solution to the following Laplacian eigenvalue problem: are the orthogonal basis for the space of functions defined on the surface of a manifold. In other words, a function f on the surface can be expressed as a sum over coefficients of these infinite bases: Furthermore, the LB operator is positive semi-definite, having non-negative eigenvalues λ i that can be sorted as follows: The sequence of eigenvalues of the LB operator is called the spectrum of the LB operator. As it is computed based on the gradient and divergence that depend on the Riemannian structure of the manifold, it possesses the isometry invariant property [17]. These significant features of the LB operator, which include the orthogonal basis and non-negative spectrum, motivated researchers to develop various local and global shape descriptors. The Shape-DNA and HKS were developed by considering the heat distribution as the function f on the surface. The WKS was obtained by solving the Schrödinger wave equation on the surface of the manifold. Also, it has been shown that the GPS descriptor is in close relation to the Green's function on the surface [20].
To approximate equation (2), despite the choice of function f , we need a discretization scheme. Different discretization schemas (e.g., Taubin [29], Mayer [30]) of the LB operator on the triangular meshes are discussed in [31].
A 3D shape, sampled from the surface of a Riemannian manifold M, is usually presented by a set of vertices V and their connectivity E in the form of the graph G = (V, E). For a surface mesh G, according to [21], the equation (2) can be discretized as: where A is the stiffness matrix and B is the lumped mass matrix. One popular approach to constructing the matrix A, is using weights: where α ij and β ij are the two angles facing the edge (i, j). Different mass normalization methods using the triangle area or the Voronoi region area are suggested to construct the matrix B. The cotangent weighting schema and its variants have been utilized in multiple FEM-based discretization methods. Another approach to constructing the matrix A is to use the heat kernel weight, also known as the exponential weighting scheme, as follows: where . denotes the Euclidean distance between two adjacent nodes i and j.
In [21], several existing discretization methods, including variants of linear FEM [32,33] and heat kernel weighting proposed in [22] are compared. According to [21] and from the discussion led by Xu in [31,34], discrete LB operator using cotangent weighting scheme may not converge in all situations, specifically when dealing with non-uniform meshes. However, the heat kernel weighting scheme proposed in [22] does not depend on the peculiarities of the triangulation and outperforms all linear approaches [18,21]. In addition, concerning the type of the function f , the cotangent scheme only converges for linear functions, while the heat kernel scheme converges well for nonlinear functions as well [22]. The proposed exponential approximation scheme provides point-wise convergence with good stability with respect to noise. It is important to note that although the method was discussed for surfaces without boundary, the results are valid for interior points of a surface with boundary [22].

Manifold Learning
To make the current contribution self-contained, we provide a brief introduction from the data analysis perspective. Dimensionality reduction of highdimensional data is a critical step in data analysis and processing. Non-linear dimensionality reduction, also known as manifold learning, is a problem of finding a low-dimensional representation for high-dimensional data. Several local and global manifold learning methods have been developed including Isomap [35,36], LLE [37,38], Laplacian Eigenmap [24], and Diffusion maps [39]. We review these next.
Consider a set x 1 , ..., x n ∈ M of n points on manifold embedded in R l . Manifold learning methods look for a set of corresponding points y 1 , ..., y n in R m (m l) as a structural representation, while respecting some local or global information. Each method attempts to minimize a cost function in this mapping.
Laplacian Eigenmap [40], proposed by Belkin and Niyogy in 2002, is a computationally efficient and mathematically well-studied manifold learning technique. It is based upon graph Laplacian and Laplace-Beltrami operator on the manifold. Accordingly, it is considered as a spectral analysis method. Laplacian Eigenmap deals with sparse, symmetric, and positive semi-definite matrices. It is in close connection to the heat flow [24,40].
Briefly speaking, for a given manifold, Laplacian Eigenmap applies the graph Laplacian operator and uses the eigenfunctions of such operator to provide the optimal embedding. Laplacian Eigenmap preserves local information by minimizing the distance between embedded points, which are mapped from adjacent points in the original high-dimensional space [24]. Aside from the locality preserving property, it provides structural equivalence and discrimination by capturing the intrinsic geometric structure of the manifold. The structural equivalence property states that two similar manifolds will have similar representation after projecting into a lower dimension space [41,42].
Some other manifold learning methods, e.g., Isomap, LLE, and Diffusion map are also based on spectral analysis of the high-dimensional manifold. In contrast to these methods that construct the orthogonal basis of their desired low-dimensional space using eigenfunctions of an LB operator, we develop our scale-invariant shape descriptor using the spectrum of the LB operator. Even though our primary focus is on the Laplacian Eigenmap, owing to its unique properties, we believe that other spectral manifold learning methods are also capable of extracting informative and discriminative shape fingerprints.

Method
In this section, we elaborate our proposed LESI global shape descriptor. A flowchart of the proposed approach is shown in Figure 1.

Laplacian Eigenmap-Based Shape Descriptor
We treat a global descriptor as a dimensionality reduction problem as it squeezes the latent information of a 3D model into a vector. Due to the fact that the Laplacian Eigenmap has two properties of structural equivalence and locality preservation [24], we propose a global shape descriptor using the spectrum of graph Laplacian. A graph Laplacian is constructed over an undirected weighted graph G = (V, E) with a set of points x i ∈ V and a set E of edges that connects nearby points ((i, j) ∈ E). The theory behind finding the optimal embedding in a Laplacian Eigenmap requires an undirected graph. Every 3D model is given in bidirectional connection and hence, we need to neither examine nor force it to the graph. However, as we will explain later, we need to remove isolated points. Considering the advantages of the heat kernel weighting scheme, which are summarized in Section 2.1 and discussed in details in [22], Laplacian Eigenmap suggests constructing the weighted graph as follows: The only parameter in equation (6) is t, which defines the extent to which distant neighbors influence the embedding of each point. The choice of parameter t is data-dependent, and there exists no unique way in the literature to select the proper value, but it can be tuned empirically. As t has neither a very high impact on the final embedding nor the convergence rate of our final derivations, we empirically recommend: Here, weights are bounded as 0.6 ≈ e −0.5 ≤ w ij ≤ 1. Laplacian Eigenmap attempts to find a low dimensional data set that preserves local information. For this purpose, it assumes two neighboring points x i and x j stay close after being mapped to y i and y j . Therefore, it minimizes the following function [24]: where L = D − W is the so called Laplacian matrix, W ij is a symmetric weight matrix, and D ii = j W ij , the degree matrix, is a diagonal matrix. The assumption that the graph is undirected yields the symmetric property of W and consequently, D and L. It plays a critical role in deriving equation (8).
By adding the orthogonality constraint y T D1 = 0 in order to eliminate the trivial solution and the constraint y T Dy = 1 for removing an arbitrary scaling factor in the embedding, the minimization problem (8) simplifies to: The matrix L is real, symmetric, and positive semi-definite. Therefore, the solution vector y (in equation (9)) is obtained by the minimum eigenvalue solution to the generalized eigenvalue problem [24]: At this point, the optimal low dimensional embedding, suggested by the Laplacian Eigenmap, is obtained by utilizing the eigenvectors. However, we focus on the spectrum of the graph Laplacian and its' properties. Eigenvalues obtained from equation (10) are real, non-negative, and sorted in increasing order as follows: As the row (or column) sum of L is equal to zero, eigenvalue λ = 0 and a corresponding eigenvector 1 are trivial solutions to equation (10). The multiplicity of eigenvalue zero is associated with the number of connected components of the graph. Eigen-solvers often obtain very small, though not precisely zero, eigenvalues due to the computational approximations. If we may know the number of connected components of L, we can discard all eigenvalues equal to zero, and form our shape fingerprint using the more informative section of the spectrum. This is easily done by Dulmage-Mendelsohn decomposition [43].
The second smallest eigenvalue, also known as the Fiedler value, is a measure of the connectivity within the graph. If the graph has c connected components, our proposed shape descriptor is a set of d eigenvalues as: The LESI descriptor is composed of the spectrum of the LB operator, and hence, it is isometric invariant, independent of the shape location, and informative. The latter, discussed in spectral graph theory, states that the spectrum of graph Laplacian contains a considerable amount of geometrical and topological information of the graph. Moreover, LESI has the similarity property, caused by the structural equivalence property of Laplacian Eigenmap, meaning that two 3D models from the same class of models have similar fingerprints. Unlike Shape-DNA and other shape descriptors that are based on the cotangent weighting scheme, LESI is not limited to triangulated mesh structures because the Laplacian Eigenmap is capable of dealing with high-dimensional data. For some applications in which scale is not a determinant factor, it is favorable to have a scale-invariant descriptor. A fast and efficient normalization method is presented in Section 3.2.
One important matter to consider is the convergence and accuracy of the proposed fingerprint, which ultimately depends on the heat kernel-based discretization of the LB operator. The cotangent weighting scheme and its variants are sensitive to the peculiarities and quality of the particular triangulation of the mesh (refer to Section 2.1). While an exponential weighting scheme has shown accurate performance in dealing with nonlinear functions over the surface and non-uniform mesh representations, it is not clear how this method can handle manifolds with boundaries [21,22]. It does, however, behave well for interior points of the surface. Therefore, we recommend removing rows and columns of L and D corresponding to isolated points, before solving equation (10). The descriptor obtained from the rest of the connected graph is an informative and discriminative descriptor of the graph.

Scale Normalization
For some applications, the size of an object is not a determinant factor in shape comparison and identification. Therefore, a scale-invariant shape descriptor with a solid normalization method is more desirable. For that purpose, some shape descriptors including Shape-DNA, have proposed multiple normalization methods. Most normalization methods of Shape-DNA focus on finding an appropriate scaling factor, such as the surface area, the volume, or coefficient of a fitting curve, which will be multiplied in the descriptor. Moreover, it is shown that eigenvalues with a higher order are more susceptible to noise. For that reason, the original Shape-DNA recommends cropping the spectrum and using no more than 100 eigenvalues [17]. In this section, we propose a simple and efficient normalization method that significantly reduces the effect of scale variations as well as noise. In this approach, we are interested in taking the scaling factor out from the descriptor in one step, rather than an extra step to find an appropriate neutralizing factor. Although the normalization seems simple, later in the experiments section, we will show its efficiency. The work presented in [18] and [17] influenced this method.
According to Weyl's law [44], when an arbitrary scale object is scaled with factor α, the eigenvalues are scaled with factor α 2 . Let Λ s = (α 2 λ 1 , . . . , α 2 λ d ) be a LESI descriptor obtained from a scaled 3D model with unknown factor α. To normalize the descriptor and eliminate the effect of the scaling factor, we recommend computing: To achieve a scale invariant shape descriptor, we first take the logarithm of the descriptor vector and then compute the difference of the new vector from its smallest element. This is equivalent to dividing the vector by its first element and taking the logarithm next. Basically, division takes the factor out, and the logarithm eliminates the influence of noise. Figure 2 illustrates the details of the proposed algorithm in an example. In Figure 2(a), two Teddy Bear models are shown. One model is in the original size whereas the other model is scaled with factor 0.7. It is clear from Figure  2(b) that the spectrum of the scaled model is almost (0.7) 2 of the spectrum of the original model. Taking the logarithm of the spectrum takes away the scaling factor from multiplicand and leaves it as augend. Therefore, subtracting one term (e.g., the first element) removes the scaling factor from all other terms. The result is a normalized and scale-free spectrum.

Algorithm
Our proposed descriptor consists of three major steps. For a given 3D polygonal model G = (V, E) with a set of vertices V = (x 1 , . . . , x n ) ∈ R 3×n and a set of neighbor connections E, the LESI descriptor is a d-dimensional vector Λ = (λ 1 , . . . , λ d ) of real and positive values.
In the first step, we compute the n × n real, symmetric, and sparse weight matrix (W ) for a 1-ring neighbor of every point as stated in equation (6) using the inner scaling factor given in equation (7). Next, we form the generalized eigenvalue problem equation (10) by constructing Laplacian and degree matrices (L and D respectively) without difficulty. The matrix L is sparse, real, symmetric and semi-positive. Utilizing the Dulmage-Mendelsohn decomposition, we find the number of connected components of L. The objective of the second step is to find the spectrum of the LB operator. For that purpose, we solve the generalized eigenvalue problem using the Lanczos method. Then, we leave out as many smallest eigenvalues as the number of connected components. Since in most cases a single 3D model is made up of one connected component, we only need to leave out one eigenvalue. The last step of the algorithm deals with scale normalization and noise reduction, in case it is required, by taking the logarithm of the spectrum and subtracting the first element from the rest of the vector. Detailed steps of the algorithm are summarized in Algorithm 1.

Algorithm 1: Laplacian Eigenmap-based scale-invariant global shape descriptor
Input: A 3D polygonal model G = (V, E) with n vertices x i and edge list E Output: A d-dimensional vector Λ = (λ 1 , . . . , λ d ) 1 Compute edge weights using equations (6) and (7); 2 Construct the sparse, real, and symmetric n × n matrices W , D, and L; 3 Find number of connected components (nConComp) from L; 4 Solve equation (10) for nConComp + d eigenvalues; 5 Sort them in increasing order and leave out nConComp smallest ones; 6 if normalization is required then

Experiments
In Section 4.1, we first present two datasets used in our experiments. Then, in Section 4.2, we qualitatively visualize and measure the competence of the proposed method in discriminating different clusters compared with candidate methods from the literature. Next, in Section 4.3, we validate the effectiveness of the LESI descriptor to distinct multiple classes by measuring the accuracy of multi-class classification. Finally, in Section 4.4, extensive experiments are carried out to study the robustness of the proposed shape descriptor with respect to noise, scale invariance, and down-sampling. All the algorithms were implemented using the MATLAB R2013a environment running on a personal computer with Intel(R) Core(TM) i3-4130 CPU @ 3.40GHz and 12GB memory.

Dataset
To validate the utility of the proposed shape descriptor, we utilized two standard, widely-used, and publicly available datasets of 3D polygon meshes. The high-resolution TOSCA dataset [45] contains 80 three-dimensional non-rigid models, including 11 cats, 6 centaurs, 9 dogs, 4 gorillas, 8 horses, 12 women poses, 3 wolves and two men with 7 and 20 poses respectively. In our experiments, we use all models except the gorilla models, as they contain isolated points. The models in each class of the TOSCA dataset are almost identical in terms of scale, the number of vertices, quality of triangulation, and structure, which all represent the same object with different poses.
The McGill dataset with articulating parts [46] is used to evaluate the ability of the descriptors to describe models with poor intra-class quality. The McGill dataset contains 3D models of 30 ants, 30 crabs, 25 glasses, 20 hands, 29 humans, 25 octopuses, 20 pliers, 25 snakes, 31 spiders, and 20 Teddy bears. The classification of the McGill dataset models is more challenging due to scale and shape variations.

Retrieval Results
In this section, we evaluate the general performance of our proposed shape descriptor and compare it with several state-of-the-art spectral-based global shape descriptors including Shape-DNA [17], cShape-DNA [27], and GPS [20] algorithms. We chose these methods because they are widely used by researchers (e.g., [47][48][49]) to develop new descriptors or applications, or to evaluate the performance of their proposed descriptors. Moreover, cShape-DNA represents the normalized version of the original Shape-DNA. Even though there are multiple ways to convert a local point descriptor to a global shape fingerprint, in this article we focus only on algorithms that have been originally introduced as global fingerprints. To this end, we take advantage of the source code made available on Dr. Kokkinos's homepage 1 [25], as well as the shape descriptor package provided by Li et al. [50] available on a GitHub repository 2 to generate the Shape-DNA and GPS descriptors, respectively. We also compare the performance of shape retrieval using the code provided for evaluation by SHREC'11 [51].
The shape descriptors are compared using the TOSCA dataset to discriminate between different classes of 3D objects. In this experiment, we use the first 33 non-zero eigenvalues (d = 33). Then, to visualize the locations of objects in the shape space, we project them onto a 2D plane using Principle Component Analysis (PCA). Figure 3 displays the effectiveness of our method compared with the fingerprints of interest.   To demonstrate the power of our method in classifying objects with low intraclass similarity compared with other shape descriptors, the same experiment is carried out on the McGill dataset. Models of the same class with articulating parts are in different scales, shape, and structure. The 2D PCA projections of 33-dimension descriptors from all four algorithms are shown in Figure 4. As illustrated in Figure 4, the original Shape-DNA is highly sensitive to scales. Multiple methods are presented in [17] to make the descriptor normalized to scale. cShape-DNA represents a normalized version of it by multiplying the descriptor with the surface area. Although cShape-DNA can separate models from each other, classes are not separated efficiently. LESI outperforms the other algorithms by providing distinct descriptors, which can separate classes. Shape descriptors offered by LESI prove superior to the other algorithms in the shape retrieval and classification tasks, as described below and in the next section respectively.
To examine the superiority of LESI quantitatively, we computed multiple standard retrieval measures including Nearest Neighbor (NN), First Tier (FT), Second Tier (ST), e-Measure (E), and Discounted Cumulative Gain (DCG). These measures represent state-of-the-art quality metrics used when evaluating matching results for shape-based search engines [52].

Multi-class Classification Results
In this section, we corroborate the findings of Section 4.2 by training a linear multi-class SVM classifier to assess the accuracy of LESI compared to other shape descriptors. For this experiment, we utilized the McGill dataset. In addition to the shape descriptors evaluated in Section 4.2, we computed another normalized version of Shape-DNA by dividing the feature vector by its first element (similar to what LESI offers) as suggested in [17]. This way we can compare the effect of the exponential weighting scheme without the influence of the normalization method or compactness (offered by cShape-DNA). Using 10-fold cross-validation and repeating the experiment 3 times, we report the average accuracy for each method in Table (2).
The new LESI approach significantly outperforms all other methods when using a two-tailed paired t-test (p < 0.05). The t-test was performed on one set of 10 folds in order to avoid violating the independence assumption of the t-test. There is a significant improvement in accuracy when comparing the Shape-DNA (Normalized) to other variants of the Shape-DNA, which is due in part to the normalization method. However, the average accuracy of the LESI descriptor is noticeably higher (95%) when compared to 90% of the Shape-DNA (Normalized). Finally, Figure 5 shows the confusion matrix obtained from the linear multiclass SVM using LESI descriptor. The number of correct classifications made for each class (indicated by the green diagonal), confirms that our method captures the discriminative features of the shapes.

Robustness
In this section, we address the robustness of the LESI shape descriptor to shape variations, including noise, scale, and down-sampling by performing another set of experiments. First, we generate the disturbed version of every model in the TOSCA dataset. Then, we test the capability of every method mentioned above in discriminating between different classes. For this purpose, besides plotting the 2D PCA projection of shape descriptors, we also compute and plot the pairwise Euclidean distance matrix, in every case. The distance matrix represents the dissimilarity between each pair of models in the set. It is often used to compute other evaluating metrics such as nearest-neighbor, and first and second tier, to name a few. The dissimilarity of descriptors increases from blue to red, and the more separate classes differ in color, the better they are discriminating from each other.
Resistance to noise. Multiple noisy versions of the TOSCA dataset are generated following the idea articulated in [53]. To this end, the surface meshes of all models are disturbed by changing the position of each point along its normal vector that is chosen randomly from an interval (−L, L) with the 0 mean, where L determines the noise level and is a fraction of the diagonal length of the model bounding box. In this experiment, three noise levels L = 0.5%, L = 1%, and L = 2% are tested, where the latter one represents a greater level of noise. Two-dimensional PCA projections of all descriptors with the presence of different levels of noise are plotted in Figure 6. Combining these with the results shown in Figure 3, where no noise is present, demonstrates that the LESI algorithm is highly noise-resistant while the performance of the Shape-DNA and cShape-DNA decreases as the level of noise increases. Moreover, GPS fails in separating different classes of models with the presence of noise.  reflects the effect of noise on the discriminative power of the descriptors. The LESI algorithm shows consistent results as the level of noise increases from 0% (top row) to 2% (bottom row).
Scale invariance. In order to validate the insensitivity of the LESI descriptor to scale variations and compare the robustness of the proposed method with other descriptors, each model of the TOSCA dataset is scaled by a factor of 0.5, 0.875, 1.25, 1.625, or 2 randomly. Figure 7 shows that the LESI algorithm surpasses other methods in discerning different classes. Comparing the result of this experiment with the results shown in Figure 3 demonstrates the consistency of the LESI algorithm with the presence of scale variation. The distance matrices in Figure 10 show that the original Shape-DNA algorithm is very susceptible to scale variations. Even though the cShape-DNA has significantly improved scale sensitivity of the original Shape-DNA, it does not provide as accurate results as the LESI algorithm does.
Resistance to the sampling rate. To investigate the effect of sampling rates on the discriminative power of the shape descriptors, Bronstein et al. [8] propose to reduce the number of vertices to 20% of its original size. Accordingly, the down-sampled version of the TOSCA dataset is generated, and shape descriptors associated with them are computed. The 2D PCA projections and distance matrices of descriptors are illustrated in Figures 8 and 11, respectively. Although the original Shape-DNA shows a more accurate result than cShape-DNA, the separation of cat, dog, and wolf models is challenging. Although  the performance of the LESI method is slightly affected, it still outperforms cShape-DNA and GPS methods.

Discussion
In this article, motivated by the unique properties of Laplacian Eigenmap (i.e., locality preservation, structural equivalence, and dimensionality reduction) and inspired by the existing spectral-based shape descriptors, we investigated the application of manifold learning in deriving a shape fingerprint in order to address the limitations tied to popular cotangent-based shape descriptors. We proposed a global descriptor (LESI) with an easy-to-compute and efficient normalization technique that facilitates applications such as shape classification and retrieval. Our method applies fewer restrictions on the class of meshes as well as improving the quality of tessellations. Analogous to other spectral descriptors, LESI uses the spectrum of the LB operator, which is independent of the shape location, is informative (contains a considerable amount of geometrical and topological information), and above all isometric invariant. We compared the discriminating power of LESI with three prominent descriptors from the literature, namely Shape-DNA, cShape-DNA, and GPS, and found it to be superior.
In the first set of experiments illustrated in Figures 3 and 4, our method sub-
Shape-DNA cShape-DNA GPS LESI Random Scale stantially outperforms the others. The superiority of LESI is more significant when the McGill dataset is used (Table 2 and Figure5). This dataset includes wide variations in mesh structure and scales, causing the failure of the other methods to generate acceptable results. However, LESI, due to utilizing a different method of discretization to form the LB operator, focuses on the vicinity rather than the quality of the triangulation. Therefore, our technique, unlike Shape-DNA cShape-DNA GPS LESI Down Sampled Figure 11: The Euclidean pairwise distance matrix of shape descriptors computed by (from left to right) Shape-DNA, cShape-DNA, GPS, and LESI algorithms from from down sampled TOSCA dataset by rate of 20%. other methods, is not affected by the low quality of polygon meshes.
The second set of experiments evaluates the reliability of our method in the presence of noise, scale variations, as well as different sampling rates. LESI shows impressive robustness against the first two sets of perturbation. Despite the negative impact of down sampling in LESI descriptor, it continues to show better performance when compared to cShape-DNA and GPS. It should be noted that the result could also be improved by increasing the size of the output vector.
In addition to the discriminating power of the descriptor, degenerate and non-uniform meshes may also cause failure of an algorithm to converge. The cotangent weight-based algorithms were not able to compute the descriptors for 2 shapes from the McGill dataset. GPS also failed to compute descriptors for 6 models of the down sampled TOSCA dataset. However, our technique converges at all times despite the quality of the polygon mesh structure.
Moreover, LESI, unlike cotangent weight-based techniques, is not confined to the triangulated meshes as it disregards the mesh geometry [54]. LESI inherits this property from the capability of manifold learning techniques in coping with high dimensional data. The discretization of the LB operator using cotangent weights on the quadrilateral meshes is not as straightforward as on triangular meshes. To compute the LB operator on a quadrilateral mesh, all rectangles need to be divided into triangles. It could be done easily, however, as for each quad there are two possible triangulations, the result is not unique.
In the original Laplacian Eigenmaps, the high dimension data requires a considerable amount of processing as the list of all connections need to be computed for the dataset. In fact, for each point in the high dimension space, a given number of nearest neighbors need to be extracted which could be challenging and unmanageable. While applying this technique to the 3D meshes, we skip this step as the neighbors are already defined and given in the mesh structure.
This work benefits from the Laplacian Eigenmap technique in a space in which the vicinities are given. LESI takes advantage of simple Laplacian computation, to form the LB operator, which provides concise and informative shape descriptors. Experimental results prove that LESI is more effective compared with the other powerful descriptors.
One limitation of LESI is its inability to separate models of different men (David and Michael). However, it was able to differentiate between the wolf and dog, as well as between women and men.
Although we investigated only the application of Laplacian Eigenmap in introducing a shape descriptor, there are some other spectral-based manifold learning methods, such as Isomap, LLE, and Diffusion map, which have not been examined. This can be considered future work.

Conclusion
This work presents LESI, a novel scale-invariant global shape descriptor based on Laplacian Eigenmap that is significantly better when compared to other shape descriptors. We conclude that manifold learning methods can be used to develop new spectral-based shape descriptors to learn the structure of manifolds despite the quality of sampled meshes.