Earthdata Forum

Dear LP DAAC Team,

I hope your day is going great.

I need vegetation index (EVI, NDVI) from Harmonized Landsat and Sentinel (HLS) time series data.

I have run this Python code.

stac_lp = [s for s in stac_response['links'] if 'LP' in s['title']] # Search for only LP-specific catalogs

# LPCLOUD is the STAC catalog we will be using and exploring today
lp_cloud = r.get([s for s in stac_lp if s['title'] == 'LPCLOUD'][0]['href']).json()
for l in lp_cloud: print(f"{l}: {lp_cloud[l]}")

lp_links = lp_cloud['links']
for l in lp_links:
try:
print(f"{l['href']} is the {l['title']}")
except:
print(f"{l['href']}")

lp_collections = [l['href'] for l in lp_links if l['rel'] == 'collections'][0] # Set collections endpoint to variable
collections_response = r.get(f"{lp_collections}").json() # Call collections endpoint
print(f"This collection contains {collections_response['description']} ({len(collections_response['collections'])} available)")

collections = collections_response['collections']
collections[1]

My Python output showed that the data referenced to 'ASTER Global Digital Elevation Model NetCDF V003'

How to fix my Python code to extract HLS NDVI/EVI time series data?

Thank you very much for your help.

Bayu

Hi @bayu.suseno I am not sure why you are receiving ASTER data as the output. Let me check with our team that wrote the script and we'll get back to you. Thanks, Danielle

Hi @bayu.suseno

There was an issue with the workflow in that tutorial due to some changes made to CMR STAC. We are currently developing a new HLS tutorial that is available here: https://github.com/nasa/HLS-Data-Resources. The new repository is still under active development, so some resources are not fully refined, but the new tutorial in the python/tutorials folder has a similar use case calculating EVI. Please follow the setup instructions for the repository to set up the necessary Python environment.

The error in the R script looks like maybe the necessary gdal configuration options were not set. If the cell setting the gdal config options ran successfully during your attempt (shown below), try making sure that the `terra` package is version is 1.5.26 or newer.

setGDALconfig("GDAL_HTTP_UNSAFESSL", value = "YES")
setGDALconfig("GDAL_HTTP_COOKIEFILE", value = ".rcookies")
setGDALconfig("GDAL_HTTP_COOKIEJAR", value = ".rcookies")
setGDALconfig("GDAL_DISABLE_READDIR_ON_OPEN", value = "EMPTY_DIR")
setGDALconfig("CPL_VSIL_CURL_ALLOWED_EXTENSIONS", value = "TIF")

Hope this helps.

Regards,

LP DAAC User Services

Dear Daniel and LP DAAC User Services,

Thank you for your reply.

I will look for the code in the new HLS tutorial link and report the result soon.

Thank you very much for your help. I appreciate it.

Best regards,
Bayu Suseno

<a href="url">text</a>
[This link](http://example.net/)
[url=http://www.dynamicdrive.com/forums]Dynamic Drive Forums[/url]

Dear LP DAAC Team,

I have run the Python code in the HLS tutorial that is available here: https://github.com/nasa/HLS-Data-Resources. I am excited that the codes are running smoothly, so I could extract EVI from HLS2.0.

Thank you very much for your help. I appreciate it.

Best regards,
Bayu

Glad to hear it @bayu.suseno, please let us know if you run into any other issues or have any additional questions. Thanks -- Danielle