1 Introduction
One of the most challenging problems in the domain of 2D image or 3D shape is to handle the nonrigid deformation, especially in the situation of anisotropy, which is universal in the real world. In the viewpoint of transformation groups, the isometric transformation is a prop subgroup of the conformal transformation, which is a prop subgroup of the diffeomorphism. Obviously, the anisotropic nonrigid transformation exceeds the boundary of isometric transformation and contains conformal transformation. Based on the Erlangen program of Klein, geometry is a discipline that studies the properties of space that remain unchanged under a particular group of transformation. In order to solve the anisotropic transformation problem, it is necessary to find the invariants under the conformal transformation.
The original motivation of conformal mapping is how to flatten the map of globe, and the Mercator projection produce an anglepreserving map that is very useful for navigation. More generally, the conformal geometry focuses on the shape in which the only measure is angle instead of usually length. The descriptions of conformal mapping contain angle preservation[12, 26, 5], metric rescaling[21, 27], preservation of circles[14, 28], etc. Some key ideas reside in the conformal surface geometry are Dirac equation[6], CauchyRiemann equation[22], Möbius transformations[27, 28], Riemann mapping[10, 9, 35, 33], Ricci flow[34], etc. The conformal geometry lies between the topology geometry and the Riemannian geometry, it studies the invariants of the conformal transformation group. The conformal structures[10, 9] based on the theories of Riemann surfaces are invariants under conformal transformation. According to conformal geometry[7], the shape factor[10] and conformal module[35] are conformal invariants. Moreover, the conformal inner product[26] defined by an inner product of function is also changeless under conformal transformation. According to the Liouville Theorem[20, 8], the Möbius transformation plays an important role in conformal mapping.
The definition of Möbius transformation[25]
shows that it is compounded by a series of simple transformations: Translation, Stretching, Rotation, Reflection and Inversion. In the domain of invariants under translation, stretching and rotation transformations, the Geometric moment invariants(GMIs)
[32] and the ShapeDNA[17] show a general method to generate the moment invariants; Hu et al[13] proposed a general construction method of surface isometric moment invariants based on the intrinsic metric. In the domain of invariants under reflection transformation, the chiral invariants[36] show the moment invariants based on the generating functions of ShapeDNA[17]. In the domain of invariants under conformal transformation, Hu[12]proposed limited conformal invariants based on geodesic tangent vectors. In the domain of invariants under Möbius transformation, the expression
proposed by Blaschke[1] is proved to be a conformal invariant by Chen[4]; based on the GaussBonnet Theorem, White[30] proposed that is a global conformal invariant if is an oriented and closed surface. The GaussBonnet Theorem associates the differential expression(Gaussian curvature) of the surface with its topological invariant (the Euler’s characteristic). This great theorem motivates us to explore the differential invariants under the Möbius transformation since the differential expressions play essential roles in some procedures of physics, mathematics, computer science and other fields. In the domain of differential invariants, rotation and affine differential invariants were proposed by Olver[23] based on the moving frame method; a special type of affine differential invariants was presented by Wang et al[29]; Li et al[19] prove the existence of projective moment invariants of images with relative projective differential invariants; the research[18] on the relationship between differential invariants and moment invariants show that they are isomorphic under affine transformation.In this article, we study invariants by combining functional map[24] and the derivatives of function(see Fig.1). In section 2, we show the background of this paper. In section 3, we propose the invariants under Möbius transformation. In section 4, we show another Möbius invariant from the functional view. Finally, we propose a conjecture about the structure of differential invariants under conformal transformation. The main contributions of this paper are as follows.

We propose two differential expressions that are invariant under 2D and 3D Möbius transformation respectively. According to the Liouville Theorem, the 3D differential invariant is a conformal invariant.

Based on the analysis on absoluteness and relativity of invariance about the two differential expressions and their components, we propose integral invariants under Möbius transformation.

We propose a conjecture about the composition of differential invariants under conformal transformation.
2 Notion and Background
2.1 Notion
The formulation in this paper is same with the functional maps framwork[24]. Assuming and are two manifolds, a bijective mapping induces the transformation of derived quantities, where is scalar function defined on manifold. It means that any function have a counterpart function and .
To make the invariants under Möbius transformation clear, we partially modify original definition and theorem in this paper with this formulation.
2.2 Theoretic Background
According to the Liouville Theorem[20], the only conformal mapping in are Möbius transformation[11, 25, 15]. Furthermore, the Generalized Liouville Theorem shows that any conformal mapping defined on () must be a restriction of Möbius transformation.
Theorem 2.1 (Generalized Liouville Theorem[8])
Suppose that , are domains in and that is a homeomorphism. If , then is 1quasiconformal if and only if or its complex conjucate is a meromorphic function of a complex variable in . If , then is 1quasiconformal if and only if is the restriction to of a Möbius transformation, i.e., the composition of a finite number of reflections in spheres and planes.
Next, we will show the common expressions of Möbius transformation in different dimensions().
In the filed of complex analysis, a Möbius transformation could be expressed as
(1) 
where , . Based on the Liouville Theorem[20], every Möbius transformation in higher dimensions could be given with the form
(2) 
where , is 0 or 2, and
is an orthogonal matrix. The choice of
decides if contains inversion transformation, and the sign of decides if contains reflection transformation.More generally, a Möbius transformation could be composed of a series of simple transformations, the definition of Möbius transformation is as below.
Definition 1 (Möbius transformation[25])
A ndimension Möbius transformation is a homemorphism of (the onepoint compactification of ), it is a mapping that is a finite composition of the following elementary transformations():
(1)Translation: , .
(2)Stretching: , and .
(3)Rotation: , and is an orthogonal matrix.
(4)Reflection about plane : , is the normal vector of , is the distance from the origin to .
(5)Inversion about sphere : , is the inversion center, is the inversion radius.
3 Möbius Invariants
3.1 Inversion Invariants
In order to derive the differential invariant under inversion transformation , in the 2D situation we assume that the map the function on domain to on domain , where and , this means that the coordinates transformations under are as follows.
(3) 
(4) 
At the same time, it means that the coordinates transformations under are as follows.
(5) 
(6) 
Based on and the equations (5)(6), we obtain the relationships between the partial derivatives of and as follows.
(7)  
(8)  
(9)  
(10)  
(11) 
Then we obtain a 2D equation under the inversion transformation, it is
(12) 
This means that
(13) 
is a differential invariant under inversion transformation. We use the same method in 3D situation and obtain a differential invariant under the inversion transformation, it is
(14) 
where
(15)  
3.2 The Boundary of Invariance
We have shown that (13) and (14) are differential invariants under inversion transformation. It is obvious that they are invariants under translation transformation. We prove that (13) and (14) are also differential invariants under rotation, stretching and reflection transformations(see Appendix A for a proof). According to the definition of Möbius transformation, we conclude that the differential expression (13) is a differential invariant under 2D Möbius transformation. Furthermore, with the Generalized Liouville Theorem we obtain that (14) is a conformal invariant.
3.3 Absoluteness and Relativity of Invariance
If expression is an invariant under transformation , the transformed expression satisfies
(16) 
where is an expression related to . In this context, is an absolute invariant if , otherwise, is a relative invariant. Base on the analysis in 3.2, (13) is an absolute invariant under Möbius transformation and (14) is an absolute invariant under conformal transformation. Next, we will show the numerator and denominator of (13) or (14) are relative invariants.
In the derivation of 2D inversion invariants, we obtain that for the numerator and denominator of (13), this means
(17) 
(18) 
where is the determinant of Jacobian matrix of transformation , is the absolute valve of . In 3D situation, we obtain for the numerator and denominator of (14). In the stretching transformation, we obtain in 2D situation, and in 3D situation. We also obtain that for the numerator and denominator of (13) or (14) under translation, rotation and reflection transformations.
3.4 Multiscale and Quantity
Assuming is a regular parameter surface defined on , if transform defined on to defined on and , based on the change of variable theorem[16] for multiple integrals and Table 1 we obtain that
(19) 
(20) 
where is the area extension factor, so we obtain that
(21) 
(22) 
are integral invariants under 2D Möbius transformation. In the same way, we obtain that
(23) 
(24) 
are integral invariants under 3D conformal transformation.
Actually a differential expression of function defined on domain accurately characterize at point of , it provides extremely wide space to describe the function .
Multiscale of Invariants
Assuming is a function of , a general method to construct descriptors in different scale is the integral of on region () with different size, and when the result is a global invariant, for example, the Willmore energy [1] applied in the theory of surfaces[31], digital geometry processing[2] and other fields.
In this view, the only difference between invariant with specifyscale and global invariant is the definition domain, the construction method of specifyscale invariant is same with global invariant. The former could be elaborately modified by selecting domain of integration in different applications.
Quantity of Invariants
A general method to construct a large number of invariants is using various functions with these functions are independent of each other[3]. We just show a simple method to construct integral invariants based on differential invariants and integral, in addition, more invariant forms can be constructed with differential invariants. Next, we give a possible form of invariants under Möbius transformation:
(25) 
(26) 
(27) 
(28) 
3.5 Another Conformal Invariant
The expression proposed by Biacchke[1] has been proved to be an invariant under Möbius transformation[4, 30]. It differs from our method in two important respects: the domain of transformation and the number of functions participated in invariants(see detailed expression at Appendix B).
4 Conjecture of Conformal Invariants
We have shown that (13) is a Möbius invariant and (14) is a conformal invariant. However, the fascinating part of (13) or (14) is that the differential expressions
(29) 
(30) 
(31) 
are differential invariants under rigid transformation. Based on this observation and the fact that the differential expressions play important roles in transformation, we have a bold conjecture about the structure of differential invariants under conformal transformation.
Conjecture:
The differential invariants under conformal transformation are composed of differential invariants under rigid transformation in a selfconsistent manner.
One of the possible selfconsistent forms in ndimensional Euclidean space may be
(32) 
where is differential invariant under rigid transformation.
5 Experimental Results
We choose a human face model from TOSCA database and treat the zcoordinate value of vertexes of the triangle mesh as a function defined on xcoordinate and ycoordinate, i.e.
. With least square method, the coordinates of a vertice and its 1ring neighbors were used to estimate parameters in Taylor expansion of
at the vertice; in order to guarantee the accuracy of descriptor calculation, we only consider vertexes that are located inside the mesh and have enough 1ring neighbors. After that, we calculate a descriptor at the vertice and the descriptor is composed by (13), (25) and (26) with different . Moreover, in integral invariants, the area around a vertice is determined by Mixed Voronoi cell.We deform the definition domain of with reflection, stretching, rotation and inversion transformation(Fig.3). In reflection transformation, and ; the in stretching transformation is 2; in rotation transformation the original data is rotated 90 degrees counterclockwise; in inversion transformation the inversion center is and inversion radius is 500(see more explanation about experiments at Appendix C).
5.1 Stability of Invariants
In this experiment we choose and the integral invariants is calculated at the local area of each vertex. After we obtain a 5dimension descriptor at vertexes of the five mesh in Fig.3, we calculate the average error of each dimension of the descriptor. In addition, we choose an isometric invariant at the vertex, the Laplacian operator, to compare with above invariants. The average error of each dimension is calculated by the following formula
(33) 
where is the value of invariant at vertex on original data, is the value of invariant at vertex on deformed data, and is the total number of vertexes participated in the calculation. The result of this experiment is in Table 2, it shows that (13) ,(25) and (26) are invariants under Möbius transformations.
5.2 Discrimination of Invariants
In this experiment we use the 5dimension descriptor of vertex at original to match its corresponding vertex in the deformed mesh with nearest neighbor rule, the metric between vertexes is standardized Euclidean distance. The error rate (percentage) of this experiment is in Table 3.
In conformal deformation scenario, this experiment shows the potential of Möbius invariants in matching task. Fig.4 shows some matchingfail situations, where the white point is the real position and the red point is the matching vertex. The reason for most matching failures is that the original white vertex and deformed red vertex have similar functional distribution environments.
6 Conclusions
In this article, we propose two differential invariants under 2D and 3D Möbius transformation respectively, in particular, the 3D expression is a conformal invariant according to the Liouville Theorem. After that, we analyze the absoluteness and relativity of invariance on the two expressions and their components, and we show an integral construction method that targets to the multiscale and quantity of invariant, the experimental results show that the invariants proposed in this paper perform well. Furthermore, we show another Möbius invariant from the functional view. Finally, we propose a conjecture about the structure of differential invariants under conformal transformation.
This article shows a method of combining functional map and derivatives of function to study conformal invariant, more research about the differential invariants under conformal transformation is necessary in the future. In addition to practical application solutions based on Möbius invariants, questing the generative structure of conformal differential invariant is also an interesting topic.
7 Acknowledgment
The authors would like to thank Dr. Antti Rrasila of Aalto University for providing help on how to distinguish Möbius invariants and conformal invariants.
This work was partly funded by National Key R&D Program of China (No. 2017YFB1002703) and National Natural Science Foundation of China (Grant No.60873164, 61227802 and 61379082).
References
 [1] Biaschke, W.: Vorlesungen über differentialgeometrie iii (1929)
 [2] Bobenko, A.I., Schröder, P.: Discrete willmore flow (2005)
 [3] Brown, A.B.: Functional dependence. Transactions of the American Mathematical Society 38(2), 379–394 (1935)
 [4] Chen, B.Y.: An invariant of conformal mappings. Proceedings of the American Mathematical Society 40(2), 563–564 (1973)
 [5] Corman, E., Solomon, J., BenChen, M., Guibas, L., Ovsjanikov, M.: Functional characterization of intrinsic and extrinsic geometry. ACM Transactions on Graphics (TOG) 36(2), 14 (2017)
 [6] Crane, K., Pinkall, U., Schröder, P.: Spin transformations of discrete surfaces. ACM Transactions on Graphics (TOG) 30(4), 104 (2011)
 [7] Farkas, H.M., Kra, I.: Riemann surfaces. In: Riemann surfaces, pp. 9–31. Springer (1992)
 [8] Gehring, F.: Topics in quasiconformal mappings. In: Quasiconformal Space Mappings, pp. 20–38. Springer (1992)
 [9] Gu, X., Wang, Y., Yau, S.T., et al.: Computing conformal invariants: Period matrices. Communications in Information & Systems 3(3), 153–170 (2003)
 [10] Gu, X., Yau, S.T.: Surface classification using conformal structures. In: null. p. 701. IEEE (2003)
 [11] Haantjes, J.: Conformal Representation of an Ndimensional Euclidean Space with a Nondefinite Fundamental Form on Itself (1937)
 [12] Hu, P.: A Class of Isometric Invariants and Their Applications(in Chinese). Ph.D. thesis, Institute of Computing Technology, Chinese Academy of Sciences (May 2011)
 [13] Hu, P., Li, H., Lin, Z.: A construction method for surface isometric invariants. Journal of Systems Science and Mathematical Sciences 9, 006 (2009)
 [14] Kharevych, L., Springborn, B., Schröder, P.: Discrete conformal mappings via circle patterns. ACM Transactions on Graphics (TOG) 25(2), 412–438 (2006)
 [15] Kühnel, W., Rademacher, H.B.: Liouville’s theorem in conformal geometry. Journal de mathématiques pures et appliquées 88(3), 251–260 (2007)
 [16] Lax, P.D.: Change of variables in multiple integrals. The American mathematical monthly 106(6), 497–501 (1999)
 [17] Li, E., Huang, Y., Xu, D., Li, H.: Shape dna: Basic generating functions for geometric moment invariants. arXiv preprint arXiv:1703.02242 (2017)
 [18] Li, E., Li, H.: Isomorphism between differential and moment invariants under affine transform. arXiv preprint arXiv:1705.08264 (2017)
 [19] Li, E., Mo, H., Xu, D., Li, H.: Image projective invariants. IEEE Transactions on Pattern Analysis and Machine Intelligence (2018)
 [20] Liouville, J.: Extension au cas des trois dimensions de la question du tracé géographique. Applications de l’analyse à la géométrie pp. 609–617 (1850)
 [21] Luo, F.: Combinatorial yamabe flow on surfaces. Communications in Contemporary Mathematics 6(05), 765–780 (2004)
 [22] Mullen, P., Tong, Y., Alliez, P., Desbrun, M.: Spectral conformal parameterization. Computer Graphics Forum 27(5), 1487–1494 (2008)
 [23] Olver, P.J.: Equivalence, invariants and symmetry. Cambridge University Press (1995)
 [24] Ovsjanikov, M., BenChen, M., Solomon, J., Butscher, A., Guibas, L.: Functional maps: a flexible representation of maps between shapes. ACM Transactions on Graphics (TOG) 31(4), 30 (2012)
 [25] Rasila, A.: Introduction to quasiconformal mappings in nspace. Proceedings of the International Workshop on Quasiconformal (2006)
 [26] Rustamov, R.M., Ovsjanikov, M., Azencot, O., BenChen, M., Chazal, F., Guibas, L.: Mapbased exploration of intrinsic shape differences and variability. ACM Transactions on Graphics (TOG) 32(4), 72 (2013)
 [27] Springborn, B., Schröder, P., Pinkall, U.: Conformal equivalence of triangle meshes. ACM Transactions on Graphics (TOG) 27(3), 77 (2008)
 [28] Vaxman, A., Müller, C., Weber, O.: Conformal mesh deformations with möbius transformations. ACM Transactions on Graphics (TOG) 34(4), 55 (2015)
 [29] Wang, Y., Wang, X., Zhang, B.: Affine differential invariants of functions on the plane. Journal of Applied Mathematics 2013 (2013)
 [30] White, J.H.: A global invariant of conformal mappings in space. Proceedings of the American Mathematical Society 38(1), 162–164 (1973)
 [31] Willmore, T.J.: Surfaces in conformal geometry. Annals of Global Analysis and Geometry 18(34), 255–264 (2000)

[32]
Xu, D., Li, H.: Geometric moment invariants. Pattern recognition
41(1), 240–249 (2008)  [33] Xu, J., Kang, H., Chen, F.: Contentaware image resizing using quasiconformal mapping. The Visual Computer 34(3), 431–442 (2018)

[34]
Yu, X., Lei, N., Wang, Y., Gu, X.: Intrinsic 3d dynamic surface tracking based on dynamic ricci flow and teichmüller map. Proceedings. IEEE International Conference on Computer Vision
2017, 5400–5408 (2017)  [35] Zeng, W., Gu, X.D.: Registration for 3d surfaces with large deformations using quasiconformal curvature flow. Computer Vision and Pattern Recognition (CVPR) (2011)
 [36] Zhang, H., Mo, H., Hao, Y., Li, S., Li, H.: Fast and efficient calculations of structural invariants of chirality. arXiv preprint arXiv:1711.05866 (2017)