写入多张dicom图片，对图片进行刚性旋转并保持tag不变，用python解决

可以使用Python中的pydicom库来对DICOM图像进行读取、编辑和写入操作。下面是一段示例代码，首先读取多个DICOM图像，进行旋转并保存为新的DICOM文件：

import pydicom
import numpy as np

def rotate_dicom(dicom_file_path, angle):
    # Load DICOM file
    ds = pydicom.dcmread(dicom_file_path)

    # Get the image data as a numpy array
    img = ds.pixel_array

    # Rotate the image using numpy
    img_rotated = np.rot90(img, k=angle//90)

    # Update the image pixel data in the DICOM dataset
    ds.PixelData = img_rotated.tobytes()

    # Adjust the image orientation and position information
    ds.ImageOrientationPatient = ds.ImageOrientationPatient[-2:] + ds.ImageOrientationPatient[:2]
    ds.ImagePositionPatient = ds.ImagePositionPatient[-1::-1]

    # Update the DICOM tag for the rotation angle
    ds[0x0020, 0x0037] = ''.join(['[', str(np.cos(np.radians(angle))), '\\', str(-np.sin(np.radians(angle))), '\\0\\',
                                  str(np.sin(np.radians(angle))), '\\', str(np.cos(np.radians(angle))), '\\0]'])

    # Save the rotated DICOM file
    ds.save_as(dicom_file_path.replace('.dcm', '_rotated.dcm'))

# Example usage
dicom_files = ['image1.dcm', 'image2.dcm', 'image3.dcm']
for file in dicom_files:
    rotate_dicom(file, 90)

这段代码中，我定义了一个rotate_dicom函数，它接受一个DICOM文件路径和旋转角度作为输入。首先，它使用pydicom库中的dcmread函数加载DICOM文件并获取像素数组。然后，它使用numpy库旋转图像数组，并更新DICOM数据集中的像素数据。接下来，它使用DICOM标签更新图像的定位和方向信息以及旋转角度。最后，它使用save_as函数将旋转后的图像保存为新的DICOM文件。在示例中，函数被循环调用以对多个DICOM文件进行操作。