Hi Markus, I have tried i.landsat.toar on different builds of WinGRASS 6.4.4svn, the latest being revision 58054-745 dated 20-Oct-2013. The problem persist.
I did the same with Linux (ubuntu) based GRASS6.4.4svn for the same Landsat scene, same metadata file (LE71710622000137SGS00). The command output (verbose) and r.info print actual values, read from the metafile or computed by i.landsat.toar, NOT ZEROS. A quick comparison - apart from the zeros, the output of i.landsat.toar (WinGRASS) looks the same as the output of i.landsat.toar (Linux based GRASS), at least for this scene. I have compared the --verbose command output and results of r.info for band 1, 5 and 61. Echoing Daniel, I think radiances/reflectances computed by i.landsat.toar in WinGRASS are correct even though r.info and command output print zeros for most of the variables. For the thermal bands the output data range is min=nan and max=nan if using the new metadata file, which was also discussed by Nikos in a different post. I think this depends on the version of GRASS being run. Please find attached some information from r.info and command output which might help to verify/refute my observations or point to the cause/remedy of this problem. Best regards, Beatrice. ________________________________________ Please try the updated version in GRASS 6.4.svn (thanks to E. Jorge Tizado). The GRASS 7 update is hopefully forthcoming, too. Markus
(Tue Oct 22 20:39:27 2013) g.list rast ---------------------------------------------- raster files available in mapset <Zone36>: le71710622000137sgs00_b1 le71710622000137sgs00_b4 le71710622000137sgs00_b62 le71710622000137sgs00_b2 le71710622000137sgs00_b5 le71710622000137sgs00_b7 le71710622000137sgs00_b3 le71710622000137sgs00_b61 le71710622000137sgs00_b8 ---------------------------------------------- (Tue Oct 22 20:39:27 2013) Command finished (0 sec) (Tue Oct 22 20:47:53 2013) r.info D4B1 +----------------------------------------------------------------------------+ | Layer: D4B1 Date: Tue Oct 22 20:41:36 2013 | | Mapset: Zone36 Login of Creator: Beatrice | | Location: Zone36 | | DataBase: C:\Beatrice\grassdata | | Title: ( D4B1 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: DCELL | | Data Units: unitless Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = 0.000623984264880047 max = 0.395979314429105 | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Reflectance of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (0.0000 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ -0.0000000 | | Sun elevation (and azimuth) ........... 0.00000 (0.00000) | | Digital number (DN) range ............. 0 to 0 | | Calibration constants (Lmin to Lmax) .. +0.00000 to +0.00000 | | DN to Radiance (gain and bias) ........ -0.00000 and +0.00000 | | Mean solar irradiance (ESUN) .......... -0.000 | | Radiance to Reflectance (divide by) ... -0.00000 | | | | Dark object (1000 pixels) DN = ........ 8 | | Mode in reflectance histogram ......... -0.00000 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "D4B" metfile="D:\Reprojected\2000\LE71710622000137SGS00\LE717106220\ | | 00137SGS00_MTL.txt" method="dos4" percent=0.01 pixel=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:47:54 2013) Command finished (0 sec) (Tue Oct 22 20:48:01 2013) r.info D4B5 +----------------------------------------------------------------------------+ | Layer: D4B5 Date: Tue Oct 22 20:43:28 2013 | | Mapset: Zone36 Login of Creator: Beatrice | | Location: Zone36 | | DataBase: C:\Beatrice\grassdata | | Title: ( D4B5 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: DCELL | | Data Units: unitless Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = 0 max = 0.553567612354257 | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Reflectance of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (0.0000 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ -0.0000000 | | Sun elevation (and azimuth) ........... 0.00000 (0.00000) | | Digital number (DN) range ............. 0 to 0 | | Calibration constants (Lmin to Lmax) .. -44836016437496836852630521603 | | 7 DN to Radiance (gain and bias) ........ -0.00000 and +0.00000 | | DN to Radiance (gain and bias) ........ -0.00000 and +0.00000 | | Mean solar irradiance (ESUN) .......... -0.000 | | Radiance to Reflectance (divide by) ... -0.00000 | | | | Dark object (1000 pixels) DN = ........ 9 | | Mode in reflectance histogram ......... -0.00000 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "D4B" metfile="D:\Reprojected\2000\LE71710622000137SGS00\LE717106220\ | | 00137SGS00_MTL.txt" method="dos4" percent=0.01 pixel=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:48:01 2013) Command finished (0 sec) (Tue Oct 22 20:48:07 2013) r.info D4B61 +----------------------------------------------------------------------------+ | Layer: D4B61 Date: Tue Oct 22 20:43:49 2013 | | Mapset: Zone36 Login of Creator: Beatrice | | Location: Zone36 | | DataBase: C:\Beatrice\grassdata | | Title: ( D4B61 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: DCELL | | Data Units: Kelvin Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = nan max = nan | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Temperature of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (0.0000 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ -0.0000000 | | Sun elevation (and azimuth) ........... 0.00000 (0.00000) | | Digital number (DN) range ............. 0 to 0 | | Calibration constants (Lmin to Lmax) .. -0.00000 to +0.00000 | | DN to Radiance (gain and bias) ........ -0.00000 and +0.00000 | | Temperature (K1 and K2) ............... 0.000 and 0.000 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "D4B" metfile="D:\Reprojected\2000\LE71710622000137SGS00\LE717106220\ | | 00137SGS00_MTL.txt" method="dos4" percent=0.01 pixel=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:48:07 2013) Command finished (0 sec)
(Tue Oct 22 20:40:27 2013) i.landsat.toar --verbose input_prefix=le71710622000137sgs00_b output_prefix=D4B metfile=D:\Reprojected\2000\LE71710622000137SGS00\LE71710622000137SGS00_MTL.txt method=dos4 Metada file is MTL file: new format RADIANCE & QUANTIZE from data of the metadata file Calculating dark pixel of <le71710622000137sgs00_b1>... ... DN = 08 [1673] : mode 63 [2886796], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b2>... ... DN = 09 [2212] : mode 56 [1606384], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b3>... ... DN = 23 [4156] : mode 55 [824147], exclusing DN > 229 ... TAUv = 0.95742, TAUz = 0.94839, Edown = 15.91089 Calculating dark pixel of <le71710622000137sgs00_b4>... ... DN = 09 [28910] : mode 11 [1699431], exclusing DN > 229 ... TAUv = 0.99984, TAUz = 0.99981, Edown = 0.04038 Calculating dark pixel of <le71710622000137sgs00_b5>... ... DN = 09 [13555] : mode 12 [1147605], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b7>... ... DN = 07 [2910] : mode 11 [1160764], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b8>... ... DN = 04 [2131] : mode 48 [1448951], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 LANDSAT: 7 SENSOR: ETM+ ACQUISITION DATE 2000-05-16 [production date 2013-01-18] earth-sun distance = 0.00000000 solar elevation angle = 0.00000000 Atmospheric correction: dos4 percent of solar irradiance in path radiance = 0.0000 ------------------- BAND 1 (code 1) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): 0.00000 to 0.00000 at-surface radiance = -0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): -0.00000 at-surface reflectance = radiance / -0.00000 the darkness DN with a least 1000 pixels is 8 the DN mode is 63 ------------------- BAND 2 (code 2) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): 0.00000 to 0.00000 at-surface radiance = 0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): -0.00000 at-surface reflectance = radiance / 0.00000 the darkness DN with a least 1000 pixels is 9 the DN mode is 56 ------------------- BAND 3 (code 3) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): -89658965690580456893483080226 372761896147515345354517380636056544393035363933305120614632 335831616473153088978045301666730651970076865416148992201999 954106797305386792565208023318320675143163555234209355804493 382474158013613963587732388823257558500751796899613154738299 245029984813637308400411049125535881285819421131323272818148 129589701836696195670037259148820901979067301010258255666604 638154157976805767114513600170929736319539546898271009988984 473324350816484348130658220517398941516662160777539356244298 973957598320873236830859814125349440968517979291387951582482 105964018515198607985862521153979174409643538279965104124325 676995862386090038589652089662036915925230910533150860808857 727821919022988556248372378673457667695968188057112835971398 410687642959179901665848833424622511605195036102882869638247 707746732651297939823383589787393064991330786427633516445450 417364891040060384688667154441148392292042577419428723717856 6321533949125733709336270751167143147474050258707442302976.0 0000 to 0.00000 at-surface radiance = 0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): -0.00000 at-surface reflectance = radiance / 0.00000 the darkness DN with a least 1000 pixels is 23 the DN mode is 55 ------------------- BAND 4 (code 4) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): 0.00000 to 0.00000 at-surface radiance = 0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): 0.00000 at-surface reflectance = radiance / 0.00000 the darkness DN with a least 1000 pixels is 9 the DN mode is 11 ------------------- BAND 5 (code 5) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): -44836016437496836852630521603 537902907990735186433048149787840711899756923193114412670906 379378072219904673704641626959523600039037900704296428553321 1741432707900020638731957174544212620458246275072.00000 to 0.00000 at-surface radiance = -0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): 0.00000 at-surface reflectance = radiance / -0.00000 the darkness DN with a least 1000 pixels is 9 the DN mode is 12 ------------------- BAND 6 thermal (code 61) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): -0.00000 to 0.00000 at-surface radiance = -0.00000000 * DN + -0.00000 at-sensor temperature = 0.00000 / log[(0.00000 / radiance) + 1.0] ------------------- BAND 6 thermal (code 62) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): -0.00000 to 0.00000 at-surface radiance = -0.00000000 * DN + -0.00000 at-sensor temperature = 0.00000 / log[(0.00000 / radiance) + 1.0] ------------------- BAND 7 (code 7) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): -0.00000 to 0.00000 at-surface radiance = 0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): 0.00000 at-surface reflectance = radiance / 0.00000 the darkness DN with a least 1000 pixels is 7 the DN mode is 11 ------------------- BAND 8 (code 8) calibrated digital number (DN): 0.0 to 0.0 calibration constants (L): 0.00000 to 0.00000 at-surface radiance = 0.00000000 * DN + -0.00000 mean solar exoatmospheric irradiance (ESUN): 0.00000 at-surface reflectance = radiance / 0.00000 the darkness DN with a least 1000 pixels is 4 the DN mode is 48 ------------------- Calculating... Writing reflectance of <le71710622000137sgs00_b1> to <D4B1>... Writing reflectance of <le71710622000137sgs00_b2> to <D4B2>... Writing reflectance of <le71710622000137sgs00_b3> to <D4B3>... Writing reflectance of <le71710622000137sgs00_b4> to <D4B4>... Writing reflectance of <le71710622000137sgs00_b5> to <D4B5>... Writing temperature of <le71710622000137sgs00_b61> to <D4B61>... Writing temperature of <le71710622000137sgs00_b62> to <D4B62>... Writing reflectance of <le71710622000137sgs00_b7> to <D4B7>... Writing reflectance of <le71710622000137sgs00_b8> to <D4B8>... (Tue Oct 22 20:45:07 2013) Command finished (4 min 39 sec)
(Tue Oct 22 19:57:17 2013) g.list rast ---------------------------------------------- raster files available in mapset <Zone36>: le71710622000137sgs00_b1 le71710622000137sgs00_b4 le71710622000137sgs00_b62 le71710622000137sgs00_b2 le71710622000137sgs00_b5 le71710622000137sgs00_b7 le71710622000137sgs00_b3 le71710622000137sgs00_b61 le71710622000137sgs00_b8 ---------------------------------------------- (Tue Oct 22 19:57:17 2013) Command finished (0 sec) (Tue Oct 22 20:03:16 2013) r.info Dos4B1 +----------------------------------------------------------------------------+ | Layer: Dos4B1 Date: Tue Oct 22 19:59:02 2013 | | Mapset: Zone36 Login of Creator: root | | Location: Zone36-64 | | DataBase: /home/beatrice/grassdata | | Title: ( Dos4B1 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 255 | | Data Type: DCELL | | Data Units: unitless Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = 0.000623984264880051 max = 0.395979314429105 | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Reflectance of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (7.9039 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ 1.0110998 | | Sun elevation (and azimuth) ........... 54.37667 (50.51917) | | Digital number (DN) range ............. 1 to 255 | | Calibration constants (Lmin to Lmax) .. -6.20000 to +191.60000 | | DN to Radiance (gain and bias) ........ +0.77874 and -1.24652 | | Mean solar irradiance (ESUN) .......... 1969.000 | | Radiance to Reflectance (divide by) ... +498.33971 | | | | Dark object (1000 pixels) DN = ........ 8 | | Mode in reflectance histogram ......... 0.09595 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "Dos4B" metfile="/media/My Passport/Reprojected/2000/LE7171062200013\ | | 7SGS00/LE71710622000137SGS00_MTL.txt" method="dos4" percent=0.01 pix\ | | el=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:03:16 2013) Command finished (0 sec) (Tue Oct 22 20:03:27 2013) r.info Dos4B5 +----------------------------------------------------------------------------+ | Layer: Dos4B5 Date: Tue Oct 22 20:00:30 2013 | | Mapset: Zone36 Login of Creator: root | | Location: Zone36-64 | | DataBase: /home/beatrice/grassdata | | Title: ( Dos4B5 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 255 | | Data Type: DCELL | | Data Units: unitless Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = 0 max = 0.553567612354257 | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Reflectance of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (7.9039 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ 1.0110998 | | Sun elevation (and azimuth) ........... 54.37667 (50.51917) | | Digital number (DN) range ............. 1 to 255 | | Calibration constants (Lmin to Lmax) .. -1.00000 to +31.06000 | | DN to Radiance (gain and bias) ........ +0.12622 and -0.56475 | | Mean solar irradiance (ESUN) .......... 225.700 | | Radiance to Reflectance (divide by) ... +57.12304 | | | | Dark object (1000 pixels) DN = ........ 9 | | Mode in reflectance histogram ......... 0.01663 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "Dos4B" metfile="/media/My Passport/Reprojected/2000/LE7171062200013\ | | 7SGS00/LE71710622000137SGS00_MTL.txt" method="dos4" percent=0.01 pix\ | | el=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:03:28 2013) Command finished (0 sec) (Tue Oct 22 20:03:33 2013) r.info Dos4B61 +----------------------------------------------------------------------------+ | Layer: Dos4B61 Date: Tue Oct 22 20:00:47 2013 | | Mapset: Zone36 Login of Creator: root | | Location: Zone36-64 | | DataBase: /home/beatrice/grassdata | | Title: ( Dos4B61 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 255 | | Data Type: DCELL | | Data Units: Kelvin Vertical datum: | | Rows: 6950 | | Columns: 7790 | | Total Cells: 54140500 | | Projection: UTM (zone -36) | | N: 9784915.00215827 S: 9576385.00215825 Res: 30.00431655 | | E: 519015.00192555 W: 285285.00192663 Res: 30.00385109 | | Range of data: min = -nan max = -nan | | | | Data Description: | | generated by i.landsat.toar | | | | Comments: | | Temperature of Landsat-7 ETM+ (method dos4) | | ----------------------------------------------------------------- | | Acquisition date (and time) ........... 2000-05-16 (7.9039 h) | | Production date ....................... 2013-01-18 | | | | Earth-sun distance (d) ................ 1.0110998 | | Sun elevation (and azimuth) ........... 54.37667 (50.51917) | | Digital number (DN) range ............. 1 to 255 | | Calibration constants (Lmin to Lmax) .. +0.00000 to +17.04000 | | DN to Radiance (gain and bias) ........ +0.06709 and -0.06709 | | Temperature (K1 and K2) ............... 0.000 and 0.000 | | ------------------------------------------------------------------ | | | | i.landsat.toar input_prefix="le71710622000137sgs00_b" output_prefix=\ | | "Dos4B" metfile="/media/My Passport/Reprojected/2000/LE7171062200013\ | | 7SGS00/LE71710622000137SGS00_MTL.txt" method="dos4" percent=0.01 pix\ | | el=1000 rayleigh=0.0 | | | +----------------------------------------------------------------------------+ (Tue Oct 22 20:03:33 2013) Command finished (0 sec)
(Tue Oct 22 19:58:16 2013) i.landsat.toar --verbose input_prefix=le71710622000137sgs00_b output_prefix=Dos4B metfile=/media/My Passport/Reprojected/2000/LE71710622000137SGS00/LE71710622000137SGS00_MTL.txt method=dos4 Metada file is MTL file: new format RADIANCE & QUANTIZE from data of the metadata file Calculating dark pixel of <le71710622000137sgs00_b1>... ... DN = 08 [1673] : mode 63 [2886796], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b2>... ... DN = 09 [2212] : mode 56 [1606384], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b3>... ... DN = 23 [4156] : mode 55 [824147], exclusing DN > 229 ... TAUv = 0.95742, TAUz = 0.94839, Edown = 15.91089 Calculating dark pixel of <le71710622000137sgs00_b4>... ... DN = 09 [28910] : mode 11 [1699431], exclusing DN > 229 ... TAUv = 0.99984, TAUz = 0.99981, Edown = 0.04038 Calculating dark pixel of <le71710622000137sgs00_b5>... ... DN = 09 [13555] : mode 12 [1147605], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b7>... ... DN = 07 [2910] : mode 11 [1160764], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 Calculating dark pixel of <le71710622000137sgs00_b8>... ... DN = 04 [2131] : mode 48 [1448951], exclusing DN > 229 ... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000 LANDSAT: 7 SENSOR: ETM+ ACQUISITION DATE 2000-05-16 [production date 2013-01-18] earth-sun distance = 1.01109982 solar elevation angle = 54.37667123 Atmospheric correction: dos4 percent of solar irradiance in path radiance = 0.0100 ------------------- BAND 1 (code 1) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -6.20000 to 191.60000 at-surface radiance = 0.77874016 * DN + -1.24652 mean solar exoatmospheric irradiance (ESUN): 1969.00000 at-surface reflectance = radiance / 498.33971 the darkness DN with a least 1000 pixels is 8 the DN mode is 63 ------------------- BAND 2 (code 2) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -6.40000 to 196.50000 at-surface radiance = 0.79881890 * DN + -2.53246 mean solar exoatmospheric irradiance (ESUN): 1840.00000 at-surface reflectance = radiance / 465.69074 the darkness DN with a least 1000 pixels is 9 the DN mode is 56 ------------------- BAND 3 (code 3) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -5.00000 to 152.90000 at-surface radiance = 0.62165354 * DN + -10.68625 mean solar exoatmospheric irradiance (ESUN): 1551.00000 at-surface reflectance = radiance / 361.17853 the darkness DN with a least 1000 pixels is 23 the DN mode is 55 ------------------- BAND 4 (code 4) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -5.10000 to 241.10000 at-surface radiance = 0.96929134 * DN + -6.08214 mean solar exoatmospheric irradiance (ESUN): 1044.00000 at-surface reflectance = radiance / 264.14782 the darkness DN with a least 1000 pixels is 9 the DN mode is 11 ------------------- BAND 5 (code 5) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -1.00000 to 31.06000 at-surface radiance = 0.12622047 * DN + -0.56475 mean solar exoatmospheric irradiance (ESUN): 225.70000 at-surface reflectance = radiance / 57.12304 the darkness DN with a least 1000 pixels is 9 the DN mode is 12 ------------------- BAND 6 thermal (code 61) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): 0.00000 to 17.04000 at-surface radiance = 0.06708661 * DN + -0.06709 at-sensor temperature = 0.00000 / log[(0.00000 / radiance) + 1.0] ------------------- BAND 6 thermal (code 62) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): 3.20000 to 12.65000 at-surface radiance = 0.03720472 * DN + 3.16280 at-sensor temperature = 0.00000 / log[(0.00000 / radiance) + 1.0] ------------------- BAND 7 (code 7) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -0.35000 to 10.80000 at-surface radiance = 0.04389764 * DN + -0.09957 mean solar exoatmospheric irradiance (ESUN): 82.07000 at-surface reflectance = radiance / 20.77133 the darkness DN with a least 1000 pixels is 7 the DN mode is 11 ------------------- BAND 8 (code 8) calibrated digital number (DN): 1.0 to 255.0 calibration constants (L): -4.70000 to 243.10000 at-surface radiance = 0.97559055 * DN + -0.44005 mean solar exoatmospheric irradiance (ESUN): 1368.00000 at-surface reflectance = radiance / 346.23094 the darkness DN with a least 1000 pixels is 4 the DN mode is 48 ------------------- Calculating... Writing reflectance of <le71710622000137sgs00_b1> to <Dos4B1>... Writing reflectance of <le71710622000137sgs00_b2> to <Dos4B2>... Writing reflectance of <le71710622000137sgs00_b3> to <Dos4B3>... Writing reflectance of <le71710622000137sgs00_b4> to <Dos4B4>... Writing reflectance of <le71710622000137sgs00_b5> to <Dos4B5>... Writing temperature of <le71710622000137sgs00_b61> to <Dos4B61>... Writing temperature of <le71710622000137sgs00_b62> to <Dos4B62>... Writing reflectance of <le71710622000137sgs00_b7> to <Dos4B7>... Writing reflectance of <le71710622000137sgs00_b8> to <Dos4B8>... (Tue Oct 22 20:01:47 2013) Command finished (3 min 30 sec)
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
