Level Set Methods and Dynamic Implicit Surfaces
Contents
Preface vii
Color Insert (facing page 146)
I Implicit Surfaces 1
1 Implicit Functions 3
1.1 Points . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4 Geometry Toolbox . . . . . . . . . . . . . . . . . . . . . 8
1.5 Calculus Toolbox . . . . . . . . . . . . . . . . . . . . . . 13
2 Signed Distance Functions 17
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.2 Distance Functions . . . . . . . . . . . . . . . . . . . . . 17
2.3 Signed Distance Functions . . . . . . . . . . . . . . . . . 18
2.4 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.5 Geometry and Calculus Toolboxes . . . . . . . . . . . . . 21
II Level Set Methods 23
3 Motion in an Externally Generated Velocity Field 25
3.1 Convection . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.2 Upwind Dierencing . . . . . . . . . . . . . . . . . . . . 29
3.3 Hamilton-Jacobi ENO . . . . . . . . . . . . . . . . . . . 31
3.4 Hamilton-Jacobi WENO . . . . . . . . . . . . . . . . . . 33
3.5 TVD Runge-Kutta . . . . . . . . . . . . . . . . . . . . . 37
4 Motion Involving Mean Curvature 41
4.1 Equation of Motion . . . . . . . . . . . . . . . . . . . . . 41
4.2 Numerical Discretization . . . . . . . . . . . . . . . . . . 44
4.3 Convection-Diusion Equations . . . . . . . . . . . . . . 45
5 Hamilton-Jacobi Equations 47
5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.2 Connection with Conservation Laws . . . . . . . . . . . . 48
5.3 Numerical Discretization . . . . . . . . . . . . . . . . . . 49
5.3.1 Lax-Friedrichs Schemes . . . . . . . . . . . . . . . 50
5.3.2 The Roe-Fix Scheme . . . . . . . . . . . . . . . . 52
5.3.3 Godunovs Scheme . . . . . . . . . . . . . . . . . 54
6 Motion in the Normal Direction 55
6.1 The Basic Equation . . . . . . . . . . . . . . . . . . . . . 55
6.2 Numerical Discretization . . . . . . . . . . . . . . . . . . 57
6.3 Adding a Curvature-Dependent Term . . . . . . . . . . . 59
6.4 Adding an External Velocity Field . . . . . . . . . . . . . 59
7 Constructing Signed Distance Functions 63
7.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 63
7.2 Reinitialization . . . . . . . . . . . . . . . . . . . . . . . 64
7.3 Crossing Times . . . . . . . . . . . . . . . . . . . . . . . 65
7.4 The Reinitialization Equation . . . . . . . . . . . . . . . 65
7.5 The Fast Marching Method . . . . . . . . . . . . . . . . 69
8 Extrapolation in the Normal Direction 75
8.1 One-Way Extrapolation . . . . . . . . . . . . . . . . . . . 75
8.2 Two-Way Extrapolation . . . . . . . . . . . . . . . . . . 76
8.3 Fast Marching Method . . . . . . . . . . . . . . . . . . . 76
9 Particle Level Set Method 79
9.1 Eulerian Versus Lagrangian Representations . . . . . . . 79
9.2 Using Particles to Preserve Characteristics . . . . . . . . 82
10 Codimension-Two Objects 87
10.1 Intersecting Two Level Set Functions . . . . . . . . . . . 87
10.2 Modeling Curves in 3 . . . . . . . . . . . . . . . . . . . 87
10.3 Open Curves and Surfaces . . . . . . . . . . . . . . . . . 90
10.4 Geometric Optics in a Phase-Space-Based Level
Set Framework . . . . . . . . . . . . . . . . . . . . . . . . 90
### 水平集方法与动态隐式表面
水平集方法(Level Set Methods)是一种用于追踪随时间变化的界面(如自由边界问题中的相界面)的数值技术。它在多个科学和工程领域中都有广泛的应用,包括图像处理、计算机视觉、流体动力学等。本文将根据提供的文件摘要来详细阐述其核心概念和技术。
#### 隐式表面
- **隐式函数**:隐式函数定义了空间中的几何结构,如点、曲线或表面。
- **点**:可以使用一个常数值作为隐式函数表示一个点。
- **曲线**:通过一维隐式函数定义,该函数的零值等高线表示曲线。
- **表面**:二维隐式函数用于表示三维空间中的表面,其零等值面即为所求表面。
- **距离函数**:
- **有符号距离函数**:一种特殊类型的隐式函数,不仅定义了零等值面,还给出了每个点到这个面的距离及方向。正值代表点位于表面外侧,负值则表示内侧。
#### 水平集方法概述
水平集方法的核心是利用水平集方程追踪界面的演化。它能够处理复杂的拓扑变化,并且可以高效地进行数值模拟。
- **运动方程**:描述界面随时间演化的偏微分方程。
- **外部速度场中的运动**:考虑外界因素对界面移动的影响。
- **涉及平均曲率的运动**:当界面的运动受到曲率控制时的情况。
- **数值离散化**:为了在计算机上模拟这些方程,需要将其离散化。
- **向上差分法**:用于处理非线性偏微分方程的一种常用技术。
- **Hamilton-Jacobi 方程**:一类非线性偏微分方程,广泛应用于控制理论和优化问题。
- **Hamilton-Jacobi ENO/WENO 方法**:基于有限差分的高精度数值方案。
- **TVD Runge-Kutta 方法**:确保数值解满足特定稳定性的高阶时间积分方法。
#### 特殊主题
- **构建有符号距离函数**:重新初始化水平集函数,使其保持为有符号距离函数。
- **快速行进方法**:一种高效算法,用于计算最近点问题。
- **正常方向的外推**:扩展水平集方法以处理更复杂的界面演化。
- **粒子水平集方法**:结合粒子方法和水平集方法的优点,提高追踪复杂界面的能力。
- **二维对象的建模**:使用两个水平集函数的交集来模拟二维空间中的曲线。
- **开放曲面和表面**:处理非封闭表面的问题,例如裂缝或孔洞。
#### 实际应用
- **图像分割**:利用水平集方法自动识别和分割图像中的物体。
- **流体动力学模拟**:追踪不混溶流体之间的界面,如水和空气。
- **医学成像**:用于分析和重建医学图像中的结构。
#### 结论
水平集方法为追踪和模拟随时间变化的界面提供了一个强大而灵活的工具。通过对隐式函数的巧妙运用,它可以有效地处理各种复杂场景下的界面演化问题。随着数值技术和计算能力的不断提高,水平集方法将在更多领域展现出其价值。
1