【python】SimpleITK 和 Nibabel 读取医学图像 nii 数据(2D显示)

SimpleITK 和 Nibabel 的区别:

SimpleITK 加载数据是channel_first,即(155,240,240);

Nibabel 是 channel_last,即(240,240,155),其中155是图像通道数,也就是155张图像,可以把nii看成二维图像,也可以看成三维。

But  nibabel加载出来的图像被旋转了90度,横过来了,不知道什么原因~亟待解决~

  • SimpleITK
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[60], cmap = 'gray')
    io.show()

path = 'F:/my_data/t1ce.nii.gz' #数据所在路径
data = read_img(path)
show_img(data)

  • Nibabel
import nibabel as nib
import matplotlib.pyplot as plt

def read_data(path):
    image_data = nib.load(path).get_data()
    return image_data

def show_img(ori_img):
    plt.imshow(ori_img[:,:,85], cmap = 'gray') #channel_last
    plt.show()

path = 'F:/my_data/t1ce.nii.gz'
data = read_data(path)
show_img(data)

欢迎批斗~~~

评论 68
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值