spreadlddzone

spreadlddzone

Shortest friction-distance path over map with friction from a source cell to cell under consideration, only paths in downstream direction from the source cell are considered

Result = spreadlddzone(ldd, points, initialfrictiondist, friction)
ldd
spatial ldd
initialfrictiondist
spatial, non spatial scalar
points
spatial nominal, ordinal, boolean
friction
spatial, non spatial scalar
Result
spatial points

Options

--unittrue or --unitcell

--unittrue
distance is measured in true distance (default)
--unitcell
distance is measured in number of cell lengths

Operation

The expression points identifies those cells from which the shortest friction-distance to every cell centre is calculated. The spreading for determination of these friction-distances starts at the centre of cells which have a non zero value on points. The initial friction distance (at the start of the spreading) is taken from the values at these point cells on initialfrictiondist. During spreading a path is followed over the consecutive neighbouring cells. While following this path the friction-distance increases. The increase of friction-distance per unit distance is specified by the cell values on friction. Using these values, increase when travelling from one cell to its neighbouring cell is calculated as follows: Let friction(sourcecell) and friction(destinationcell) be the friction values at the cell where is moved from and where is moved to, respectively. While moving from the source cell to the destination cell the increase of friction- distance is:

distance x {(friction(sourcecell)+friction(destinationcell)}/2

where distance is the distance between the sourcecell and the destination cell. This distance equals the cell length if the source cell and the destination cell are neighbours in horizontal or vertical directions; it equals sqrt(2) multiplied by the cell length if the cells are neighbours in diagonal directions.

During operation of the command, the spreading is executed from all non zero cells on points, over all possible paths. For determination of the friction-distance cell values on Result, for each cell the path from a non zero cell on points is chosen with the shortest friction-distance. So during the execution of the spreadlddzone operation, for each cell, the friction-distance for each possible path from the non zero cells on points to the cell under consideration is calculated and then the path with the shortest friction-distance is chosen. On Result each cell is assigned the points cell value of the cell where the shortest path to the cell begins.

Notes

The values on friction must be larger than zero.

Missing value cells on points, initialfrictiondist and friction are assigned a missing value on Result. Additionally potential shortest paths that cross missing value cells on points, initialfrictiondist or friction are ignored.

Group

This operation belongs to the group of Neighbourhood operators; spread operators

Examples

  1. • pcrcalc
    binding
    Result1 = Result1.map;
    Ldd2 = Ldd2.map;
    Points1 = Points1.map;
    initial
    report Result1 = spreadlddzone(Ldd2, Points1,0,1);

    • python
    Ldd2 = readmap(“Ldd2.map”)
    Points1 = readmap(“Points1.map”)
    Result1 = spreadlddzone(Ldd2, Points1,0,1)
    Result1.map Ldd2.map Points1.map
    _images/spreadlddzone_Result1.png _images/accu_Ldd2.png _images/spreadldd_Points1.png

  2. • pcrcalc
    binding
    Result2 = Result2.map;
    Ldd2 = Ldd2.map;
    Points2 = Points2.map;
    Initial = Initial.map;
    FrictMat = FrictMat.map;
    initial
    report Result2 = spreadlddzone(Ldd2, Points2,Initial,FrictMat);

    • python
    Ldd2 = readmap(“Ldd2.map”)
    Points2 = readmap(“Points2.map”)
    Initial = readmap(“Initial.map”)
    FrictMat = readmap(“FrictMat.map”)
    Result2 = spreadlddzone(Ldd2, Points2,Initial,FrictMat)
    Result2.map Ldd2.map Points2.map Initial.map FrictMat.map
    _images/spreadlddzone_Result2.png _images/accu_Ldd2.png _images/spreadldd_Points2.png _images/spreadldd_Initial.png _images/spreadldd_FrictMat.png

Table Of Contents

Previous topic

spreadldd

Next topic

spreadmax