在Python中怎样用NumPy对三维医疗影像矩阵进行重采样

云敏小哥_5366

发布时间：2026-10-29 09:47:32

280人浏览过

来源于php中文网

原创

医疗影像重采样需同步更新体素间距与affine矩阵：先用nibabel读取原始zooms，按(z,y,x)顺序计算各向异性缩放因子，scipy.ndimage.zoom重采样时图像用order=3、标签用order=0，并用新affine=nibabel.affine×diag(zoom_x,zoom_y,zoom_z,1)更新空间坐标系。

在python中怎样用numpy对三维医疗影像矩阵进行重采样

重采样前必须确认原始体素间距和目标体素间距

医疗影像（如CT、MRI的NIfTI或DICOM序列）的重采样不是单纯缩放数组，而是物理空间对齐。NumPy本身不处理空间元数据，affine矩阵或pixdim信息必须从原始文件中读取（例如用nibabel或pydicom），否则重采样结果在空间上是错位的。

常见错误：直接对numpy.ndarray调用scipy.ndimage.zoom或cv2.resize，却忽略Z轴与XY轴体素尺寸通常不同（比如Z轴层厚2.5mm，XY为0.5mm），导致各向异性形变。

用nibabel.load(path).header.get_zooms()获取原始(dx, dy, dz)
目标间距通常由临床需求指定，例如统一到(1.0, 1.0, 1.0) mm 或匹配某参考序列
缩放因子不是固定值，而是逐轴计算：zoom_factors = (orig_dz / target_dz, orig_dy / target_dy, orig_dx / target_dx) —— 注意顺序：NumPy数组索引是[z,y,x]，而NIfTI的zooms是(x,y,z)，需翻转

用`scipy.ndimage.zoom`做各向异性重采样要设`order`和`mode`

scipy.ndimage.zoom是NumPy生态中最常用、可控性强的重采样工具，但它默认order=3（三次样条）在医学图像中可能引入虚假纹理，尤其对二值掩膜（segmentation mask）会破坏边缘锐度。

对原始图像（灰度值连续）：用order=3保持平滑，mode='constant'（边界填0）或'reflect'（更自然）
对分割标签图（整数类别）：必须用order=0（最近邻），否则插值会产生非整数值，导致标签污染
务必传入zoom为3元组，且顺序对应数组维度：zoom=(zoom_z, zoom_y, zoom_x)，别写成(x,y,z)
若内存受限，避免一次性重采样整个3D块；可沿Z轴分块处理，再拼接

import numpy as np
from scipy.ndimage import zoom
<h1>假设img_3d.shape == (64, 256, 256)，原始体素(x,y,z) = (0.8, 0.8, 5.0)，目标(1.0,1.0,1.0)</h1><p>orig_zooms = (0.8, 0.8, 5.0)  # nibabel.get_zooms()返回顺序
target_zooms = (1.0, 1.0, 1.0)
zoom_factors = (orig_zooms[2]/target_zooms[2],  # z轴缩放
orig_zooms[1]/target_zooms[1],  # y轴缩放
orig_zooms[0]/target_zooms[0])  # x轴缩放</p><p>resampled = zoom(img_3d, zoom=zoom_factors, order=3, mode='reflect')</p><div class="aritcle_card flexRow">
                                                        <div class="artcardd flexRow">
                                                                <a class="aritcle_card_img" href="/xiazai/skill9965" title="macOS UI 自动化技能包 —— 用 Python 控制你的 Mac"><img
                                                                                src="https://img.php.cn/upload/skill/001/246/273/6adaf20dbae7c507.jpg" alt="macOS UI 自动化技能包 —— 用 Python 控制你的 Mac" ></a>
                                                                <div class="aritcle_card_info flexColumn">
                                                                        <a href="/xiazai/skill9965" title="macOS UI 自动化技能包 —— 用 Python 控制你的 Mac">macOS UI 自动化技能包 —— 用 Python 控制你的 Mac</a>
                                                                        <p>让 macOS 自动化像说话一样简单 —— 用代码控制你的 Mac，无需学习 AppleScript

系统要求
macOS 14+ (Sonoma)

Python 3.9+

需要开启：辅助功能 + 屏幕录制权限</p>
                                                                </div>
                                                                <a href="/xiazai/skill9965" title="macOS UI 自动化技能包 —— 用 Python 控制你的 Mac" class="aritcle_card_btn flexRow flexcenter"><b></b><span>下载</span> </a>
                                                        </div>
                                                </div><p><span>立即学习</span>“<a href="https://pan.quark.cn/s/00968c3c2c15" style="text-decoration: underline !important; color: blue; font-weight: bolder;" rel="nofollow" target="_blank">Python免费学习笔记（深入）</a>”；</p>

重采样后必须更新空间坐标系（affine）

只改数组尺寸不改affine矩阵，会导致后续用FSL、ITK或SimpleITK加载时定位错误——图像内容被拉伸/压缩，但头文件仍声称它在原空间里。

affine是一个4×4齐次变换矩阵，前三列前三行描述体素尺寸和方向，第四列表示原点偏移
重采样后，新affine = 原affine × 缩放矩阵，其中缩放矩阵为np.diag([zoom_x, zoom_y, zoom_z, 1.0])
注意：如果用了zoom函数，缩放是相对于图像中心做的，但affine更新只需考虑体素尺寸变化，**不涉及平移校正**（除非你显式做了中心对齐）
推荐用nibabel保存时传入新affine，不要手动改header字段