AfriSAR Search and Visualize

Description

This dataset provides gridded estimates of aboveground biomass (AGB) for four sites in Gabon at 0.25 ha (50 m) resolution derived with field measurements and airborne LiDAR data collected from 2010 to 2016. The sites represent a mix of forested, savannah, and some agricultural and disturbed landcover types: Lope site, within Lope National Park; Mabounie, mostly forested site; Mondah Forest, protected area; and the Rabi forest site, part of the Smithsonian Institution of Global Earth Observatories world-wide network of forest plots. Plot-level biophysical measurements of tree diameter and tree height (or estimated by allometry) were performed at 1 ha and 0.25 ha scales on multiple plots at each site and used to derive AGB for each tree and then summed for each plot. Aerial LiDAR scans were used to construct digital elevation models (DEM) and digital surface models (DSM), and then the DEM and DSM were used to construct a canopy height model (CHM) at 1 m resolution. After checking site-plot locations against the CHM, mean canopy height (MCH) was computed over each 0.25 ha. A single regression model relating MCH and AGB estimates, incorporating local height based on the trunk DBH (HD) relationships, was produced for all sites and combined with the CHM layer to construct biomass maps at 0.25 ha resolution.

Source: https://daac.ornl.gov/AFRISAR/guides/AfriSAR_AGB_Maps.html

[1]:

# Uncomment and run the line below if missing packages
# !pip install pystac-client rioxarray hvplot

[ ]:


import rioxarray
import rasterio as rio
from pystac_client import Client
import hvplot.xarray

import warnings
warnings.filterwarnings("ignore")

Declare your collection of interest

You can discover available collections the follwoing ways:

  1. Use the {STAC_API_URL}/collections API endpoint (JSON response)
  2. Programmatically using pystac (see example in the list-collections.ipynb notebook)
  3. In the STAC Browser : http://delta-staging-stac-browser.s3-website-us-east-1.amazonaws.com/
[2]:

STAC_API_URL = 'https://stac.maap-project.org/'
collection = 'AfriSAR_AGB_Maps_1681'

[3]:

# custom headers
headers = []

cat = Client.open(STAC_API_URL, headers=headers)

Listing all the available collections

[4]:

for collection in cat.get_collections():
    print(collection)

<CollectionClient id=AfriSAR_UAVSAR_Coreg_SLC>
<CollectionClient id=Landsat8_SurfaceReflectance>
<CollectionClient id=AfriSAR_AGB_Maps_1681>
<CollectionClient id=ABLVIS1B>
<CollectionClient id=GEDI02_B>
<CollectionClient id=AFLVIS2>
<CollectionClient id=BIOSAR1>
<CollectionClient id=AFRISAR_DLR>
<CollectionClient id=GEDI02_A>

Getting the collection of intrest ‘AfriSAR_AGB_Maps_1681’

[5]:

collection = cat.get_collection('AfriSAR_AGB_Maps_1681')
collection

[5]:

Use satsearch to discover items in the ‘AfriSAR_AGB_Maps_1681’ collection

[6]:

search = cat.search(
    max_items = 15,
    limit = 5,
    collections = collection,
)

Adding items to an iterable

[7]:

items = [i.to_dict() for i in search.get_all_items()]
items[0]

[7]:

{'type': 'Feature',
 'stac_version': '1.0.0',
 'id': 'Rabi_AGB_50m',
 'properties': {'boxes': ['-1.94602 9.85914 -1.90031 9.90636'],
  'links': [{'rel': 'http://esipfed.org/ns/fedsearch/1.1/s3#',
    'href': 's3://nasa-maap-data-store/file-staging/nasa-map/AfriSAR_AGB_Maps_1681___1/Rabi_AGB_50m.tif',
    'title': 'File to download',
    'hreflang': 'en-US'},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/metadata#',
    'href': 'https://daac.ornl.gov/AFRISAR/guides/AfriSAR_AGB_Maps.html',
    'title': 'ORNL DAAC Data Set Documentation (USERS GUIDE)',
    'hreflang': 'en-US'},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/metadata#',
    'href': 'https://dx.doi.org/10.3334/ORNLDAAC/1681',
    'title': 'Data set Landing Page DOI URL (DATA SET LANDING PAGE)',
    'hreflang': 'en-US'},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/metadata#',
    'href': 'https:/daac.ornl.gov/daacdata/afrisar/AfriSAR_AGB_Maps/comp/AfriSAR_AGB_Maps.pdf',
    'title': 'Data Set Documentation (GENERAL DOCUMENTATION)',
    'hreflang': 'en-US'},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/metadata#',
    'href': 'https:/daac.ornl.gov/daacdata/afrisar/AfriSAR_AGB_Maps/comp/AfriSAR_AGB_Maps_PlotDetails.csv',
    'title': 'Data Set Documentation (GENERAL DOCUMENTATION)',
    'hreflang': 'en-US'},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/data#',
    'href': 's3://nasa-maap-data-store/file-staging/nasa-map/AfriSAR_AGB_Maps_1681___1',
    'hreflang': 'en-US',
    'inherited': True},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/documentation#',
    'href': 'https://daac.ornl.gov/AFRISAR/guides/AfriSAR_Mondah_Field_Data.html',
    'hreflang': 'en-US',
    'inherited': True},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/metadata#',
    'href': 'https://doi.org/10.3334/ORNLDAAC/1580',
    'hreflang': 'en-US',
    'inherited': True},
   {'rel': 'http://esipfed.org/ns/fedsearch/1.1/documentation#',
    'href': 'https://daac.ornl.gov/daacdata/afrisar/AfriSAR_Mondah_Field_Data/comp/AfriSAR_Mondah_Field_Data.pdf',
    'hreflang': 'en-US',
    'inherited': True}],
  'title': 'AfriSAR_AGB_Maps.Rabi_AGB_50m.tif',
  'updated': '2019-05-13T12:45:07+00:00',
  'datetime': '2016-02-01T00:00:00Z',
  'time_end': '2016-03-31T23:59:59.000Z',
  'proj:bbox': [595552.0, 9784881.0, 600802.0, 9789931.0],
  'proj:epsg': 32732.0,
  'concept_id': 'G1200115789-NASA_MAAP',
  'dataset_id': 'AfriSAR: Aboveground Biomass for Lope, Mabounie, Mondah, and Rabi Sites, Gabon',
  'proj:shape': [101.0, 105.0],
  'time_start': '2016-02-01T00:00:00.000Z',
  'browse_flag': False,
  'data_center': 'NASA_MAAP',
  'granule_size': '0.014109',
  'proj:geometry': {'type': 'Polygon',
   'coordinates': [[[595552.0, 9784881.0],
     [600802.0, 9784881.0],
     [600802.0, 9789931.0],
     [595552.0, 9789931.0],
     [595552.0, 9784881.0]]]},
  'day_night_flag': 'BOTH',
  'proj:transform': [50.0,
   0.0,
   595552.0,
   0.0,
   -50.0,
   9789931.0,
   0.0,
   0.0,
   1.0],
  'original_format': 'ECHO10',
  'coordinate_system': 'CARTESIAN',
  'online_access_flag': False,
  'collection_concept_id': 'C1200115768-NASA_MAAP'},
 'geometry': {'type': 'Polygon',
  'coordinates': [[[9.859161573944016, -1.9460214484122889],
    [9.906363400458433, -1.9459965989464993],
    [9.906339295667408, -1.9003139362380275],
    [9.859138724101811, -1.900338201941798],
    [9.859161573944016, -1.9460214484122889]]]},
 'links': [{'rel': 'collection',
   'href': 'https://stac.maap-project.org/collections/AfriSAR_AGB_Maps_1681',
   'type': 'application/json'},
  {'rel': 'parent',
   'href': 'https://stac.maap-project.org/collections/AfriSAR_AGB_Maps_1681',
   'type': 'application/json'},
  {'rel': <RelType.ROOT: 'root'>,
   'href': 'https://stac.maap-project.org/',
   'type': <MediaType.JSON: 'application/json'>,
   'title': 'maap-stac'},
  {'rel': 'self',
   'href': 'https://stac.maap-project.org/collections/AfriSAR_AGB_Maps_1681/items/Rabi_AGB_50m',
   'type': 'application/geo+json'}],
 'assets': {'data': {'href': 's3://nasa-maap-data-store/file-staging/nasa-map/AfriSAR_AGB_Maps_1681___1/Rabi_AGB_50m.tif',
   'type': 'image/tiff; application=geotiff; profile=cloud-optimized',
   'roles': ['data']},
  'metadata': {'href': 'https://daac.ornl.gov/AFRISAR/guides/AfriSAR_AGB_Maps.html',
   'roles': ['metadata']},
  'documentation': {'href': 'https://daac.ornl.gov/AFRISAR/guides/AfriSAR_Mondah_Field_Data.html',
   'roles': ['documentation']}},
 'bbox': [9.859138724101811,
  -1.9460214484122889,
  9.906363400458433,
  -1.9003139362380275],
 'stac_extensions': ['https://stac-extensions.github.io/projection/v1.0.0/schema.json',
  'https://stac-extensions.github.io/raster/v1.1.0/schema.json'],
 'collection': 'AfriSAR_AGB_Maps_1681'}

Visualizing the read data with hvplot

[9]:

ds_masked = da.where(da != da._FillValue)

ds_masked.hvplot(
    'x', 'y',
    cmap='viridis',
    frame_height=400,
    frame_width=400
).redim.range(value=(0,da.max().values))

[9]:

[ ]:



[ ]:



[ ]: