Calculates and compares the radiance of a thermal LANSAT scene to the "ground truth" radiance as measured by a NOAA buoy. Based on work by Frank Padula and Monica Cook at RIT.
If you want to use this code, you should have a basic knowledge of python and/or basic coding. No warranty. Use it on armstrong or related servers for best results. Developed on Fedora x64 by Nathan Dileas (nid4986@g.rit.edu). Copyright RIT 2015-2018
This code essentially has two funtions: calculating the radiance from the landsat image provided, and calculating the corresponding ground truth radiance from outside data, atmospheric (NARR or MERRA-2), NOAA buoy data, and MODTRAN. If atmospheric data or landsat images need to be downloaded, it will take between 5-7 minutes for NARR, and 2-3 for MERRA. Use the file forward_model.py as a convinient command line interface:
$ python forward_model.py -h
usage: forward_model.py [-h] [-a {merra,narr}] [-v] [-b BANDS [BANDS ...]]
scene_id buoy_id
Compute and compare the radiance values of a landsat image to the propogated
radiance of a NOAA buoy, using atmospheric data and MODTRAN.
positional arguments:
scene_id LANDSAT or MODIS scene ID. Examples:
LC08_L1TP_017030_20170703_20170715_01_T1,
MOD021KM.A2011154.1650.006.2014224075807.hdf
buoy_id NOAA Buoy ID. Example: 45012
optional arguments:
-h, --help show this help message and exit
-a {merra,narr}, --atmo {merra,narr}
Choose atmospheric data source, choices:[narr, merra].
-v, --verbose
-b BANDS [BANDS ...], --bands BANDS [BANDS ...]
See INSTALL.txt
- landsat test scene
python forward_model.py LC08_L1TP_017030_20170703_20170715_01_T1 45012
- modis test scene
python forward_model.py MOD021KM.A2011154.1650.006.2014224075807.hdf 45012
- all the downloaded data will be in a directory: buoy_calib/Landsat-Buoy-Calibration/downloaded_data/
- the settings and paths can be changed in: buoy_calib/Landsat-Buoy-Calibration/buoycalib/settings.py
- these websites will help you search to find your own scenes / buoys
- for buoys: http://www.ndbc.noaa.gov/
- for landsat scenes: https://earthexplorer.usgs.gov/, Collection 1 Level 1 - Landsat 8 OLI/TIRS C1 Level 1
- for modis scenes: https://ladsweb.modaps.eosdis.nasa.gov/search/, look for Terra-MODIS, then MOD021KM
- landsat scenes before 2017 do not work currently
- http://scholarworks.rit.edu/theses/2961/ - Padula 08 Thesis
- http://scholarworks.rit.edu/theses/8513/ - Cook 14 Thesis
- tools/to_csv.py: used to compile results quickly and easily.
- tools/generate_atmo_figure.py : generate a figure using information from a already processed scene.
- test/functional/run_all_scenes.bash: run a batch of scenes. Move it to this directory before use.