8.3 DICOM数据处理：医疗影像预处理全流程
在医疗影像AI应用中，DICOM（Digital Imaging and Communications in Medicine）是最常用的医学影像存储和传输标准。掌握DICOM数据的处理技术对于构建医疗AI系统至关重要。本章将详细介绍DICOM数据格式、处理方法以及完整的预处理流程，帮助你构建高效、标准化的医疗影像处理管道。
DICOM格式概述
DICOM是医学数字成像和通信的国际标准，它不仅定义了医学影像的文件格式，还定义了网络通信协议。DICOM文件包含两个主要部分：
DICOM文件结构
- 文件元信息（File Meta Information）：包含文件的基本信息，如传输语法、媒体存储SOP类等
- 数据集（Dataset）：包含详细的医学信息和图像数据
重要DICOM标签
DICOM使用标签（Tag）来组织信息，每个标签由组号和元素号组成。以下是一些重要的DICOM标签：
| 标签 | 名称 | 描述 |
|---|---|---|
| (0010,0010) | Patient’s Name | 患者姓名 |
| (0010,0020) | Patient ID | 患者ID |
| (0010,0030) | Patient’s Birth Date | 患者出生日期 |
| (0010,0040) | Patient’s Sex | 患者性别 |
| (0008,0020) | Study Date | 研究日期 |
| (0008,0030) | Study Time | 研究时间 |
| (0008,0060) | Modality | 成像模态（CT、MRI等） |
| (0028,0002) | Samples per Pixel | 每像素样本数 |
| (0028,0004) | Photometric Interpretation | 光度解释 |
| (0028,0010) | Rows | 图像行数 |
| (0028,0011) | Columns | 图像列数 |
| (0028,0030) | Pixel Spacing | 像素间距 |
| (0028,1050) | Window Center | 窗位 |
| (0028,1051) | Window Width | 窗宽 |
DICOM数据读取与处理
1. 使用pydicom库处理DICOM文件
importpydicomimportnumpyasnpimportmatplotlib.pyplotaspltfrompydicom.pixel_data_handlers.utilimportapply_modality_lut,apply_voi_lutclassDICOMProcessor:"""DICOM处理器"""def__init__(self):passdefread_dicom(self,file_path):""" 读取DICOM文件 Args: file_path: DICOM文件路径 Returns: DICOM数据对象 """try:dicom_data=pydicom.dcmread(file_path)returndicom_dataexceptExceptionase:print(f"读取DICOM文件失败:{e}")returnNonedefget_dicom_info(self,dicom_data):""" 获取DICOM文件信息 Args: dicom_data: DICOM数据对象 Returns: 包含重要信息的字典 """ifdicom_dataisNone:returnNoneinfo={'PatientName':str(dicom_data.get('PatientName','N/A')),'PatientID':str(dicom_data.get('PatientID','N/A')),'PatientBirthDate':str(dicom_data.get('PatientBirthDate','N/A')),'PatientSex':str(dicom_data.get('PatientSex','N/A')),'StudyDate':str(dicom_data.get('StudyDate','N/A')),'Modality':str(dicom_data.get('Modality','N/A')),'Rows':dicom_data.get('Rows','N/A'),'Columns':dicom_data.get('Columns','N/A'),'PixelSpacing':dicom_data.get('PixelSpacing','N/A'),'WindowCenter':dicom_data.get('WindowCenter','N/A'),'WindowWidth':dicom_data.get('WindowWidth','N/A')}returninfodefextract_pixel_array(self,dicom_data,apply_transforms=True):""" 提取像素数组 Args: dicom_data: DICOM数据对象 apply_transforms: 是否应用变换（模态LUT和VOI LUT） Returns: 像素数组 """ifdicom_dataisNone:returnNonetry:# 获取原始像素数据pixel_array=dicom_data.pixel_arrayifapply_transforms:# 应用模态LUT（将存储值转换为物理值）pixel_array=apply_modality_lut(pixel_array,dicom_data)# 应用VOI LUT（应用窗位窗宽）pixel_array=apply_voi_lut(pixel_array,dicom_data)returnpixel_array