[论文解读] A class of second-order geometric quasilinear hyperbolic PDEs and their application in imaging science

Guozhi Dong, Michael Hintermüller+1 more|arXiv (Cornell University)|Jan 1, 2019

Numerical methods in inverse problems参考文献 31被引用 1

一句话总结

本文提出了一类用于图像处理的二阶几何拟线性双曲型偏微分方程，重点研究阻尼二阶全变差流与平均曲率流。该文建立了解的存在性与唯一性，为简单初值提供了解析解，并通过数值比较展示了其在去噪与去抖动应用中相较于一阶流的更优收敛性与边缘保持性能。

ABSTRACT

In this paper, we study damped second-order dynamics, which are quasilinear hyperbolic partial differential equations (PDEs). This is inspired by the recent development of second-order damping systems for accelerating energy decay of gradient flows. We concentrate on two equations: one is a damped second-order total variation flow, which is primarily motivated by the application of image denoising; the other is a damped second-order mean curvature flow for level sets of scalar functions, which is related to a non-convex variational model capable of correcting displacement errors in image data (e.g. dejittering). For the former equation, we prove the existence and uniqueness of the solution. For the latter, we draw a connection between the equation and some second-order geometric PDEs evolving the hypersurfaces which are described by level sets of scalar functions, and show the existence and uniqueness of the solution for a regularized version of the equation. The latter is used in our algorithmic development. A general algorithm for numerical discretization of the two nonlinear PDEs is proposed and analyzed. Its efficiency is demonstrated by various numerical examples, where simulations on the behavior of solutions of the new equations and comparisons with first-order flows are also documented.

研究动机与目标

• 开发用于改进图像恢复的二阶几何拟线性双曲型PDE，动机源于优化中二阶动力学的优越数值表现。
• 将阻尼二阶流扩展至全变差流与平均曲率流，以实现图像去噪与位移误差校正（例如去抖动）。
• 为这些新型PDE建立解析与数值基础，包括解的存在性、唯一性及收敛性。
• 通过数值实验表明，二阶流在保持边缘与减少振荡方面优于一阶梯度流。

提出的方法

• 提出一种阻尼二阶全变差流（TVF）作为双曲型PDE：$\ddot{w} + \eta(t)\dot{w} = -\partial\Phi(w)$，其中$\Phi(w)$为全变差泛函。
• 为水平集引入阻尼二阶平均曲率流，用于建模超曲面演化以最小化周长并校正位移误差。
• 采用辛Euler格式进行时间半离散化，确保能量稳定与数值收敛。
• 推导出矩阵形式的半离散系统$\mathbf{z}_{k+1} = \mathbf{A}_k \mathbf{z}_k$，其中$\mathbf{A}_k$表示时间推进算子。
• 利用Gershgorin圆定理与特征值分析，证明在CFL型条件下，离散格式具有压缩性。
• 实现平均曲率流PDE的正则化版本，以确保适定性与数值稳定性。

实验结果

研究问题

• RQ1 能否为去噪与去抖动等图像恢复任务构建并分析二阶几何PDE？
• RQ2 这些新型二阶双曲型PDE的解是否存在且唯一？
• RQ3 与一阶梯度流相比，这些二阶流的解在收敛速度与边缘保持方面表现如何？
• RQ4 在二阶TVF框架下，能否为简单初值推导出解析解？
• RQ5 时间步长需满足何种条件，才能保证半离散数值格式的稳定性与收敛性？

主要发现

• 阻尼二阶全变差流存在唯一解，且对简单初值可显式导出解析解。
• 对正则化后的二阶平均曲率流，严格证明了解的存在性与唯一性。
• 基于辛Euler方法的数值格式在涉及阻尼参数$\eta$与系统矩阵最大特征值的时间步长条件下，被证明具有收敛性。
• 证明了离散迭代矩阵$\mathbf{A}_k$的特征值满足$|\mu_{i,\pm}^k| \leq 1$，确保压缩性与收敛性。
• 数值比较表明，二阶流收敛速度更快，且在边缘丰富区域更有效地保持边缘特征。
• 该方法通过平均曲率流演化水平集，成功校正图像中的位移误差（如去抖动），恢复原始图像结构。