Hello, I have a massive list (~1million) of locations at given date and need to find matchup raw spectral bands of MODIS, MERIS, sentinel-3 OLCI sensors at the given date/location if available. some example locations below:
I can do this task by downloading the imagery that covers the location/date, and then use ESA SNAP software, however, it is impossible with such a massive dataset and the download will take months or years, also I don’t have space to store that much image. I tried the CREODIAS platform, but it is very slow to access its eodata currently and don’t have access to MODIS data yet.
After reading around I feel there maybe luck with SentinelHub, I probably need to use the Statistical API, but I did not find some example to start with. I am familiar with python, so if I can setup a local environment to do that would be great.
You are correct with the Statistical API option. If you already had a look at the provided examples, they may not be specific to your workflow, but they are a good starting point to start adjusting to fit your own workflow.
In this case you will need to construct Bounding Boxes from the coordinates, for example using a certain buffer size. Check this example of how to do this with geopandas.
Then you can start with this example to extract raw band values for a single band and single date. Also here are other examples in python . Remember to adjust the config to fit the data source i.e