如何使用 vtkImageReSlice？

Question

我正在尝试使用 vtkImageReSlicer 从 3d 中提取 2d 切片 vtkImageData 对象。但我似乎无法得到正确的食谱。我做对了吗？
我对 ResliceAxes Matrix 也有点困惑。它代表切割平面吗？如果我移动 ReSliceAxes 原点是否也会移动切割平面？当我在 vtkImageReSlicer 上调用 Update，程序崩溃。但是当我不调用它，输出为空。

这是我目前所拥有的。

#my input is any vtkactor that contains a closed curve of type vtkPolyData
ShapePolyData = actor.GetMapper().GetInput()
boundingBox = ShapePolyData.GetBounds()
for i in range(0,6,2):
     delta = boundingBox[i+1]-boundingBox[i]
     newBoundingBox.append(boundingBox[i]-0.5*delta)
     newBoundingBox.append(boundingBox[i+1]+0.5*delta)

voxelizer = vtk.vtkVoxelModeller()
voxelizer.SetInputData(ShapePolyData)
voxelizer.SetModelBounds(newBoundingBox)
voxelizer.SetScalarTypeToBit()
voxelizer.SetForegroundValue(1)
voxelizer.SetBackgroundValue(0)
voxelizer.Update()

VoxelModel =voxelizer.GetOutput()

ImageOrigin = VoxelModel.GetOrigin()
slicer = vtk.vtkImageReslice()
#Am I setting the cutting axis here. x axis set at 1,0,0 , y axis at 0,1,0 and z axis at 0,0,1
slicer.SetResliceAxesDirectionCosines(1,0,0,0,1,0,0,0,1)
#if I increase the z value, will the cutting plane move up?
slicer.SetResliceAxesOrigin(ImageOrigin[0],ImageOrigin[1],ImageOrigin[2])
slicer.SetInputData(VoxelModel)
slicer.SetInterpolationModeToLinear()
slicer.SetOutputDimensionality(2)
slicer.Update() #this makes the code crash

voxelSurface = vtk.vtkContourFilter()
voxelSurface.SetInputConnection(slicer.GetOutputPort())
voxelSurface.SetValue(0, .999)

voxelMapper = vtk.vtkPolyDataMapper()
voxelMapper.SetInputConnection(voxelSurface.GetOutputPort())
voxelActor = vtk.vtkActor()
voxelActor.SetMapper(voxelMapper)
Renderer.AddActor(voxelActor)

Answer 1

我从来没有用过vtkImageReslice，但是我用vtkExtractVOI代替vtkImageData，我想这可以让你达到类似的结果。这是用后者修改的示例，而不是：

ImageOrigin = VoxelModel.GetOrigin()
slicer = vtk.vtkExtractVOI()
slicer.SetInputData(VoxelModel)
#With the setVOI method you can define which slice you want to extract
slicer.SetVOI(xmin, xmax, ymin, ymax, zslice, zslice)
slicer.SetSampleRate(1, 1, 1)
slicer.Update()

voxelSurface = vtk.vtkContourFilter()
voxelSurface.SetInputConnection(slicer.GetOutputPort())
voxelSurface.SetValue(0, .999)

voxelMapper = vtk.vtkPolyDataMapper()
voxelMapper.SetInputConnection(voxelSurface.GetOutputPort())
voxelActor = vtk.vtkActor()
voxelActor.SetMapper(voxelMapper)
Renderer.AddActor(voxelActor)

如何使用 vtkImageReSlice？

how to use vtkImageReSlice?

vtk

python-3.x