def test_apply_sunearth_corr(self):
"""Test the correction of reflectances with sun-earth distance."""
out_refl = hf.apply_earthsun_distance_correction(self.raw_refl)
np.testing.assert_allclose(out_refl, self.corr_refl)
self.assertTrue(
out_refl.attrs['sun_earth_distance_correction_applied'])
assert isinstance(out_refl.data, da.Array)
def vis_calibrate(self, data, solar_irradiance):
"""Calibrate to reflectance.
This uses the method described in Conversion from radiances to
reflectances for SEVIRI warm channels: https://tinyurl.com/y67zhphm
"""
reflectance = np.pi * data * 100.0 / solar_irradiance
return apply_earthsun_distance_correction(reflectance, self._scan_time)
def vis_calibrate(self, data, solar_irradiance):
"""Calibrate to reflectance.
This uses the method described in Conversion from radiances to
reflectances for SEVIRI warm channels: https://www-cdn.eumetsat.int/files/2020-04/pdf_msg_seviri_rad2refl.pdf
"""
reflectance = np.pi * data * 100.0 / solar_irradiance
return apply_earthsun_distance_correction(reflectance, self._scan_time)