[GRASS-user] i.landsat.toar comments have only 0 value
Pekka Hurskainen
pekka.hurskainen at gmail.com
Mon Dec 16 03:39:30 PST 2013
Dear Markus,
Thanks for the reply. I have just now tested with Landsat8 OLI scene
(27-08-2013) 168/062. I installed the latest 6.4.4svn (58463) and tried
i.landsat.toar with DOS2 and DOS3 methods. I run it for a smaller subset,
with clouds and image edges masked out. Below is the command output, first
for dos3 and then for dos2. Note that for dos2, I manually entered solar
elev angle, but the output is still 0 ? Everything else I left default
values.
BTW, what is the difference between DOS2 and DOS2b methods?
-------clip------
(Mon Dec 16 12:33:37
2013)
i.landsat.toar --verbose input_prefix=crop_b output_prefix=dos3_3_b
metfile=C:\HY-Data\pjhurska\Kilimanjaro\Landsat\Landsat8_2013-27-08_168_062_OLI\LC81680622013239LGN00_MTL.txt
sensor=oli8 method=dos3 pixel=5
Metada file is MTL file: new format
RADIANCE & QUANTIZE from MIN_MAX_(RADIANCE|PIXEL_VALUE)
WARNING: ESUN evaluated from REFLECTANCE_MAXIMUM_BAND
Calculating dark pixel of <crop_b1>...
... DN = 7840 [5] : mode 9627 [2891], exclusing DN > 58981
... TAUv = 0.78835, TAUz = 0.76206, Edown = 0.00000
Calculating dark pixel of <crop_b2>...
... DN = 7135 [5] : mode 9083 [2409], exclusing DN > 58981
... TAUv = 0.84575, TAUz = 0.82578, Edown = 0.00000
Calculating dark pixel of <crop_b3>...
... DN = 6200 [5] : mode 8799 [1880], exclusing DN > 58981
... TAUv = 0.91432, TAUz = 0.90271, Edown = 0.00000
Calculating dark pixel of <crop_b4>...
... DN = 5785 [6] : mode 6663 [1183], exclusing DN > 58981
... TAUv = 0.95287, TAUz = 0.94633, Edown = 0.00000
Calculating dark pixel of <crop_b5>...
... DN = 6024 [7] : mode 14327 [927], exclusing DN > 58981
... TAUv = 0.98442, TAUz = 0.98222, Edown = 0.00000
Calculating dark pixel of <crop_b6>...
... DN = 5377 [5] : mode 11683 [721], exclusing DN > 58981
... TAUv = 0.99871, TAUz = 0.99852, Edown = 0.00000
Calculating dark pixel of <crop_b7>...
... DN = 5166 [5] : mode 10961 [843], exclusing DN > 58981
... TAUv = 0.99963, TAUz = 0.99958, Edown = 0.00000
Calculating dark pixel of <crop_b8>...
... DN = 5906 [6] : mode 9018 [5640], exclusing DN > 58981
... TAUv = 0.92888, TAUz = 0.91916, Edown = 0.00000
Calculating dark pixel of <crop_b9>...
... DN = 4992 [5] : mode 5039 [133488], exclusing DN > 58981
... TAUv = 0.99757, TAUz = 0.99722, Edown = 0.00000
LANDSAT: 8 SENSOR: OLI/TIRS
ACQUISITION DATE 2013-08-27 [production date 2013-08-27]
Earth-sun distance = 0.00000000
Solar elevation angle = 0.00000000
Atmospheric correction: dos3
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): 968367649887747415141875521950
925908889125765915809137356521975470948435722379319114916517
109681257194579805986272450667354678861757535842162547084653
462145348007752249138970238330548827196232199351249578581263
248711758443327663076681220844948160752094949553961518385593
288738424407009196744807773476057507004619157871083200622000
320486055830136042443031065077914185459331507512811830598511
152020147139378105907351040013136335438675174223655826910682
669147134563183589280887783585976043716156080824364735851118
183768026326159743796155841003601516760621076483451614820141
371685937002372918255194596159954067917411389793697522459155
015399457565770890778803961998816054585345162915220244370597
148522387936504415391544888346434689270315995072771028338997
995727094089743160130590886592801301509071858609333431506640
587793968989583037059469778982726604216980071724204307021521
635177260328502996043398383316279662847376171701723950431400
757812176165083327750696905811530775310780096413788963894683
905077061031539383276262428122121528157514609481673317570684
435219591138433952053577693737258063060125527179863820119388
770600138288962007966562469949923301233492925582540474476621
929036185589173077700588636992204749746454624011817064126344
632245678976560320724068426895929872763359600665707344187778
891947245102403001861675618641623023833022185852578363120467
361109134430306902862293966757899915186924084339575936028017
743399931822767869305933232821423387100836970448169974752735
074832657312708225553213351832437155472265900799102082698841
521263830766549846516292538579807583602895447499728905211008
406860112609801193185282209285068919448526126257634719443776
306035623536657971781547855260622470703165809089946123219856
394296744397591189893603472912483484277176386264413758344436
893155901820405208007799764626033880161122800644723691809874
031293719596978952957939025121263874432512217233471764941837
349859301579149060534601918858766260236011294832113712139455
698498970564770118735817734482120560350401517034380439852381
686780485705510656529346796692015942241400805762289636824034
380012173568901522372612940698148203094590690817052813619697
272631456298510647656983077748357220329595173179533274877623
690664824083139848333834479764874064891277168954263595672454
741560283118661607390700810458152871418035457554976493129153
6739339299897124622932302102528.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 5 pixels is 7840
the DN mode is 9627
-------------------
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 5 pixels is 7135
the DN mode is 9083
-------------------
BAND 3 (code 3)
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 5 pixels is 6200
the DN mode is 8799
-------------------
BAND 4 (code 4)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 183658085783589171116602975780
576110409539519711468540243914721655904198177980762961939111
367605627562671716937127859554390051130688152018870728220925
072842007885913864783619283774492980062937015097153316113130
872437715945381379739519113601595608897104545933411437714766
964222032417753040188158845816964370334214805486222793513041
392415948524253835628431799076337651616691066413747194706092
588365827695512535302789689209508312728090357719245288559169
261388427087636514887083888570441112579969923313176980814385
626616270691354605624826831921938182203368446708429571503396
673727480549490161133959472989086466065342736246984472386035
984295478176953900049240996927707494912239369479071549791691
006705649149921302160640916967785780445041977901000120712458
693774631060343641966598542706661151042560211066330115225274
882082359272986006411181268975110156015073466360964844002614
923902881591929598593433704046771611849396994917616009430564
636178472140317666513101744160796308166097582587535447676030
395599188338919046557681998905510895192924869482845313899438
243168242481424950199879384022305364926654416732220470906145
003035888135318340747367027552221587693565002548093465052373
146810072538472500364217293690318091228705920528403356637468
515206939874589256642932818842704726527146610272815374838100
196715703591772856115604922172749774983079205445635592246176
519699396484278000996109577040304620492915499240374708941194
899957749080944402792983491485854033303690426232983954252982
089745184902031081081879977151639415526477184204943465088423
197754901126349507777941665413693263267204404932798742686446
175425099611417857336481903514883801243132058814981214699218
223830439709555518391135039601061824428651849472616761261741
947438101584718397168382882672948322751691531336445977140261
8847363072.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 5 pixels is 5785
the DN mode is 6663
-------------------
BAND 5 (code 5)
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 5 pixels is 6024
the DN mode is 14327
-------------------
BAND 6 (code 6)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 339445567057160696751644954879
970406447096914339798998362968503888855957504.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 5 pixels is 5377
the DN mode is 11683
-------------------
BAND 7 (code 7)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): -12134434465588505055597499511
824959727712116382507612325036502948239426558532389801597855
534050652218225524835841065680524520096072520033585910067094
864419565702430636692147519347440048706033789010796831449675
609397487722631716840845517015635835853269211738830091674017
125859807797852579371667606442129288215138677555089849122233
496258894611200291018407686794712556602863468629113156241991
764561853870347050462715860062969867544293321888080259438748
975116227397144258398582860683531866628235263927302795102064
876102210094737464694328482051327435929619847514195837359477
453721352442867031673331032270809238607484683630583314906933
393968396774392452257227334055169647686113972735180171551034
094817132579151468969733950921315554583854644271192962176107
852279850198497591348543860021633007244464057151870155597981
968744623129348330236016103587706033254928248063411627431203
939460766333363995335790489660112468340889170867819092108790
845331032250512936792412018170190148253165743211304227984791
418825798114108108620478601673212231026083304899473015984957
797069887035981967325806061173567396266179556816201804747668
19328.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 5 pixels is 5166
the DN mode is 10961
-------------------
BAND 8 (code 8)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 177790655472633694188941980802
047503613813833685602824969586860483988457303292419157292935
380388916793723522118532246597515650465972468098640750969175
323067583642284054472809092057610201799293465748504101117625
584350450442885820883877261635394365497242394779959210685927
483970642556853937908745511045113355423116353373754214638700
488299349599339888965342105384720125168911889987302964445071
690804617411665559673197765683461268733759690639062267875652
078576707772030813055249319922273769990973469604361776775922
848041675652772916924480281074566297232981629176185283191324
907861265119304021792937302365364770562183880430197029735849
971627803188571242704480618796802428300307738635105278630447
639039111579815598101840646206149312815075326886519428941714
704663858921998805414402540381763432286269864844293075921085
556261175910325106160771636615881978320283927409653766398503
610130923687169102643128269846432862140570867543535243138038
049809462751772897752234749314486560174042187499283473193614
10650569707755492771627008.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 5 pixels is 5906
the DN mode is 9018
-------------------
BAND 9 (code 9)
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 5 pixels is 4992
the DN mode is 5039
-------------------
BAND 10 thermal (code 10)
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 11 thermal (code 11)
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]
-------------------
Calculating...
Writing reflectance of <crop_b1> to <dos3_3_b1>...
Writing reflectance of <crop_b2> to <dos3_3_b2>...
Writing reflectance of <crop_b3> to <dos3_3_b3>...
Writing reflectance of <crop_b4> to <dos3_3_b4>...
Writing reflectance of <crop_b5> to <dos3_3_b5>...
Writing reflectance of <crop_b6> to <dos3_3_b6>...
Writing reflectance of <crop_b7> to <dos3_3_b7>...
Writing reflectance of <crop_b8> to <dos3_3_b8>...
Writing reflectance of <crop_b9> to <dos3_3_b9>...
Writing temperature of <crop_b10> to <dos3_3_b10>...
Writing temperature of <crop_b11> to <dos3_3_b11>...
(Mon Dec 16 12:34:44 2013) Command finished (1 min 7
sec)
(Mon Dec 16 12:43:05
2013)
i.landsat.toar --verbose input_prefix=crop_b output_prefix=cost_b
metfile=C:\HY-Data\pjhurska\Kilimanjaro\Landsat\Landsat8_2013-27-08_168_062_OLI\LC81680622013239LGN00_MTL.txt
sensor=oli8 method=dos2 sun_elevation=60.02214600 pixel=5
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
Calculating dark pixel of <crop_b1>...
... DN = 7840 [5] : mode 9627 [2891], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b2>...
... DN = 7135 [5] : mode 9083 [2409], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b3>...
... DN = 6200 [5] : mode 8799 [1880], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b4>...
... DN = 5785 [6] : mode 6663 [1183], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b5>...
... DN = 6024 [7] : mode 14327 [927], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b6>...
... DN = 5377 [5] : mode 11683 [721], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000
Calculating dark pixel of <crop_b7>...
... DN = 5166 [5] : mode 10961 [843], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000
Calculating dark pixel of <crop_b8>...
... DN = 5906 [6] : mode 9018 [5640], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000
Calculating dark pixel of <crop_b9>...
... DN = 4992 [5] : mode 5039 [133488], exclusing DN > 58981
... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000
LANDSAT: 8 SENSOR: OLI/TIRS
ACQUISITION DATE 2013-08-27 [production date 2013-08-27]
Earth-sun distance = 0.00000000
Solar elevation angle = 0.00000000
Atmospheric correction: dos2
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): 968367649887747415141875521950
925908889125765915809137356521975470948435722379319114916517
109681257194579805986272450667354678861757535842162547084653
462145348007752249138970238330548827196232199351249578581263
248711758443327663076681220844948160752094949553961518385593
288738424407009196744807773476057507004619157871083200622000
320486055830136042443031065077914185459331507512811830598511
152020147139378105907351040013136335438675174223655826910682
669147134563183589280887783585976043716156080824364735851118
183768026326159743796155841003601516760621076483451614820141
371685937002372918255194596159954067917411389793697522459155
015399457565770890778803961998816054585345162915220244370597
148522387936504415391544888346434689270315995072771028338997
995727094089743160130590886592801301509071858609333431506640
587793968989583037059469778982726604216980071724204307021521
635177260328502996043398383316279662847376171701723950431400
757812176165083327750696905811530775310780096413788963894683
905077061031539383276262428122121528157514609481673317570684
435219591138433952053577693737258063060125527179863820119388
770600138288962007966562469949923301233492925582540474476621
929036185589173077700588636992204749746454624011817064126344
632245678976560320724068426895929872763359600665707344187778
891947245102403001861675618641623023833022185852578363120467
361109134430306902862293966757899915186924084339575936028017
743399931822767869305933232821423387100836970448169974752735
074832657312708225553213351832437155472265900799102082698841
521263830766549846516292538579807583602895447499728905211008
406860112609801193185282209285068919448526126257634719443776
306035623536657971781547855260622470703165809089946123219856
394296744397591189893603472912483484277176386264413758344436
893155901820405208007799764626033880161122800644723691809874
031293719596978952957939025121263874432512217233471764941837
349859301579149060534601918858766260236011294832113712139455
698498970564770118735817734482120560350401517034380439852381
686780485705510656529346796692015942241400805762289636824034
380012173568901522372612940698148203094590690817052813619697
272631456298510647656983077748357220329595173179533274877623
690664824083139848333834479764874064891277168954263595672454
741560283118661607390700810458152871418035457554976493129153
6739339299897124622932302102528.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 5 pixels is 7840
the DN mode is 9627
-------------------
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 5 pixels is 7135
the DN mode is 9083
-------------------
BAND 3 (code 3)
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 5 pixels is 6200
the DN mode is 8799
-------------------
BAND 4 (code 4)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 183658085783589171116602975780
576110409539519711468540243914721655904198177980762961939111
367605627562671716937127859554390051130688152018870728220925
072842007885913864783619283774492980062937015097153316113130
872437715945381379739519113601595608897104545933411437714766
964222032417753040188158845816964370334214805486222793513041
392415948524253835628431799076337651616691066413747194706092
588365827695512535302789689209508312728090357719245288559169
261388427087636514887083888570441112579969923313176980814385
626616270691354605624826831921938182203368446708429571503396
673727480549490161133959472989086466065342736246984472386035
984295478176953900049240996927707494912239369479071549791691
006705649149921302160640916967785780445041977901000120712458
693774631060343641966598542706661151042560211066330115225274
882082359272986006411181268975110156015073466360964844002614
923902881591929598593433704046771611849396994917616009430564
636178472140317666513101744160796308166097582587535447676030
395599188338919046557681998905510895192924869482845313899438
243168242481424950199879384022305364926654416732220470906145
003035888135318340747367027552221587693565002548093465052373
146810072538472500364217293690318091228705920528403356637468
515206939874589256642932818842704726527146610272815374838100
196715703591772856115604922172749774983079205445635592246176
519699396484278000996109577040304620492915499240374708941194
899957749080944402792983491485854033303690426232983954252982
089745184902031081081879977151639415526477184204943465088423
197754901126349507777941665413693263267204404932798742686446
175425099611417857336481903514883801243132058814981214699218
223830439709555518391135039601061824428651849472616761261741
947438101584718397168382882672948322751691531336445977140261
8847363072.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 5 pixels is 5785
the DN mode is 6663
-------------------
BAND 5 (code 5)
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 5 pixels is 6024
the DN mode is 14327
-------------------
BAND 6 (code 6)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 339445567057160696751644954879
970406447096914339798998362968503888855957504.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 5 pixels is 5377
the DN mode is 11683
-------------------
BAND 7 (code 7)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): -12134434465588505055597499511
824959727712116382507612325036502948239426558532389801597855
534050652218225524835841065680524520096072520033585910067094
864419565702430636692147519347440048706033789010796831449675
609397487722631716840845517015635835853269211738830091674017
125859807797852579371667606442129288215138677555089849122233
496258894611200291018407686794712556602863468629113156241991
764561853870347050462715860062969867544293321888080259438748
975116227397144258398582860683531866628235263927302795102064
876102210094737464694328482051327435929619847514195837359477
453721352442867031673331032270809238607484683630583314906933
393968396774392452257227334055169647686113972735180171551034
094817132579151468969733950921315554583854644271192962176107
852279850198497591348543860021633007244464057151870155597981
968744623129348330236016103587706033254928248063411627431203
939460766333363995335790489660112468340889170867819092108790
845331032250512936792412018170190148253165743211304227984791
418825798114108108620478601673212231026083304899473015984957
797069887035981967325806061173567396266179556816201804747668
19328.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 5 pixels is 5166
the DN mode is 10961
-------------------
BAND 8 (code 8)
calibrated digital number (DN): 0.0 to 0.0
calibration constants (L): 177790655472633694188941980802
047503613813833685602824969586860483988457303292419157292935
380388916793723522118532246597515650465972468098640750969175
323067583642284054472809092057610201799293465748504101117625
584350450442885820883877261635394365497242394779959210685927
483970642556853937908745511045113355423116353373754214638700
488299349599339888965342105384720125168911889987302964445071
690804617411665559673197765683461268733759690639062267875652
078576707772030813055249319922273769990973469604361776775922
848041675652772916924480281074566297232981629176185283191324
907861265119304021792937302365364770562183880430197029735849
971627803188571242704480618796802428300307738635105278630447
639039111579815598101840646206149312815075326886519428941714
704663858921998805414402540381763432286269864844293075921085
556261175910325106160771636615881978320283927409653766398503
610130923687169102643128269846432862140570867543535243138038
049809462751772897752234749314486560174042187499283473193614
10650569707755492771627008.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 5 pixels is 5906
the DN mode is 9018
-------------------
BAND 9 (code 9)
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 5 pixels is 4992
the DN mode is 5039
-------------------
BAND 10 thermal (code 10)
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 11 thermal (code 11)
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]
-------------------
Calculating...
Writing reflectance of <crop_b1> to <cost_b1>...
Writing reflectance of <crop_b2> to <cost_b2>...
Writing reflectance of <crop_b3> to <cost_b3>...
Writing reflectance of <crop_b4> to <cost_b4>...
Writing reflectance of <crop_b5> to <cost_b5>...
Writing reflectance of <crop_b6> to <cost_b6>...
Writing reflectance of <crop_b7> to <cost_b7>...
Writing reflectance of <crop_b8> to <cost_b8>...
Writing reflectance of <crop_b9> to <cost_b9>...
Writing temperature of <crop_b10> to <cost_b10>...
Writing temperature of <crop_b11> to <cost_b11>...
(Mon Dec 16 12:43:54 2013) Command finished (49
sec)
-------clip------
Pekka
On Mon, Dec 9, 2013 at 10:43 PM, Markus Neteler <neteler at osgeo.org> wrote:
> On Mon, Nov 25, 2013 at 11:05 AM, Pekka Hurskainen
> <pekka.hurskainen at gmail.com> wrote:
> > Hi,
> >
> > Unless I have missed something obvious, the problem described by Beatrice
> > with i.landsat.toar still exists?
> >
> > I've tried both WinGrass 7.0.svn build 2013-11-24 and WinGrass 6.4.4svn,
> > revision 58275 (2013-11-24).
>
> I uploaded a fix on 21 Nov 2013:
>
> http://trac.osgeo.org/grass/changeset/58274/grass/branches/releasebranch_6_4/imagery/i.landsat.toar/landsat_met.c
> (likewise for GRASS 6.5 and 7)
>
>
> > I've tested with a L8 and L5 scene (both with the new metadata file),
> using
> > 'uncorrected' atmospheric correction method. I just chose the sensor type
> > and left everything else default.
> >
> > Both command output (verbose) and r.info print zeros for Earth-sun
> distance,
> > Sun elevation, DN range, Calibration constants, DN to Radiance, ESUN and
> > Radiance to Reflectance values.
>
> This sounds strange and should be verified.
>
> > The output files give sensible reflectance values, but can I trust them
> > (i.e. the metadata file was read correctly and the right parameters were
> > used)? I haven't compared with Linux version.
>
> Which scene do you work with?
>
> Markus
>
<pekka.hurskainen at gmail.com>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20131216/a02fc058/attachment-0001.html>
More information about the grass-user
mailing list