thanks for advicing me with the SH Process API tutorial. I followed along and managed to obtain the same results. In fact, now I have a TIFF file containing all the 13 bands (I checked for this with QGIS, and every georeferenced pixel contains information on all the bands).
I have a few other questions. Hope not to bother you.
Is more practical having a file for each band? If so, how can I save each band of a different TIFF file?
How can I extract information from a TIFF file? I’m reading more about the GDAL Python library. Am I on the right way? If not, could you advise me the best way to manage this?
In the tutorial we are " downloading raw digital numbers in the
INT16 format instead of reflectances in the
FLOAT32 format" for performances purpouses. Moreover, “the digital numbers are in the range from 0-10000, so we have to scale the downloaded data appropriately”.
- Does using integers affect accuracy so much? Is it better using FLOAT32 for research purposes?
In the tutorial, there is a parameter “resolution” used to dimension the output bounding box (right?).
- How does it relates to the resolution of the different bands?
I mean bands 2, 3, 4, 8 have a resolution of 10m, bands 1, 9, 10 have a resolution of 60m and the remaing a resolution of 20m. I was reading also this topic in the forum (10 m/px resolution for all Sentinel-2 bands : how?), but I cannot get the full sense of this.
It’s a lot of stuff, I know. Sorry!!
Thanks in advance!