blurred_grayscale = scipy.ndimage.filters.gaussian_filter(original_grayscale,
    blur_parameter)
difference_image = original_grayscale - multiplier * blurred_grayscale
Counting particles using image processing in python