Technical documentation

SIAC_S2

SIAC.SIAC_S2.SIAC_S2(s2_t, send_back=False, mcd43='/home/docs/MCD43/', vrt_dir='/home/docs/MCD43_VRT/', aoi=None, global_dem=None, cams_dir=None, jasmin=False)[source]
if not os.path.exists(file_path + ‘/emus/’):

os.mkdir(file_path + ‘/emus/’)

if len(glob(file_path + ‘/emus/’ + ‘isotropic_MSI_emulators_*_x?p_S2?.pkl’)) < 12:

url = ‘http://www2.geog.ucl.ac.uk/~ucfafyi/emus/’ req = requests.get(url) to_down = [] for line in req.text.split():

if ‘MSI’ in line:

fname = line.split(‘”’)[1].split(‘<’)[0] if ‘MSI’ in fname:

to_down.append([fname, url])

f = lambda fname_url: downloader(fname_url[0], fname_url[1], file_path + ‘/emus/’) parmap(f, to_down)

SIAC_L8

SIAC.SIAC_L8.SIAC_L8(l8_dir, send_back=False, mcd43='/home/docs/MCD43/', vrt_dir='/home/docs/MCD43_VRT/', aoi=None, global_dem=None, cams_dir=None, jasmin=False)[source]

Atmospheric parameters Inversion

class SIAC.the_aerosol.solve_aerosol(sensor_sat, toa_bands, band_wv, band_index, view_angles, sun_angles, obs_time, cloud_mask, gamma=10.0, a_z_order=1, pixel_res=None, aoi=None, aot_prior=None, wv_prior=None, o3_prior=None, aot_unc=None, wv_unc=None, o3_unc=None, log_file=None, ref_scale=0.0001, ref_off=0, ang_scale=0.01, ele_scale=0.001, prior_scale=[1.0, 0.1, 46.698, 1.0, 1.0, 1.0], emus_dir='SIAC/emus/', mcd43_dir='~/DATA/Multiply/MCD43/', global_dem='/vsicurl/http://www2.geog.ucl.ac.uk/~ucfafyi/eles/global_dem.vrt', cams_dir='/vsicurl/http://www2.geog.ucl.ac.uk/~ucfafyi/cams/', spec_m_dir='SIAC/spectral_mapping', aero_res=1000)[source]

Bases: object

Prepareing modis data to be able to pass into atmo_cor for the retrieval of atmospheric parameters.

Atmospheric parameters Solver

class SIAC.atmo_solver.solving_atmo_paras(boa, toa, sza, vza, saa, vaa, aot_prior, tcwv_prior, tco3_prior, elevation, aot_unc, tcwv_unc, tco3_unc, boa_unc, Hx, Hy, mask, full_res, aero_res, emulators, band_indexs, band_wavelength, pix_res=10.0, gamma=0.5, alpha=- 1.6, subsample=1, subsample_start=0, log_file=None)[source]

Bases: object

A simple implementation of dark dense vegitation method for the restieval of prior aot.

Atmospheric correction

class SIAC.the_correction.atmospheric_correction(sensor_sat, toa_bands, band_index, view_angles, sun_angles, aoi=None, aot=None, tcwv=None, tco3=None, aot_unc=None, tcwv_unc=None, tco3_unc=None, cloud_mask=None, obs_time=None, log_file=None, a_z_order=1, ref_scale=0.0001, ref_off=0, ang_scale=0.01, ele_scale=0.001, atmo_scale=[1.0, 1.0, 1.0, 1.0, 1.0, 1.0], global_dem='/vsicurl/http://www2.geog.ucl.ac.uk/~ucfafyi/eles/global_dem.vrt', cams_dir='/vsicurl/http://www2.geog.ucl.ac.uk/~ucfafyi/cams/', emus_dir='SIAC/emus/', cams_scale=[1.0, 0.1, 46.698, 1.0, 1.0, 1.0], block_size=600, rgb=[None, None, None])[source]

Bases: object

A class doing the atmospheric coprrection with the input of TOA reflectance angles, elevation and emulators of 6S from TOA to surface reflectance.