Matlab中niftiwrite保存后MRI图像方向异常如何修正？

NIfTI图像的空间方向由affine变换矩阵和元数据中的Orientation字段共同控制，单纯使用permute/rot/flip调整数组存储维度，不会同步更新空间定位信息，这是导致保存后方向异常的核心原因。以下是具体解决方案：

步骤1：读取图像时完整获取元数据 #

使用niftiread时必须同时获取图像数据和元信息，避免丢失空间定位参数：

[img3D, info3D] = niftiread('your_3d_image.nii');
[binary_img, info_binary] = niftiread('your_binary_image.nii');

步骤2：确保两幅图像空间维度匹配 #

合并前必须验证两幅图像的尺寸完全一致，否则会破坏空间结构：

assert(all(size(img3D) == size(binary_img)), '两幅图像尺寸不匹配，请先配准');

如果尺寸不同，需用imresize3或imregister3将二值图像配准到3D图像的空间网格上。

步骤3：合并图像并同步元数据 #

根据合并需求选择对应方式，同时复用正确的元数据：

• 若将二值图像作为新通道生成4D图像：

  merged_img = cat(4, img3D, binary_img);
  merged_info = info3D;
  merged_info.Dimensions(4) = 2; % 更新第4维度的通道数
  

• 若将二值图像作为掩码与3D图像融合：

  merged_img = img3D .* binary_img;
  merged_info = info3D; % 直接复用3D图像的元数据
  

步骤4：传入完整元数据保存图像 #

必须将包含affine矩阵的元数据传入niftiwrite，避免生成默认affine导致方向错误：

niftiwrite(merged_img, 'merged_result.nii', merged_info);

常见误区修正 #

如果确实需要调整数组维度，必须同步更新affine矩阵。例如交换X/Y轴时：

img_permuted = permute(img3D, [2 1 3]);
% 同步交换affine矩阵对应行列
affine_permuted = info3D.Affine([2 1 3 4], :);
info_permuted = info3D;
info_permuted.Affine = affine_permuted;
niftiwrite(img_permuted, 'adjusted_image.nii', info_permuted);

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