from pyspectral import rayleigh
import numpy as np
sunz = np.array([[50., 60.], [50., 60.]])
satz = np.array([[40., 50.], [40., 50.]])
azidiff = np.array([[160, 160], [20, 20]])
corr = rayleigh.Rayleigh('GOES-16', 'abi')
print corr.get_reflectance(sunz, satz, azidiff, 'ch1')
[[ 7.66551598 11.84851174] [ 12.30975667 17.83044933]]
print corr.get_reflectance(sunz, satz, azidiff, 'ch2')
[[ 2.01927932 3.20415785] [ 3.52058364 5.18308456]]
print corr.get_reflectance(sunz, satz, azidiff, 'ch3')
[[ 0. 0.] [ 0. 0.]]
No handlers could be found for logger "pyspectral.rayleigh"
corr = rayleigh.Rayleigh('GOES-16', 'abi', atmosphere='midlatitude summer', aerosol_type='rural_aerosol')
print corr.get_reflectance(sunz, satz, azidiff, 'ch1')
[[ 11.25806044 19.05767067] [ 14.21482674 20.60287676]]