[GRASS-user] Issues related to lack of Landsat 8 Metadata usage in i.landsat.toar run

Date: November 30, 2015

In my application of i.landsat.toar, I am coming across an issue for which there is previous precedence.

http://gis.stackexchange.com/questions/72259/landsat-8-i-landsat-toar-error-grass-7/72789#72789

Computer Basics: Data Used: Landsat 8; Operating System: Windows 7 32 bit; Grass version: 7.0.0.1 (Stable version) I have tried with 7.02.RC1-1 version, landsat 7 data, but there was no difference)

My metadata file records the following values:
SUN_ELEVATION = 49.52284282
EARTH_SUN_DISTANCE = 0.9939432

Also,

RADIANCE_MULT_BAND_1 = 1.2709E-02
RADIANCE_ADD_BAND_1 = -63.54622
Similarly, none of the other band’s radiance parameter is zero

  1. Issue # 01: The code is not picking up the Earth-sun distance or the Solar elevation angle from the metadata file

  2. Issue # 02: As par USGS reference: http://landsat.usgs.gov/Landsat8_Using_Product.php , at-sensor radiance = RADIANCE_MULT_BAND_X * DN + RADIANCE_ADD_BAND_X However, the output is not using the metadata values. Instead using 0. Any explanation?

I am coming across identical situation like the previous case, evident from the attached command output. Will appreciate some guideline for solution (if required), explanation for the circumstance (if this is how it is supposed to be).

Thanks and Best Wishes,

Imtiaz

Command Output:

Mon Nov 30 02:30:07 2015)
i.landsat.toar --verbose input=LC81370442015300LGN00_B output=TOAR_B metfile=E:\LandCover_Classification\LC81370442015300LGN00\LC81370442015300LGN00_MTL.txt sensor=oli8 date=2015-10-27 sun_elevation=49.52284282 product_date=2015-10-27
Metada file is MTL file: new format
RADIANCE & QUANTIZE from MIN_MAX_(RADIANCE|PIXEL_VALUE)
WARNING: ESUN evaluated from REFLECTANCE_MAXIMUM_BAND
WARNING: Overwriting solar elevation of metadata file
LANDSAT: 8 SENSOR: OLI/TIRS
ACQUISITION DATE 2015-10-27 [production date 2015-10-27]
Earth-sun distance = 0.00000000 /* Not using data from metedata file /
Solar elevation angle = 0.00000000 /
Not using data from metedata file */
Atmospheric correction: UNCORRECTED

BAND 1 (code 1)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): -12921358752688410020721955688
115920665908972319060407734713394360286772747253243677217490
460072444863477515873024485329406587544102380733964553773209
470174482838252720913256857159594704835548262628428491099806
824007743633806603939891682505986476500036813488471710794679
812718339464991426221557174273397692968656770569470784372007
509468405212310749439605762528891325540967113062645506423667
906579067093107701745245095315224335459710698106660563473926
87291100325217965687908018786027399961269512990865888378880.
00000 to 0.00000
at-sensor radiance = -0.00000000 * DN + -0.00000
mean solar exoatmospheric irradiance (ESUN): 0.00000
at-sensor reflectance = radiance / -0.00000

BAND 2 (code 2)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): -0.00000 to 0.00000
at-sensor radiance = -0.00000000 * DN + -0.00000
mean solar exoatmospheric irradiance (ESUN): -0.00000
at-sensor reflectance = radiance / -0.00000