nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Map Algebra in Python Using Variable

RustyRex

I am trying to move from raster calculator to the python window for my spatial analysis and moving from ERDAS Model Maker altogether. I am already struggling on only my first task. Here I am just trying to 'scale' my raster from 0-1

import arcpy from arcpy import env from arcpy.sa import * #Get Max and Min gridMaxResult = arcpy.GetRasterProperties_management("RASTER1", "MAXIMUM") gridMinResult = arcpy.GetRasterProperties_management("RASTER1", "MINIMUM") #Assign output gridMax = gridMaxResult.getOutput(0) gridMin = gridMinResult.getOutput(0)  #Rescale or 'normalize' OutRas = (((Raster("RASTER1") - gridMin) * (1 - 0)) / (gridMax - gridMin)) + 0 OutRas.save("NormalizedRas")

I get:
Runtime error <type 'exceptions.RuntimeError'>: ERROR 000732: Input Raster: Dataset 0 does not exist or is not supported

It is reading the variable correctly (gridMin is equal to 0) but it seems like it doesnt know it is a number. Any advice would help me with this problem and my syntax going forward.

Thanks

Find more posts tagged with

Accepted answers

curtvprice

What would the syntax of this code look like if i were to make it a python script to import as a script tool with input and output parameters?

The normal way to provide information to an arcpy script tool is to use GetParameterAsText(), setting up the input parameter as a Raster Layer, and the output as a Raster Dataset. Raster() is used to convert the input string to a raster data arcpy map algebra can work with.

import arcpy from arcpy.sa import * inRaster = arcpy.GetParameterAsText(0) # input raster layer name or raster dataset path outRaster = arcpy.GetParameterAsText(1) # output: raster dataset  ras = Raster(inRaster) range = ras.maximum - ras.minimum mappedras = (ras - ras.minimum) / range mappedras.save(outRaster)

Beyond that, you're kind of out of the scope of this thread.

All comments

Luke_Pinner

It is reading the variable correctly (gridMin is equal to 0) but it seems like it doesnt know it is a number

That's because it isn't a number.

The getOutput method returns a unicode string, see the help for the result object.

Cast your gridMax/Min variables to float before using them and the map algebra expression will work.

>>> type(gridMax)
<type 'unicode'>
>>> gridMax = float(gridMaxResult.getOutput(0))
>>> gridMin = float(gridMinResult.getOutput(0))

curtvprice

I am just trying to 'scale' my raster from 0-1

Three thoughts:

1) The Raster object has the properties you want without having to use that tool:

from arcpy.sa import *
rasPath = "RASTER1"
ras = Raster(rasPath)
print ras.minimum
print ras.maximum

2) I think your mapping expression needs some tweaking.

range = ras.maximum - ras.minimum
mappedras = (ras - ras.minimum) / range
mappedras.save("RASTER1_MAP") # save as permanent raster

3) The Slice tool can be very handy for this task depending on whether you need the resolution (it outputs integers so there may be some accuracy issues with this approach).

mappedras = Slice(ras, 1001, "EQUAL_INTERVAL", 0) * .001

RustyRex

Perfect, guys. Thanks also for the other tips. I kept my expression more complex than it needed to be so I could keep everything straight and in case I wanted to scale within a range not quite as simple as the 0-1 I have now. But for the code I will def. construct it the way you laid it out.

RustyRex

Thanks for the help. What would the syntax of this code look like if i were to make it a python script to import as a script tool with input and output parameters?

Thanks.

curtvprice

What would the syntax of this code look like if i were to make it a python script to import as a script tool with input and output parameters?

import arcpy from arcpy.sa import * inRaster = arcpy.GetParameterAsText(0) # input raster layer name or raster dataset path outRaster = arcpy.GetParameterAsText(1) # output: raster dataset  ras = Raster(inRaster) range = ras.maximum - ras.minimum mappedras = (ras - ras.minimum) / range mappedras.save(outRaster)

Beyond that, you're kind of out of the scope of this thread.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.