AnsweredAssumed Answered

Raster Calculator Expression used during iteration loop

Question asked by wtgeographer on Apr 22, 2014
Latest reply on May 23, 2014 by wtgeographer

Like • Show 0 Likes0
Comment • 2

I was wondering if someone can help with the following Raster Calculator Expression. I am trying to calculate the output least cost path raster to "null" where the end points overlap full least cost path raster. Additionally I am trying to use in -line variable substitution in the expression as the following snippet iterates 3 times.

I hope I have made sense. I appreciate any suggestions?

Entire function:

# Process: Cost Distance           arcpy.gp.CostDistance_sa(destValue, weighted_Raster, "distance_Raster", "", 'backlink_Raster')      # Process: Cost Path     Least_Cost_Path = arcpy.gp.CostPath_sa(originValue,"distance_Raster", "backlink_Raster", os.path.join(savepath + "\\Results.gdb" + "\\lcp_" + str(count)), "EACH_CELL", "OBJECTID")      global count      while count < 3:          # Local variables:         End_Pts = 'End_Pts'         Output_polyline_features = 'Output_polyline_features'         End_Pts_Raster = 'End_Pts_Raster'         Least_Cost_Route = 'Least_Cost_Route'         Reclass_rast1 = 'Reclass_rast1'          # Process: Raster to Polyline         arcpy.RasterToPolyline_conversion(Least_Cost_Path, Output_polyline_features, "ZERO", "0", "SIMPLIFY", "")          # Process: Feature Vertices To Points         arcpy.FeatureVerticesToPoints_management(Output_polyline_features, End_Pts, "BOTH_ENDS")          # Process: Point to Raster         arcpy.PointToRaster_conversion(End_Pts, "", End_Pts_Raster, "MOST_FREQUENT", "NONE", Least_Cost_Path)          # Process: Raster Calculator         arcpy.gp.RasterCalculator_sa("SetNull(~IsNull("%End_Pts_Raster%"), \"%Least_Cost_Path%\")", Least_Cost_Route)

Set null statement:

       arcpy.gp.RasterCalculator_sa("SetNull(~IsNull("%End_Pts_Raster%"), \"%Least_Cost_Path%\")", Least_Cost_Route)

curtvprice

Apr 24, 2014 6:33 PM

Correct Answer

Are you using ModelBuilder and just exporting the script for us to see?

If you're working Python, arcpy map algebra is a lot easier to work with.

For example instead of

arcpy.gp.RasterCalculator_sa("SetNull(~IsNull("%End_Pts_Raster%"), \"%Least_Cost_Path%\")", Least_Cost_Route)

the arcpy map algebra equivalent is something like this. Note how easy it is to save your rasters using the count variable to uniquely named grids with the .save() method.

from arcpy.sa import * while ....     ...     eptras = Raster(End_Pts_Raster)     LCostPth = Raster(Least_Cost_Path)     LCR = SetNull( ~ IsNull(eptras), LCostPth)     LCR.save("lcr{0}".format(count)     count -= 1

Desktop 10.1 Help: What is Map Algebra

See the reply in context

No one else had this question

Outcomes

from arcpy.sa import * while .... ... eptras = Raster(End_Pts_Raster) LCostPth = Raster(Least_Cost_Path) LCR = SetNull( ~ IsNull(eptras), LCostPth) LCR.save("lcr{0}".format(count) count -= 1

# Returns a value of 1 where End_Raster is NoData and 0 End_Raster points. isNull1 = IsNull(End_Raster) isNull1.save("isNull1_{0}".format(count)) #Returns LCP minus the end points # setNull1 = SetNull(isNull1,get_LCP(count),"Value = 0") setNull1.save("setNull1_{0}".format(count))

2 Replies

curtvprice Apr 24, 2014 6:33 PM
Unmark Correct
Correct Answer
Are you using ModelBuilder and just exporting the script for us to see?

If you're working Python, arcpy map algebra is a lot easier to work with.

For example instead of
```
arcpy.gp.RasterCalculator_sa("SetNull(~IsNull("%End_Pts_Raster%"), \"%Least_Cost_Path%\")", Least_Cost_Route)
```
the arcpy map algebra equivalent is something like this. Note how easy it is to save your rasters using the count variable to uniquely named grids with the .save() method.
```
from arcpy.sa import * while ....     ...     eptras = Raster(End_Pts_Raster)     LCostPth = Raster(Least_Cost_Path)     LCR = SetNull( ~ IsNull(eptras), LCostPth)     LCR.save("lcr{0}".format(count)     count -= 1
```
Desktop 10.1 Help: What is Map Algebra
Like • Show 0 Likes0
Actions

wtgeographer May 23, 2014 7:41 AM

Map algebra tools were used in my solution. Thanks Curtis! Here is what worked for me...

Actions

GeoNet

The Esri Community

Raster Calculator Expression used during iteration loop

Outcomes

Related Content

Recommended Content