在Python中读取 `.nii` 格式的医学图像文件，通常使用 `nibabel` 库。该库提供了对NIFTI格式的完整支持，并且可以轻松访问图像的元数据和像素数据。 ### 1. 安装 `nibabel` 库 如果你尚未安装 `nibabel`，可以使用以下命令进行安装： ```bash pip install nibabel ``` ### 2. 读取 `.nii` 文件 以下是一个完整的示例代码，展示如何使用 `nibabel` 读取 `.nii` 文件并提取图像数据和元数据： ```python import nibabel as nib # 加载 .nii 文件 file_path = 'path/to/your/file.nii' img = nib.load(file_path) # 输出图像对象的基本信息 print(img) # 获取图像的头部信息 header = img.header print("Header information:") print(header) # 获取图像的维度信息 dimensions = img.shape print("Image dimensions (x, y, z):", dimensions) # 获取图像的体素数据（像素值） data = img.get_fdata() print("Data shape:", data.shape) print("Data type:", data.dtype) ``` ### 3. 显示 `.nii` 文件的切片图像 如果你想可视化 `.nii` 文件中的图像，可以使用 `matplotlib` 来显示特定切片： ```python import matplotlib.pyplot as plt import nibabel as nib # 加载 .nii 文件 file_path = 'path/to/your/file.nii' img = nib.load(file_path) # 获取图像数据 data = img.get_fdata() # 显示某一特定切片（例如中间的 Z 轴切片） slice_index = data.shape[2] // 2 slice_data = data[:, :, slice_index] # 显示图像 plt.imshow(slice_data, cmap='gray') plt.title(f'Slice at Z={slice_index}') plt.axis('off') plt.show() ``` ### 4. 使用 `OrthoSlicer3D` 显示三维视图 如果需要查看三维图像的整体结构，可以使用 `nibabel.viewers.OrthoSlicer3D` 来交互式地查看图像的三个正交切片视图： ```python from nibabel.viewers import OrthoSlicer3D # 显示三维图像的正交切片视图 OrthoSlicer3D(img.dataobj).show() ``` ### 5. 处理压缩的 `.nii.gz` 文件 `nibabel` 也支持直接读取压缩的 `.nii.gz` 文件，无需手动解压。只需将文件路径指定为 `.nii.gz` 文件即可，例如： ```python file_path = 'path/to/your/file.nii.gz' img = nib.load(file_path) ``` ### 6. 提取特定帧或时间点的数据 如果 `.nii` 文件包含多维数据（如4D图像），可以通过索引提取特定时间点的图像数据： ```python # 假设 data 是一个 4D 数组 (x, y, z, t) time_point = 0 # 提取第一个时间点 time_data = data[:, :, :, time_point] # 显示该时间点的某一 Z 轴切片 slice_index = time_data.shape[2] // 2 slice_data = time_data[:, :, slice_index] plt.imshow(slice_data, cmap='gray') plt.title(f'Slice at Z={slice_index}, Time={time_point}') plt.axis('off') plt.show() ``` 通过上述方法，可以高效地读取和处理 `.nii` 格式的医学图像文件，并结合 `matplotlib` 进行可视化分析。 ---