def main():
"""Main"""
handler = logging.StreamHandler(sys.stderr)
handler.setLevel(logging.DEBUG)
formatter = logging.Formatter(fmt=_DEFAULT_LOG_FORMAT,
datefmt=_DEFAULT_TIME_FORMAT)
handler.setFormatter(formatter)
logging.getLogger('').addHandler(handler)
logging.getLogger('').setLevel(logging.DEBUG)
modis = RelativeSpectralResponse('EOS-Aqua', 'modis')
modis.read(channel='20', scale=0.001)
solar_irr = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM,
dlambda=0.005)
solar_irr.read()
# Calculate the solar-flux:
sflux = solar_irr.solar_flux_over_band(modis.rsr)
print("Solar flux over Band: ", sflux)
from pyspectral.nir_reflectance import Calculator
#refl37 = Calculator(modis.rsr, solar_flux=sflux)
refl37 = Calculator(modis.rsr)
sunz = 80.
tb3 = 290
tb4 = 282
refl = refl37.reflectance_from_tbs(sunz, tb3, tb4)
print refl
refl37.make_tb2rad_lut('./modis_EOS-Aqua_band20_tb2rad_lut.npz')
def _get_solarflux(self):
"""Derive the in-band solar flux from rsr over the Near IR band (3.7 or 3.9 microns)."""
solar_spectrum = \
SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM,
dlambda=0.0005,
wavespace=self.wavespace)
self.solar_flux = solar_spectrum.inband_solarflux(self.rsr[self.bandname])
def test_solar_flux(self):
"""Calculate the solar-flux"""
self.solar_irr = SolarIrradianceSpectrum(
TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=0.005)
# rsr function (se above) is given in micronsm therefore the scale
# factor is 1.0 and not 1e+6 (default)!
sflux = self.solar_irr.inband_solarflux(self.rsr, scale=1.0)
self.assertAlmostEqual(sflux, 2.0029277645144234)
def test_read(self):
"""Test that solar irradiance spctrum"""
self.solar_irr = \
SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM,
dlambda=0.005)
self.assertTrue(os.path.exists(self.solar_irr.filename))
self.assertEqual(self.solar_irr.wavelength.shape[0], 1697)
self.assertEqual(self.solar_irr.irradiance.shape[0], 1697)
def srf_exatmospheric_irradiance(srf, dlambda_nm=0.5, start_nm=200.0, end_nm=2000.0):
"""
Compute the exatmospheric solar irradiance of some sensor (band) at 1a.u.
:param srf: SRF
`dlambda_nm`, `start_nm`, `end_nm`: interpolation parameters for the solar spectrum
:return: band irradiance at 1au, in [W/(m^2 nm)]
"""
srr = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=dlambda_nm/1000)
srr.interpolate(ival_wavelength=(start_nm/1000, end_nm/1000))
lambdas = srr.ipol_wavelength * 1000
values = srr.ipol_irradiance / 1000.0
response = srf(lambdas)
avg = np.dot(response, values) / np.sum(response)
return avg
def exatmospheric_irradiance(srf, dlambda_nm=0.5, start_nm=200.0, end_nm=2000.0):
"""
Compute the exatmospheric solar irradiance of some sensor (band) at 1a.u.
:param srf: SRF
`dlambda_nm`, `start_nm`, `end_nm`: interpolation parameters for the solar spectrum
:return: band irradiance at 1au, in [W/(m^2 um)]
"""
warnings.warn(DeprecationWarning(
'`exatmospheric_irradiance` is scheduled to be removed, '
'use `srf_exatmospheric_irradiance` (nm based units) instead.'))
srr = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=dlambda_nm/1000)
srr.interpolate(ival_wavelength=(start_nm/1000, end_nm/1000))
lambdas = srr.ipol_wavelength * 1000
response = srf(lambdas)
avg = np.dot(response, srr.ipol_irradiance) / np.sum(response)
return avg
def __init__(self, wave1=200*u.nm, wave2=1200*u.nm, dlambda=1*u.nm,
i=75*u.deg, d=1.5):
self.wave1 = wave1
self.wave2 = wave2
self.dlambda = dlambda
self.i = i # incidence angle
self.d = d # Mars distance in AU (scaling the solar flux)
sol = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM)
sol.interpolate(dlambda=dlambda.to(u.micron).value,
ival_wavelength=(wave1.to(u.micron).value,
wave2.to(u.micron).value))
self.waves = (sol.ipol_wavelength*u.micron).to(u.nm)
self.E_w = (sol.ipol_irradiance*self.E_w_unit_in).to(self.E_w_unit_out)
self.read_reflectance()
self.read_QE()
for k, v in dic.items():
setattr(self, k, v)
def setUp(self):
"""Set up."""
self.solar_irr = SolarIrradianceSpectrum(
TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=0.005)
self.rsr = TEST_RSR
print("*** plot day microphysics RGB for ", str(tslot))
#glbd = GeostationaryFactory.create_scene("meteosat", "09", "seviri", tslot)
glbd = GeostationaryFactory.create_scene("Meteosat-9", "", "seviri", tslot)
print("... load sat data")
glbd.load(['VIS006','VIS008','IR_016','IR_039','IR_108','IR_134'], area_extent=europe.area_extent)
#area="EuropeCanaryS95"
area="EuroMercator" # blitzortung projection
local_data = glbd.project(area, precompute=True)
print("... read responce functions")
from pyspectral.near_infrared_reflectance import Calculator
from pyspectral.solar import (SolarIrradianceSpectrum, TOTAL_IRRADIANCE_SPECTRUM_2000ASTM)
solar_irr = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=0.0005)
#from pyspectral.seviri_rsr import load
#seviri = load()
#rsr = {'wavelength': seviri['IR3.9']['wavelength'],
# 'response': seviri['IR3.9']['met10']['95']}
#sflux = solar_irr.solar_flux_over_band(rsr)
solar_irr = SolarIrradianceSpectrum(TOTAL_IRRADIANCE_SPECTRUM_2000ASTM, dlambda=0.0005)
from pyspectral.rsr_reader import RelativeSpectralResponse
seviri = RelativeSpectralResponse('Meteosat-10', 'seviri')
sflux = solar_irr.inband_solarflux(seviri.rsr['IR3.9'])
ch39 = local_data['IR_039']
ch11 = local_data['IR_108']
ch13 = local_data['IR_134']
