directional

Result = directional(expression)
expression

spatial, non spatial boolean, nominal, ordinal, scalar, ldd

Result

dimension of expression directional

Options

--degrees or --radians

--degrees

values on expression are interpreted as degrees (default)

--radians

values on expression are interpreted as radians

Operation

If expression is a PCRaster map or an calculation resulting in a PCRaster map, it is converted: if expression is of one of the data types boolean, nominal, ordinal, scalar the cell values on expression are converted to the circular scale of Result, on a cell-by-cell basis. If expression is of the data type ldd, the codes on expression representing local drain directions are converted to real directions of drainage and saved as Result. The directions are converted to the circular scale of Result clockwise, assigning 0 degrees (or radians if the option --radians is set) to cells with a local drainage towards the top of the map (ldd code 8). A flat cell (cell value 5 on expression) is assigned a -1 on Result. The command can also generate a map of directional data type with one constant value.

If expression has no PCRaster data type, a Result with data type directional is generated. This is the case if expression is a number or a calculation with numbers. The value of expression must be in the domain of the directional data type, i.e. if the option --degrees is set: equal to 0 or between 0 and 360; if the option radians is set equal to 0 or between 0 and 2pi. Result will be a map with the same location attributes as the global clone map; all cells will have the value of expression.

Notes

A cell with missing value on expression is assigned a missing value on Result.

Group

This operation belongs to the group of Conversion and assignment

See Also

Data types, Directional data type, Numbers and data types

Examples

  1. • pcrcalc
    #! –degrees
    binding
    Result2 = Result2.map;
    Expr = Expr.map;
    initial
    report Result2 = directional(Expr);

    • python
    setglobaloption(“degrees”)
    Expr = readmap(“Expr.map”)
    Result2 = directional(Expr)

    Result2.map

    Expr.map

    ../../_images/directional_Result2.png ../../_images/boolean_Expr.png

  2. • pcrcalc
    binding
    Result1 = Result1.map;
    Expr = Expr.map;
    initial
    report Result1 = directional(Expr);

    • python
    Expr = readmap(“Expr.map”)
    Result1 = directional(Expr)

    Result1.map

    Expr.map

    ../../_images/directional_Result1.png ../../_images/directional_Expr.png