code: def shape_based_slice_insertation_object_wise(sl1, sl2, dim, nslices, order=3):
      
      out = None
      sandwich = numpy.concatenate((sl1[numpy.newaxis], sl2[numpy.newaxis]), 0)
      label_image, n_labels = label(sandwich)
      for lid in range(1, n_labels + 1):
          _sl1, _sl2 = label_image == lid
          _out = shape_based_slice_insertation(_sl1, _sl2, dim, nslices, order=3)
          if out is None:
              out = _out
          else:
              out |= _out
      return out
docs: Wrapper to apply `shape_based_slice_insertation()` for each binary object
      separately to ensure correct extrapolation behaviour.