Should possibly go in segysak, if best done from volume. Or possibly in xarray-spatial, if not really an I/o thing. Could do any contour in principle. (Part of volumetric workflow.)