nicholas_smale

Roads in Least Cost Analysis

Discussion created by nicholas_smale on May 28, 2011
Latest reply on May 30, 2011 by curtvprice
Hi!

I would like to include the roads layer in my least cost grid by creating regions that the path would be dissuaded from crossing, to be done by making all the cells above rasterised road have their cost value changed to '10' (on a 10 point scale).
I created this thread a while ago to ask this question, and it was answered but we got sidetracked and now the thread is dead. The solution provided is this:

The overwriting can be done with a combination of con and isnull, as in
Con( Isnull( [roads] ), [impedance], 10 )
Here, 10 is the value to put in the [impedance] grid at all road cells.

A similar approach is to add a large impedance at all road cells, as in
Not( Isnull( [roads] ) ) * 10 + [impedance]
This exploits the (documented) convention that false values are 0 and true values are 1.


Here is my problem:

I can't make the statements in your first post work, for instance 'Con( Isnull( [RasterRoads] ), [Cost], 10 )' with 'RasterRoads' being a rasterised road layer that is nodata for everything except road cells, and 'cost' being the cost layer already made up of slope and landuse.
I'm supposed to be using Single Output Map Algebra, right? Why does putting that statement into it not work? (the model thing remains white, it won't let me run it).


Could someone please help me, I am desperate :( :( :(

Outcomes