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 PIL和CV对 图片的读取,显示,裁剪,保存实现方法
今天小编就为大家分享一篇python PIL和CV对 图片的读取,显示,裁剪,保存实现方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2019-08-08
这篇文章主要介绍了django foreignkey(外键)的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2019-07-07
这篇文章主要介绍了Python合并多个Excel数据的方法也就是说将多个excel中的数据合并到另一个表中，本文通过实例代码相结合的形式给大家介绍的非常详细，具有一定的参考借鉴价值，需要的朋友参考下吧2018-07-07
这篇文章主要介绍了Python计算程序运行时间的方法,分别记录起始时间与结束时间,计算两者之间的差值来获得程序的运行时间,需要的朋友可以参考下2014-12-12
这篇文章主要为大家详细介绍了python实现微信打飞机游戏，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-03-03
这篇文章主要介绍了python单元测试框架pytest的使用示例，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2020-10-10
这篇文章主要为大家详细介绍了python扫描线填充算法，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-02-02
pytorch 实现cross entropy损失函数计算方式
今天小编就为大家分享一篇pytorch 实现cross entropy损失函数计算方式，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2020-01-01
python通过exifread模块获得图片exif信息的方法
这篇文章主要介绍了python通过exifread模块获得图片exif信息的方法,实例分析了Python操作exifread模块的技巧,需要的朋友可以参考下2015-03-03
浏览器引擎是用来处理、渲染和显示网页内容的核心组件,其主要任务是将用户输入的URL所代表的网页资源加载并呈现出来,通常包括HTML、CSS、JavaScript以及各种多媒体内容,本文给大家介绍了如何基于python实现一个简单的浏览器引擎,需要的朋友可以参考下2024-10-10
最新评论