I have been using a simple gpt script (see below) using BandSelect to write S2 (level 2A products) bands as geotiff files. My output is two multiband geotiff files for 10m and 20m bands at their native resolution for analysis in R. I assume I do not need to resample as I am keeping these original resolutions? extract_to_tif_files.xml (1.5 KB)
The resulting 10m band file loads perfectly as a R raster. However the 20m band file fails to load properly - showing as a double image. Running gdalinfo suggests the tif metadata is incorrect, the 20m band file having the same dimensions and resolution as the 10m band file. gdalinfo_geotif_20m_output.txt (2.1 KB)
I have tried exporting as Geotif-Bigtif with slightly different results - correct dimensions in terms of total pixels for the 20m file, but incorrect resolution (10m) resulting in the geography being wrong.
The 20m tif files do reload in SNAP but geography/projection is altered from original.
This seems such a simple task, I can’t understand what I’m doing wrong?? Is there a more reliable way of accomplishing the export? Many thanks.