Apollonius diagram 简介

Apollonius diagram (or the Additively weighted Voronoi diagram) 其距离函数定义为点到球表面的最短距离(可以是负值,若点在球内):

δ(x,Pi)=xciwi\delta(x, P_i)=||x-c_i||-w_i
阅读更多

Delaunay 三角剖分笔记

delaunay三角剖分 是一种根据点集获得满足一定良好性质的三角剖分方式,等价于 上举映射(lifting map) 后,其 下凸壳(Lower Convex Hull) 在原平面的投影。如果点集中不存在四点共圆的情况(退化情况),Delaunay 三角剖分是唯一的。

阅读更多

Alpha Shape 与复形

抽象单纯形 σV\sigma \subset V ,是有限点集,其维度定义为 dimσ=σ1\mathrm{dim}\,\sigma =|\sigma |-1 。(这里的 |\cdot| 和下面的不表示同一个意思)

几何单纯形 σ|\sigma| 是抽象单纯形的实现,是仿射无关的一组顶点的凸包:

σ:=conv{xv:vσ}Rd|\sigma|:=\mathrm{conv}\{x_v:v\in\sigma\}\subset \mathbb{R}^d

凸包定义为:

conv{xi}i=1n:={Σiλixiλi0,Σiλi=1}\mathrm{conv}\{x_i\}_{i=1}^n:=\{\Sigma_i\lambda_ix_i|\lambda_i\ge0,\Sigma_i\lambda_i=1\}

单纯性复形 KK 为抽象单纯形的集合,要求 σK,τσ    τK\forall\sigma\in K,\forall \tau\subset\sigma \implies \tau\in K ,比如

{\{ {1},{2},{3},{4},{1,2},{2,3},{1,3},{1,2,3}\{1\},\{2\},\{3\},\{4\},\{1,2\},\{2,3\},\{1,3\},\{1,2,3\} }\}

单纯形复形的实现 K|K| ,就是把 KK 中的单纯形实现了。

给定点集 PRdP\subset \mathbb{R}^d ,

Delaunay complex D(P)D(P) ,是指按Delaunay三角剖分得到的实现对应的抽象单纯形集合。

Alpha complex Aα(P)A_\alpha(P) ,是由α-空球条件选出的 D(P)D(P) 的子复形,带有自然的几何实现 Alpha Shape Sα(P)=Aα(P)S_\alpha(P) = |A_\alpha(P)|

Rd\Aα\mathbb{R}^d\backslash|A_\alpha| 分成: 一个 无界连通分支 → 外部;若干 有界连通分支 → 内部空腔(voids)。在一点紧致化下无限连通分支连通到 \infty 点。在同胚意义下, D(P{})\Aα(P)(外部)(所有空腔)D(P\cup\{\infty\})\backslash A_\alpha(P) \cong (外部)\cup(所有空腔) (因此在工程实现中,可以在包围盒上取一点作为无穷远点)。

Alpha shape Sα(P)S_\alpha(P) 同胚于球并集 Uα=iB(pi,α)U_\alpha =\bigcup_i B(p_i, \alpha)

半径不同的球 → Weighted Alpha Complex / Power Diagram。构造方式:空球条件改为加权空球条件,Delaunay → Power complex

三球交点的简单求法

三球交点实际上是在求幂轴在任意球上的交点,幂轴用一个不定方程组 Ax=bAx=b 表示(球方程两两相减),而幂点是该轴上距离球心最近的点,因此以某一球心为参考点 OO ,用拉格朗日乘子法( L=OX2+λT(Axb)\mathcal{L}=|OX|^2+\lambda^T(Ax-b) )得到以下方程,求解就可以直接得到幂点 xx ,知道幂点求交点就显然了:

x=AT(AAT)1b其中,A=[OO2OO3],b=12[R12R22+OO22R12R32+OO32]\begin{aligned} x&=A^T(AA^T)^{-1}b \\ \text{其中,}A&=\begin{bmatrix} \overrightarrow{OO_2} \\ \overrightarrow{OO_3} \end{bmatrix} ,\\ b&=\frac{1}{2} \begin{bmatrix} R_1^2-R_2^2+|OO_2|^2 \\ R_1^2-R_3^2+|OO_3|^2 \end{bmatrix} \end{aligned}
Your browser is out-of-date!

Update your browser to view this website correctly.&npsb;Update my browser now

×