MUSCATE SENTINEL-2 L2A

[joalofe]

To whom may be interested:
I'm currently using the pyGeodes library to download some Sentinel-2 L2A (Muscate format) tiles. I'm trying to build a query for the "collection MUSCATE SENTINEL2 L2A". The query is as follows:
query = {
"grid:code": {"eq": "T33TVF"},
"start_datetime": {"gte": complete_datetime_from_str("2025-01-01")},
}
items,dataframe = geodes.search_items(query=query)
filtered_dataframe = dataframe[dataframe["collection"] == "MUSCATE_SENTINEL2_SENTINEL2_L2A"]
filtered_dataframe

The point is that I'm not able to establish a further filter based on the cloud cover percentage (even though I've seen in some other examples of your documentation that parameter applied to other collections, such as PEPS_S2_L1C). Is that parameter present in the collection I'm interested about? If not, how can I do this 'advanced' filter?

Moreover, regarding the date, I have two questions:

• Is it possible to do a query with a start_datetime and a end_datetime? I've tried to do it but I wasn't able to, that's why I decided to do a 'simpler' query like the one abovementioned.
• If not, is it possible to establish a logic in which I can search for the tiles based on the exact date? Something like:
  query = {
  "grid:code": {"eq": "T33TVF"},
  "start_datetime": {"eq": complete_datetime_from_str("2025-01-01")},
  }

Thank you in advance,
Jose Antonio.

[hfrcnes]

Hello @joalofe ,

You can use the following query to filter the data as needed (simply replace the values with your own):

query = {
    "grid:code": {"eq": "T33TVF"}, 
    "start_datetime": {"gte": complete_datetime_from_str("2025-01-01")}, 
    "end_datetime": {"lte": complete_datetime_from_str("2025-01-31")}, 
    "eo:cloud_cover": {"lte": 50}
}
items, dataframe = geodes.search_items(query=query, collections=["MUSCATE_SENTINEL2_SENTINEL2_L2A"])

We don't support query intervals at the moment, so you need to use both start_datetime and end_datetime.
I just checked, and the cloud cover percentage search works for me on this collection. If it's still not working for you, please share your code so I can take a further look.
We also added support for collection search, so you can simply add collections=["MUSCATE_SENTINEL2_SENTINEL2_L2A"] in geodes.search_items.

Regards,
Hugo

[joalofe]

Thank you very much, you have been very helpful. It seems to work.
If there were any other problems, I'll let you know

[joalofe]

By the way, is there any way to build the query with a .shp like it is done in the GEODES Portal? Or the only ways are by using 'geometry' or 'bbox'?

[hfrcnes]

Hello again,
No, we don't support that at the moment but you could extract the bounding box of your shapefile using the tools/framework of your choice and use it to search in 'bbox' with pyGeodes.