批量提取高光谱图像平均反射率的代码

以下是一种可能的Python代码样例，它可以根据给定的高光谱图像，提取所有波段的平均反射率，并返回每个波段的平均反射率向量：

import numpy as np
import spectral # 需要安装spectral包，可使用pip安装

def get_average_reflectance(image):
    """
    获取高光谱图像的平均反射率向量
    :param image: 高光谱图像，numpy数组，shape为(height, width, bands)
    :return: 平均反射率向量，numpy数组，shape为(bands,)
    """
    # 提取图像的高度、宽度和波段数
    height, width, bands = image.shape
    
    # 将图像重新排列为二维数组，每个样本表示一个像素的所有波段反射率
    pixels = image.reshape((height * width, bands))
    
    # 计算每个波段的平均反射率
    average_reflectance = np.mean(pixels, axis=0)
    
    return average_reflectance

使用示例：

# 加载高光谱图像
image = spectral.open_image('path/to/image.hdr')

# 提取平均反射率向量
average_reflectance = get_average_reflectance(image.load())

# 打印平均反射率向量
print(average_reflectance)