描述：

MSFM2D/MSFM3D函数使用多模板快速行进方法(MSFM)计算从一个点列表到2D或3D图像中所有其他像素的最短距离。该方法利用二阶导数和交叉邻域给出了更精确的距离。

-功能骨架将使用快速行进距离变换计算由二值图像/体积表示的对象的精确骨架（中心线）。

-函数Shortestpath使用二维或三维距离图中的Euler或Runge Kutta 4追踪从起点到源点的最短路径。

实施：
二维快速行进法用Matlab和c代码实现。c代码使用自定义的构建未排序二叉树最小搜索，它比普通的二叉排序树执行得更快。c代码的速度是matlab代码（使用Microsoft Visual编译器编译）的500多倍。

文献：
我们使用了两篇论文：
-J.Andreas Baerentzen“关于3D晶格快速推进方法的实现”
-M.Sabry Hassouna等人，“多阶快速行进法：笛卡尔域上Eikonal方程的高精度解”
-R.van Uitert等人，“基于快速行进法的体数据亚体素精确骨架”

我们将实施结果与论文中的结果进行了比较：
-正常快速行军第1步，完全相同的结果。
-2阶，显著小于论文中的误差。
-Multistencil 1阶，比论文中的误差更大
- Multistencil排名第2，明显比在论文中发表的结果差。(注意:我们的结果与其他现有实现一致)

我们的代码的最后一个版本产生了比论文或其他文献中更好的结果。这是通过使用所有可用信息求解多项式根来实现的，如下面Olivier Roy的评论所述。

示例：
用mex-msfm2d.c编译c代码；mex-msfm3d.c；mex-rk4.c；

在msfm2d、最短路径和骨架的帮助下尝试这些示例