忠实的从地震震源重建构造俯冲带的 REEB 图-计算机科学-机器学习-拓扑数据分析-图形重构技术.pdf

忠实的从地震震源重建构造俯冲带的REEB图

HALLEYFRITZE,SUSHOVANMAJHI,MARISSAMASDEN,ATISHMITRA,ANDMICHAELSTICKNEY

摘要.拓扑数据分析（TDA）中的一个重要问题—既有理论意义也有实际应用价值—是从可能带有噪

声的数据中重建底层的（通常是未知的）度量图的拓扑和几何结构。Reeb图最近成功地用于在Gromov–Hausdorff

噪声下抽象度量图的重建：假设样本与真实值在度量上接近。然而，这种强有力的全局密度保证通常

不可用，使得现有的基于Reeb图的方法无法使用。另一个截然不同但更相关的范式关注的是从仅在

Hausdorff距离内接近的真实样本中重构嵌入到欧几里得空间中的度量图。我们放宽了对密度的假设，

给出了可证明的几何重建方案，即使样本只在局部上是度量上的接近，但仍为在Hausdorff噪声模型

下成功地进行欧几里得图的几何重建提供了可证明的保证。

我们将图形重构技术应用于从全球地震目录中重建地震板块构造边界。SLAB2.0模型是对地球上

本所有已知俯冲带板块的综合空间总结。我们从可能含有噪声的地震震源数据中重构SLAB2.0模型的部

分内容。

译

中

2

v

0

1

4

9

1

.

0

1

4

2

:

v

i

x

r

a

Keywordsandphrases.里布图，地壳俯冲带，几何重建，Vietoris–Rips复形.

1

2

1.介绍

现代技术使得大规模几何数据的获取变得容易快捷。在许多实际应用中，数据来自一个度量图

（在生物科学和天文科学中被称为丝状结构）。虽然其中一些图形是摘要—例如，传感器网络—其他

常见应用涉及嵌入某些欧几里得空间中的图表样本；我们称这种图为嵌入式。示例包括围绕道

路网络或地图[1]的GPS轨迹、沿断层和板块构造边界[2]的地震位置等。数据中的噪声使重建问

题成为拓扑和几何数据分析中一个计算上具有挑战性的问题。

非正式地说，一个度量图是一个组合图，其顶点集为，边集为，并带有一

个长度度量使得在拓扑上等价于一个维单纯复形；更正式的定义见[3]。一个围绕的样

本通常被建模为一个有限度量空间。然后，在适当的比例尺上计算一个邻域图作为重建

算法的输入。我们选择在原始数据上使用Vietoris–Rips图，比例尺为，如第3节中所

定义。

度量图的重构——无论是抽象图还是嵌入图——在文献中已有深入研究；参见例如[4,5,6,7,8,

9,10]。尽管在Gromov-Hausdorff框架下对抽象图进行了大量重建工作，但从仅与真实值Hausdorff

接近的欧几里得样本中成功地拓扑和几何重构嵌入图的发展并不多。对于嵌入图，我们特别提到

Majhi[9]使用Vietoris-Rips复形的技术。所得抽象单纯复形仅仅保持了拓扑保真度，未能生成在相

同空间中的几何图如输出所示。本文主要目标是缓解这种缺陷，开发从噪声样本中基于Hausdorff

距离重构嵌入图的可验证方案。

动机与贡献.我们受到最近Reeb图在几何重建抽象度量图方面关键成功[11,5,7]的启发。这项工

作的动机是双重的。(1)首先，在抽象设置中，我们将过于严格且不切实际的Gromov–Hausdorff假

设放松，以提供定理4.2中的重建保证，从局部接近的-近似样本或其邻域图（定义A.3）。(2)

其次，我们考虑嵌入度量图的重构问题，使得欧几里得距离诱导出内在距离在上，

并且原始数据仅在中的豪斯多夫距离接近于，使用图的受限失真（定义3.1）作为嵌

入图的新采样参数——参见定理4.3。除了我们的理论结果，我们还展示了从传感器数据中重建地

震板块构造边界的技术。