conformal mapping的理解

先来复习下高数里面的第二类积分变换

以前对conformal mapping的理解都是盲人摸象, 片面局部的. 因为看的资料都是从不同的角度对解释它的. 随着阅历的丰富和理解的深入, 现已能将原来零星的知识点串连起来并找出其中的关系, 于是开篇博客讲讲自己对conformal mapping的学习历程以及深化后的理解, 好了, 开始正题.

以前看的资料呢对conformal mapping有不同的解释.

最始接触其实是在复变函数的共形映射, 可惜只停留在概念, 没有深入

对映射讲得比较多,比较透彻的是”As-Conformal-As-Possible Surface Registration”

然后就很奇怪为啥conformal是harmonic的子集

另外一篇论文是”As-Conformal-As-Possible Discrete Volumetric Mapping”

然后就想, cauchy-riemann条件是个什么鬼, 为什么满足了这个条件就能保角?

再来看看http://mathworld.wolfram.com/ConformalMapping.html

By letting w=f(z), the real and imaginary parts of w(z) must satisfy the Cauchy-Riemann equations and Laplace’s equation.

又多了一个laplace’s equation,它又是什么鬼?跟cauchy-Riemann equations有啥关系.
继续看http://mathworld.wolfram.com/Cauchy-RiemannEquations.html

然后, 我们画张图来形象的描述一下曲面映射吧

对于这张图有几点说明, 在$xOy$平面内向量[1;0]被映射到$[u_x; v_x]$, 向量[0;1]映射到$[u_y; v_y]$, 也就是矩阵$J$的两个列向量

现在回答上面提出的问题

cauchy-riemann条件是个什么鬼, 为什么满足了这个条件就能保角?
其实图中的jacobian矩阵$J$就是来测量面积变化的, 就是我们所讲的deformation gradient, 如何满足了cauchy-riemann条件, 那么 $u_x =v_y$, $u_y = -v_y$, 其实就代表了$J$可以表示为一个标量乘一个旋转矩阵$J = sR$, 注意旋转矩阵可以表示为

$$R_\theta = \begin{bmatrix} \cos\theta & -\sin\theta\\ \sin\theta & \cos\theta\end{bmatrix}$$

当$J = R$的时候, 就是rigid变换, 也就是上面说的isometric mapping

注意这是连续的情况, 离散的情况, 以mesh来讲, 可以像As-Conformal-As-Possible Discrete Volumetric Mapping一样可以离散到点, 也可以离散到每个三角形, 使他们的形变矩阵满足sR条件.

然后就很奇怪为啥conformal是harmonic的子集?
如果映射f满足cauchy-riemann方程, 对其求二阶导, 则可得出它满足laplace方程, 而满足laplace方程的映射称为harmonic mapping, 其实就是最小化Dirichlet energe $\int_\Omega |\nabla f|^2$, 关于Dirichlet energe 和laplace equation的关系, 可以参看: http://blog.csdn.net/seamanj/article/details/50898781.但是反之却不一定成立, 即一个映射满足laplace方程,却不一定满足cauchy-riemann方程.

最后看下”computing Discrete minimal surfaces and their conjugates”里面用的conformal mapping

这里我们需要证明 $jacobian(f) \le \frac{1}{2}\nabla f^2$

$$det(J) = \begin{vmatrix} u_x & u_y\\ v_x & v_y\end{vmatrix} = u_xv_y - u_yv_x$$
$$\frac{1}{2}\nabla f^2 = \frac{1}{2}trace({\partial f}^T\partial f)= \frac{1}{2}\begin{vmatrix} u_x & u_y\\ v_x & v_y\end{vmatrix}:\begin{vmatrix} u_x & u_y\\ v_x & v_y\end{vmatrix}= \frac{1}{2}(u_x^2+u_y^2+v_x^2+v_y^2)$$
$$\frac{1}{2}\nabla f^2 - det(J) = \frac{1}{2}((u_x - v_y)^2+(u_y+v_x)^2)$$

可知, 当其为0时,只有f满足cauchy-riemann方程才成立, 即映射f为conformal mapping.