image - 如何将 DICOM 从 Monochrome 1 转换为 Monochrome 2？

Question

我正在处理一个包含 DICOM 图像的项目，我需要比较两个 DICOM 图像。问题是，一个是单色 1，另一个是单色 2（零分别表示白色和黑色）。如何转换这些像素强度以进行比较？我正在使用“pydicom”工具包。

Answer 1

如果有帮助，在使用 matplotlib.pyplot 进行可视化时，您可以使用

plt.imshow(image, cmap='gray_r')

在不改变像素值的情况下将像素反转回 Monochrome2 以进行视觉比较。

还，

np.invert(image)

可能是一种解决方法。

Answer 2

您的主要问题不是光度解释（MONO1/2）。

您无法比较两个 DICOM 图像的像素强度，除非它们使用相同的比例（例如 Hounsfield 单位）。

如果你有

(0028,1052) RescaleIntercept - present with any value
(0028,1053) RescaleSlope - present with any value
(0028,1054) RescaleType - present with value "OD" or "HU"

然后很容易：应用线性变换：

<measured value> = <pixel value> * RescaleSlope + RescaleIntercept

可以比较测量值。

如果您将非线性模态 LUT 作为查找表存储在标题中，情况也是如此，但同样的限制适用于重新缩放类型。

否则我会避免比较像素值。当然，将两张图像中的一张反转似乎很容易，但它们具有不同的光度解释这一事实告诉我，它们是通过不同的设备或技术获得的。这意味着，像素数据在视觉上是正确且可比较的，但在数学上不相关。