论文阅读：Mean Value Coordinates

均值坐标–Mean Value Coordinates(MVC) 是一种重要的插值技术，广泛应用于参数化、网格变形等几何处理方向。本文主要讲解Floater等人2003年发表在CAGD上的论文Mean value coordinates，主要探讨2D多边形中均值坐标的计算方法。

1. 重心坐标

对于三角形$△v_1{v}_2{v}_3$ 中任一点 $p$ 表示为三角形三个顶点的凸组合, 即满足如下的等式：

$$p={\lambda }_1{v}_1+{\lambda }_2{v}_2+{\lambda }_3{v}_{3}and{\lambda }_1+{\lambda }_2+{\lambda }_3=1\text{\hspace{1em}(1)}$$

对于三角形来讲，如下图所示，权重很好计算（三角形面积之比）：
$$\{\begin{array}{r}{\lambda }_1=S△p{v}_2{v}_3/S△v_1{v}_2{v}_3\\ {\lambda }_2=S△p{v}_1{v}_3/S△v_1{v}_2{v}_3\\ {\lambda }_3=S△p{v}_1{v}_2/S△v_1{v}_2{v}_3\end{array}\text{\hspace{1em}(2)}$$

为什么叫重心坐标呢？ 以下内容摘自知乎：重心坐标（Barycentric coordinates）

等式(1)可以转化为：

$${\lambda }_1(v_1-p)+{\lambda }_2(v_2-p)+{\lambda }_3(v_3-p)=0$$
进一步可以改写为：

$${\lambda }_1\overrightarrow{p{v}_1}+{\lambda }_2\overrightarrow{p{v}_2}+{\lambda }_3\overrightarrow{p{v}_3}=0$$

相当于我们在三角形的三个顶点 $v_1$处挂重量为${\lambda }_1$, $v_2$处挂重量为${\lambda }_2$, $v_3$处挂重量为${\lambda }_3$的重物时，如果把 $p$ 点放避雷针上，能保持平衡，则这种情况下 $p$ 是重心。

2. 广义重心坐标

如上所述，对于三角形的重心坐标是比较容易计算的，但我们不满足于此。我们希望求解多边形内部顶点的广义重心坐标， 即将多边形内部的顶点表示为多边形顶点的凸组合：

$$p=\sum _{i=1}^n{\lambda }_i{v}_{i}and\sum _{i=1}^n{\lambda }_i=1\text{\hspace{1em}(3)}$$

广义重心坐标有很多求法比如：均值坐标(MVC)，格林坐标(GC)，调和坐标(HC)等，这里我们介绍一下广泛使用的均值坐标（Mean Value coordinates）的计算方法。

3. 2D 均值坐标

对于一般的2D的多边形，Floater (Floater M S . Mean value coordinates[M]. Elsevier Science Publishers B. V. 2003.) 给出的均值坐标为：

$$p=\sum _{i=1}^n{\lambda }_i{v}_{i}and{\lambda }_i=w_i/\sum _{j=1}^n{w}_j\text{\hspace{1em}(4)}$$
$$w_i=[tan({\alpha }_{i-1}/2)+tan({\alpha }_i/2)]/||v_i-p||\text{\hspace{1em}(5)}$$

可以看出这种形式表示的坐标始终为正，下面我把这个证明过程简单搬运一下，

我们的目标是将多边形内部的点 $p$ 表示为多边形顶点的凸组合：
$$p=\sum _{i=1}^n{\lambda }_i{v}_{i}and{\lambda }_i=w_i/\sum _{j=1}^n{w}_j$$

等价于： ${\sum }_{i=1}^n(w_i/{\sum }_{j=1}^n{w}_j)(v_i-p)=0$, 左右两边分别乘以 ${\sum }_{j=1}^n{w}_j$:

$$\sum _{i=1}^n{w}_i(v_i-p)=0\text{\hspace{1em}(6)}$$

接下来，使用 $p$ 处的极坐标表示 $v_i$：
$$v_i=p+||v_i-p||(cos{\theta }_i,sin{\theta }_i)$$

那么有：
$$(v_i-p)/||v_i-p||=(cos{\theta }_i,sin{\theta }_i)\text{\hspace{1em}(7)}$$

将（5）代入（6）：

$$\sum _{i=1}^n[tan({\alpha }_{i-1}/2)+tan({\alpha }_i/2)]\cdot (v_i-p)/||v_i-p||=0\text{\hspace{1em}(8)}$$

并将（7）代入（8）：

$$\sum _{i=1}^n[tan({\alpha }_{i-1}/2)+tan({\alpha }_i/2)](cos{\theta }_i,sin{\theta }_i)=0\text{\hspace{1em}(9)}$$

上式进一步等效于（注意类似于：${\sum }_{i=2}^{n=1}$表示2,3,4,…n, 0）：

$$\sum _{i=2}^{n=1}tan({\alpha }_i/2)[(cos{\theta }_i,sin{\theta }_i)+(cos{\theta }_{i+1},sin{\theta }_{i+1})]=0\text{\hspace{1em}(10)}$$

如果能证明等式（10）成立，则多边形内部的顶点就可以用等式（5）所示的坐标表示（插值）。

首先，我们知道，从圆心指向圆周的所有单位向量之和为 0 向量，因为任意一个向量都能找出一个与之大小相等方向相反的向量。

$$0={\int }_0^{2\pi }(cos\theta ,sin\theta )=\sum _{i=1}^k{\int }_{{\theta }_i}^{{\theta }_{i+1}}(cos\theta ,sin\theta )=\sum _{i=1}^k{\int }_{{\theta }_i}^{{\theta }_{i+1}}\frac{sin({\theta }_{i+1}-\theta )}{sin{\alpha }_i}(cos{\theta }_i,sin{\theta }_i)+\frac{sin(\theta -{\theta }_i)}{sin{\alpha }_i}(cos{\theta }_{i+1},sin{\theta }_{i+1})$$

这里的等式构造非常巧妙。对上式积分等于：

$$\sum _{i=1}^k\frac{(cos{\theta }_i,sin{\theta }_i)}{sin{\alpha }_i}{\int }_{{\theta }_i}^{{\theta }_{i+1}}sin({\theta }_{i+1}-\theta )d\theta +\frac{(cos{\theta }_{i+1},sin{\theta }_{i+1})}{sin{\alpha }_i}{\int }_{{\theta }_i}^{{\theta }_{i+1}}sin(\theta -{\theta }_{i+1})d\theta \text{\hspace{1em}(11)}$$

而 ${\int }_{{\theta }_i}^{{\theta }_{i+1}}sin({\theta }_{i+1}-\theta )d\theta ={\int }_{{\theta }_i}^{{\theta }_{i+1}}sin(\theta -{\theta }_{i+1})d\theta =1-cos{\alpha }_i$， 因此（11）计算结果为：

$$\sum _{i=1}^k\frac{(1-cos{\alpha }_i)(cos{\theta }_i,sin{\theta }_i)}{sin{\alpha }_i}+\frac{(1-cos{\alpha }_i)(cos{\theta }_{i+1},sin{\theta }_{i+1})}{sin{\alpha }_i}\text{\hspace{1em}(12)}$$
因为： $\frac{1-cos{\alpha }_i}{sin{\alpha }_i}=tan({\alpha }_i/2)$， 所以（12）等价于：

$$\sum _{i=1}^{k}tan({\alpha }_i/2)[(cos{\theta }_i,sin{\theta }_i)+(cos{\theta }_{i+1},sin{\theta }_{i+1})]\text{\hspace{1em}(13)}$$

因此，等式（10）得证。