I am trying to recreate gdal_calc for simple calculation using below code. import numpy import rasterio def merge_bands(calculationType, band1Image, band2Image, outputFilePath): # Read raster bands di

Hi Sean, Sure, I will try those changes you have mentioned. here is the gdal_calc expression I used : gdal_calc.py -A T1_B5.TIF --A_band=1 -B T1_B4.TIF --B_band=1 --outfile=ndvi_gal_1.TIF --calc='(A.a