Subset Sentinel-2 files?


I mean to try subsetting with the Sentinel-2 files (.SAFE) in python (not with GUI) to automate some of the work. Does anybody have any advice what to use or what not to use?

All the best

In addition of subsetting I would like to keep the original number of bands after the subsetting is done. With gdal_translate for example only 3 bands remain after subsetting.

You don’t necessarily need to use python. Maybe GPF/gpt can already do the job for you. Have a look in the SNAP help.
Either way, gpt or python you would need to use the Resample and the Subset operator.

In Handling sentinel data in python ( working with snappy) I show how to call the Resample operator in a next step you would need to call the Subset operator.