-
Notifications
You must be signed in to change notification settings - Fork 12
Open
Description
I’ve noticed a discrepancy when querying the landsat-c2-l2 collection between Planetary Computer and AWS EarthSearch for a given area and time range. Is there any known reason for the missing scenes in the Planetary Computer catalog?
Please find the code bellow to reproduce this.
from pystac_client import Client
from shapely.geometry import box
import datetime as dt
aoi_bounds = (12.413501113658578, 55.657770784144326, 12.583960803598425, 55.71951501284911)
date_start = dt.datetime(2022, 8, 1)
date_end = dt.datetime(2022, 9, 2)
# 1/ planetarycomputer
import planetary_computer
pc_url = "https://planetarycomputer.microsoft.com/api/stac/v1"
pc_catalog = Client.open(pc_url, modifier=planetary_computer.sign_inplace)
pc_search = pc_catalog.search(
collections=["landsat-c2-l2"],
bbox=aoi_bounds,
datetime=f"{date_start:%Y-%m-%d}/{date_end:%Y-%m-%d}",
)
pc_items = pc_search.item_collection()
print(f"Planetary Computer Landsat scenes: {len(pc_items)}")
#Planetary Computer Landsat scenes: 7
# 2/ AWS element84
e84_url = "https://earth-search.aws.element84.com/v1"
e84_catalog = Client.open(e84_url)
e84_search = e84_catalog.search(
collections=["landsat-c2-l2"],
bbox=aoi_bounds,
datetime=f"{date_start:%Y-%m-%d}/{date_end:%Y-%m-%d}",
)
e84_items = e84_search.item_collection()
print(f"EarthSearch Landsat scenes: {len(e84_items)}")
# EarthSearch Landsat scenes: 14Metadata
Metadata
Assignees
Labels
No labels