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

GRASS GRASS user
1

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.png

I entered the metadata

image.png

I left these two following tabs as default settings

image.png

image.png

and got this error:

i.landsat.toar input=B4@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@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

2

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@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.png

I entered the metadata

image.png

I left these two following tabs as default settings

image.png

image.png

and got this error:

i.landsat.toar input=B4@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@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

3

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@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@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.png

I entered the metadata

image.png

I left these two following tabs as default settings

image.png

image.png

and got this error:

i.landsat.toar input=B4@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@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

4

On Sun, Jun 16, 2019 at 12:25 PM Gabriel Cotlier <gabiklm01@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@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.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@gmail.com> wrote:

Thanks a lot Vero,
Rregards

On Fri, Jun 14, 2019 at 1:10 PM Veronica Andreo <veroandreo@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@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@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.png

I entered the metadata

image.png

I left these two following tabs as default settings

image.png

image.png

and got this error:

i.landsat.toar input=B4@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@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user