Anisotropic stretching

一般的曲面形变时的延伸更容易展现出非匀质的特性(anisotropic)， 比如一块布更容易向某一个方向拉伸。所以在讨论曲面形变的能量时，需要考虑其不同方向所产生的能量不同。

Xianfeng Gu等人的《A Generic Deformation Model for Dense Non-Rigid Surface Registration: a Higher-Order MRF-based Approach》讲的是用MRF模型来解决曲面配准的问题。其中谈到了局部形变的能量可以由两个主方向的拉伸来衡量。

-----------------------

需要一些黎曼几何的背景，但如果没有应该也能了解大致思想。用(M,g)来表示一个曲面M与其黎曼度量。根据拓扑流形，曲面是局部开曲面的集合，而每个开曲面与一个R2欧式空间同构，即可以有局部参数化表示. 曲面每个点p的度量张量(metric  tensor) 是一个2*2正定矩阵。这点p的局部肯定可以有不同的参数化表示，其对应的2*2正定矩阵也不用，但是可以由以下关系表示

其中alpha和beta表示不同的参数化。J表示Jacobian矩阵，也就是两个参数坐标系之间的线性变换。同理（此处paper中没提到），如果两曲面间有一个微分同胚映射f，那么对应点的参数坐标任然可以用以上关系表示。由于metric tensor正定，必然可以找到一坐标系使其变为单位矩阵，那么如果当对应点的metric tensor都为单位矩阵时，的两个特征值就表示了两个主方向的拉伸。

几何意义上，如图，这两个数描述了一个微小的局部里，一个圆是如何形变成一个椭圆的。

----------------------

我们可以想象形变的能量与这两个特征值有关，那么在实际应用时，需要对其求导。特征值的函数是不可微的，但是应该可以求次导数(subdifferential)

但单单用这个作为形变能量是不够的，我们知道一般来说thin shell的形变一般分为拉伸和弯曲（bending）。 后者应该可以通过类似的方法对shape operator(对应second fundamental form) 的变化进行描述，即两个主曲率的变化。