Skip to content

Accessing CCMEO STAC API

Examples on how to use python libraries to discover the collections from the CCMEO datacube STAC API.

Using pystac-client

pip install pystac-client

The developers of the pystac-client library provide multiple examples for interacting with STAC API. Please refer to those and consult the pystac-client documentation for additional information.

For certificate error while accessing the STAC API, see using custom certificate documentation

List collections

In this example you will learn how to query a STAC API with pystac-client to get the available collections

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
# Import the library needed
import pystac_client

# Link to ccmeo datacube stac-api
stac_root = "https://datacube.services.geo.ca/stac/api"
# Initialize the STAC client
catalog = pystac_client.Client.open(stac_root)

# Get the list of available collection id
for collection in catalog.get_collections():
    print(collection.id)

List items

In this example you will learn how to :

  • Query a STAC API with pystac-client with a portion of filter;
  • Access the data href for the specific filter.

Info

This specific example uses the collection hrdem-mosaic-2m from CCMEO's datacube

For more details on search parameters available with pystac-client, see items search documentation. 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
# Import the library needed
import pystac_client

# Link to ccmeo datacube stac-api
stac_root = "https://datacube.services.geo.ca/stac/api"
# Initialize the STAC client
catalog = pystac_client.Client.open(stac_root)

# Define search parameters
# In this example, we filter on the mosaic collection of HRDEM tiles
collections = ['hrdem-mosaic-2m']
# We define a bounding box extent
extent = [-79.28229773059192, 44.31501485755303, 
          -79.1702187573089, 44.3929540402247]

# Build the search query with the filtering parameters
search = catalog.search(
    collections=collections, 
    bbox=extent
    ) 

# Use the pagination to improve efficiency
links=[]
# Iterate over each returned page and get the dtm COG link
for page in search.pages():
    for item in page:
        links.append(item.assets['dtm'].href)