<div dir="ltr"><div>Thanks a lot I will try all bands.</div><div>Best regards,</div><div>Gabriel<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 18, 2019 at 9:36 AM Carlos Henrique Grohmann de Carvalho <<a href="mailto:carlos.grohmann@gmail.com">carlos.grohmann@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Indeed it needs ALL bands. <div><br></div><div>I used this module in class a couple of weeks ago. You need to import and use all bands, including the panchromatic.<div><br></div><div>Carlos</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 18, 2019 at 9:32 AM Veronica Andreo <<a href="mailto:veroandreo@gmail.com" target="_blank">veroandreo@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Seems to me that the module needs all Landsat bands. Try importing all of them</div><div><br></div><div>Vero<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun., 17 jun. 2019 a las 15:54, Gabriel Cotlier (<<a href="mailto:gabiklm01@gmail.com" target="_blank">gabiklm01@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Vero, Nikos and Markus,</div><div><br></div><div>I'm trying to get an NDVI from different preprocessing atmospheric correction methods DOS1, DOS2, DOS 2b, DOS3, DOS4. However, unsuccessfully, I'm not being able to run my supposed workflow on GRASS GIS which is:</div><div><br></div><div>Step 1. DN to TOA Radiance (
<em><b>i.landsat.toar</b></em>)
</div><div>Step 2. TOA Radiance to Land Surface Reflectance (<em><b>i.atcorr</b></em>)</div><div><br></div><div>that thought maybe is the correct workflow for Landdsat 8 Level 1 products, maybe is not?</div><div><br></div><div>Since I need it only for Band 4 and Band 5 I proceeded as follows:</div><div><br></div><div>1. in a grassdata directory named with the correspondent Landsat 8 scene "LC82270832018074LGN00" I put Band 4 and Band 5 together with the MTL metadata txt file.</div><div><br></div><div>2. set a "Location" for this data set using Band 4 data by reading projection and Datum from Band 4.</div><div><br></div><div>3.started grass gis and d.rast for displaying the already automatically loaded Band 4<br></div><div><br></div><div>5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4 loaded "LC82270832018074LGN00_B5"</div><div><br></div><div><div><div><div><img src="cid:ii_jx0q1n4k2" alt="image.png" style="margin-right: 0px;" width="781" height="352"><br></div></div></div> </div><div>6. I changed working directory to where the data is:</div><div><br></div><div><div><img src="cid:ii_jx0q4gcf3" alt="image.png" style="margin-right: 0px;" width="786" height="328"><br></div> </div><div>7. in command line (the withe one in grass interface) I tried to
<em><b>i.landsat.toar</b></em> for obtaining TOA Radiance as follows :</div><div><br></div><div>i.landsat.toar -r input=LC82270832018074LGN00_B output=TOAR_LC82270832018074LGN00_B metfile=LC82270832018074LGN00_MTL.txt sensor=oli8 date=2018-03-15 sun_elevation=45.06637424 product_date=2018-03-20</div><div><br></div><div>and the resulting error was:</div><div><br></div><div><div><img src="cid:ii_jx0qc6yx4" alt="image.png" style="margin-right: 25px;" width="967" height="139"><br></div></div><div><br></div><div>I do not really understand what could be happening related to the bands names?</div><div><br></div><div>Thanks a lot I will appreciate a way to solve this problem. May target is to get finally both bands 4 and 5 radiometrically callibrated and further atmospherically corrected with
<em><b>i.atcorr</b></em>
using different DOS methos. </div><div><br></div><div>Thanks a lot again for the help and guidance.</div><div><br></div><div>Gabriel<br></div><div><br></div></div>
</blockquote></div>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_2067763056458882436gmail_signature"><div dir="ltr">Prof. Carlos Henrique Grohmann<br>Institute of Energy and Environment - Univ. of São Paulo, Brazil<div>- Digital Terrain Analysis | GIS | Remote Sensing - </div><div><br></div><div><a href="http://carlosgrohmann.com/" target="_blank">http://carlosgrohmann.com</a></div><div><a href="http://orcid.org/0000-0001-5073-5572" style="font-size:13px;color:rgb(17,85,204)" target="_blank">http://orcid.org/0000-0001-5073-5572</a><br><div>________________<br>Can’t stop the signal.</div></div></div></div>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div>