python code to exclude pixel islands of a given size, requires mahotas for building a label matrix (could use pymorph instead)

gistfile1.py

import mahotas

def excludePixels(image, size_cutoff=1):
    labeled_image = mahotas.label(image)[0]

    for label_id in range(labeled_image.max()+1):
        label_id_index = labeled_image == label_id
        if label_id_index.sum() <= size_cutoff:
            labeled_image[label_id_index] = 0

    return labeled_image>0