[GRASS-user] Issues related to lack of Landsat 8 Metadata usage in i.landsat.toar run
iat bd
imtiazataher at gmail.com
Mon Nov 30 07:03:22 PST 2015
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20151130/60b2fe69/attachment.html>
More information about the grass-user
mailing list