Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Sentinel-3 projection #353

Open
FedericOldani opened this issue May 30, 2024 · 5 comments
Open

Sentinel-3 projection #353

FedericOldani opened this issue May 30, 2024 · 5 comments

Comments

@FedericOldani
Copy link

Hi there,
I am trying to download Sentinel-3 LST images but I am stuck with the following error:
Cannot automatically compute the resolution, since asset 'lst-in' on item 0 'S3B_SL_2_LST_20230101T091518_20230101T091818_0179_074_264_2160' doesn't provide enough metadata to determine its native resolution.

This happens with any item in my tests. Should I compute the params in some way?
Thank you in advance

@FedericOldani FedericOldani changed the title Sentinel-3 resolution Sentinel-3 projection May 30, 2024
@TomAugspurger
Copy link

Can you share a code snippet that reproduces the probelm?

@FedericOldani
Copy link
Author

Here the code:

import pystac_client
import planetary_computer
import stackstac

catalog = pystac_client.Client.open(
    'https://planetarycomputer.microsoft.com/api/stac/v1',
    modifier=planetary_computer.sign_inplace,
)

bbox = (7.566238975059995, 44.323814339090575, 7.568329074631334, 44.326032959203175)

search = catalog.search(
    collections=['sentinel-3-slstr-lst-l2-netcdf'],
    bbox=bbox,
    datetime='2023-04-01/2024-05-12'
)
items = search.item_collection()

data = (
    stackstac.stack(
        items,
        assets=['lst-in'],
        bounds_latlon=bbox,
        rescale=False,
        xy_coords='center',
        properties=True,
        epsg=4326
        )
)

Note that without epsg param it returns the error:
ValueError: Cannot pick a common CRS, since asset 'lst-in' of item 0 'S3A_SL_2_LST_20230401T092051_20230401T092351_0179_097_150_2160' does not have one.

@TomAugspurger
Copy link

I'm not too familiar with sentinel-3 data, but do you know if those assets have a CRS? Looking at gdalinfo on the href from https://planetarycomputer.microsoft.com/api/sas/v1/sign?href=https://sentinel3euwest.blob.core.windows.net/sentinel-3/SLSTR/SL_2_LST___/2024/05/29/S3B_SL_2_LST____20240529T213622_20240529T213922_20240530T112857_0180_093_285_5760_PS2_O_NT_004.SEN3/LST_in.nc, I don't see anything that looks too promising.

@neoporteria
Copy link

The georeference of Sentinel-3 data comes in a different file (geodetic-in and cartesian-in), so it's necessary to extract the coordinates from there and georeference them with snap or using them as gcps

@777arc
Copy link
Contributor

777arc commented Nov 8, 2024

@FedericOldani would you consider this still an open issue or can we close?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants