Landsat 9 script for surface area calculation is missing from TMSOS.
Closed this issue · 0 comments
SanchitMinocha commented
Please add landsat9 also in run_sarea_for_res() in run_sarea.py
# Sentinel-2
log.debug(f"Reservoir: {reservoir_name}; Downloading Sentinel-2 data from {start_date} to {end_date}")
sarea_s2(reservoir_name, reservoir_polygon, start_date, end_date, os.path.join(datadir, 's2'))
s2_dfpath = os.path.join(datadir, 's2', reservoir_name+'.csv')
# Landsat-8
log.debug(f"Reservoir: {reservoir_name}; Downloading Landsat-8 data from {start_date} to {end_date}")
sarea_l8(reservoir_name, reservoir_polygon, start_date, end_date, os.path.join(datadir, 'l8'))
l8_dfpath = os.path.join(datadir, 'l8', reservoir_name+'.csv')
# Sentinel-1
log.debug(f"Reservoir: {reservoir_name}; Downloading Sentinel-1 data from {start_date} to {end_date}")
s1_dfpath = sarea_s1(reservoir_name, reservoir_polygon, start_date, end_date, os.path.join(datadir, 'sar'))
s1_dfpath = os.path.join(datadir, 'sar', reservoir_name+'_12d_sar.csv')
tmsos = TMS(reservoir_name, reservoir_area)
result = tmsos.tms_os(l8_dfpath=l8_dfpath, s2_dfpath=s2_dfpath, s1_dfpath=s1_dfpath)
Also add sarea_cli_l9.py script which got missed in some previous commit.