image = sitk.ReadImage(image_path)
image_array = sitk.GetArrayFromImage(image)
seg = sitk.ReadImage(seg_path)
seg_array = sitk.GetArrayFromImage(seg)
# 输出图像数据和元信息检验是否读取成功
print("图像数据类型:", image_array.dtype)
print("图像形状:", image.GetSize())
print("像素间距:", image.GetSpacing())
print("原点坐标:", image.GetOrigin())
print("方向矩阵:", image.GetDirection())
# 绘制切片查看图像是否读取成功 在同一个画布上画两个主题的图像
plt.subplot(1, 2, 1)
plt.imshow(image_array[34, :, :], cmap="gray")
plt.title("Image")
plt.subplot(1, 2, 2)
plt.imshow(seg_array[34, :, :], cmap="gray")
plt.title("Segmentation")
plt.show()
nifti格式数据读取及验证
最新推荐文章于 2024-07-24 23:49:50 发布