I am creating NDVI from Sentinel2 bands 8 and 4 on QGIS, using the raster calculator, but I need to convert the format or clip to a polygon beforehand, because my computer crashes if I try to perform on the raw bands.
If I convert the bands to 32bit beforehand (needed for the end product anyway) or if I clip to polygon to make the raw band a smaller extent, I am finding the raster values change, and not in a consistent manner. These are some of the differences found between my 'processed' bands and the original bands, comparing my processed NDVI to the calculated NDVI from the raw bands.
I used GDAL warp (near resampling, LZW compression) to change the data type from UInt16 to 32b Float. I also trialed using GDAL translate (no compression) to do the same, and still found the values changed significantly and inconsistently.
Is this all to do with resampling? Is there any way to do this sort of processing on QGIS without altering my raster values?
Any advice would be greatly appreciated, I have searched stackexchange for similar questions, but haven't found anything which could solve my problem.