参数化的扭曲度量

参数化(Parameterization)的扭曲(Distortion)常常使用雅可比矩阵(Jacobian Matrix)SVD分解得到的奇异值来量化。

扭曲(Distortion)在参数化中无处不在,但对扭曲进行度量是有必要的,它能帮助我们衡量不同参数化之间的利弊,也能够帮助我们得到一个满足特定最优条件的参数化。其中,最普遍的方法是使用雅可比矩阵(Jacobian Matrix)的奇异值来度量扭曲。在不少能量优化的方法中,构建一个以奇异值为自变量的能量函数是其最基本的一步。下面就二维曲面三角参数化来介绍上述度量的定义和计算。


微分几何定义

我们假设参数域上点\mathbf{u}(u,v)对应原曲面上一点\mathbf{p}=f(u,v)。我们在uv方向的移动一个小变量\Delta u\Delta v,那么对应曲面上新的点f(u+\Delta u, v + \Delta v)可以通过ff(u,v)附近的的一阶泰勒展开\tilde{f}来近似:

\tilde{f}(u+\Delta u, v + \Delta v)=f(u,v)+f_u(u,v)\Delta u + f_v(u,v)\Delta v

上述式子用雅可比矩阵形式即可写作:

\tilde{f}(u+\Delta u, v + \Delta v)=\mathbf{p}+J_f(\mathbf{u})\begin{pmatrix} \Delta u\\ \Delta v \end{pmatrix}\bigr

其中J_f = \begin{pmatrix} f_u & f_v \end{pmatrix}。对雅可比矩阵进行SVD分解,即可得到奇异值\sigma_1\sigma_2

J_f=U\Sigma V^T=U\begin{pmatrix} \sigma_1 & 0 \\ 0 & \sigma_2 \\ 0 & 0 \end{pmatrix}V^T


离散定义

一个比较容易理解和计算的雅可比矩阵的定义就是,以原网格和参数域上对应三角形之间2×2的线性变换矩阵作为其雅可比矩阵(对于不在同一平面的情形,可将其旋转+平移变换到相同平面上再作比较)。令原网格上的点为P(x,y,0),参数域上的点为\Omega(u,v),则三角形T_{abc}对应的雅可比矩阵J_{abc}满足:

J_{abc}\begin{bmatrix} u_a & u_b & u_c\\ v_a & v_b & v_c \end{bmatrix}=\begin{bmatrix} x_a & x_b & x_c\\ y_a & y_b & y_c \end{bmatrix}

之后,再对J_{abc}进行SVD分解,即可得到离散情形下对应的奇异值\sigma_1\sigma_2


奇异值的几何意义

要考察奇异值\sigma_1\sigma_2的几何意义,我们可以考察J_f=U\Sigma V^T是如何将参数域上的点\mathbf{u}及其领域变换到原曲面上点\mathbf{p}及其领域的:

由图可以非常直观地看出,\sigma_1\sigma_2相当于参数域某点\mathbf{u}及其领域在拉回到原曲面时,在两个正交方向上的拉伸比例。我们可以根据\sigma_1\sigma_2来定义一些比较常见的保测度的参数化:

并且也可以明显地看出,保长度=保角度+保面积。


常见能量定义

名称 定义
As-Similar-As-Possible \sum_{t \in T}^{ } (\rho_t(\sigma_1-\sigma_2)^2)
As-Rigid-As-Possible \sum_{t \in T}^{ } (\rho_t(\sigma_1-1)^2+(\sigma_2-1)^2)
Green-Lagrange \sum_{t \in T}^{ } (\rho_t(\sigma_1^2-1)^2+(\sigma_2^2-1)^2)
Symmetric Dirichlet \sum_{t \in T}^{ } (\sigma_1^2+\sigma_1^{-2}+\sigma_2^2+\sigma_2^{-2})

 

称谓(*)
邮箱
留言(*)