### 计算辐射定标影像中每个波段的平均地表反射率 要计算已辐射定标的影像中每个波段的平均地表反射率，通常需要使用遥感图像处理库（如GDAL或Rasterio）来读取影像数据，并使用NumPy进行数值计算。 以下是一个完整的Python代码示例，用于读取多波段遥感影像并计算每个波段的平均地表反射率： ```python import rasterio import numpy as np def calculate_average_reflectance(image_path): """ 计算已辐射定标影像中每个波段的平均地表反射率。 参数: image_path (str): 辐射定标后的影像文件路径（如GeoTIFF格式）。 返回: list: 每个波段的平均地表反射率列表。 """ with rasterio.open(image_path) as src: num_bands = src.count average_reflectances = [] for band_idx in range(1, num_bands + 1): # 读取单个波段数据 band_data = src.read(band_idx).astype(np.float32) # 掩膜无效值（如0或NaN） valid_data = band_data[(band_data != 0) & (~np.isnan(band_data))] if valid_data.size == 0: average_reflectance = 0.0 else: # 计算平均反射率 average_reflectance = np.mean(valid_data) average_reflectances.append(average_reflectance) return average_reflectances # 示例用法 image_path = "path/to/radiometrically_calibrated_image.tif" avg_reflectances = calculate_average_reflectance(image_path) for i, avg_ref in enumerate(avg_reflectances, start=1): print(f"Band {i} average surface reflectance: {avg_ref:.4f}") ``` #### 代码说明： - **rasterio** 用于读取GeoTIFF等格式的遥感影像。 - **NumPy** 用于处理图像数据并计算平均值。 - 在读取每个波段后，使用掩膜去除无效值（如0或NaN），以避免影响平均值计算。 - 输出每个波段的平均地表反射率。 #### 注意事项： - 输入的影像应为已辐射定标的格式，反射率值通常在0到1之间或以其他比例表示（如0-10000表示0-1）。 - 如果影像使用数字量化值（DN），则需要在计算前将其转换为物理反射率值，通常需要应用增益、偏移等定标参数。 ### 示例扩展：应用增益和偏移进行定标 如果影像数据为数字量化值（Digital Numbers, DN），需要应用增益和偏移转换为地表反射率： ```python # 示例：将DN转换为反射率 def convert_dn_to_reflectance(band_data, gain=0.0001, offset=0): return (band_data * gain) + offset ``` 在读取波段数据后调用此函数，再进行平均值计算。