I noticed this composite appear in the list of available_composite_names for abi_l1b and ahi_hsd readers in a recent release of SatPy. Unfortunately, I haven't been able to get it to work.
I've downloaded the BlackMarble_2016_3km_geo.tif from here:
However, when I try to run:
scene = Scene(filenames=filenames, reader=reader)
scene.load(['true_color_with_night_ir_hires'])
resampled = scene.resample(
scene.max_area(),
resampler='native',
)
I get the following exception:
Traceback (most recent call last):
File "./geocolor.py", line 55, in <module>
scene.max_area(),
File "/usr/local/lib/python3.8/site-packages/satpy/scene.py", line 284, in max_area
return self._compare_areas(datasets=datasets, compare_func=max)
File "/usr/local/lib/python3.8/site-packages/satpy/scene.py", line 262, in _compare_areas
raise ValueError("Can't compare areas with different "
ValueError: Can't compare areas with different projections.
Any tips on how to get this seemingly useful DayNightComposite to work?