Snappy subsetting problems for S3A data

Hello Marie,

If found the cause for this problem. OLCI RR Orbits can contain the north pole. And this wasn’t properly handled for all cases in our code (SNAP-759). This is fixed and will be released with the next update. Not sure when will be published, maybe in two week or so.

As a temporary workaround, you can do a pixel based subset first and afterwards a geo-region subset.
The pixel based subset might look like this:

from snappy import Rectangle
op1 = SubsetOp()
op1.setSourceProduct(p)
size = p.getSceneRasterSize()
op1.setRegion(Rectangle(0, 500, size.width, size,height))
pixelSubset = op1.getTargetProduct()

and in your second op you use pixelSubset instead of p

op.setSourceProduct(pixelSubset)

1 Like