3D图像侵蚀 - Python

Question

我想使用python对3D CT卷进行侵蚀。我已经尝试过侵蚀功能舞会SciPy和Scikit-image，但它们似乎没有正常工作（两者都没有扩张）。这是我试过的：

import numpy as np
from scipy import ndimage
from skimage import morphology

np_image_data = sitk.GetArrayFromImage(imageData) #Numpy array with CT data
boneMask = np_image_data>=1000
struct = ndimage.generate_binary_structure(3, 1)

# Scipy erosion
erodedMask1 = ndimage.binary_erosion(boneMask.astype(uint), structure=struct, iterations=1)

# Scikit-image erosion
erodedMask2 = morphology.binary_erosion(boneMask.astype(uint), struct)

它们都侵蚀了整个图像（只剩下一些体素）。有什么建议吗？我尝试使用SimpleITK的BinaryErodeImageFilter，但它很慢并且占用了大量内存。

提前谢谢大家！