python - 通过 Nibabel 加载 Nifti 并使用形状函数

Question

我有一个 nifti 文件1.nii.gz

现在，我从来没有处理过 nifti 文件。

所以，只要用这个软件打开它，我就意识到 nii.gz 是一种包含3 个二维图片数组的容器。事实上，如果我滚动鼠标，我可以看到在图片中标记为 1 的“方向”的 448 张 2d 图片，“方向”2 的 448 张 2d 图片和“方向”3 的 25 张 2d 图片。

在此之后，我打开了 shell 并尝试将这个 nii.gz 与 Nibabel 库一起使用

import nibabel as nib
img = nib.load(1.nii.gz)

但是，如果我输入

img.shape

结果我得到 (448,448,25)，所以这个 .nii.gz 似乎是一个 3d 矩阵，而不是一个包含 3 个 2d 图片数组的容器。你能给我解释一下吗？

Answer 1

Nifti是一种医学图像格式，用于存储图像和伴随数据，图像通常是灰度的，并且它们被视为切片，每个切片具有不同的身体横截面。

它们将所有切片存储在同一个数组中，有时它们会在不同时间获取切片，因此有时它们会向数组添加第四维。

因此，要显示图像或操作它们，您可以对它们进行切片并查看其中的图像。

在您的情况下，数据的形状 (448,448,25) 表明：

有25张图像（切片），尺寸为 448 x 448

import nibabel as nib
import matplotlib.pyplot as plt

# Change the path to your path
path = 'path to img.nii.gz'
Nifti_img  = nib.load(path)
nii_data = my_img.get_fdata()
nii_aff  = my_img.affine
nii_hdr  = my_img.header
print(nii_aff ,'\n',nii_hdr)
print(nii_data.shape)
if(len(nii_data.shape)==3):
   for slice_Number in range(nii_data.shape[2]):
       plt.imshow(nii_data[:,:,slice_Number ])
       plt.show()
if(len(nii_data.shape)==4):
   for frame in range(nii_data.shape[3]):
       for slice_Number in range(nii_data.shape[2]):
           plt.imshow(nii_data[:,:,slice_Number,frame])
           plt.show()