Does anybody know how to process landsat 8 data with POLYMER-v4.11?

Post Reply
southocean
Posts: 11
Joined: Fri Jun 07, 2019 12:05 pm
company / institution: Sun Yat-Sen UNIVERSITY
Location: guangdong province China

Does anybody know how to process landsat 8 data with POLYMER-v4.11?

Post by southocean »

Hello!

I am new using POLYMER. And I want to use POLYMER-v4.11 to processlandsat 8 data. But I don't know how to start.
Because in the README.txt file and example.py I can't find any instructions on how to process landsat8 data. Does the README.txt file not updated? Can someone tell me?

Many thanks!
southocean
User avatar
fsteinmetz
Site Admin
Posts: 306
Joined: Fri Sep 07, 2018 1:34 pm
company / institution: Hygeos
Location: Lille, France
Contact:

Re: Does anybody know how to process landsat 8 data with POLYMER-v4.11?

Post by fsteinmetz »

Hello,

Landsat8 processing is quite new and experimental and the documentation is not ideal, but all sensors are processed in the same way.
For Landsat8, you can write a python script in polymer's root directory, which looks like :

Code: Select all

from polymer.main import run_atm_corr
from polymer.level1_landsat8 import Level1_OLI
from polymer.level2_nc import Level2_NETCDF

run_atm_corr(
	Level1_OLI(<landsat8_dirname>),
	Level2_NETCDF()
)
But please note that Level1_OLI, unlike the other sensors, requires the angles to be pre-calculated using a program (l8_angles, using the following command: l8_angles LC08_..._ANG.txt BOTH 1 -b 1) that is not distributed with Polymer, and which is available at https://github.com/DHI-GRAS/l8-angles

Kind regards,
François
southocean
Posts: 11
Joined: Fri Jun 07, 2019 12:05 pm
company / institution: Sun Yat-Sen UNIVERSITY
Location: guangdong province China

Re: Does anybody know how to process landsat 8 data with POLYMER-v4.11?

Post by southocean »

Hello François,

Thank you for your reply.But it took me a long time to try to install l8angles and I failed. I browse the website you provide(https://github.com/DHI-GRAS/l8-angles) ,and I found that the two installation methods provided in the readme file are based on the windows platform.
The l8angles installation package of anaconda is based on the windows platform(https://anaconda.org/search?q=l8angles), while I use the Linux platform to run this command

Code: Select all

conda config --add channels DHI-GRAS
conda install py-l8angles
it showed that the package could not be found.

Do you have any good suggestions for me?

Many thanks!
southocean
User avatar
fsteinmetz
Site Admin
Posts: 306
Joined: Fri Sep 07, 2018 1:34 pm
company / institution: Hygeos
Location: Lille, France
Contact:

Re: Does anybody know how to process landsat 8 data with POLYMER-v4.11?

Post by fsteinmetz »

Hi

Apologies, I was not pointing you to the proper code repository. Please find the landsat angles generation code and documentation at the following address :
https://www.usgs.gov/land-resources/nli ... ient-files
You can download the source code at the bottom of this page, see section Tools and Documentation, Landsat 8 Operational Land Imager (OLI)/Thermal Infrared Sensor (TIRS) (.tgz) - note that this is a command line tool that needs to be run prior to Polymer execution, not a Python module.

With kind regards,
François
Post Reply