topic Re: how to replace RasterCalculator with arcpy in ArcGIS Spatial Analyst Questions

how to replace RasterCalculator with arcpy

Gisbert61 — Sun, 24 Nov 2013 09:00:51 GMT

Hello, newbie here.

I'm combining two rasters into an output using the Rastercalculator tool and that work fine, got the tool into modelbuilder and that works great as well.
Now I need to do additional steps like looping through more files to repeat the process and prompting the user for some choices so I'm looking into python using arcpy to do the rastermath. I cannot find the syntax to do this:

[INDENT]arcpy.gp.RasterCalculator_sa("""(Int("2010cl3" / 10) * 10) + (Int("2013cl5" / 10))""","Y:/case13/first/scrap_db.gdb/outputras")[/INDENT]

in arcpy. I understand I should not use rastercalculator in python but instead use arcpy.sa, can someone give an example please?

Thanks in advance

Re: how to replace RasterCalculator with arcpy

XanderBakker — Mon, 25 Nov 2013 08:34:43 GMT

Hi Bert,

Try something like this:

import arcpy from arcpy.sa import *   # settings (paths to input and output data) rasname1 = "Y:/Path/To/Input/Data.gdb/2010cl3" # edit this rasname2 = "Y:/Path/To/Input/Data.gdb/2013cl5" # edit this rasnameout = "Y:/case13/first/scrap_db.gdb/outputras"   # grab a SA license if arcpy.CheckExtension("Spatial") == "Available":     arcpy.CheckOutExtension("Spatial") else:     exit("No SA license available...")   # make Raster objects for easy processing ras1 = arcpy.Raster(rasname1) ras2 = arcpy.Raster(rasname2)   # apply math rasout = (Int(ras1 / 10) * 10) + Int(ras2 / 10)   # save the result to disk rasout.save(rasnameout)

Kind regards,

Xander

Re: how to replace RasterCalculator with arcpy

Gisbert61 — Tue, 26 Nov 2013 00:24:05 GMT

Xander,

Got a result! thanks!