📚 课程目录第3章医学影像数据基础与DICOM工程化理解

3. 医学影像数据基础与DICOM工程化理解

在3D Slicer中处理医学影像，第一步不是点按钮，而是理解数据本身。CT、MR、PET、US四类模态的物理成像机制不同，决定了它们的灰度含义、空间分辨率和临床用途完全不同。CT的像素值直接反映X射线衰减系数，经过Hounsfield单位（HU）标定后，水为0 HU、空气约-1000 HU、致密骨可达+1000 HU以上，因此CT适合做阈值分割和骨窗观察。MR没有统一的绝对灰度标准，信号强度取决于质子密度、T1/T2弛豫时间及序列参数，同一组织在不同序列中可能呈现高信号或低信号，所以MR分割更依赖多序列对比和相对关系。PET记录的是正电子湮灭产生的符合事件，像素值反映放射性示踪剂浓度，常以SUV（标准化摄取值）表达代谢活性，空间分辨率低但功能信息强。US则是回声强度成像，受探头频率、增益和声束角度影响大，图像常带斑点噪声，且没有固定的解剖坐标系。

DICOM是这些模态共同的工程容器。一个DICOM文件由文件头和数据集组成，文件头包含128字节前导码和"DICM"魔数，数据集则由一系列数据元素构成，每个元素用（组号，元素号）唯一标识，例如(0010,0010)是患者姓名，(0020,0032)是图像位置Patient Position，(0028,0030)是像素间距Pixel Spacing。理解标签的关键在于区分"患者相关""设备相关""图像相关""几何相关"四类信息。工程化处理时，我们通常先读取Series Instance UID判断序列归属，再读取Image Orientation (Patient)和Image Position (Patient)重建空间方向，最后用Pixel Spacing和Slice Thickness/Slice Location计算体素尺寸。

体素是三维空间中的最小采样单元，它同时携带灰度值和空间位置。一个体素在物理空间中的尺寸由Pixel Spacing（行、列方向）和层间距共同决定。例如某CT序列Pixel Spacing为0.7mm×0.7mm，层厚1.25mm，则体素体积为0.7×0.7×1.25=0.6125mm³。若层间距与层厚不一致，说明存在层间插值或重叠，重建时要特别小心。

空间坐标是3D Slicer中最容易出错的部分。DICOM使用患者坐标系（LPS：Left-Posterior-Superior），即X轴指向患者左侧，Y轴指向后侧，Z轴指向头侧。而3D Slicer内部使用RAS坐标系（Right-Anterior-Superior），X轴指向患者右侧，Y轴指向前侧，Z轴指向头侧。两者X、Y方向相反，Z方向相同。因此从DICOM读取的Image Position (Patient)不能直接当作Slicer中的坐标使用，必须经过LPS→RAS转换：X_RAS = -X_LPS，Y_RAS = -Y_LPS，Z_RAS = Z_LPS。这个符号约定必须在整个处理流程中保持一致，否则会出现左右翻转或前后颠倒。

窗宽窗位是CT阅片中调节对比度的核心参数。窗宽（Window Width, WW）定义显示的HU范围，窗位（Window Level, WL）定义该范围的中心。显示范围 = [WL - WW/2, WL + WW/2]。例如WW=400、WL=40时，显示范围为-160~240 HU，低于-160的像素显示为纯黑，高于240的显示为纯白，中间线性映射到0~255灰度。肺窗常用WW=1500、WL=-600，骨窗常用WW=2000、WL=400，脑窗常用WW=80、WL=40。在3D Slicer中，Volume模块的Display节点可以实时调整WW/WL，但要注意这只影响显示，不改变原始像素值。做阈值分割时，应基于原始HU值设定阈值，而不是基于显示灰度。

多帧与序列组织是工程化处理的另一个关键。一个DICOM序列可能包含数百个切片，每个切片有独立的SOP Instance UID，但共享同一个Series Instance UID。3D Slicer的DICOM模块会自动按Series分组，但有时会出现同一序列被拆成多个Series的情况，通常是因为Slice Thickness或Image Position不一致。此时需要手动检查Instance Number和Slice Location，必要时用"Add Series"合并。对于多帧DICOM（如超声或动态PET），一个文件内包含多帧图像，需要用NumberOfFrames标签识别，并在Slicer中按帧展开。

匿名化与伦理合规是医学影像处理不可回避的环节。DICOM标签中直接标识患者身份的有(0010,0010)姓名、(0010,0020)患者ID、(0010,0030)出生日期、(0008,0080)机构名称等。间接标识包括检查日期、设备序列号、医生姓名等。工程化匿名化不能只删姓名，还要考虑日期偏移、ID重映射和私有标签清理。3D Slicer的DICOM模块提供匿名化导出功能，但建议在导出后再次用工具检查残留PHI。伦理合规方面，即使数据已脱敏，使用前仍需获得机构审查委员会（IRB）批准或数据使用协议，并保留处理日志以备审计。

实际工程中，建议建立一份DICOM检查清单：确认Modality、检查Pixel Spacing与Slice Thickness是否一致、核对Image Orientation方向余弦、验证LPS→RAS转换符号、检查窗宽窗位是否适合目标组织、确认匿名化字段已清除。这份清单能显著降低后续分割和配准中的空间错误。

把上述内容串起来，一个典型的CT处理流程是：读取DICOM序列→按Series分组并排序→检查体素尺寸和方向→转换到RAS坐标→设置窗宽窗位观察→基于HU阈值分割→三维重建→导出前匿名化。每一步都依赖对DICOM标签和空间坐标的正确理解。只有把数据基础打牢，后续的配准、分割和可视化才不会因为坐标错位或灰度误解而返工。