Matlab 基于Laplacian约束的收缩算法（骨架提取）

一、简介

骨骼提取一直是广泛研究的焦点，可以追溯到20世纪60年代。现有算法针对的模型特征和应用多种多样，提取过程大致可分为基于体素方法和几何方法。其中，几何方法可以直接作用于三角形网格或点云，就有学者直接利用点云中每个点的邻域信息，通过Laplacian算子迭代式的对点云进行收缩实现了对点云骨架的粗提取，感觉这个过程非常有意思，具体流程如下所述：

1、首先，根据之前的博客Matlab 基于局部三角网的曲面重建的做法，我们先生成一个初始的表面（它由各个局部的多边形组成）。
2、点云收缩，直至收缩到“零容积”。值得注意的是，收缩过程通过应用拉普拉斯平滑去除网格表面的细节和噪声，该平滑将点云中的点沿其近似曲率法线方向进行移动，这是一个迭代过程，会直至点云收缩到很小很细的状态。在真正了解其过程前，我们需要一些基础知识，如下迭代公式所示：