如何将DICOM的pixel_array像素值缩放到[0,1]范围 pydicom有无内置函数

pydicom没有提供直接将像素值缩放到[0,1]这类自定义区间的单一内置函数，但提供了DICOM像素值预处理的相关内置工具，可以配合简单运算完成自定义区间的缩放，不需要手动查找.pixel_array的取值范围。

具体实现步骤 #

• 第一步：调用pydicom内置的apply_modality_lut函数，将.pixel_array返回的原始存储值转换为真实物理值（比如CT的HU值），自动适配DICOM文件中的rescale slope、rescale intercept参数
• 第二步：如果是用于图像显示场景，可调用内置的apply_voi_lut函数，根据DICOM内置的窗宽窗位参数，将真实物理值转换为适合显示的灰度值
• 第三步：通过简单的数值运算，将处理后的像素值缩放到你需要的自定义区间

代码示例 #

from pydicom import dcmread
from pydicom.pixel_data_handlers.util import apply_modality_lut, apply_voi_lut
import numpy as np

# 读取目标DICOM文件
ds = dcmread("sample.dcm")
# 转换为真实物理值
real_pixel = apply_modality_lut(ds.pixel_array, ds)
# 应用窗宽窗位得到显示灰度值（如果不需要窗位处理可跳过这一步）
windowed_pixel = apply_voi_lut(real_pixel, ds)
# 缩放到[0,1]区间
normalized_pixel = windowed_pixel / np.max(windowed_pixel)

如果不需要窗宽窗位处理，直接基于全局取值范围缩放可以用以下逻辑：

min_val = np.min(real_pixel)
max_val = np.max(real_pixel)
normalized_pixel = (real_pixel - min_val) / (max_val - min_val)

注意事项 #

• 部分DICOM文件可能存在无效填充像素（比如边缘的0值），建议缩放前先过滤无效值，避免最大最小值计算偏差
• 以上逻辑对单帧、多帧DICOM文件均适用，无需额外适配

内容的提问来源于stack exchange，提问作者deepthought