【医疗影像处理】DICOM Rescale Intercept / Rescale Slope

What is Rescale Intercept / Rescale Slope?

DICOM tags :

• Rescale intercept (0028|1052)
• Rescale slope (0028|1053)

缩放斜率和截距由硬件制造商决定。
它指定从存储在磁盘表示中的像素到存储在内存表示中的像素的线性转换。磁盘存储的值定义为SV。而转化到内存中的像素值uints就需要两个dicom tag : Rescale intercept和Rescale slope。

$$OutputUnits=m\ast SV+b$$

$$RescaleIntercept:b$$

$$RescaleSlope:m$$

为什么磁盘上的表示与内存中不同?
存储在磁盘上的表示可能与存储在内存中的表示有不同的值范围。
例如:CT图像的像素值是用Hounsfield单位来测量的，而Hounsfield单位可以是负值，因此通常存储为无符号整数。
因此，CT DICOM文件通常具有负截距。线性缩放也适用于这样的情况，一个像素可能有一个大范围的值，而存储的值尽可能少的比特，同时避免量化错误。这通常应用于PET成像，在PET成像中，所测量的活度值的范围可能超过小整数。此外，PET的最大活性可能因切片而异;例如，与健康组织切片相比，带有肿瘤的切片可能具有很高的值。由于这个原因，Rescale intercept/Rescale slope在一个PET图像系列的切片基础上经常变化。

转化为Hu

def rescale_pixelarray(dataset):
    image = dataset.pixel_array
    rescaled_image = image * dataset.RescaleSlope + dataset.RescaleIntercept
    return rescaled_image