python使用nibabel和sitk读取保存nii.gz文件实例

更新时间：2020年07月01日 11:00:50 作者：小鹿学长

这篇文章主要介绍了python使用nibabel和sitk读取保存nii.gz文件实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧

nii.gz格式是医学图像常用的压缩格式，python中可用nibabel和sitk来读取保存。

使用nibabel

由于使用nibabel图像会旋转90度，所以读取保存的时候还得保存映射信息，3维图像格式为（z, y, x）

读取nii.gz文件

img = nib.load('xxxxx.nii.gz')
img_affine = img.affine
img = img.get_data()

保存nii.gz文件

nib.Nifti1Image(img,img_affine).to_filename('xxxxx.nii.gz')

使用sitk

使用sitk读取nii时，读取出来的还是图片格式，可以使用他自带的函数进行处理，不过速度比较慢，建议使用GetArrayFromImage转换成numpy格式再处理，3维图像格式为（x, y, z）

读取nii.gz文件

img = sitk.ReadImage('xxxxx.nii.gz')

img = sitk.GetArrayFromImage(img)

保存nii.gz文件

out = sitk.GetImageFromArray(img)

sitk.WriteImage(out,'xxxxx.nii.gz')

在numpy数组和nibabel或sitk中相互转换时，要注意数据的格式，一般保存为int或uint类型。比如输入nii为16位有符号整型时，我们可能需要转换成0~255灰度图，可用如下代码：

img = sitk.ReadImage('xxxxx.nii.gz')
img = sitk.Cast(sitk.RescaleIntensity(img),sitk.sitkUInt8)
img = sitk.GetArrayFromImage(img)

补充知识：SimpleITK保存Nii文件与错误处理方式

Reason:

把处理好的分割结果保存为nii文件，用ITKsnap读取时出现了如下错误。

SimpleITK读取和保存Nii文件

1. 读取

import SimpleITK as sitk
filename = './xxx.nii'
ct = sitk.ReadImage(filename)
ct_array = sitk.GetArrayFromImage(ct)

origin =ct.GetOrigin()
direction = ct.GetDirection()
space = ct.GetSpacing()

2. 保存

savedImg = sitk.GetImageFromArray(ct_array)
savedImg.SetOrigin(origin)
savedImg.SetDirection(direction)
savedImg.SetSpacing(space)
sitk.WriteImage(savedImg, saved_name)

Note：被保存的ct_array数组一定是ndarray，float类型的才能被ITKsnap正确读取，如果是int类型的，就会出现上图中的错误。

以上这篇python使用nibabel和sitk读取保存nii.gz文件实例就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持脚本之家。

您可能感兴趣的文章:

python人工智能tensorflow函数tensorboard使用方法
这篇文章主要为大家介绍了python人工智能tensorflow函数tensorboard使用方法，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2022-05-05
python实现简单聊天功能
这篇文章主要为大家详细介绍了python实现简单聊天功能，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2021-07-07
python3 使用OpenCV计算滑块拼图验证码缺口位置(场景示例)
这篇文章主要介绍了python3 使用OpenCV计算滑块拼图验证码缺口位置,本文通过场景示例给大家详细介绍，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2021-08-08
python 利用jieba.analyse进行关键词提取
这篇文章主要介绍了python 利用jieba.analyse进行关键词提取的方法，帮助大家更好的利用python，感兴趣的朋友可以了解下
2020-12-12
Python数据分析之缺失值检测与处理详解
在实际的数据处理中，缺失值是普遍存在的，如何使用 Python 检测和处理缺失值，就是本文要讲的主要内容。感兴趣的同学可以关注一下
2021-12-12
python PIL中ImageFilter模块图片滤波处理和使用方法
这篇文章主要介绍PIL中ImageFilter模块几种图片滤波处理和使用方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
2023-11-11
Python assert语句的简单使用示例
这篇文章主要给大家介绍了关于Python assert语句的简单使用，文中通过示例代码介绍的非常详细，对大家学习或者使用Python具有一定的参考学习价值，需要的朋友们下面来一起学习学习吧
2019-07-07
python中redis的安装和使用
本文给大家介绍的是在Python中安装和使用redis数据库的方法以及简单示例，有需要的小伙伴可以参考下
2016-12-12
Python Requests安装与简单运用
requests是python的一个HTTP客户端库，跟urllib，urllib2类似，那为什么要用requests而不用urllib2呢？带着这个问题来一起学习本教程吧
2016-04-04
浅析Python自带性能强悍的标准库itertools
itertools是python内置的模块，使用简单且功能强大。这篇文章就主要介绍了通过itertools实现可迭代对象的无限迭代、有限迭代和排列组合。感兴趣的同学可以关注一下
2021-12-12