Reprojecting S3B-OLCI Level-2 Reduced Resolution to GeoTIFF

[jvaldezch]

I'm downloading some IOP from different sensor (Aqua, Terra, Virrs and Sentinel L2) data from your website and so far I managed to properly reproject Aqua, Terra and Virrs raster information and in those cases I have use GDAL (GDAL Warp) to do it. But when I try to do it wth Sentinel-3 Level-2 Reduced Resolution A and B it fails. It seems that metadata is different that data contained.

I'm missing a file or which is the projection coordinates used for this files?

Any help will be appreciated.

Code: Select all

file_in = "S3B_OLCI_ERRNT.20250101T164657.L2.IOP.NRT.nc"

band = xr.open_rasterio(
    f"NETCDF:{file_in}:/geophysical_data/bbp_443"
)

printr(band.rio.crs)

band.rio.write_crs("EPSG:4326 ", inplace=True)
band_reprojected = band.rio.reproject("EPSG:4326")
band_reprojected.rio.to_raster(file_tmp)

Take a look at coordinates in QGIS 3.4 (image below).

qgis_3_4.png (110.06 KiB) Not viewed yet

[OB SeaDAS - dshea]

Nothing about this file jumps out as a problem. Seadas had no problem displaying this L2 file.

S3B_OLCI_ERRNT_20250101T164657_L2_IOP_NRT_a_443.png (26.42 KiB) Not viewed yet

Maybe GDAL has a problem with full orbit files that cross the dateline.

[jvaldezch]

Got it! So subsetting the image first might help. Thanks!