Geolocation issue with L2 Sentinel2 and Landsat8 products
Posted: Wed Feb 26, 2020 11:38 am America/New_York
Hi All,
I'm working on a process that requires,
1. L2 processing Sentinel2 and Landsat8 scenes through l2gen
2. Post-processing the Rrs outputs from step 1 outside seadas
3. create geolocated geotiffs from step 2 results
In order to create geolocated geotiffs in step 3 I use python Rasterio package. For that Rasterio takes the 4 boundary values (eastmost_lon, northmost_lat, ...). Now, the issue is that some of the scenes processed in step 1 through l2gen in the output have invalid boundary coordinates under
However,
I'd love to know if anyone has suggestion on how can I correctly create geolocated geotiffs from output of step 2.
Thanks,
Armin
I'm working on a process that requires,
1. L2 processing Sentinel2 and Landsat8 scenes through l2gen
2. Post-processing the Rrs outputs from step 1 outside seadas
3. create geolocated geotiffs from step 2 results
In order to create geolocated geotiffs in step 3 I use python Rasterio package. For that Rasterio takes the 4 boundary values (eastmost_lon, northmost_lat, ...). Now, the issue is that some of the scenes processed in step 1 through l2gen in the output have invalid boundary coordinates under
/Metadata/Global_Attributes
like northernmost_latitude 90.0 float32
southernmost_latitude -90.0 float32
easternmost_longitude 180.0 float32
westernmost_longitude -180.0 float32
geospatial_lat_units degrees_north ascii
geospatial_lon_units degrees_east ascii
geospatial_lat_max 90.0 float32
geospatial_lat_min -90.0 float32
geospatial_lon_max 180.0 float32
geospatial_lon_min -180.0 float32
However,
Metadata/Scan_Line_Attributes
seems to have array valid scan line coordinates under slat, elat, slon, elon, clat, clon
. If you pick the 4 extreme coordinates from here they do not form a rectangular bounding box. As a side note when I load the l2gen processed output in QGIS, the Extend
values, which seems to correspond to the 4 bounding coordinates does not agree with any of the Metadata/Scan_Line_Attributes
coordinates.I'd love to know if anyone has suggestion on how can I correctly create geolocated geotiffs from output of step 2.
Thanks,
Armin