使用MATLAB创建DICOM-RT轮廓数据3D掩码时遇到问题

问题场景 #

从DICOM CT图像中利用RT结构的GTV提取肿瘤掩码，读取CT和RT结构的代码可正常运行，但通过createMask生成掩码后，用volshow(rtMask)无法看到掩码痕迹，怀疑referenceInfo设置有误。

排查与解决步骤 #

1. 修正referenceInfo的空间参数匹配 #

手动设置3,3,3可能与CT图像的实际空间分辨率不匹配，导致掩码与CT空间错位。正确做法是从CT的DICOM元数据中读取真实参数：

% 读取CT文件夹中任意一个DICOM文件获取元数据
ctSampleInfo = dicominfo(fullfile('DicomCT', dir('DicomCT/*.dcm').name));
% 获取x/y方向像素间距、z方向层厚
pixelSpacing = ctSampleInfo.PixelSpacing;
sliceThickness = ctSampleInfo.SliceThickness;
% 创建匹配CT的referenceInfo
referenceInfo = imref3d(size(V), pixelSpacing(1), pixelSpacing(2), sliceThickness);

2. 确认GTV轮廓的索引正确性 #

createMask的第二个参数是轮廓索引，需确认该索引对应GTV：

% 遍历所有轮廓，查看名称与索引对应关系
for idx = 1:length(rtContours)
    fprintf('%d: %s\n', idx, rtContours(idx).ROIName);
end

若GTV不是第一个轮廓，替换createMask中的索引值即可。

3. 验证掩码本身的有效性 #

可视化前先检查掩码是否为空或尺寸匹配：

% 统计掩码中非零像素数量（肿瘤区域）
fprintf('掩码非零像素数：%d\n', sum(rtMask(:)));
% 对比CT与掩码的尺寸
fprintf('CT尺寸：%s\n', num2str(size(V)));
fprintf('掩码尺寸：%s\n', num2str(size(rtMask)));

如果非零像素数为0，说明轮廓索引错误或空间匹配失败；尺寸不一致则需检查referenceInfo的尺寸设置。

4. 调整可视化方式增强辨识度 #

默认volshow可能难以区分掩码，可尝试叠加显示或调整阈值：

% 叠加掩码到CT图像上（红色半透明）
volshow(V);
hold on;
volshow(rtMask, 'DisplayStyle', 'Volume', 'Opacity', 0.5, 'Color', [1 0 0]);
hold off;

% 单独显示掩码时设置阈值
volshow(rtMask, 'Threshold', [0.5 1]); % 仅显示掩码区域

5. 排查RT与CT的切片位置匹配 #

RT轮廓的切片位置需与CT切片位置对应，否则掩码会生成在错误层位：

% 获取所有CT切片的位置并排序
ctSlicePositions = [];
ctFiles = dir('DicomCT/*.dcm');
for i = 1:length(ctFiles)
    info = dicominfo(fullfile('DicomCT', ctFiles(i).name));
    ctSlicePositions = [ctSlicePositions info.SlicePosition];
end
ctSlicePositions = sort(ctSlicePositions);

% 获取RT轮廓覆盖的切片位置
rtSlicePositions = unique([rtContours.ContourData(:,3)]);

% 检查匹配的切片数量
fprintf('CT切片数：%d，RT覆盖切片数：%d，匹配切片数：%d\n',...
    length(ctSlicePositions), length(rtSlicePositions),...
    length(intersect(ctSlicePositions, rtSlicePositions)));

若匹配数为0，说明RT与CT的空间坐标系不匹配，需进一步检查DICOM中的ImageOrientationPatient等参数。

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