图像文件加载与保存2（包括自然图像，灰度图像，JP2图像）

1.JP2文件打开为数组

import imageio
image_arr=imageio.imreade("./file.jp2")

2.使用nib/np.load加载图像

##np 加载数据
import numpy as np
import nibabel as nib
data_arr=np.load('./file')
#TODO nib 加载、保存数据
#加载
img_nii=nib.load(fliepath+name)
affine=img_nii.affine#读取图像信息：方向等
img_arr = img_nii.get_fdata()#加载为数组
#保存
####赋予原图的图像相关信息，类似下文中sitk赋予Spacing
new_image = nib.Nifti1Image(img_arr, affine)###把affine赋给新图像
#new_image = nib.Nifti1Image(data_arr, np.eye(4))###np.eye(4)指指定某个方向，对角为1，其余为0，4维数组
nib.save(new_image, 'nifti.nii.gz')

3.使用sitk加载、保存数据

import SimpleITK as sitk
data= sitk.ReadImage(labelimage)####data:order=3label:order=0#####
data_np = sitk.GetArrayFromImage(data)  ########转化为数组##########
# 获取图像的信息：spacing等
new_data.SetDirection(data.GetDirection())  ###########图像方向不变
new_data.SetOrigin(data.GetOrigin())  ###########图像原点不变
new_data.SetSpacing((data.GetSpacing()[0] ,data.GetSpacing()[1] data.GetSpacing()[2]))#####层间距变为()
## TODO 保存数据
sitk.WriteImage(new_data,r'.\output\newdatas\\'+'data\\'+dataimage.split('N')[-1]+'.nii')#保存格式（保存文件，路径+文件名+文件格式）

4.调试中显示2D图像

import matplotlib.pyplot as plt
plt.imshow()
plt.show()