# Differential and integral invariants under Mobius transformation

One of the most challenging problems in the domain of 2-D image or 3-D shape is to handle the non-rigid deformation. From the perspective of transformation groups, the conformal transformation is a key part of the diffeomorphism. According to the Liouville Theorem, an important part of the conformal transformation is the Mobius transformation, so we focus on Mobius transformation and propose two differential expressions that are invariable under 2-D and 3-D Mobius transformation respectively. Next, we analyze the absoluteness and relativity of invariance on them and their components. After that, we propose integral invariants under Mobius transformation based on the two differential expressions. Finally, we propose a conjecture about the structure of differential invariants under conformal transformation according to our observation on the composition of the above two differential invariants.

## Authors

• 41 publications
• 12 publications
• 8 publications
• 67 publications
• 31 publications
• ### Image Projective Invariants

In this paper, we propose relative projective differential invariants (R...
07/19/2017 ∙ by Erbo Li, et al. ∙ 0

• ### Image Differential Invariants

Inspired by the methods of systematic derivation of image moment invaria...
11/13/2019 ∙ by Hanlin Mo, et al. ∙ 22

• ### Differential Invariants

Validation is a major challenge in differentiable programming. The state...
01/09/2021 ∙ by Uwe Naumann, et al. ∙ 0

• ### Isomorphism between Differential and Moment Invariants under Affine Transform

The invariant is one of central topics in science, technology and engine...
05/20/2017 ∙ by Erbo Li, et al. ∙ 0

• ### Finding models through graph saturation

We give a procedure that can be used to automatically satisfy invariants...
06/25/2018 ∙ by Sebastiaan J. C. Joosten, et al. ∙ 0

• ### A Root-Free Splitting-Lemma for Systems of Linear Differential Equations

We consider the formal reduction of a system of linear differential equa...
11/02/2019 ∙ by Eckhard Pflügel, et al. ∙ 0

• ### Linear Time Recognition Algorithms for Topological Invariants in 3D

In this paper, we design linear time algorithms to recognize and determi...
04/12/2008 ∙ by Li Chen, et al. ∙ 0

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Introduction

One of the most challenging problems in the domain of 2-D image or 3-D shape is to handle the non-rigid 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 non-rigid 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 angle-preserving 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], Cauchy-Riemann 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 Gauss-Bonnet Theorem, White[30] proposed that is a global conformal invariant if is an oriented and closed surface. The Gauss-Bonnet 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 2-D and 3-D Möbius transformation respectively. According to the Liouville Theorem, the 3-D 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 1-quasiconformal if and only if or its complex conjucate is a meromorphic function of a complex variable in . If , then is 1-quasiconformal 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

 T(z)=az+bcz+d, (1)

where , . Based on the Liouville Theorem[20], every Möbius transformation in higher dimensions could be given with the form

 T(x)=b+γA(x−a)∥x−a∥ϵ2, (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 n-dimension Möbius transformation is a homemorphism of (the one-point 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 2-D 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.

 u=ax+r2(x−ax)(x−ax)2+(y−ay)2 (3)
 v=ay+r2(y−ay)(x−ax)2+(y−ay)2 (4)

At the same time, it means that the coordinates transformations under are as follows.

 x=ax+r2(u−ax)(u−ax)2+(v−ay)2 (5)
 y=ay+r2(v−ay)(u−ax)2+(v−ay)2 (6)

Based on and the equations (5)(6), we obtain the relationships between the partial derivatives of and as follows.

 gu = fxxu+fyyu (7) gv = fxxv+fyyv (8) guu = (fxxxu+fxyyu)xu+fxxuu+(fyxxu+fyyyu)yu+fyyuu (9) guv = (fxxxv+fxyyv)xu+fxxuv+(fyxxv+fyyyv)yu+fyyuv (10) gvv = (fxxxv+fxyyv)xv+fxxvv+(fyxxv+fyyyv)yv+fyyvv (11)

Then we obtain a 2-D equation under the inversion transformation, it is

 guu+gvvg2u+g2v=fxx+fyyf2x+f2y (12)

This means that

 fxx+fyyf2x+f2y (13)

is a differential invariant under inversion transformation. We use the same method in 3-D situation and obtain a differential invariant under the inversion transformation, it is

 fA+fB(f2x+f2y+f2z)2 (14)

where

 fA=(fxx+fyy+fzz)(f2x+f2y+f2z) (15) fB=f2xfxx+f2yfyy+f2zfzz+2fxfxyfy+2fxfxzfz+2fyfyzfz

### 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 2-D 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

 Inv′T=WT⋅InvT (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 2-D inversion invariants, we obtain that for the numerator and denominator of (13), this means

 guu+gvv=||J||−1(fxx+fyy) (17)
 g2u+g2v=||J||−1(f2x+f2y) (18)

where is the determinant of Jacobian matrix of transformation , is the absolute valve of . In 3-D situation, we obtain for the numerator and denominator of (14). In the stretching transformation, we obtain in 2-D situation, and in 3-D situation. We also obtain that for the numerator and denominator of (13) or (14) under translation, rotation and reflection transformations.

The result of absoluteness and relativity of invariance on (13) and (14) is shown in Table 1.

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

 ∬D′(guu+gvv)dudv=∬DWT(fxx+fyy)||JT||dxdy=∬D(fxx+fyy)dxdy (19)
 ∬D′(g2u+g2v)dudv=∬DWT(f2x+f2y)||JT||dxdy=∬D(f2x+f2y)dxdy (20)

where is the area extension factor, so we obtain that

 ∬D(fxx+fyy)dxdy (21)
 ∬D(f2x+f2y)dxdy (22)

are integral invariants under 2-D Möbius transformation. In the same way, we obtain that

 ∭D(f2x+f2y+f2z)32dxdydz (23)
 ∭D(fA+fB)34dxdydz (24)

are integral invariants under 3-D 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 specify-scale and global invariant is the definition domain, the construction method of specify-scale 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:

 ∬D(fxx+fyy)n+1(f2x+f2y)ndxdy (25)
 ∬D(f2x+f2y)n+1(fxx+fyy)ndxdy (26)
 ∭D(fA+fB)34(n+1)(f2x+f2y+f2z)32ndxdydz (27)
 ∭D(f2x+f2y+f2z)32(n+1)(fA+fB)34ndxdydz (28)

if the denominators of (25), (26), (27), (28) are not zero.

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

 f2x+f2yorf2x+f2y+f2z (29)
 fxx+fyyorfxx+fyy+fzz (30)
 f2xfxx+f2yfyy+f2zfzz+2fxfxyfy+2fxfxzfz+2fyfyzfz (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 self-consistent manner.

One of the possible self-consistent forms in n-dimensional Euclidean space may be

 n−1∑i=1∏aij=1DRIj(f2x1+f2x2+⋯+f2xn)n−1 (32)

where is differential invariant under rigid transformation.

## 5 Experimental Results

We choose a human face model from TOSCA database and treat the z-coordinate value of vertexes of the triangle mesh as a function defined on x-coordinate and y-coordinate, i.e.

. With least square method, the coordinates of a vertice and its 1-ring 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 1-ring 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 5-dimension 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

 Err=1N∑i|InvT;i−InvO;i||InvT;i|+|InvO;i|×100% (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 5-dimension 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 matching-fail 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 2-D and 3-D Möbius transformation respectively, in particular, the 3-D 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., Ben-Chen, 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 N-dimensional Euclidean Space with a Non-definite 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., Ben-Chen, 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 n-space. Proceedings of the International Workshop on Quasiconformal (2006)
• [26] Rustamov, R.M., Ovsjanikov, M., Azencot, O., Ben-Chen, M., Chazal, F., Guibas, L.: Map-based 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(3-4), 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.: Content-aware image resizing using quasi-conformal 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 quasi-conformal 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)