Search code examples
pythonopencvimage-processingscikit-imageimage-morphology

Failed to remove noise by remove_small_objects


I have a white and black image. I try to remove noise by remove_small_objects.

import cv2 as cv
import numpy as np
from skimage import morphology

img = np.array([[255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255],
                [255, 255,   0, 255,   0,   0,   0,   0, 255, 255, 255],
                [255, 255, 255, 255,   0,   0,   0,   0, 255,   0,   0],
                [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                [255, 255,   0,   0,   0,   0,   0, 255,   0,   0,   0],
                [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                [255, 255, 255,   0,   0,   0,   0,   0,   0,   0,   0],
                [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0]])

cleaned = morphology.remove_small_objects(img, min_size=10, connectivity=1)
print(cleaned)

while True:
    cv.imshow('Demo', cleaned.astype(np.uint8))
    if cv.waitKey(1) & 0xFF == 27:
        break

cv.destroyAllWindows()

However, it didn't work as I expected. The white pixel 255 in the middle is still there.

Did I do something wrong? Thanks

enter image description here


Solution

  • From the docs (emphasis mine):

    skimage.morphology.remove_small_objects(ar, min_size=64, connectivity=1, in_place=False)

    Remove objects smaller than the specified size.

    Expects ar to be an array with labeled objects, and removes objects smaller than min_size. If ar is bool, the image is first labeled. This leads to potentially different behavior for bool and 0-and-1 arrays.

    import numpy as np
    from skimage import io, morphology
    import matplotlib.pyplot as plt
    
    img = np.array([[255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255],
                    [255, 255,   0, 255,   0,   0,   0,   0, 255, 255, 255],
                    [255, 255, 255, 255,   0,   0,   0,   0, 255,   0,   0],
                    [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                    [255, 255,   0,   0,   0,   0,   0, 255,   0,   0,   0],
                    [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                    [255, 255, 255,   0,   0,   0,   0,   0,   0,   0,   0],
                    [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0],
                    [255, 255,   0,   0,   0,   0,   0,   0,   0,   0,   0]])
    
    arr = img > 0
    cleaned = morphology.remove_small_objects(arr, min_size=2)
    cleaned = morphology.remove_small_holes(cleaned, min_size=2)
    
    fig, axs = plt.subplots(1, 2)
    axs[0].imshow(img, cmap='gray')
    axs[0].set_title('img')
    axs[1].imshow(cleaned, cmap='gray')
    axs[1].set_title('cleaned')
    plt.show(fig)
    

    plot