1

我正在尝试T2从我的文件制作地图dicom，但我的代码遇到以下错误：

RuntimeWarning: invalid value encountered in ushort_scalars
T2 map:   if math.log((I8[j][k]) / (I1[j][k])) == 0:

ExampleName.py:91: RuntimeWarning: invalid value encountered in ushort_scalars
  if (3.03 - (17.85 + 14.87)) / (math.log((I8[j][k]) / (I1[j][k]))) < 0:

ExampleName.py:93: RuntimeWarning: invalid value encountered in ushort_scalars
  elif (3.03-(17.85+14.87))/(math.log((I8[j][k]) / (I1[j][k]))) > 100:

ExampleName.py:96: RuntimeWarning: invalid value encountered in ushort_scalars
  T2_map[j,k,i] = ((3.03-(17.85+14.87))/(math.log((I8[j][k]) / (I1[j][k]))))
[[[ 0. nan nan ...  0.  0.  0.]

这是我正在使用的代码示例：

#start of code snippet
T2_map = np.empty((number_of_segments, number_of_segments, number_of_segments))
for i in range(1, number_of_slices + 1):
    if 3 + (i - 1) <= number_of_segments:
        I1 = (dcm_pixel_array["array" + str(3 + (i - 1))]) * 2
    else:
        I1 = (dcm_pixel_array["array" + str(number_of_segments)]) * 2
    if 3 + ((7 * number_of_slices) + (i - 1)) <= number_of_segments:
        I8 = (dcm_pixel_array["array" + str(3 + ((7*number_of_slices) + (i-1)))]) * 2
    else:
        I8 = (dcm_pixel_array["array" + str(number_of_segments)]) * 2
    for j in range(I1.ndim):
        for k in range(I1[0].ndim):
            if math.log((I8[j][k]) / (I1[j][k])) == 0:
                T2_map[j,k,i] = "NAN"
            else:
                if (3.03 - (17.85 + 14.87)) / (math.log((I8[j][k]) / (I1[j][k]))) < 0:
                    T2_map[j,k,i] = 0
                elif (3.03-(17.85+14.87))/(math.log((I8[j][k]) / (I1[j][k]))) > 100:
                    T2_map[j,k,i] = 100
                else:
                    T2_map[j,k,i] = ((3.03-(17.85+14.87))/(math.log((I8[j][k]) / (I1[j][k]))))
#end of code snippet

澄清：number_of_segments = 288

我不确定这些错误是什么意思或如何纠正它们。欢迎任何建议。

感谢您的帮助！

4

1 回答 1

1

问题是我的程序在某个点被零除，所以我必须放入一个 if 语句来检查 I1 和 I8/I1 是否为零值。如果其中一个为 0，我设置 T2[j,k,i] = "NAN"。

于 2020-08-06T21:22:57.020 回答