Added some output validation after AC.

Signed-off-by: Daniel Scheffler <danschef@gfz-potsdam.de>

enpt/options/config.py

@@ -93,7 +93,7 @@ config_for_testing_dlr = dict(
     n_lines_to_append=50,
     disable_progress_bars=True,
     is_dummy_dataformat=False,
-    enable_ac=False,
+    enable_ac=True,
     enable_ice_retrieval=False,
     CPUs=32,
     ortho_resampAlg='gauss',

enpt/processors/atmospheric_correction/atmospheric_correction.py

@@ -93,6 +93,11 @@ class AtmosphericCorrector(object):
         enmap_l2a_vnir, enmap_l2a_swir, cwv_model, cwc_model, ice_model, toa_model, se, scem, srem = \
             sicor_ac_enmap(enmap_l1b=enmap_ImageL1, options=options, logger=enmap_ImageL1.logger)
 
+        # validate results
+        for detectordata, detectorname in zip([enmap_l2a_vnir, enmap_l2a_swir]):
+            if np.mean(detectordata[:, :, 0]) == 0 or np.std(detectordata[:, :, 0]) == 0:
+                enmap_ImageL1.logger.warning('The atmospheric correction returned empty %s bands!' % detectorname)
+
         # join results
         enmap_ImageL1.logger.info('Joining results of atmospheric correction.')