Question

我正在从此链接运行simpleitk的python示例 https://simpleitk.readthedocs.io/en/master/Examples/ImageRegistrationMethod1/Documentation.html

在间谍水蟒中。

我已经通过anaconda Prompt安装了软件包simpleitk，system，os。

from __future__ import print_function

import SimpleITK as sitk
import sys
import os


def command_iteration(method) :
    print("{0:3} = {1:10.5f} : {2}".format(method.GetOptimizerIteration(),
                                   method.GetMetricValue(),
                                   method.GetOptimizerPosition()))

if len ( sys.argv ) < 4:
    print( "Usage: {0} <fixedImageFilter> <movingImageFile> <outputTransformFile>".format(sys.argv[0]))
cc


fixed = sitk.ReadImage(sys.argv[1], sitk.sitkFloat32)
print(sys.argv[1])
moving = sitk.ReadImage(sys.argv[2], sitk.sitkFloat32)

R = sitk.ImageRegistrationMethod()
R.SetMetricAsMeanSquares()
R.SetOptimizerAsRegularStepGradientDescent(4.0, .01, 200 )
R.SetInitialTransform(sitk.TranslationTransform(fixed.GetDimension()))
R.SetInterpolator(sitk.sitkLinear)

R.AddCommand( sitk.sitkIterationEvent, lambda: command_iteration(R) )

outTx = R.Execute(fixed, moving)

print("-------")
print(outTx)
print("Optimizer stop condition: {0}".format(R.GetOptimizerStopConditionDescription()))
print(" Iteration: {0}".format(R.GetOptimizerIteration()))
print(" Metric value: {0}".format(R.GetMetricValue()))

sitk.WriteTransform(outTx,  sys.argv[3])

if ( not "SITK_NOSHOW" in os.environ ):

    resampler = sitk.ResampleImageFilter()
    resampler.SetReferenceImage(fixed);
    resampler.SetInterpolator(sitk.sitkLinear)
    resampler.SetDefaultPixelValue(100)
    resampler.SetTransform(outTx)

    out = resampler.Execute(moving)
    simg1 = sitk.Cast(sitk.RescaleIntensity(fixed), sitk.sitkUInt8)
    simg2 = sitk.Cast(sitk.RescaleIntensity(out), sitk.sitkUInt8)
    cimg = sitk.Compose(simg1, simg2, simg1//2.+simg2//2.)
    sitk.Show( cimg, "ImageRegistration1 Composition" )

在运行此代码时，出现错误后出现错误用法：E：/registration/simpleitk.py 发生异常，请使用％tb查看完整的追溯。

SystemExit：1

C：\ Users \ aBC \ AppData \ Local \ Continuum \ anaconda3 \ lib \ site-packages \ IPython \ core \ interactiveshell.py：3275：   UserWarning：要退出：请使用“退出”，“退出”或Ctrl-D。警告（“退出：   使用“退出”，“退出”或Ctrl-D。”，stacklevel = 1）

我如何纠正此错误？

Answer 1

我怀疑脚本无法在对sitk.ReadImage的调用中加载图像。期望在sys.argv [1]变量中提供文件名。该脚本被设计为在命令行上运行，并以命令行参数的形式提供固定图像文件名，运动图像文件名和输出转换文件名。

如果您无法在Spyder中提供命令行参数（我不知道；我不熟悉），则可以在脚本中对文件名进行硬编码。

如果您对SimpleITK还有其他疑问，建议您查看ITK Discourse，discourse.itk.org。

Python在“ sys.exit（1）”上引发异常“ SystemExit：1”

1 个答案: