I'm trying to do Map Algebra on a bunch of rasters by iterating through a for loop and saving them in an empty geodatabase. However, they aren't saving properly. Only the last raster is created and it is named "rescaledRaster" rather than the outName. Here's the code:
import arcpy from arcpy.sa import * arcpy.env.workspace = r"[abspath]" arcpy.env.scratchWorkspace = r"[abspath]\rescale.gdb" aprx = arcpy.mp.ArcGISProject(r"[abspath]\filename.aprx") # Select correct interpolated rasters for lyr in aprx.listMaps().listLayers(): if lyr.isRasterLayer: if 'Temp' in lyr.name: if 'Int_Raster' in lyr.name: # Set local variables raster = Raster(lyr.dataSource) thresholds = [14, 23] bounds = [14, 23] scale = [1, 0.496714] transformation = TfLinear( bounds, bounds, thresholds, 1, thresholds, 0 ) outName = arcpy.env.scratchWorkspace + "\\" + raster.name.replace('Int_Raster','Resc') # Rescale the raster rescaledRaster = RescaleByFunction( raster, transformation, scale, scale ) # Save the rescaled raster rescaledRaster.save(outName) print(rescaledRaster.isTemporary) print(rescaledRaster.name)
The print statements seem correct each time; .isTemporary comes out False and .name returns the outName. Despite this, the output rasters keep getting overwritten.