在SimpleITK中,各术语对应如下:
Width: 宽度,矢状面(Sagittal)
Height: 高度,冠状面(Coronal)
Depth: 深度, 横断面(Axial)
使用GetSize()的输出为:(矢状面Width, 冠状面Height, 横断面Depth),也即原始SimpleITK数据的存储形式
使用GetArrayFromImage()方法后,X轴与Z轴发生了对调,输出形状为:(横断面Depth, 冠状面Height, 矢状面Width)
import SimpleITK as sitk
path = './data/patient001/patient001_frame01.nii.gz'
image = sitk.ReadImage(path)
# 直接获取尺寸
shape_img = image.GetSize()
print(f'GetSize of image: {shape_img}')
# 转换为numpy矩阵后输出矩阵维度
data = sitk.GetArrayFromImage(image)
shape_data = data.shape
print(f'numpy of image: {shape_data}')
'''
输出信息
GetSize of image: (216, 256, 10)
numpy of image: (10, 256, 216)
'''