python显示和查看nii文件

# encoding=utf8
'''
查看和显示nii文件
'''
 
import matplotlib
matplotlib.use('TkAgg')
 
from matplotlib import pylab as plt
import nibabel as nib
from nibabel import nifti1
from nibabel.viewers import OrthoSlicer3D
 
example_filename = "C:/Users/94496/Desktop/sub-verse402_split-verse251_ct.nii"
 
img = nib.load(example_filename)
print (img)
print (img.header['db_name'])   #输出头信息
 
width,height,queue=img.dataobj.shape
 
OrthoSlicer3D(img.dataobj).show()
 
num = 1
for i in range(0,queue,10):
 
    img_arr = img.dataobj[:,:,i]
    plt.subplot(5,4,num)
    plt.imshow(img_arr,cmap='gray')
    num +=1
 
 
plt.show()