<div dir="ltr"><div><div>Dear Markus,<br><br></div>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.<br>
<br></div><div></div>BTW, what is the difference between DOS2 and DOS2b methods?<br><div><br>-------clip------<br><br>(Mon Dec 16 12:33:37 2013) <br>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<br>
Metada file is MTL file: new format<br>RADIANCE & QUANTIZE from MIN_MAX_(RADIANCE|PIXEL_VALUE)<br>WARNING: ESUN evaluated from REFLECTANCE_MAXIMUM_BAND<br>Calculating dark pixel of <crop_b1>...<br>... DN = 7840 [5] : mode 9627 [2891], exclusing DN > 58981<br>
... TAUv = 0.78835, TAUz = 0.76206, Edown = 0.00000<br>Calculating dark pixel of <crop_b2>...<br>... DN = 7135 [5] : mode 9083 [2409], exclusing DN > 58981<br>... TAUv = 0.84575, TAUz = 0.82578, Edown = 0.00000<br>
Calculating dark pixel of <crop_b3>...<br>... DN = 6200 [5] : mode 8799 [1880], exclusing DN > 58981<br>... TAUv = 0.91432, TAUz = 0.90271, Edown = 0.00000<br>Calculating dark pixel of <crop_b4>...<br>... DN = 5785 [6] : mode 6663 [1183], exclusing DN > 58981<br>
... TAUv = 0.95287, TAUz = 0.94633, Edown = 0.00000<br>Calculating dark pixel of <crop_b5>...<br>... DN = 6024 [7] : mode 14327 [927], exclusing DN > 58981<br>... TAUv = 0.98442, TAUz = 0.98222, Edown = 0.00000<br>
Calculating dark pixel of <crop_b6>...<br>... DN = 5377 [5] : mode 11683 [721], exclusing DN > 58981<br>... TAUv = 0.99871, TAUz = 0.99852, Edown = 0.00000<br>Calculating dark pixel of <crop_b7>...<br>... DN = 5166 [5] : mode 10961 [843], exclusing DN > 58981<br>
... TAUv = 0.99963, TAUz = 0.99958, Edown = 0.00000<br>Calculating dark pixel of <crop_b8>...<br>... DN = 5906 [6] : mode 9018 [5640], exclusing DN > 58981<br>... TAUv = 0.92888, TAUz = 0.91916, Edown = 0.00000<br>
Calculating dark pixel of <crop_b9>...<br>... DN = 4992 [5] : mode 5039 [133488], exclusing DN > 58981<br>... TAUv = 0.99757, TAUz = 0.99722, Edown = 0.00000<br> LANDSAT: 8 SENSOR: OLI/TIRS<br> ACQUISITION DATE 2013-08-27 [production date 2013-08-27]<br>
Earth-sun distance = 0.00000000<br> Solar elevation angle = 0.00000000<br> Atmospheric correction: dos3<br> Percent of solar irradiance in path radiance = 0.0000<br>-------------------<br> BAND 1 (code 1)<br>
calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): 968367649887747415141875521950<br>925908889125765915809137356521975470948435722379319114916517<br>109681257194579805986272450667354678861757535842162547084653<br>462145348007752249138970238330548827196232199351249578581263<br>
248711758443327663076681220844948160752094949553961518385593<br>288738424407009196744807773476057507004619157871083200622000<br>320486055830136042443031065077914185459331507512811830598511<br>152020147139378105907351040013136335438675174223655826910682<br>
669147134563183589280887783585976043716156080824364735851118<br>183768026326159743796155841003601516760621076483451614820141<br>371685937002372918255194596159954067917411389793697522459155<br>015399457565770890778803961998816054585345162915220244370597<br>
148522387936504415391544888346434689270315995072771028338997<br>995727094089743160130590886592801301509071858609333431506640<br>587793968989583037059469778982726604216980071724204307021521<br>635177260328502996043398383316279662847376171701723950431400<br>
757812176165083327750696905811530775310780096413788963894683<br>905077061031539383276262428122121528157514609481673317570684<br>435219591138433952053577693737258063060125527179863820119388<br>770600138288962007966562469949923301233492925582540474476621<br>
929036185589173077700588636992204749746454624011817064126344<br>632245678976560320724068426895929872763359600665707344187778<br>891947245102403001861675618641623023833022185852578363120467<br>361109134430306902862293966757899915186924084339575936028017<br>
743399931822767869305933232821423387100836970448169974752735<br>074832657312708225553213351832437155472265900799102082698841<br>521263830766549846516292538579807583602895447499728905211008<br>406860112609801193185282209285068919448526126257634719443776<br>
306035623536657971781547855260622470703165809089946123219856<br>394296744397591189893603472912483484277176386264413758344436<br>893155901820405208007799764626033880161122800644723691809874<br>031293719596978952957939025121263874432512217233471764941837<br>
349859301579149060534601918858766260236011294832113712139455<br>698498970564770118735817734482120560350401517034380439852381<br>686780485705510656529346796692015942241400805762289636824034<br>380012173568901522372612940698148203094590690817052813619697<br>
272631456298510647656983077748357220329595173179533274877623<br>690664824083139848333834479764874064891277168954263595672454<br>741560283118661607390700810458152871418035457554976493129153<br>6739339299897124622932302102528.00000 to 0.00000<br>
at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): -0.00000<br> at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 7840<br> the DN mode is 9627<br>
-------------------<br> BAND 2 (code 2)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): -0.00000 to 0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br>
at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 7135<br> the DN mode is 9083<br>-------------------<br> BAND 3 (code 3)<br> calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): -0.00000 to 0.00000<br>
at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 6200<br> the DN mode is 8799<br>
-------------------<br> BAND 4 (code 4)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 183658085783589171116602975780<br>576110409539519711468540243914721655904198177980762961939111<br>
367605627562671716937127859554390051130688152018870728220925<br>
072842007885913864783619283774492980062937015097153316113130<br>872437715945381379739519113601595608897104545933411437714766<br>964222032417753040188158845816964370334214805486222793513041<br>392415948524253835628431799076337651616691066413747194706092<br>
588365827695512535302789689209508312728090357719245288559169<br>261388427087636514887083888570441112579969923313176980814385<br>626616270691354605624826831921938182203368446708429571503396<br>673727480549490161133959472989086466065342736246984472386035<br>
984295478176953900049240996927707494912239369479071549791691<br>006705649149921302160640916967785780445041977901000120712458<br>693774631060343641966598542706661151042560211066330115225274<br>882082359272986006411181268975110156015073466360964844002614<br>
923902881591929598593433704046771611849396994917616009430564<br>636178472140317666513101744160796308166097582587535447676030<br>395599188338919046557681998905510895192924869482845313899438<br>243168242481424950199879384022305364926654416732220470906145<br>
003035888135318340747367027552221587693565002548093465052373<br>146810072538472500364217293690318091228705920528403356637468<br>515206939874589256642932818842704726527146610272815374838100<br>196715703591772856115604922172749774983079205445635592246176<br>
519699396484278000996109577040304620492915499240374708941194<br>899957749080944402792983491485854033303690426232983954252982<br>089745184902031081081879977151639415526477184204943465088423<br>197754901126349507777941665413693263267204404932798742686446<br>
175425099611417857336481903514883801243132058814981214699218<br>223830439709555518391135039601061824428651849472616761261741<br>947438101584718397168382882672948322751691531336445977140261<br>8847363072.00000 to 0.00000<br>
at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 5785<br> the DN mode is 6663<br>
-------------------<br> BAND 5 (code 5)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): -0.00000<br>
at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 6024<br> the DN mode is 14327<br>-------------------<br> BAND 6 (code 6)<br> calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): 339445567057160696751644954879<br>970406447096914339798998362968503888855957504.00000 to<br>0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br>
at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 5377<br> the DN mode is 11683<br>-------------------<br> BAND 7 (code 7)<br> calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): -12134434465588505055597499511<br>824959727712116382507612325036502948239426558532389801597855<br>534050652218225524835841065680524520096072520033585910067094<br>864419565702430636692147519347440048706033789010796831449675<br>
609397487722631716840845517015635835853269211738830091674017<br>125859807797852579371667606442129288215138677555089849122233<br>496258894611200291018407686794712556602863468629113156241991<br>764561853870347050462715860062969867544293321888080259438748<br>
975116227397144258398582860683531866628235263927302795102064<br>876102210094737464694328482051327435929619847514195837359477<br>453721352442867031673331032270809238607484683630583314906933<br>393968396774392452257227334055169647686113972735180171551034<br>
094817132579151468969733950921315554583854644271192962176107<br>852279850198497591348543860021633007244464057151870155597981<br>968744623129348330236016103587706033254928248063411627431203<br>939460766333363995335790489660112468340889170867819092108790<br>
845331032250512936792412018170190148253165743211304227984791<br>418825798114108108620478601673212231026083304899473015984957<br>797069887035981967325806061173567396266179556816201804747668<br>19328.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br>
mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 5166<br> the DN mode is 10961<br>-------------------<br> BAND 8 (code 8)<br>
calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 177790655472633694188941980802<br>047503613813833685602824969586860483988457303292419157292935<br>380388916793723522118532246597515650465972468098640750969175<br>
323067583642284054472809092057610201799293465748504101117625<br>584350450442885820883877261635394365497242394779959210685927<br>483970642556853937908745511045113355423116353373754214638700<br>488299349599339888965342105384720125168911889987302964445071<br>
690804617411665559673197765683461268733759690639062267875652<br>078576707772030813055249319922273769990973469604361776775922<br>848041675652772916924480281074566297232981629176185283191324<br>907861265119304021792937302365364770562183880430197029735849<br>
971627803188571242704480618796802428300307738635105278630447<br>639039111579815598101840646206149312815075326886519428941714<br>704663858921998805414402540381763432286269864844293075921085<br>556261175910325106160771636615881978320283927409653766398503<br>
610130923687169102643128269846432862140570867543535243138038<br>049809462751772897752234749314486560174042187499283473193614<br>10650569707755492771627008.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br>
mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 5906<br> the DN mode is 9018<br>-------------------<br> BAND 9 (code 9)<br>
calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br>
the darkness DN with a least 5 pixels is 4992<br> the DN mode is 5039<br>-------------------<br> BAND 10 thermal (code 10)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br>
at-surface radiance = -0.00000000 * DN + -0.00000<br> at-sensor temperature = -0.00000 / log[(0.00000 /<br>radiance) + 1.0]<br>-------------------<br> BAND 11 thermal (code 11)<br> calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> at-sensor temperature = -0.00000 / log[(0.00000 /<br>radiance) + 1.0]<br>-------------------<br>Calculating...<br>
Writing reflectance of <crop_b1> to <dos3_3_b1>...<br>Writing reflectance of <crop_b2> to <dos3_3_b2>...<br>Writing reflectance of <crop_b3> to <dos3_3_b3>...<br>Writing reflectance of <crop_b4> to <dos3_3_b4>...<br>
Writing reflectance of <crop_b5> to <dos3_3_b5>...<br>Writing reflectance of <crop_b6> to <dos3_3_b6>...<br>Writing reflectance of <crop_b7> to <dos3_3_b7>...<br>Writing reflectance of <crop_b8> to <dos3_3_b8>...<br>
Writing reflectance of <crop_b9> to <dos3_3_b9>...<br>Writing temperature of <crop_b10> to <dos3_3_b10>...<br>Writing temperature of <crop_b11> to <dos3_3_b11>...<br>(Mon Dec 16 12:34:44 2013) Command finished (1 min 7 sec) <br>
(Mon Dec 16 12:43:05 2013) <br>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=<a href="tel:60.02214600" value="+16002214600" target="_blank">60.02214600</a> pixel=5<br>
Metada file is MTL file: new format<br>RADIANCE & QUANTIZE from MIN_MAX_(RADIANCE|PIXEL_VALUE)<br>WARNING: ESUN evaluated from REFLECTANCE_MAXIMUM_BAND<br>WARNING: Overwriting solar elevation of metadata file<br>Calculating dark pixel of <crop_b1>...<br>
... DN = 7840 [5] : mode 9627 [2891], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>Calculating dark pixel of <crop_b2>...<br>... DN = 7135 [5] : mode 9083 [2409], exclusing DN > 58981<br>
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>Calculating dark pixel of <crop_b3>...<br>... DN = 6200 [5] : mode 8799 [1880], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>
Calculating dark pixel of <crop_b4>...<br>... DN = 5785 [6] : mode 6663 [1183], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>Calculating dark pixel of <crop_b5>...<br>... DN = 6024 [7] : mode 14327 [927], exclusing DN > 58981<br>
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>Calculating dark pixel of <crop_b6>...<br>... DN = 5377 [5] : mode 11683 [721], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000<br>
Calculating dark pixel of <crop_b7>...<br>... DN = 5166 [5] : mode 10961 [843], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000<br>Calculating dark pixel of <crop_b8>...<br>... DN = 5906 [6] : mode 9018 [5640], exclusing DN > 58981<br>
... TAUv = 1.00000, TAUz = 0.86622, Edown = 0.00000<br>Calculating dark pixel of <crop_b9>...<br>... DN = 4992 [5] : mode 5039 [133488], exclusing DN > 58981<br>... TAUv = 1.00000, TAUz = 1.00000, Edown = 0.00000<br>
LANDSAT: 8 SENSOR: OLI/TIRS<br> ACQUISITION DATE 2013-08-27 [production date 2013-08-27]<br> Earth-sun distance = 0.00000000<br> Solar elevation angle = 0.00000000<br> Atmospheric correction: dos2<br> Percent of solar irradiance in path radiance = 0.0000<br>
-------------------<br> BAND 1 (code 1)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 968367649887747415141875521950<br>925908889125765915809137356521975470948435722379319114916517<br>
109681257194579805986272450667354678861757535842162547084653<br>
462145348007752249138970238330548827196232199351249578581263<br>248711758443327663076681220844948160752094949553961518385593<br>288738424407009196744807773476057507004619157871083200622000<br>320486055830136042443031065077914185459331507512811830598511<br>
152020147139378105907351040013136335438675174223655826910682<br>669147134563183589280887783585976043716156080824364735851118<br>183768026326159743796155841003601516760621076483451614820141<br>371685937002372918255194596159954067917411389793697522459155<br>
015399457565770890778803961998816054585345162915220244370597<br>148522387936504415391544888346434689270315995072771028338997<br>995727094089743160130590886592801301509071858609333431506640<br>587793968989583037059469778982726604216980071724204307021521<br>
635177260328502996043398383316279662847376171701723950431400<br>757812176165083327750696905811530775310780096413788963894683<br>905077061031539383276262428122121528157514609481673317570684<br>435219591138433952053577693737258063060125527179863820119388<br>
770600138288962007966562469949923301233492925582540474476621<br>929036185589173077700588636992204749746454624011817064126344<br>632245678976560320724068426895929872763359600665707344187778<br>891947245102403001861675618641623023833022185852578363120467<br>
361109134430306902862293966757899915186924084339575936028017<br>743399931822767869305933232821423387100836970448169974752735<br>074832657312708225553213351832437155472265900799102082698841<br>521263830766549846516292538579807583602895447499728905211008<br>
406860112609801193185282209285068919448526126257634719443776<br>306035623536657971781547855260622470703165809089946123219856<br>394296744397591189893603472912483484277176386264413758344436<br>893155901820405208007799764626033880161122800644723691809874<br>
031293719596978952957939025121263874432512217233471764941837<br>349859301579149060534601918858766260236011294832113712139455<br>698498970564770118735817734482120560350401517034380439852381<br>686780485705510656529346796692015942241400805762289636824034<br>
380012173568901522372612940698148203094590690817052813619697<br>272631456298510647656983077748357220329595173179533274877623<br>690664824083139848333834479764874064891277168954263595672454<br>741560283118661607390700810458152871418035457554976493129153<br>
6739339299897124622932302102528.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): -0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 7840<br>
the DN mode is 9627<br>-------------------<br> BAND 2 (code 2)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): -0.00000 to 0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br>
mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 7135<br> the DN mode is 9083<br>-------------------<br> BAND 3 (code 3)<br>
calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): -0.00000 to 0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / 0.00000<br>
the darkness DN with a least 5 pixels is 6200<br> the DN mode is 8799<br>-------------------<br> BAND 4 (code 4)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 183658085783589171116602975780<br>
576110409539519711468540243914721655904198177980762961939111<br>367605627562671716937127859554390051130688152018870728220925<br>072842007885913864783619283774492980062937015097153316113130<br>872437715945381379739519113601595608897104545933411437714766<br>
964222032417753040188158845816964370334214805486222793513041<br>392415948524253835628431799076337651616691066413747194706092<br>588365827695512535302789689209508312728090357719245288559169<br>261388427087636514887083888570441112579969923313176980814385<br>
626616270691354605624826831921938182203368446708429571503396<br>673727480549490161133959472989086466065342736246984472386035<br>984295478176953900049240996927707494912239369479071549791691<br>006705649149921302160640916967785780445041977901000120712458<br>
693774631060343641966598542706661151042560211066330115225274<br>882082359272986006411181268975110156015073466360964844002614<br>923902881591929598593433704046771611849396994917616009430564<br>636178472140317666513101744160796308166097582587535447676030<br>
395599188338919046557681998905510895192924869482845313899438<br>243168242481424950199879384022305364926654416732220470906145<br>003035888135318340747367027552221587693565002548093465052373<br>146810072538472500364217293690318091228705920528403356637468<br>
515206939874589256642932818842704726527146610272815374838100<br>196715703591772856115604922172749774983079205445635592246176<br>519699396484278000996109577040304620492915499240374708941194<br>899957749080944402792983491485854033303690426232983954252982<br>
089745184902031081081879977151639415526477184204943465088423<br>197754901126349507777941665413693263267204404932798742686446<br>175425099611417857336481903514883801243132058814981214699218<br>223830439709555518391135039601061824428651849472616761261741<br>
947438101584718397168382882672948322751691531336445977140261<br>8847363072.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br>
the darkness DN with a least 5 pixels is 5785<br> the DN mode is 6663<br>-------------------<br> BAND 5 (code 5)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br>
at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): -0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 6024<br> the DN mode is 14327<br>
-------------------<br> BAND 6 (code 6)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 339445567057160696751644954879<br>970406447096914339798998362968503888855957504.00000 to<br>0.00000<br>
at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 5377<br> the DN mode is 11683<br>
-------------------<br> BAND 7 (code 7)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): -12134434465588505055597499511<br>824959727712116382507612325036502948239426558532389801597855<br>
534050652218225524835841065680524520096072520033585910067094<br>
864419565702430636692147519347440048706033789010796831449675<br>609397487722631716840845517015635835853269211738830091674017<br>125859807797852579371667606442129288215138677555089849122233<br>496258894611200291018407686794712556602863468629113156241991<br>
764561853870347050462715860062969867544293321888080259438748<br>975116227397144258398582860683531866628235263927302795102064<br>876102210094737464694328482051327435929619847514195837359477<br>453721352442867031673331032270809238607484683630583314906933<br>
393968396774392452257227334055169647686113972735180171551034<br>094817132579151468969733950921315554583854644271192962176107<br>852279850198497591348543860021633007244464057151870155597981<br>968744623129348330236016103587706033254928248063411627431203<br>
939460766333363995335790489660112468340889170867819092108790<br>845331032250512936792412018170190148253165743211304227984791<br>418825798114108108620478601673212231026083304899473015984957<br>797069887035981967325806061173567396266179556816201804747668<br>
19328.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / -0.00000<br> the darkness DN with a least 5 pixels is 5166<br>
the DN mode is 10961<br>-------------------<br> BAND 8 (code 8)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 177790655472633694188941980802<br>047503613813833685602824969586860483988457303292419157292935<br>
380388916793723522118532246597515650465972468098640750969175<br>323067583642284054472809092057610201799293465748504101117625<br>584350450442885820883877261635394365497242394779959210685927<br>483970642556853937908745511045113355423116353373754214638700<br>
488299349599339888965342105384720125168911889987302964445071<br>690804617411665559673197765683461268733759690639062267875652<br>078576707772030813055249319922273769990973469604361776775922<br>848041675652772916924480281074566297232981629176185283191324<br>
907861265119304021792937302365364770562183880430197029735849<br>971627803188571242704480618796802428300307738635105278630447<br>639039111579815598101840646206149312815075326886519428941714<br>704663858921998805414402540381763432286269864844293075921085<br>
556261175910325106160771636615881978320283927409653766398503<br>610130923687169102643128269846432862140570867543535243138038<br>049809462751772897752234749314486560174042187499283473193614<br>10650569707755492771627008.00000 to 0.00000<br>
at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br> at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 5906<br> the DN mode is 9018<br>
-------------------<br> BAND 9 (code 9)<br> calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = 0.00000000 * DN + -0.00000<br> mean solar exoatmospheric irradiance (ESUN): 0.00000<br>
at-surface reflectance = radiance / 0.00000<br> the darkness DN with a least 5 pixels is 4992<br> the DN mode is 5039<br>-------------------<br> BAND 10 thermal (code 10)<br> calibrated digital number (DN): 0.0 to 0.0<br>
calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> at-sensor temperature = -0.00000 / log[(0.00000 /<br>radiance) + 1.0]<br>-------------------<br> BAND 11 thermal (code 11)<br>
calibrated digital number (DN): 0.0 to 0.0<br> calibration constants (L): 0.00000 to 0.00000<br> at-surface radiance = -0.00000000 * DN + -0.00000<br> at-sensor temperature = -0.00000 / log[(0.00000 /<br>radiance) + 1.0]<br>
-------------------<br>Calculating...<br>Writing reflectance of <crop_b1> to <cost_b1>...<br>Writing reflectance of <crop_b2> to <cost_b2>...<br>Writing reflectance of <crop_b3> to <cost_b3>...<br>
Writing reflectance of <crop_b4> to <cost_b4>...<br>Writing reflectance of <crop_b5> to <cost_b5>...<br>Writing reflectance of <crop_b6> to <cost_b6>...<br>Writing reflectance of <crop_b7> to <cost_b7>...<br>
Writing reflectance of <crop_b8> to <cost_b8>...<br>Writing reflectance of <crop_b9> to <cost_b9>...<br>Writing temperature of <crop_b10> to <cost_b10>...<br>Writing temperature of <crop_b11> to <cost_b11>...<br>
(Mon Dec 16 12:43:54 2013) Command finished (49 sec) <br><br></div><div class="gmail_extra">-------clip------<br><br></div><div class="gmail_extra">Pekka<br></div><div class="gmail_extra"><br><div class="gmail_quote">
On Mon, Dec 9, 2013 at 10:43 PM, Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>On Mon, Nov 25, 2013 at 11:05 AM, Pekka Hurskainen<br>
<<a href="mailto:pekka.hurskainen@gmail.com" target="_blank">pekka.hurskainen@gmail.com</a>> wrote:<br>
> Hi,<br>
><br>
> Unless I have missed something obvious, the problem described by Beatrice<br>
> with i.landsat.toar still exists?<br>
><br>
> I've tried both WinGrass 7.0.svn build 2013-11-24 and WinGrass 6.4.4svn,<br>
> revision 58275 (2013-11-24).<br>
<br>
</div>I uploaded a fix on 21 Nov 2013:<br>
<a href="http://trac.osgeo.org/grass/changeset/58274/grass/branches/releasebranch_6_4/imagery/i.landsat.toar/landsat_met.c" target="_blank">http://trac.osgeo.org/grass/changeset/58274/grass/branches/releasebranch_6_4/imagery/i.landsat.toar/landsat_met.c</a><br>
(likewise for GRASS 6.5 and 7)<br>
<div><br>
<br>
> I've tested with a L8 and L5 scene (both with the new metadata file), using<br>
> 'uncorrected' atmospheric correction method. I just chose the sensor type<br>
> and left everything else default.<br>
><br>
> Both command output (verbose) and <a href="http://r.info" target="_blank">r.info</a> print zeros for Earth-sun distance,<br>
> Sun elevation, DN range, Calibration constants, DN to Radiance, ESUN and<br>
> Radiance to Reflectance values.<br>
<br>
</div>This sounds strange and should be verified.<br>
<div><br>
> The output files give sensible reflectance values, but can I trust them<br>
> (i.e. the metadata file was read correctly and the right parameters were<br>
> used)? I haven't compared with Linux version.<br>
<br>
</div>Which scene do you work with?<br>
<span><font color="#888888"><br>
Markus<br>
</font></span></blockquote></div><br><br clear="all"><br><a href="mailto:pekka.hurskainen@gmail.com" target="_blank"></a>
</div></div>