python 将多张 dicom 图像融合为一张 dicom 图像

可以使用Python的Pydicom和Numpy库来将多张DICOM图像融合为一张DICOM图像。以下是一个简单的示例代码：

import numpy as np
import pydicom

# 读取多张DICOM图像
ds_list = []
for file_path in file_path_list:
    ds = pydicom.dcmread(file_path)
    ds_list.append(ds)

# 将多张DICOM图像融合为一张
merged_ds = pydicom.dataset.Dataset()
merged_ds.update(ds_list[0]) # 使用第一张图像的元数据

pixel_array_list = [ds.pixel_array for ds in ds_list]
merged_pixel_array = np.dstack(pixel_array_list)

merged_ds.PixelData = merged_pixel_array.tobytes()
merged_ds.Rows, merged_ds.Columns, _ = merged_pixel_array.shape

# 保存融合后的DICOM图像
pydicom.filewriter.dcmwrite(output_file_path, merged_ds, write_like_original=False)

其中，`file_path_list`是多张DICOM图像的文件路径列表，`output_file_path`是融合后的DICOM图像的保存路径。在代码中，我们首先使用Pydicom读取多张DICOM图像，然后使用Numpy库将它们融合为一张图像。最后，我们使用Pydicom将融合后的图像保存为一个新的DICOM文件。