对地震道分别进行30o、45o、60o常相位旋转的代码？

我可以回答这个问题。对于地震道进行常相位旋转，可以使用下面的代码：

import numpy as np

def rotate_phase(tr, angle):
    """
    对地震道进行常相位旋转
    :param tr: 地震道
    :param angle: 旋转角度，单位为度
    :return: 旋转后的地震道
    """
    # 将角度转换为弧度
    angle = np.radians(angle)
    # 构造旋转矩阵
    rot_matrix = np.array([[np.cos(angle), np.sin(angle)],
                           [-np.sin(angle), np.cos(angle)]])
    # 将地震道转换为矩阵
    data_matrix = np.array([tr.data, tr.data])
    # 进行矩阵乘法，得到旋转后的矩阵
    rotated_matrix = np.dot(rot_matrix, data_matrix)
    # 将矩阵转换为地震道
    rotated_tr = tr.copy()
    rotated_tr.data = rotated_matrix[0]
    return rotated_tr

使用这个函数，可以对地震道进行不同角度的常相位旋转。例如，对于一个地震道 tr，可以使用下面的代码将其旋转 30 度：

rotated_tr = rotate_phase(tr, 30)