[GRASS-user] question on Landsat 8 radiometric and atmospheric corr.

[Veronica Andreo]

Sorry, I missed one part of your email.

Of course, you can run it from the terminal, the manual page has proper
examples: https://grass.osgeo.org/grass76/manuals/i.landsat.toar.html
and you can also follow the examples in this presentation:
https://gitpitch.com/veroandreo/grass-gis-conae/master?p=slides/04_imagery&grs=gitlab#/
(feedback is welcome, btw)

cheers,
Vero

El vie., 14 jun. 2019 a las 13:06, Veronica Andreo (<veroandreo at gmail.com>)
escribió:

> Hi Gabriel,
>
> AFAIU, you have to put the basename only, i.e., only B if your files are
> named B1, B2, B3 and so forth
>
> HTH,
> Vero
>
> El vie., 14 jun. 2019 a las 11:01, Gabriel Cotlier (<gabiklm01 at gmail.com>)
> escribió:
>
>> Dear grass users,
>>
>> My question is regrading the preprocessing of Landsat 8 data.
>> I have LC08_L1TP, level 1 Landsat 8 data product loaded in grass project
>> and want to perform the preprocessing of these data sets.
>> The procedure of turning DN to Radiance (radiometric calibration) by
>> means of *i.landsat.toar *grass function, and thereafter turning
>> obtained Radiance to Land Surface Reflectance (Atmospheric correction) by
>> means of grass function *i.atcorr* to get final reflectance results by
>> means different DOS methods to be further compared.
>>
>> I enter the data band
>> [image: image.png]
>> I entered the metadata
>>
>> [image: image.png]
>>
>>   I left these two following tabs as default settings
>> [image: image.png]
>>
>> [image: image.png]
>>
>> and got this error:
>>
>> i.landsat.toar input=B4 at PERMANENT output=B4_radiance
>> metfile=C:\Users\Gabriel\Documents\grassdata\LC08_L1TP_227083_20180315_20180320_01_T1\LC08_L1TP_227083_20180315_20180320_01_T1_MTL.txt
>> sensor=oli8 date=2018-03-15 sun_elevation=45.06637424
>> product_date=2018-03-20
>> WARNING: ESUN evaluated from REFLECTANCE_MAXIMUM_BAND
>> WARNING: Overwriting solar elevation of metadata file
>> Calculating...
>> ERROR: Unable to open header file for raster map <B4 at PERMANENT1@>
>> (Fri Jun 14 10:52:44 2019) Command finished (0 sec)
>>
>> I would like to load and run this process for many bands at once but the
>> GUI for this functions do not permit to input more than one band at once.
>> Could be this done more efficiently by command line? Any suggestion for
>> this preprocessing by command line or either by python scrip ?
>>
>> Thanks a lot.
>> Gabriel
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190614/2807766f/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 27831 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190614/2807766f/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 52585 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190614/2807766f/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 44350 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190614/2807766f/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 42666 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190614/2807766f/attachment-0007.png>