Vegetation Health Index - 231 m 8 days
FAIR Overall Score
Abstract
The Vegetation Health Index (VHI) is based on a combination of products extracted from vegetation signals, namely the Normalized Difference Vegetation Index (NDVI) and the land surface temperature, both derived from MODIS satellite data. The NDVI is based on 8 day maximum value composite MOD09Q1 (v006) reflectance and the land surface temperature (LST) on 8 day MOD11A2 (v006) LST products. The spatial resolution is 231 m, therefore the original 1000 m resolution of the MOD11A2 LST is downscaled to 231 m of the MOD09Q1 reflectance. Both products are masked to the highest quality standards using the provided quality layers. Missing pixel values in the time series are linearly interpolated. Non-vegetated areas are masked using the most recent Corine Land Cover product version for the according year. The final product is regridded to the LAEA Projection (EPSG:3035). The VHI relies on a strong inverse correlation between NDVI and land surface temperature, since increasing land temperatures are assumed to act negatively on vegetation vigour and consequently to cause stress. The data is provided as 8 day measures. The time series is starting from 2001. The VHI values range from 0-100, whereas high values correspond to healthy vegetation and low values indicate stressed vegetation.
Keywords
Legal constraints
CC-BY-4.0
FAIR Overall Score: 77%
| Principle | Score | Earned | Level |
|---|---|---|---|
| Findable | 57% | 4 of 7 | moderate |
| Accessible | 100% | 7 of 7 | advanced |
| Interoperable | 67% | 4 of 6 | moderate |
| Reusable | 83% | 5 of 6 | moderate |
Evaluated by F-UJI web service: Anusuriya Devaraju, & Robert Huber. (2020). F-UJI - An Automated FAIR Data Assessment Tool. Zenodo. https://doi.org/10.5281/zenodo.6361400
Snippet code
##### ----Explore and download STAC data with Python ----- #####
from pystac_client import Client
## Read the example catalog
URL = 'https://stac.eurac.edu/'
catalog = Client.open(URL)
## List the Collections in the given Catalog
stac_collections = list(catalog.get_collections())
print(f"Number of collections: {len(stac_collections)}")
## Print collection IDs
print("Collections IDs:")
for collection in stac_collections:
print(f"- {collection.id}")
print("-------------------------")
## Retrieve a specific collection
collection = catalog.get_collection("ID_OfTheCollection")
## Search for items in the collection
collection_items = list(catalog.search(collections=['ID_OfTheCollection'], max_items=10).items())
print(collection_items)
## Retrieve a list of the first 10 items belonging to a specific collection
item = collection.get_item("NameOfTheItem")
#print(list(item.assets.items())[0:10])
##Print the band name and the href of a specific item retrieved from the previous list
print(item.assets["ID_OfTheCollection"].title)
print(item.assets["ID_OfTheCollection"].href)
## Print Item’s assets through the assets attribute, which is a dictionary
for asset_key in item.assets:
asset = item.assets[asset_key]
print("{}: {} ({})".format(asset_key, asset.href, asset.media_type))
## ------------- DONWLOAD COG files from a specific collection ------- ####
from pystac_client import Client
import requests
## Read the example catalog
URL = 'https://stac.eurac.edu/'
catalog = Client.open(URL)
## List the Collections in the given Catalog
stac_collections = list(catalog.get_collections())
# print(f"Number of collections: {len(stac_collections)}")
## Print collection IDs
print("Collections IDs:")
for collection in stac_collections:
print(f"- {collection.id}")
print("-------------------------")
## Retrieve a specific collection
collection = catalog.get_collection("ID_OfTheCollection")
## Search for items in the collection
collection_items = list(catalog.search(collections=['ID_OfTheCollection'], max_items=10).items())
print(collection_items)
url_basepath = "https://eurac-eo.s3-eu-west-1.amazonaws.com/"
collection_url = url_basepath + "ID_OfTheCollection" + "/"
print(collection_url)
for item_id in collection_items:
# item = item.assets[asset_key]
print(item_id.id)
cog_file = str(item_id.id + '.tif')
url = collection_url + cog_file
response = requests.get(url)
with open(cog_file, "wb") as f:
f.write(response.content)
Related docs
| # | Name | Description | Link | Date published | Category |
|---|---|---|---|---|---|
| 1 | MOOC Cubes and Clouds | Free Online Course teaching the concepts of data cubes, cloud platforms and open science in geospatial and EO. | Link | March 8, 2024 | OpenEO, STAC |
| 2 | STAC guidelines | Documentation to browse and download items of the STAC catalog | Link | Feb. 22, 2024 | STAC |