BinaryReconstructionByErosion works on 2D but not 3D images

[blowekamp (Bradley Lowekamp)]

The filter is working correctly.

The segmented thighs are open at either end, so the erosion of the mask image occurs from those open ends.

I don’t see a problem doing this operation on a series of 2d slices.

Here is some example code of doing slice-by-slice operations in python with SimpleITK:

ResampleImageFilter 4D Images

Looks like a good solution to me. I have done similar things when working with 4D XYZC images. This approach is really asking for a reusable solution something like a python function decorator to operate slice_by_slice. Here is what I came up with: def slicing_decorator(func): """ A function decorator which extracts image slices of N-1 dimensions and calls func on each slice. The resulting images are then concatenated together with JoinSeries. :param func: A function whi…