itk 快速行进输出

Question

我正在使用 itk 进行一些预处理，并且我想使用 Fast Marching 过滤器和 Geodesic Active Contour 过滤器进行测试。 我遵循 ITK 软件指南第 9.3.3 节 (http://www.itk.org/ItkSoftwareGuide.pdf) 中描述的算法。 但是我没有得到预期的结果。我正在处理 3D 图像。

你可以看到我的代码here

我正在使用this paper 中描述的指标，目标与我的相同。

我有几个问题：

1. 快速行进过滤器应该输出距离图吧？相反，当我将音量输出到一系列 png（值介于 0 和 4095 之间）时，我有一个二进制图像（像素为 0 或 4095）。我想我应该得到一个灰度体积，指示从种子中获得每个像素所需的时间。
2. 按照 Suzuki 描述的过程，我成功地使算法或多或少地工作，但是我更改了测地线过滤器的参数值。我不记得确切的值，但它与论文中描述的值不接近。当我们处理在 0 和 1 之间归一化的 sigmoid 输入时，发生了什么？
3. 我应该为快速行进滤波器还是 sigmoid 图像使用恒速函数？什么时候应该首选这两种方法？
4. 我正在使用重新缩放器来输出我的浮动图像（过滤器的输出）。这可能是我看到的不一致的原因吗？
5. 关于我可能做错的任何建议？

谢谢。

Answer 1

好的，我发现了我的问题。 Fast Marching 过滤器确实输出了一个时间交叉图（距离图），但是当我在算法中指定了一个停止值时，所有未访问的像素都有一个高值（1.7e+38，因为它是最大值的一半用于输出图像的类型，在我的情况下是浮动的 3.4e+38)。因此，当我使用重新缩放滤镜时，它压缩了我所有的图像动态，结果是一个二值图像。 我认为使用 sigmoid 图像作为快速行进滤波器的输入可以获得更好的结果。 感谢@nav 的建议。