使用饱和，混合多个卫星图像

Question

我要混合多个卫星图像。但是发生了一个错误。我学习了satpy文档上的示例。

下面是代码，netcdf文件在这里：https://drive.google.com/drive/folders/1zp6EBVfjuh41LDRRZo4PJoeGGGn13AKy?usp=sharing

from glob import glob
from satpy import Scene, MultiScene, DataQuery
from satpy.utils import debug_on
debug_on()


areaid = 'worldeqc3km70'

eumetsat = glob('E:/Global/combine_test/MSG4-SEVI-MSG15-0100-NA-20210801000010.306000000Z-20210801001259-4774254.nat')
goes17 = glob('E:/Global/combine_test/OR_ABI-L1b-RadF-M6C13_G17_s20212130000319_e20212130009396_c20212130009445.nc')
gk2a = glob('E:/Global/combine_test/gk2a_ami_le1b_ir105_fd020ge_202108010000.nc')
goes17_scene = Scene(reader="abi_l1b", filenames=goes17)
eumetsat_scene =  Scene(reader="seviri_l1b_native", filenames=eumetsat)
gk2a_scene = Scene(reader="ami_l1b", filenames=gk2a)

goes17_scene.load(["C13"])
eumetsat_scene.load(['IR_108'])
gk2a_scene.load(["IR105"])

mscn = MultiScene([goes17_scene, eumetsat_scene, gk2a_scene])

#groups = {DataQuery(name='IR_group', wavelength=(9.8, 10.8, 11.8)): ['C13', 'IR105', 'IR_108']}
groups = {DataQuery(name="IR_group", wavelength=(10, 11, 12)): ['C13', 'IR_108', 'IR105']}

mscn.group(groups)
print(mscn.loaded_dataset_ids)



resampled = mscn.resample(areaid, reduce_data=False)
blended = resampled.blend()

blended.save_datasets(filename='./test_{area}.png'.format(area=areaid))

错误消息: RuntimeError:没有生成或无法加载所请求的数据集。请求的复合输入可能需要有匹配的尺寸(例如。通过重采样)。

Answer 1

正如注释中提到的，这是一个已知的bug，希望在未来几周内修复。请跟随第2089期获取更多信息。

短期的解决办法是制定自己的“混合”方法，以你期望的方式处理事情：

from satpy.multiscene import stack

def my_blend(mscn, common_datasets, blend_function=stack):
    new_scn = Scene()
    for ds_id in common_datasets:
        datasets = [scn[ds_id] for scn in mscn.scenes if ds_id in scn]
        new_scn[ds_id] = blend_function(datasets)

    return new_scn

blended = my_blend(resampled, ["ir_group"])