windowdiversity

windowdiversity

Number of unique values within a specified square neighbourhood

Result = windowdiversity(expression, windowlength)
expression
spatial boolean, nominal, ordinal
windowlength
spatial, non spatial scalar
Result
spatial scalar

Options

--unittrue or --unitcell

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

Operation

For each cell a square window with the cell in its centre is defined by windowlength. The windowlength is the length of the window in horizontal and vertical directions. For each cell on expression the number of unique values within its window is counted. This number is assigned to the corresponding cell on Result. Both cells on expression which are entirely in the window and cells which are partly in the window are considered.

Notes

The cell value on windowlength should be greater than 0, else a missing value is assigned to the corresponding cell on Result.

A cell on windowlength with a missing value results in a missing value on Result at the corresponding cell. However, if a missing value on windowlength occurs in the window in a cell which is not the centre cell of the window the expression value in that cell is considered for determination of the number of unique values in the window.

Group

This operation belongs to the group of Neigbourhood operators; window operators

Examples

  1. • pcrcalc
    binding
    Result2 = Result2.map;
    Expr = Expr.map;
    WinLen2 = WinLen2.map;
    initial
    report Result2 = windowdiversity( Expr, WinLen2);

    • python
    Expr = readmap(“Expr.map”)
    WinLen2 = readmap(“WinLen2.map”)
    Result2 = windowdiversity( Expr, WinLen2)
    Result2.map Expr.map WinLen2.map
    _images/windowdiversity_Result2.png _images/windowdiversity_Expr.png _images/windowaverage_WinLen2.png

  2. • pcrcalc
    binding
    Result3 = Result3.map;
    Expr2 = Expr2.map;
    initial
    report Result3 = (windowdiversity(Expr2,celllength() * 1.1)) > 1;

    • python
    Expr2 = readmap(“Expr2.map”)
    Result3 = (windowdiversity(Expr2,celllength() * 1.1)) > 1
    Result3.map Expr2.map
    _images/windowdiversity_Result3.png _images/windowdiversity_Expr2.png

  3. • pcrcalc
    binding
    Result1 = Result1.map;
    Expr = Expr.map;
    initial
    report Result1 = windowdiversity( Expr, 6);

    • python
    Expr = readmap(“Expr.map”)
    Result1 = windowdiversity( Expr, 6)
    Result1.map Expr.map
    _images/windowdiversity_Result1.png _images/windowdiversity_Expr.png

Table Of Contents

Previous topic

windowaverage

Next topic

windowhighpass