import pickle
import numpy as np
import SimpleITK as sitk
img = sitk.ReadImage(r’C:\Users\320133514\Desktop\2-9\needle_newdata206.nii.gz’)
img = sitk.GetArrayFromImage(img)
#print(img)
print(img.shape)
print(type(img))
print(img.ndim)
print(img[0:24, 50:370,0:320])
C = img
print©
print(“img.shape:”,C.shape)
D = np.vstack((C,img))
print(D)
print(“img.shape:”,D.shape)
out = sitk.GetImageFromArray(D)
sitk.WriteImage(out,r’C:\Users\320133514\Desktop\2-9\D48.nii.gz’)
print(“ok”)