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

保存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读取时出现了如下错误。