question on Landsat 8 radiometric and atmospheric corr.

classic Classic list List threaded Threaded
4 messages Options
Reply | Threaded
Open this post in threaded view
|

question on Landsat 8 radiometric and atmospheric corr.

Gabriel Cotlier
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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on Landsat 8 radiometric and atmospheric corr.

Veronica Andreo
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 (<[hidden email]>) 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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: 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 (<[hidden email]>) 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 (<[hidden email]>) 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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on Landsat 8 radiometric and atmospheric corr.

Gabriel Cotlier


On Sun, Jun 16, 2019 at 12:25 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
Thanks a lot Vero,
Rregards


On Fri, Jun 14, 2019 at 1:10 PM Veronica Andreo <[hidden email]> 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 (<[hidden email]>) 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 (<[hidden email]>) 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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user