python光谱角代码

以下是一个使用Python计算光谱角度的简单示例代码：

import numpy as np

def calculate_spectral_angle(band1, band2):
    """
    计算光谱角度
    参数:
        - band1: 第一个波段的像素值数组
        - band2: 第二个波段的像素值数组
    返回值:
        - spectral_angle: 光谱角度值（弧度）
    """
    # 将像素值归一化为单位向量
    band1_norm = band1 / np.linalg.norm(band1)
    band2_norm = band2 / np.linalg.norm(band2)

    # 计算光谱角度
    spectral_angle = np.arccos(np.dot(band1_norm, band2_norm))

    return spectral_angle

# 示例数据
band1 = np.array([0.2, 0.4, 0.6, 0.8])
band2 = np.array([0.3, 0.5, 0.7, 0.9])

# 计算光谱角度
angle = calculate_spectral_angle(band1, band2)
print("光谱角度：", angle)

在这个示例中，我们使用`numpy`库来进行向量计算。首先，我们将两个波段的像素值归一化为单位向量。然后，使用`numpy.dot`函数计算两个向量的点积，并使用`numpy.arccos`函数计算其反余弦值，得到光谱角度。最后，打印出计算得到的光谱角度值。

请注意，示例代码中的数据是简单的一维数组，你可以根据实际情况修改代码以适应多维数组或其他数据结构。