1. 程式人生 > 程式設計 >讀取nii或nii.gz檔案中的資訊即輸出影象操作

讀取nii或nii.gz檔案中的資訊即輸出影象操作

讀取nii或者nii.gz檔案中的資訊,並且輸出影象。

import matplotlib
from matplotlib import pylab as plt
import nibabel as nib
from nibabel.viewers import OrthoSlicer3D
file = '' #你的nii或者nii.gz檔案路徑
img = nib.load(file) 
 
print(img)
print(img.header['db_name']) #輸出nii的標頭檔案
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()

補充知識:SimpleITK讀取醫學影象 .nii 資料(2D顯示)

【環境】win10 + python3.6 + SimpleITK

nii檔案是NIFTI格式的檔案,出現的原因是原來一種影象格式是ANALYZE 7.5 format,但是這個影象格式缺少一些資訊,比如沒有方向資訊,病人的左右方位等,如果需要包括額外的資訊,就需要一個額外的檔案,比如ANALYZE7.5就需要一對<.hdr,.img>檔案來儲存影象的完整資訊。

因此,解決這個問題Data Format Working Group (DFWG) 將影象格式完整的定義為NIFTI(Neuroimaging Informatics Technology Initiative)格式

import SimpleITK as sitk
import skimage.io as io

def read_img(path):
 img = sitk.ReadImage(path)
 data = sitk.GetArrayFromImage(img)
 return data
#顯示一個系列圖
def show_img(data):
 for i in range(data.shape[0]):
  io.imshow(data[i,:],cmap = 'gray')
  print(i)
  io.show()
 
#單張顯示
def show_img(ori_img):
 io.imshow(ori_img[100],cmap = 'gray')
 io.show()

#window下的資料夾路徑 
path = 'D:\\datasets\\Naso_GTV\\1\\data.nii.gz'
data = read_img(path)
show_img(data)

讀取nii或nii.gz檔案中的資訊即輸出影象操作

img = sitk.ReadImage(path)
#檢視圖片深度
print(img.GetDepth())
#144 共144張圖
#檢視Size
print(img.GetSize())
#(512,512,144) 畫素:512*512, 144張圖片

更多的函式自己去發現

讀取nii或nii.gz檔案中的資訊即輸出影象操作

以上這篇讀取nii或nii.gz檔案中的資訊即輸出影象操作就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支援我們。