使用pydicom从轴位视图提取矢状和冠状切面

import pydicom
import numpy as np
import matplotlib.pyplot as plt

# 加载DICOM文件
ds = pydicom.dcmread("path/to/dicom/file.dcm")

# 提取图像数据
image = ds.pixel_array

# 获取切片厚度和间距
slice_thickness = ds.SliceThickness
slice_spacing = ds.SpacingBetweenSlices

# 计算切片数量和方向
num_slices = ds.NumberOfFrames if "NumberOfFrames" in ds else 1
slice_direction = ds.ImageOrientationPatient[0:3]

# 计算切片间隔向量
slice_spacing_vector = np.multiply(slice_direction, slice_spacing)

# 提取矢状切面
sagittal_slice = image[int(num_slices / 2), :, :]

# 提取冠状切面
coronal_slice = image[:, int(image.shape[1] / 2), :]

# 显示矢状切面
plt.imshow(sagittal_slice, cmap="gray")
plt.title("Sagittal Slice")
plt.show()

# 显示冠状切面
plt.imshow(coronal_slice, cmap="gray")
plt.title("Coronal Slice")
plt.show()