Python project that handles Landsat data processing

Requires Functions.py from the GDAL_Python3 repository, make sure its in the PYTHONPATH

Uses original *.tar.gz packed Landsat 4,5,7 or 8 files (e.g. from https://earthexplorer.usgs.gov/) in an input folder, creates a subfolder with the Landsat name and:

• extracts bands
• calculates radiation bands
• calculates TOA reflection bands
• calculates surface temperature in Kelvin
• calculates one or all indices (NDVI, EVI, SAVI, MSAVI, NDMI, NBR, NBR2, NDSI) according to the Landsat Surface Reflectance High Level Data Products

The script uses flags to determine if intermediate outputs should be saved as rasters to HDD or removed after the process is finished