question on Landsat 8 processing

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

question on Landsat 8 processing

Gabriel Cotlier
Dear Vero, Nikos and Markus,

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:

Step 1. DN to TOA Radiance ( i.landsat.toar)
Step 2. TOA Radiance to Land Surface Reflectance  (i.atcorr)

that thought maybe is the correct workflow for Landdsat 8 Level 1 products, maybe is not?

Since I need it only for Band 4 and Band 5 I proceeded as follows:

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.

2. set a "Location" for this data set using Band 4 data by reading projection and Datum from Band 4.

3.started grass gis and d.rast for displaying the already automatically loaded Band 4

5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4 loaded "LC82270832018074LGN00_B5"

image.png
6. I changed working directory to where the data is:

image.png
7. in command line (the withe one in grass interface) I tried to  i.landsat.toar for obtaining TOA Radiance as follows :

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

and the resulting error was:

image.png

I do not really understand what could be happening related to the bands names?

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 i.atcorr  using different DOS methos.  

Thanks a lot again for the help and guidance.

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 processing

Veronica Andreo
Seems to me that the module needs all Landsat bands. Try importing all of them

Vero

El lun., 17 jun. 2019 a las 15:54, Gabriel Cotlier (<[hidden email]>) escribió:
Dear Vero, Nikos and Markus,

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:

Step 1. DN to TOA Radiance ( i.landsat.toar)
Step 2. TOA Radiance to Land Surface Reflectance  (i.atcorr)

that thought maybe is the correct workflow for Landdsat 8 Level 1 products, maybe is not?

Since I need it only for Band 4 and Band 5 I proceeded as follows:

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.

2. set a "Location" for this data set using Band 4 data by reading projection and Datum from Band 4.

3.started grass gis and d.rast for displaying the already automatically loaded Band 4

5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4 loaded "LC82270832018074LGN00_B5"

image.png
6. I changed working directory to where the data is:

image.png
7. in command line (the withe one in grass interface) I tried to  i.landsat.toar for obtaining TOA Radiance as follows :

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

and the resulting error was:

image.png

I do not really understand what could be happening related to the bands names?

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 i.atcorr  using different DOS methos.  

Thanks a lot again for the help and guidance.

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 processing

Carlos Henrique Grohmann de Carvalho-2
Indeed it needs ALL bands. 

I used this module in class a couple of weeks ago. You need to import and use all bands, including the panchromatic.

Carlos



On Tue, Jun 18, 2019 at 9:32 AM Veronica Andreo <[hidden email]> wrote:
Seems to me that the module needs all Landsat bands. Try importing all of them

Vero

El lun., 17 jun. 2019 a las 15:54, Gabriel Cotlier (<[hidden email]>) escribió:
Dear Vero, Nikos and Markus,

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:

Step 1. DN to TOA Radiance ( i.landsat.toar)
Step 2. TOA Radiance to Land Surface Reflectance  (i.atcorr)

that thought maybe is the correct workflow for Landdsat 8 Level 1 products, maybe is not?

Since I need it only for Band 4 and Band 5 I proceeded as follows:

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.

2. set a "Location" for this data set using Band 4 data by reading projection and Datum from Band 4.

3.started grass gis and d.rast for displaying the already automatically loaded Band 4

5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4 loaded "LC82270832018074LGN00_B5"

image.png
6. I changed working directory to where the data is:

image.png
7. in command line (the withe one in grass interface) I tried to  i.landsat.toar for obtaining TOA Radiance as follows :

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

and the resulting error was:

image.png

I do not really understand what could be happening related to the bands names?

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 i.atcorr  using different DOS methos.  

Thanks a lot again for the help and guidance.

Gabriel

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


--
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing - 

http://orcid.org/0000-0001-5073-5572
________________
Can’t stop the signal.

_______________________________________________
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 processing

Gabriel Cotlier
Thanks a lot I will try all bands.
Best regards,
Gabriel

On Tue, Jun 18, 2019 at 9:36 AM Carlos Henrique Grohmann de Carvalho <[hidden email]> wrote:
Indeed it needs ALL bands. 

I used this module in class a couple of weeks ago. You need to import and use all bands, including the panchromatic.

Carlos



On Tue, Jun 18, 2019 at 9:32 AM Veronica Andreo <[hidden email]> wrote:
Seems to me that the module needs all Landsat bands. Try importing all of them

Vero

El lun., 17 jun. 2019 a las 15:54, Gabriel Cotlier (<[hidden email]>) escribió:
Dear Vero, Nikos and Markus,

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:

Step 1. DN to TOA Radiance ( i.landsat.toar)
Step 2. TOA Radiance to Land Surface Reflectance  (i.atcorr)

that thought maybe is the correct workflow for Landdsat 8 Level 1 products, maybe is not?

Since I need it only for Band 4 and Band 5 I proceeded as follows:

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.

2. set a "Location" for this data set using Band 4 data by reading projection and Datum from Band 4.

3.started grass gis and d.rast for displaying the already automatically loaded Band 4

5. import with r.in.gdal Band 5 using Band 5 complete name as was Band 4 loaded "LC82270832018074LGN00_B5"

image.png
6. I changed working directory to where the data is:

image.png
7. in command line (the withe one in grass interface) I tried to  i.landsat.toar for obtaining TOA Radiance as follows :

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

and the resulting error was:

image.png

I do not really understand what could be happening related to the bands names?

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 i.atcorr  using different DOS methos.  

Thanks a lot again for the help and guidance.

Gabriel

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


--
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing - 

http://orcid.org/0000-0001-5073-5572
________________
Can’t stop the signal.
_______________________________________________
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 processing

Markus Neteler
On Tue, Jun 18, 2019 at 2:39 PM Gabriel Cotlier <[hidden email]> wrote:
>
> Thanks a lot I will try all bands.

If there is anything to be improved in the manual page, please suggest
through (ideally) a pull request:

https://github.com/OSGeo/grass/blob/master/imagery/i.landsat.toar/i.landsat.toar.html

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