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

   
   
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

   
   
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