Multiscene blend does not work well.

[jinkooY]

Describe the bug
I want to stack (or merge) multiple satellite images. But when I tried Blend(), it only produced each satellite image, not a blended image.

netcdf file can download from here : https://drive.google.com/drive/folders/1zp6EBVfjuh41LDRRZo4PJoeGGGn13AKy?usp=sharing

To Reproduce

from glob import glob
import os
import satpy
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']}
mscn.group(groups)


resampled = mscn.resample(areaid, reduce_data=False)
blended = resampled.blend()
blended.save_datasets(filename="./png/{name}.png",writer='simple_image')

Expected behavior
A clear and concise description of what you expected to happen.

Actual results
Text output of actual results or error messages including full tracebacks if applicable.

Screenshots
If applicable, add screenshots to help explain your problem.

Environment Info:

• OS: [ Windows]
• Satpy Version: [0.35]

Additional context
Add any other context about the problem here.

[djhoese]

I see you closed this issue. If it was a mistake that you created it or you fixed the problem yourself then 👍. If this is still a problem comment here and we can reopen the issue.