MATLAB读取并显示多幅MRI图像时出现维度不匹配错误的原因排查

这个Subscripted assignment dimension mismatch错误的核心原因很明确：你硬编码预分配了256×256×1×20的数组，但实际通过dicomread读取的MRI图像，其维度和你预设的不匹配。

具体可能的情况： #

• 你的DICOM图像并不是256×256分辨率的（比如常见的MRI图像可能是512×512）；
• 部分图像可能包含额外的通道（虽然MRI大多是单通道，但少数情况下可能有衍生的多通道数据）；
• 图像的维度顺序和你预设的不一致（比如有些DICOM数据可能返回[行数, 列数, 切片数]的三维数组，而你只预留了单切片的位置）。

修复方案：动态匹配图像维度，避免硬编码 #

不要直接预设图像尺寸，而是先读取第一幅图像获取真实维度，再基于这个维度预分配数组。如果存在少数尺寸不一致的图像，还可以选择性地调整尺寸（需要Image Processing Toolbox支持）。

修改后的代码如下：

% 先读取第一幅图像，获取真实的图像维度
first_img = dicomread('brain_001.dcm');
img_dims = size(first_img);

% 基于真实维度预分配数组，最后一维为图像数量（20）
X = repmat(int16(0), [img_dims, 20]);

% 循环读取所有图像
for p = 1:20
    filename = sprintf('brain_%03d.dcm', p);
    current_img = dicomread(filename);
    
    % 可选：如果当前图像维度和基准维度不一致，调整尺寸（按需选择）
    if ~isequal(size(current_img), img_dims)
        % 使用imresize统一尺寸，需要Image Processing Toolbox
        current_img = imresize(current_img, img_dims);
        % 如果你没有工具箱，也可以手动裁剪/填充到目标尺寸
    end
    
    % 将图像赋值到预分配的数组中
    X(:,:,:,p) = current_img;
end

% 显示图像堆栈
montage(X, [])

调试小技巧 #

如果想快速定位问题，可以在循环中加入一句disp(size(current_img))，打印每幅图像的维度，这样就能直观看到哪幅图像的尺寸和你预设的256×256不匹配。

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