我有一个表示为体素化图像的薄结构。我正在使用来自 scikit-image 的行进立方体,但我注意到它不能很好地处理薄结构。
相同结构的示例移动了半个体素:
import numpy as np
from skimage.measure import marching_cubes
arr1 = np.zeros((4,4,4))
arr2 = np.zeros((4,4,4))
# Perfect match of structure and voxel layout
arr1[:,1,:] = 0.5
# Mismatch of structure and voxel layout
arr2[:,1:3,:] = 0.25
# Finds the surface
v1, s2, n1, val1 = marching_cubes(arr1, level=0.4)
# Does not find the surface
v2, s2, n2, val2 = marching_cubes(arr2, level=0.4)
所以第二个失败的原因是没有高于的值level
。但是,知道我的结构是二进制结构(不存在中间值),我想知道是否可以通过以下方式补偿这些插值伪影
- 使用不同的行进立方体算法?
- 使用分辨率增加?