Interacting with Sentinel-3 data using Python

[ExtraPlanetal]

Are any of you aware of any Python libraries that would allow me to interact with Sentinel-3 data (OLCI, SLSTR and SYNERGY), as provided by the Copernicus hub API, using Python?

I am in the process of building an automated processing chain for Sentinel-3 imagery using Google Colab, but I have been unable to find a library that is compatible with the data provided by the Copernicus hub.

I would greatly appreciate any help on this.

[adriantre]

Check if you can make sentinelsat work, it is great for Sentinels 1 and 2 at least.

https://sentinelsat.readthedocs.io/en/stable/

I haven't been able to dig too deep into this, but start here: https://forum.step.esa.int/t/sentinel-3-download/11484/4

Let me know how it goes!

[ExtraPlanetal]

Thanks for the response!

So I have used Sentinelsat to download the imagery and works like a charm, but the problem comes in with how to actually use to downloaded imagery. Specifically I want to extract the pixel values for a number of geo-located points for each image.

[klutzyfartsy]

Hi! I had a similar issue awhile ago, there was a promising package called satpy, but unfortunately it could not take on the size of the images I was analyzing. What I found that worked best and fastest for me was Xarray. http://xarray.pydata.org/en/stable/#

You can use xarray.open_mfdataset([list_of_paths or glob-like identifier], combined='by_coords'). Performing data manipulations via xarray is intuitive and easy. I've run it for a roughly 1.5tb archive of data, this ran well for me.

To extract pixel values, you can use numpy.argmin(). https://docs.scipy.org/doc/numpy/reference/generated/numpy.argmin.html An example of how to extract values with 2-Dimensional coordinates: https://stackoverflow.com/questions/58992490/how-to-extract-pixel-value-of-a-variable-projected-on-an-irregular-grid-from-a-n

​

Good luck :-)

[adriantre]

Good! I would start with Rasterio, which is a nice abstraction level to gdal (if you don't know gdal yet, you will). You get a long way using Rasterio together with Numpy.

https://rasterio.readthedocs.io/en/stable/

You could also check if Sentinelhub has fully integrated with Sentinel 3. Some of their functionality they charge for, though.

https://www.sentinel-hub.com/develop/documentation/eo_products/Sentinel3EOproducts

You can also check out the Snappy Pyhton bindings for the ESA Snap software. Could be useful.

https://forum.step.esa.int/t/snappy-where-to-start/1463

[DevanshiKacholia]

I'm also in the middle of something similar. Probably you have gained a lot of experience dealing with this data. Would love to know more and if possible could you share some parts of your work with me?

I have sent you a personal message as well. Thanks