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

Gabriel Cotlier gabiklm01 at gmail.com
Sun Jun 16 13:07:02 PDT 2019


On Sun, Jun 16, 2019 at 12:25 PM Gabriel Cotlier <gabiklm01 at gmail.com>
wrote:

> Hello Vero,
> The course material is excellent and careful explained very useful and
> interesting presentation and exercises I'm learning a lot with it.
> Thank you very much for sharing these important material with me.
>
> I have a questions I could construct the Location using already
> georeferenced data and now I'm running the following code:
>
> i.landsat.toar -r input=B4 at PERMANENT output=B4_toar
> metfile=C:\Users\Gabriel\Documents\grassdata\LC08_L1TP_227083_20180315_20180320_01_T1\LC08_L1TP_227083_20180315_20180320_01_T1_MTL.txt
> sensor=oli8 method=dos2 date=2018-03-15 sun_elevation=45.06637424
> product_date=2018-03-20
> and I'm getting the following error:
> [image: image.png]
>
> what could be happening?
> Thanks a lot again for your help
> Regards,
> Gabriel
>
>
> On Fri, Jun 14, 2019 at 1:44 PM Gabriel Cotlier <gabiklm01 at gmail.com>
> wrote:
>
>> Thanks a lot Vero,
>> Rregards
>>
>>
>> On Fri, Jun 14, 2019 at 1:10 PM Veronica Andreo <veroandreo at gmail.com>
>> wrote:
>>
>>> 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/20190616/98f09b91/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/20190616/98f09b91/attachment-0005.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/20190616/98f09b91/attachment-0006.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/20190616/98f09b91/attachment-0007.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/20190616/98f09b91/attachment-0008.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 8175 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190616/98f09b91/attachment-0009.png>


More information about the grass-user mailing list