在醫(yī)學(xué)圖像處理中，我們經(jīng)常使用一種NIFTI格式圖像（.nii文件），現(xiàn)在我們來看看

什么是.nii文件？

該如何讀取.nii文件？

1. NIFTI格式圖像

什么是NIFTI（Neuroimaging Informatics Technology Initiative）格式圖像？
在講解什么是NIFTI格式之前，得先了解一下Analyze格式。Analyze格式儲存的每組數(shù)據(jù)組包含2個(gè)文件，一個(gè)為數(shù)據(jù)文件，其擴(kuò)展名為.img，包含二進(jìn)制的圖像資料；另外一個(gè)為頭文件，擴(kuò)展名為.hdr，包含圖像的元數(shù)據(jù)。在fMRI的早期，Analyze格式最常用的格式，但現(xiàn)在逐漸被NIFTI格式所取代。Analyze格式主要不足就是頭文件不能真正反映元數(shù)據(jù)。
標(biāo)準(zhǔn)NIFTI圖像的擴(kuò)展名是.nii，也包含了頭文件及圖像資料。由于NIFTI格式和Analyze格式的關(guān)系，因此NIFTI格式也可使用獨(dú)立的圖像文件（.img）和頭文件（.hdr）。單獨(dú)的.nii格式文件的優(yōu)勢就是可以使用標(biāo)準(zhǔn)的壓縮軟件（如gzip）進(jìn)行壓縮，而且一些分析軟件包（比如FSL）可以直接讀取和寫入壓縮的.nii文件（擴(kuò)展名為.nii.gz）。

為什么會出現(xiàn)NIFTI格式圖像？
原來的ANALYZE 7.5 format圖像格式缺少一些信息，比如沒有方向信息，病人的左右方位等，如果需要包括額外的信息，就需要一個(gè)額外的文件，比如ANALYZE 7.5就需要一對<.hdr, .img>文件來保存圖像的完整信息。因此，解決這個(gè)問題Data Format Working Group (DFWG) 將圖像格式完整的定義為NIFTI格式。

詳細(xì)了解NIFTI格式請參見：?https://brainder.org/2012/09/23/the-nifti-file-format/

2. 讀取NIFTI格式圖像

2.1 ITK-SNAP

推薦ITK-SNAP下載地址：http://www.itksnap.org/pmwiki/pmwiki.php?n=Downloads.SNAP3

可以看到：itk-snap獲得了分別來自三個(gè)角度的視圖（水平面、矢狀面、冠狀面）

2.2? Python讀取npy文件代碼

import os image = np.load('/content/drive/MyDrive/coursework1-data/2D_projection_R_sub-0.npy', encoding='bytes') print(image.shape) A = image[:, :, 0] B = image[:, :, 1] C = image[:, :, 2] plt.imshow(A) plt.show() plt.imshow(B) plt.show() plt.imshow(C) plt.show()

?2.3? Python讀取pkl文件代碼

file_path='/content/drive/MyDrive/Colab Notebooks/Training_meta.pkl' meta_path=os.path.join(file_path) meta=pd.read_pickle(meta_path) print(meta)

?

?2.4? Python讀取nii文件代碼并保存成png格式圖片

import numpy as np import os #遍歷文件夾 import nibabel as nib #nii格式一般都會用到這個(gè)包 import imageio #轉(zhuǎn)換成圖像 from PIL import Image np.set_printoptions(threshold=np.inf) def nii_to_image(niifile):return 0filepath = './' # 讀取本代碼同個(gè)文件夾下所有的nii格式的文件 filenames = os.listdir(filepath) imgfile = './'slice_trans = []for f in filenames: #開始讀取nii文件s = f[-3:]print(s)if s != '.nii':continues1 = f[:-4]print(s1)imgfile_path = imgfile + s1print("imgfile_path:"+imgfile_path)img_path = os.path.join(filepath, f)img = nib.load(img_path) #讀取niiprint("img:")print(img)img_fdata = img.get_fdata()fname = f.replace('.nii','') #去掉nii的后綴名img_f_path = os.path.join(imgfile, fname)if not os.path.exists(img_f_path):os.mkdir(img_f_path)#創(chuàng)建nii對應(yīng)的圖像的文件夾 # # if not os.path.exists(img_f_path): # os.mkdir(img_f_path) #新建文件夾 # #開始轉(zhuǎn)換為圖像if '.gz' in s1:(x, y, z, _) = img.shapeprint("img2:")print(img.shape)else:(x,y,z) = img.shapeprint("img3:")print(img.shape)for i in range(z): #z是圖像的序列silce = img_fdata[:, :, i] #選擇哪個(gè)方向的切片都可以imageio.imwrite(os.path.join(img_f_path,'{}_mask.png'.format(i)), silce)img = Image.open(os.path.join(img_f_path,'{}_mask.png'.format(i)))img.save(os.path.join(img_f_path,'{}_mask.png'.format(i)))

任何程序錯(cuò)誤，以及技術(shù)疑問或需要解答的，請?zhí)砑幼髡遃X：1755337994

?

?

總結(jié)

以上是生活随笔為你收集整理的全网最详细 Python如何读取NIFTI格式图像（.nii文件）和 .npy格式文件和pkl标签文件内容的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。