python中的matlab bwmorph（图像，'spur'）答案

【问题标题】：matlab's bwmorph(image, 'spur') in pythonpython中的matlab bwmorph（图像，'spur'）
【发布时间】：2018-10-24 18:59:12
【问题描述】：

我正在将一个 matlab 图像处理脚本移植到 python/skimage 中，但无法找到 Matlab 的 bwmorph 函数，特别是 skimage 中的 'spur' 操作。 matlab 文档对spur 操作这样说：

Removes spur pixels. For example:
0  0  0  0           0  0  0  0
0  0  0  0           0  0  0  0
0  0  1  0  becomes  0  0  0  0
0  1  0  0           0  1  0  0
1  1  0  0           1  1  0  0

我已经在 python 中实现了一个版本，可以很好地处理上述情况：

def _neighbors_conv(image):
    image = image.astype(np.int)
    k = np.array([[1,1,1],[1,0,1],[1,1,1]])
    neighborhood_count = ndimage.convolve(image,k, mode='constant', cval=1)
    neighborhood_count[~image.astype(np.bool)] = 0
    return neighborhood_count

def spur(image):
    return _neighbors_conv(image) > 1

def bwmorph(image, fn, n=1):
    for _ in range(n):
        image = fn(image)
    return image

t= [[0, 0, 0, 0],
    [0, 0, 1, 0],
    [0, 1, 0, 0],
    [1, 1, 0, 0]]
t = np.array(t)
print('neighbor count:')
print(_neighbors_conv(t))
print('after spur:')
print(bwmorph(t,spur).astype(np.int))

neighbor count:
[[0 0 0 0]
 [0 0 1 0]
 [0 3 0 0]
 [7 5 0 0]]
after spur:
[[0 0 0 0]
 [0 0 0 0]
 [0 1 0 0]
 [1 1 0 0]]

上述方法通过删除任何只有一个相邻像素的像素来工作。

我注意到上述实现的行为与 matlab 的 spur 操作不同。以matlab为例：

0     0     0     0     0
0     0     1     0     0
0     1     1     1     1
0     0     1     0     0
0     0     0     0     0

becomes, via bwmorph(t,'spur',1):

0      0     0     0     0
0      0     0     0     0
0      0     1     1     1
0      0     0     0     0
0      0     0     0     0

支线操作比查看 8 个邻居计数要复杂一些。我不清楚如何扩展我的实现以满足这种情况而又不使其过于激进（即删除有效像素）。

matlab 的spur 的底层逻辑是什么，或者有没有我可以使用的python 实现？

更新：我发现 Octave 使用 LUT 实现的支线：

case('spur')
      ## lut=makelut(inline("xor(x(2,2),(sum((x&[0,1,0;1,0,1;0,1,0])(:))==0)&&(sum((x&[1,0,1;0,0,0;1,0,1])(:))==1)&&x(2,2))","x"),3);
      ## which is the same as
      lut=repmat([zeros(16,1);ones(16,1)],16,1); ## identity
      lut([18,21,81,273])=0; ## 4 qualifying patterns
      lut=logical(lut);
      cmd="BW2=applylut(BW, lut);";

（来自https://searchcode.com/codesearch/view/9585333/）假设这是正确的，我只需要能够在 python 中创建这个 LUT 并应用它......

【问题讨论】：

我认为 MATLAB 的所有 bwmorph 选项都使用 LUT。创建 LUT 时，请考虑所需形状的旋转！

标签： python matlab scikit-image

【解决方案1】：

我最终实现了我自己的spur 版本和bwmorph 的其他操作。对于有相同需求的未来互联网旅行者，这里是我最终使用的一个方便的要点：

https://gist.github.com/bmabey/4dd36d9938b83742a88b6f68ac1901a6

【讨论】：