小编B. *_*her的帖子

我从NIfTI文件中读取了一个3D图像到SimpleITK(使用python),获取每个轴向切片,用它做一些事情并将新的2D切片重新插入到具有(希望)适当尺寸的3D体积中.例如,

output = sitk.Image(original.GetSize(), sitk.sitkFloat32)
output.CopyInformation(original)
for z in numpy.arange(original.GetDepth()):
    image = original[:,:,z]
    << Do Something in SimpleITK>>
    << Produce a new 2D image = newimage >>
    output[:,:,z] = newimage

最后一步是抛出一个错误

In [???]: (executing line ??? of "code.py")
Traceback (most recent call last):
  File "code.py", line ???, in <module>
    output[:,:,z] = newimage
  File "/Library/Python/2.7/site-packages/SimpleITK-0.8.1-py2.7-macosx-10.10-intel.egg/SimpleITK/SimpleITK.py", line 3894, in __setitem__
    raise IndexError("invalid index")
IndexError: invalid index

在for循环中完成最后一步的正确语法(或命令集)是什么？