[论文解读] A convergent point integral method for isotropic elliptic equations on point cloud

Zhen Li, Zuoqiang Shi|arXiv (Cornell University)|Jun 11, 2015

Numerical methods in engineering参考文献 21被引用 4

一句话总结

本文提出了一种收敛的点积分方法，用于在点云上求解各向同性椭圆方程，通过核函数的积分公式近似微分算子。该方法保持了强制性并实现了收敛，从而在无需网格或几何先验信息的情况下，能够准确求解流形上的PDE，适用于图像处理和数据分析。

ABSTRACT

In this paper, we propose a numerical method to solve isotropic elliptic equations on point cloud by generalizing the point integral method. The idea of the point integral method is to approximate the differential operators by integral operators and discretize the corresponding integral equation on point cloud. The key step is to get the integral approximation. In this paper, with proper kernel function, we get an integral approximation for the elliptic operators with isotropic coefficients. Moreover, the integral approximation has been proved to keep the coercivity of the original elliptic operator. The convergence of the point integral method is also proved.

研究动机与目标

• 开发一种无需网格或额外几何信息的数值方法，用于在点云上求解各向同性椭圆方程。
• 将点积分方法从泊松方程推广到具有变系数的一般各向同性椭圆方程。
• 在光滑流形的点云采样下，建立该方法的理论收敛性。
• 在离散化形式中保持原始椭圆算子的强制性，以确保稳定性与收敛性。
• 在图像处理和数据分析中实现应用，例如基于流形模型的图像重建。

提出的方法

• 该方法通过基于非局部扩散理论推导出的核函数 $ R_t(\boldsymbol{x}, \boldsymbol{y}) $，利用积分公式近似椭圆微分算子。
• 基于逐点权重和核函数评估，采用基于求积规则的方法在点云上离散化积分方程。
• 核函数的选择确保对各向同性椭圆算子的积分近似，同时保持强制性。
• 采用自适应核宽 $ t(\boldsymbol{x}_i) = \rho(\boldsymbol{x}_i)^2 $，其中 $ \rho(\boldsymbol{x}_i) $ 为第20个最近邻的半径，以处理非均匀采样。
• 通过涉及差值 $ u(\boldsymbol{x}) - u(\boldsymbol{y}) $ 的积分公式，近似流形上的梯度与散度算子，确保与底层几何的一致性。
• 通过正则化避免系数矩阵退化，迭代求解离散系统。

实验结果

研究问题

• RQ1 点积分方法能否推广至具有变系数的各向同性椭圆方程在点云上的求解？
• RQ2 积分近似是否在离散设置下保持原始椭圆算子的强制性？
• RQ3 当点云趋于稠密时，所得离散格式是否收敛？
• RQ4 该方法能否应用于图像处理中的PDE，例如在图块流形上的总变差最小化？
• RQ5 从点积分方法导出的离散特征值问题的谱收敛行为如何？

主要发现

• 所提出的点积分方法在从光滑流形采样的点云上，对各向同性椭圆方程实现了收敛。
• 积分近似保持了原始椭圆算子的强制性，这对稳定性和收敛性至关重要。
• 该方法适用于带边界或无边界流形，且无需网格、水平集或最近点函数。
• 在图像重建的数值实验中，采用 $ L_1 $ 扩展方法取得了高质量结果，证明了该方法在实际应用中的有效性。
• 在均匀采样下，离散特征值问题的谱收敛性已建立，特征值与特征函数均实现收敛。
• 通过基于局部邻域密度的自适应核宽选择，该方法对非均匀采样具有鲁棒性。