topic Re: raster operations in python in Python Questions

raster operations in python

RA1 — Sun, 25 Mar 2012 11:56:25 GMT

How can I sum two rasters and after reclassify the value field of the result according to the value of both input rasters ?

I can do the outPlus = Plus(inRaster1, inRaster2)
after I can do a loop for each row in outPlus

rows = arcpy.SearchCursor(outPlus)
fields=arcpy.ListFields(outPlus)
for row in rows:
for field in fields:
(...)

How do I grab the row value that are in inRaster1 and inRaster2 from inside the loop ?
I want to do something like: if (inRaster=1 and inRaster2=8): new.row.value="area under development"

Re: raster operations in python

FabianBlau — Mon, 26 Mar 2012 06:24:59 GMT

Hi!

Read this:

http://help.arcgis.com/de/arcgisdesktop/10.0/help/index.html#//000v0000012z000000

You can adress the cells via array-indices with the numpy-array.

Re: raster operations in python

curtvprice — Sun, 12 Dec 2021 06:34:06 GMT

How do I grab the row value that are in inRaster1 and inRaster2 from inside the loop ?
I want to do something like: if (inRaster=1 and inRaster2=8): new.row.value="area under development"

If you meant cell by cell, this is easily done using map algebra, no need to resort to numpy.

The following code creates a new raster "new_area" for all cells where inRaster1 = 1 and inRaster2 = 8:

from arcpy.sa import *
inRaster1 = "e:\\work\\myraster1.tif"
inRaster2 = "e:\\work\\myraster2.tif"
outRaster = Con(inRaster1 == 1 and inRaster2 == 8,1)
outRaster.save("e:\\work\\new_area")