我的目标是将辐射参数转换为来自哨兵三个 slstr 数据的反射率。
我应用的公式是：
反射率 = pi*radiance/solar_radiance/cos(solar_zenith_angle)
数据包含多个 netCDF 格式的文件。它包含 S1-S6 辐射和质量数据以及 S7-S9 BT 数据。它还包含包含太阳天顶角的几何文件。
我正在使用 xarray 来读取这些文件。
例如，我使用 S6 辐射数据来存储辐射： radiance = S6_radiance['S6_radiance_an']
存储太阳辐照度的 S6 质量数据：solar_irradiance = S6_quality['solar_irradiance_an']
cosine_zenith = np.cos(几何['solar_zenith_tn']
当我替换公式中的数据时：反射率 = pi*radiance/solar_irradiance/cosine_zenith
它给了我一个错误。它说：
ValueError：沿维度“行”没有标签的参数无法对齐，因为它们具有不同的维度大小：{2400, 1200}
尺寸为：
radiance.dims = ('rows', 'columns')
radiance.shape = (2400, 3000)
solar_irradiance.dims = ('detectors',)
solar_irradiance.shape = (4,)
cosine_zenith.dims = ('rows', 'columns')
cosine_zenith.shape = (1200, 130)