itk::LevelTracingImageFilter工作原理及实现机制咨询

这个滤波器是3D Slicer在vtkITK适配层实现的定制ITK类，不属于ITK官方标准发布模块，所以官方ITK文档站没有收录对应的原理说明，核心逻辑可以直接从头文件声明梳理清楚，本质是一个快速等值线/等值面追踪分割工具，具体运行机制如下：

• 初始化阶段：读取用户输入的种子点坐标，提取种子点位置的灰度值作为固定等值基准，整个追踪过程不会调整这个基准值。
• 邻域遍历阶段：从种子点出发，按照二维4连通/三维6连通规则扫描相邻体素，通过线性插值判断两个相邻体素之间是否穿过目标等值面——当两个相邻体素的灰度值分别落在基准值的两侧时，就判定两点间存在等值面交点。
• 轮廓追踪阶段：每找到一个等值交点，就计算该位置的图像梯度方向，沿着垂直于梯度的等值切线方向移动到下一个邻域位置，持续延伸轮廓；追踪过程中会记录所有已经访问过的轮廓点，当路径回到已访问点、形成闭合轮廓时就停止生长。
• 结果生成阶段：最终输出和输入图像尺寸一致的二值标记图，既可以配置为只标记轮廓线本身为前景，也可以标记闭合轮廓包围的整个内部区域为前景，同时支持输出所有轮廓点的物理坐标集合。

实用提示：这个滤波器和ITK中常见的迭代演化类水平集滤波器逻辑完全不同，它没有数值迭代求解过程，运行速度非常快，但对图像噪声和灰度不均的抗性很差，实际使用前建议先对输入图像做高斯平滑预处理，避免追踪出的轮廓出现锯齿、断裂或者偏移。

内容的提问来源于stack exchange，提问作者Eleonora Tiribilli