将 pixel_array 交换为 DICOM RT Dose 文件 - 值更改 - pydicom

Question

我有一个 DICOM RT Dose 文件，我想用一个新文件更改原始文件 pixel_array。问题是当我给某个像素赋予新值时，它会在过程中的某个地方随机变化。下面给出了交换的代码块。我一次交换一个像素只是为了检查问题出在哪里，所以我知道可以更快地完成：

print pixel_arrayFLK[86][85][78]
print dosem.pixel_array[86][85][78]
for a in range(zBinsTPS):
    for b in range(yBinsTPS):
        for c in range(xBinsTPS):
            dosem.pixel_array[a][b][c] = pixel_arrayFLK[a][b][c]
            if a == 86 and b == 85 and c == 78:
                print dosem.pixel_array[a][b][c]
                print pixel_arrayFLK[a][b][c]

原来的pixel_array是dosem.pixel_array，新的是pixel_arrayFLK。作为发生情况的示例：

pixel_arrayFLK[86][85][78]的值为65813，原来的dosem.pixel_array[86][85][78]的值为62947，交换后pixel_arrayFLK[86][85][78] 仍然具有值 65813 而 dosem.pixel_array[86][85][78] 突然现在具有值 277.

打印函数产量：

但是，如果我将 pixel_arrayFLK 的所有值除以 2，交换工作正常。我真的不明白发生了什么。它可能与 DICOM 文件中的数据类型有关吗？ DICOM 中包含的更多信息：

(0028, 0002) Samples per Pixel                   US: 1
(0028, 0004) Photometric Interpretation          CS: 'MONOCHROME2'
(0028, 0100) Bits Allocated                      US: 16
(0028, 0101) Bits Stored                         US: 16
(0028, 0102) High Bit                            US: 15
(0028, 0103) Pixel Representation                US: 0
(7fe0, 0010) Pixel Data                          OW or OB: Array of 6750000 bytes

我有另一个 DICOM 文件，其中分配的位 = 32，存储的位 = 32，高位 = 31。此文件中的 pixel_array 交换工作正常。

Answer 1

想通了。由于位数为 16，因此最大值为 65536，其中 65813 覆盖此值 277，因此从 0 开始。

将 pixel_array 交换为 DICOM RT Dose 文件 - 值更改 - pydicom

Exchanging the pixel_array for a DICOM RT Dose file - values change - pydicom

python

dicom

pydicom