如何用插值为NIFTI格式MRI图像添加切片？（Matlab报错咨询）

一、修复interp3报错的核心原因 #

你当前代码的问题有两个：

• 直接将文件名'subject.nii'传入size()和interp3()，这是错误的——必须先读取NIfTI文件的数值数据数组，而非文件名。
• 原始MRI数据通常是uint8/uint16等整数类型，而interp3要求输入的样本值必须是double或single类型，这就是报错的直接原因。

二、修复后的interp3实现代码 #

% 1. 读取NIfTI数据，获取图像数组和头文件信息
nii = niftiinfo('subject.nii');
img_data = niftiread(nii);
[rows, cols, slices] = size(img_data); % 此时slices=25

% 2. 将数据转换为interp3支持的类型（double）
img_data_double = double(img_data);

% 3. 定义原始网格坐标
[x1, y1, z1] = meshgrid(1:cols, 1:rows, 1:slices);

% 4. 定义目标网格：z轴生成30个均匀分布的采样点（从1到25）
[x2, y2, z2] = meshgrid(1:cols, 1:rows, linspace(1, slices, 30));

% 5. 执行3D线性插值，超出范围补0
img_resampled = interp3(x1, y1, z1, img_data_double, x2, y2, z2, 'linear', 0);

% 6. 将数据转回原始类型（可选，保持和原图像一致的位深）
img_resampled = cast(img_resampled, class(img_data));

% 7. 保存为新的NIfTI文件
niftiwrite(img_resampled, 'subject_30slices.nii', nii);

三、更简便的替代方案：使用imresize #

imresize支持3D图像的缩放，只需指定z轴的目标尺寸，代码更简洁：

% 读取数据
nii = niftiinfo('subject.nii');
img_data = niftiread(nii);

% 按z轴缩放至30切片，使用线性插值（默认也是线性）
img_resampled = imresize(img_data, [size(img_data,1), size(img_data,2), 30], 'linear');

% 保存结果
img_resampled = cast(img_resampled, class(img_data));
niftiwrite(img_resampled, 'subject_30slices_imresize.nii', nii);

关键注意事项 #

• 插值后建议将数据转回原图像的类型（如uint16），避免存储冗余数据或不符合医疗工具的格式要求。
• 如果需要在末尾/开头添加切片而非全局均匀插值，可以调整z2的采样点：比如在原25层后添加5层，可将linspace(1, slices+5, 30)，但这种方式会在末尾产生外插数据，建议优先使用全局均匀插值保证图像连续性。

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