SimpleITK保存DICOM时BitsAllocated/BitsStored标签异常求助
嘿，我来帮你拆解下这个问题：你的核心困扰是手动设置的BitsStored=12被SimpleITK覆盖成了16，对吧？这其实是因为SimpleITK的DICOM Writer会自动根据图像的像素数据类型生成并覆盖这类像素相关标签——你的图像是uint16类型，Writer默认会认为所有16位都在被使用，所以直接把BitsStored改成了16，无视了你手动设置的元数据。
一步步解决问题
1. 先确保数据符合12位的范围要求
首先，BitsStored=12意味着图像数据的最大值不能超过2^12 - 1 = 4095，这是DICOM标准的硬性要求。如果你的修改后数据超出了这个范围，得先做截断或缩放：
# 截断数据到0-4095之间，确保符合12位有效范围 sitk_stack = sitk.Clamp(sitk_stack, lowerBound=0, upperBound=4095) # 如果原数据范围更大，也可以用缩放把数据映射到0-4095： # sitk_stack = sitk.RescaleIntensity(sitk_stack, outputMinimum=0, outputMaximum=4095)
2. 配置Writer，阻止自动覆盖标签
默认情况下，SimpleITK的DICOM Writer会自动生成和像素类型相关的标签，我们需要显式配置DICOMImageIO来让它优先使用我们手动设置的元数据：
writer = sitk.ImageFileWriter() writer.KeepOriginalImageUIDOn() # 创建并配置DICOM IO对象 dicom_io = sitk.DICOMImageIO() dicom_io.KeepOriginalImageUIDOn() dicom_io.CopyPrivateTagsOn() # 如果需要保留私有标签的话 # 关键：让IO不要自动根据像素类型生成BitsStored/BitsAllocated这类标签 # 当我们手动设置了这些元数据后，IO会优先使用我们的值 # 把配置好的IO绑定到Writer writer.SetImageIO(dicom_io)
3. 修正代码里的小错误
我注意到你代码里的Instance Number标签键写错了："0020,0013"应该是"0020|0013"（用竖线分隔标签组和元素，而不是逗号），这个小问题可能导致Instance Number标签设置失败，改过来就行：
# 修正Instance Number的标签键 image_slice.SetMetaData("0020|0013", str(i)) # Instance Number
另外，你设置的HighBit(0028|0102)=11是对的——因为BitsStored=12，HighBit是从0开始计数的最高有效位，所以11是正确值。
关于“能不能仅用12位存储图像”的疑问
DICOM标准里并没有原生的12位无符号整数存储格式，行业内通用的做法是用16位的容器（BitsAllocated=16）来存储12位的有效数据（BitsStored=12）。这样既符合标准，又能在使用压缩（比如JPEG 2000）时达到近似12位存储的空间效率。所以你没办法直接存成纯12位的DICOM，只能通过设置BitsStored=12并约束数据范围来实现等效的效果。
内容的提问来源于stack exchange，提问作者F.Lotte
