更多数学原理小文请关注公众号:未名方略
首先安装
conda install -c simpleitk simpleitk
如果你有dicom文件或npy文件,单层显示如下
dicom可以查体素信息,将体素信息spacing“插入”三维数据,存为.mha文件,那么三维文件才完整。
python 代码
import SimpleITK as itk
import numpy as np
import nibabel as nib
image_path = ‘’
ds = nib.load(image_path)
img = ds.get_fdata()
header = ds.header
pixdim=header['pixdim']
px=round(pixdim[1],2)
py=round(pixdim[2],2)
pz=round(pixdim[3],2)
outpath='.mha'
img=np.uint8(img)
image:itk.Image = itk.GetImageFromArray(img)
spacing = (px,py,pz)
image.SetSpacing(spacing)
save_mha = itk.ImageFileWriter()
save_mha.SetFileName(outpath)
save_mha.Execute(image)