SENTINEL1_GRD missing pixels

hi, just a heads up: I’m trying to reproduce, but there seems to be a problem with SentinelHub at the moment, so I can’t even get to download something

There was a short interruption of our services on that day for 1h. Everything is up again and should be working normally again. If you notice anything strange, please let us know.

It seems that problem is still there.

I can at most reproduce a line artifact around the -74 deg longitude meridian (north-west of that lake), but not the other lines you are seeing, e.g. in QGis:

Are we sure these lines you see aren’t just some kind of visualization issue? have you tried loading the results in another tool?

My spatial extent is bigger than yours so you can not see it.

Additionally,I can see the artifacts in your output too:

Yes that was what I was trying to explain:

  • I can only reproduce that artifact near the lake, which only spans half the vertical extent
  • I can not reproduce the other artifacts that span the whole vertical extent (I have the same spatial extent in my file, but zoomed in on the only artifact I could find)

I suspect these two types of artifacts have a different source: e.g. the first is a DEM issue and the second is just a visualization issue.

Can you share the netCDF you get, e.g. through email?

Please find a file here: worldwater-toolbox/openEO.nc at main · DHI-GRAS/worldwater-toolbox · GitHub

FYI - I was downloading another site and I can see the similar vertical line

{'west': -99.51913833618165, 'east': -99.37459945678712, 'south': 24.89967211146524, 'north': 24.961938324196538, 'crs': 4326}

I can see that my jobs are waiting in the line for hours. Is the server down?

Thanks for your version of the netcdf download, I can now also see the other artifact.

One think I noticed is that the the size in pixels is different for the same spatial extent:

  • your file 945x400 pixels
  • my result (based on your original snippet) is 1800x382 pixels

You also have 10 timestamps, while I only get 3 with the time range of you original post:

can you post an update of your code to make sure we are talking about the same use case?

Both production and dev instance of the VITO backend look fine at my end. What connection URL are you using in this case?

  • I just increase the spatial extent due to checking the consistency in data:
    Here the updated code:
start_date = '2021-12-01'
end_date = '2021-12-31'

s1_cube = connection.load_collection(
   'SENTINEL1_GRD',
    spatial_extent = {'west': -74.23921464628793, 'east': -73.59376786894418, 'south': 4.63822888776919, 'north': 4.775094862425593, 'crs': 4326},
    temporal_extent = [start_date, end_date],
    bands = ['VH','VV'],
    properties = {"polarization": lambda p: p == "DV"})


s1_cube.download("s1_cube.nc")
  • This is another place which I tested and results seems witht the same problem:
start_date = '2021-12-01'
end_date = '2021-12-31'


s1_cube = connection.load_collection(
    'SENTINEL1_GRD',
     spatial_extent = {'west': -74.23921464628793, 'east': -73.59376786894418, 'south': 4.63822888776919, 'north': 4.775094862425593, 'crs': 4326},
     temporal_extent = [start_date, end_date],
     bands = ['VH','VV'],
     properties = {"polarization": lambda p: p == "DV"})


s1_cube.download("s1_cuba.nc")

I think the spatial (and temporal) extent are the same as the other example

There you go: {'west': -99.5945352421805, 'east': -99.28125322924105, 'south': 24.888086962860328, 'north': 24.950359025232224, 'crs': 4326}

I keep failing to reproduce these artifacts you see.
For example, see this full notebook where I try to reproduce your last example.

I also can not find these line artifacts when I inspect the netCDF file in QGIS

@sulova.andrea Can you let us know how in which way you are visualizing those results when you get the artefacts. It looks to me very similar to this issue connected to the leaflet library.

I am using xarray to show data.
Here is full [notebook]

When we use download it works fine but when we use `job, it provides missing pixels in results (vertical lines)

Do you have access to the noteebook? @daniel.thiex @stefaan.lippens

I think @jeroen.dries recently fixed this issue, but I can’t seem to find a pointer to it (github commit or ticket) at the moment

1 Like

I can agree the problem was fixed based on my outputs :slight_smile:

1 Like