From 074cf294e566c7f543042467c186d7895d36ade1 Mon Sep 17 00:00:00 2001 From: Paolo Veglio <paolo.veglio@ssec.wisc.edu> Date: Tue, 7 Feb 2023 22:51:34 +0000 Subject: [PATCH] more spectral tests implemented --- tests/test_spectral_tests.py | 60 ++++++++++++++++++++++++++++++++++++ 1 file changed, 60 insertions(+) diff --git a/tests/test_spectral_tests.py b/tests/test_spectral_tests.py index 814f539..44cfeff 100644 --- a/tests/test_spectral_tests.py +++ b/tests/test_spectral_tests.py @@ -119,6 +119,14 @@ def test_11_12um_diff(data, thresholds, ref_confidence): assert np.allclose(cmin, ref_confidence.diff11_12um_confidence.values) +def test_bt_difference_11_4um_test_ocean(): + pass + + +def test_bt_difference_11_4um_test_land(): + pass + + def test_variability_11um_test(data, thresholds, ref_confidence): cmin = np.ones(data.latitude.shape) for scene_name in ['Polar_Day_Ocean', 'Polar_Night_Ocean']: @@ -129,3 +137,55 @@ def test_variability_11um_test(data, thresholds, ref_confidence): cmin, bit = SceneType.variability_11um_test('M15', cmin) assert np.allclose(cmin, ref_confidence.var11um_confidence.values) + + +def test_oceanic_stratus_11_4um_test(data, thresholds, ref_confidence): + cmin = np.ones(data.latitude.shape) + for scene_name in ['Land_Day', 'Land_Day_Coast', 'Land_Day_Desert', 'Land_Day_Desert_Coast', + 'Ocean_Day', 'Ocean_Night', 'Polar_Day_Land', 'Polar_Day_Coast', + 'Polar_Day_Desert', 'Polar_Day_Desert_Coast', 'Polar_Day_Ocean', + 'Polar_Night_Ocean']: + SceneType = tst.CloudTests(data=data, + scene_name=scene_name, + thresholds=thresholds) + + cmin, bit = SceneType.oceanic_stratus_11_4um_test('M15-M12', cmin) + + assert np.allclose(cmin, ref_confidence.oceanic_stratus_11_4um_confidence.values) + + +def test_nir_reflectance_test(data, thresholds, ref_confidence): + cmin = np.ones(data.latitude.shape) + for scene_name in ['Ocean_Day', 'Polar_Day_Ocean']: + SceneType = tst.CloudTests(data=data, + scene_name=scene_name, + thresholds=thresholds) + + cmin, bit = SceneType.nir_reflectance_test('M07', cmin) + + assert np.allclose(cmin, ref_confidence.NIR_reflectance_confidence.values) + + +def test_vis_nir_ratio_test(data, thresholds, ref_confidence): + cmin = np.ones(data.latitude.shape) + for scene_name in ['Ocean_Day', 'Polar_Day_Ocean']: + SceneType = tst.CloudTests(data=data, + scene_name=scene_name, + thresholds=thresholds) + + cmin, bit = SceneType.vis_nir_ratio_test('M07-M05ratio', cmin) + + assert np.allclose(cmin, ref_confidence.Vis_NIR_ratio_confidence.values) + + +def test_16_21um_reflectance_test(data, thresholds, ref_confidence): + cmin = np.ones(data.latitude.shape) + for scene_name in ['Ocean_Day', 'Polar_Day_Ocean']: + SceneType = tst.CloudTests(data=data, + scene_name=scene_name, + thresholds=thresholds) + + cmin, bit = SceneType.test_16_21um_reflectance('M10', cmin) + + assert np.allclose(cmin, ref_confidence.refl_16_21um_confidence.values) + -- GitLab