I am using two GeoTIFFs to calculate NDVI.
Relating to bands 5 and 4 of Landsat 8, and using the raster calculator in QGIS, I import the TIFFs and use the formula (B5-B4)/(B5+B4).
The projection for the source TIFFs is EPSG:32611 - WGS 84 / UTM zone 11N The output file shares the same CRS.
The output GeoTIFF does not seem to reflect the same CRS, though, and calculates values beyond the bounds of the original data which creates extreme values (see screenshot).