python - SimpleITK 有选择地改变像素/切片

Question

我已在 SimpleITK 中加载 CT 扫描。我想在 NumPy 中做一些非常简单的事情，但还没有弄清楚如何在 SimpleITK 中完成它们。为了速度，我想在 SimpleITK 中完成它们。

# NumPy: Changes all values of 100 to now become 500
nparr = nparr[nparr == 100] = 500

# SimpleITK:
???

SimpleITK image==100 将生成相同尺寸的二值图像，其中所有强度==100 均为 1/True。这是所希望的。但不幸的是，我不相信 SimpleITK 支持 bool 索引。实现这一目标最有效的方法是什么？

我想出了这个看起来很时髦的东西；但我希望找到这样做的预期方法/最佳实践手段:

# Cast because data type returned is uint8 otherwise
difference = 500 - 100
offset = SimpleITK.Cast( image == 100), sitk.sitkInt32 ) * difference
image += offset

Answer 1

您可以使用 BinaryTheshold 过滤器。

结果 = satk.BinaryThreshold( 图像, 100, 101, 500, 0 )

这应该只选择强度为 100 的像素。