如何使用VTK体素索引获得世界坐标？

import vtk

# 创建一个vtkImageData对象
imageData = vtk.vtkImageData()

# 定义体素的大小和数量
voxelSize = [1.0, 1.0, 1.0]  # 体素的尺寸
voxelCount = [100, 100, 100]  # 体素的数量

# 设置vtkImageData的尺寸和原点
imageData.SetDimensions(voxelCount)
imageData.SetSpacing(voxelSize)
imageData.SetOrigin(0.0, 0.0, 0.0)

# 填充三维数据到vtkImageData中
# 这里假设已经将数据填充到一个名为data的numpy数组中
data = ...  # 填充三维数据的过程

imageData.AllocateScalars(vtk.VTK_DOUBLE, 1)
for z in range(voxelCount[2]):
    for y in range(voxelCount[1]):
        for x in range(voxelCount[0]):
            voxelIndex = [x, y, z]
            worldCoord = imageData.GetPoint(voxelIndex)
            dataValue = data[z, y, x]
            imageData.SetScalarComponentFromDouble(voxelIndex, 0, dataValue)

# 使用体素索引获取世界坐标
voxelIndex = [50, 50, 50]  # 体素索引
worldCoord = imageData.GetPoint(voxelIndex)

print("Voxel Index:", voxelIndex)
print("World Coordinate:", worldCoord)