derive characteristic shape from a set of grain boundaries
cshape = characteristicShapeN(grains.boundary('b','b'))