How to find the diameter of objects using image processing in Python?

Question

Given an image with some irregular objects in it, I want to find their individual diameter.

Thanks to this answer, I know how to identify the objects. Howeve

Answer 1

I would propose using a distance transform. So once you've got your labeled image you do:

dt = ndimage.distance_transform_edt(blobs)
slices = ndimage.find_objects(input=labels)
radii = [np.amax(dt[s]) for s in slices]

This gives the largest inscribed circle (or sphere in 3D). The find_objects function is quite handy. It returns a list of Python slice objects, which you can use to index into the image at the specific locations containing the blobs. These slices can of course be used to index into the distance transform image. Thus the largest value of the distance transform inside the slice is the radius you're looking for.

There is one potential gothcha of the above code: the slice is a square (or cubic) section so might contain small pieces of other blobs if they are close together. You can get around this with a bit more complicated logic as follows:

radii = [np.amax(dt[slices[i]]*(labels[slices[i]] == (i+1))) for i in range(nlabels)]

The above version of the list comprehension masks the distance transform with the blob that is supposed to be indexed by the slice, thereby removing any unwanted interference from neighboring blobs.