I'm trying to do a simple subtraction of 2 raster layers using QgsRasterCalculator in PyQGIS. I've searched around and thought I nailed it but the code below didn't yield anything, but there wasn't even a crash.
What could have gone wrong?
Each raster only contains one band.
target_ras_input = 'D:/Users/abc/Desktop/target.tif'target_ras = QgsRasterLayer(target_ras_input,"target_raster")QgsMapLayerRegistry.instance().addMapLayer(target_ras)initial_ras_input = 'D:/Users/abc/Desktop/initial.tif'initial_ras=QgsRasterLayer(initial_bathy_input,"initial_raster")QgsMapLayerRegistry.instance().addMapLayer(initial_ras)target=QgsRasterCalculatorEntry()target.raster = target_rastarget.bandNumber = 1target.ref = 'target_raster@1'initial=QgsRasterCalculatorEntry()initial.raster = initial_rasinitial.bandNumber = 1initial.ref = 'initial_raster@1'entries = [ target , initial ]final_output = 'D:/Users/abc/Desktop/final_raster.tif'calc=QgsRasterCalculator ( 'target.ref - initial.ref', final_output , 'GTiff', target_ras.extent(), target_ras.width(), target_ras.height(), entries )calc.processCalculation()final_raster=QgsRasterLayer(final_output, "please_work")QgsMapLayerRegistry.instance().addMapLayer(final_raster)