Best way to get raw data in Python

Apps / Services Python SDK
3

Hi @william.ray,

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.

  1. Is more practical having a file for each band? If so, how can I save each band of a different TIFF file?

  2. 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”.

  1. 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?).

  1. 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 (https://shforum.sinergise.com/t/10-m-px-resolution-for-all-sentinel-2-bands-how/2472), but I cannot get the full sense of this.

It’s a lot of stuff, I know. Sorry!!

Thanks in advance!