Hi Luca,
From the explanation, I understand that you have a geotiff file corresponding to a Sentinel-2 (S2) observation. You want to read this file, modify some bands, and save the modified version as a new geotiff file. However, when you open the modified file in SNAP (ESA software), the organization of the optical observation is different, causing problems because SNAP doesn't recognize the bands with their original names.
To address this issue, you would like to find a way to read a geotiff file, make changes to it while preserving all the metadata (including band names) of the original observation, and then save it as a new geotiff file.
To achieve your goal of modifying the geotiff file while preserving all the metadata, including band names, you can follow these steps:
[A, R] = geotiffread('trial_subset.tif');
info = geotiffinfo('trial_subset.tif');
Make the desired modifications to the `A` array, which represents the pixel values of the bands.
Write the modified data to a new geotiff file using the geotiffwrite function. To preserve the metadata, including band names, you need to pass the `info` structure obtained from `geotiffinfo` as an argument:
geotiffwrite('trial_subset_modified.tif', A, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);
By following the above-mentioned steps, you can modify bands with the metadata preserved.
Thanks
