Most occurring cell value within a specified square neighbourhood
Result = windowmajority(expression, windowlength)
spatial, non spatial scalar
spatial boolean, nominal, ordinal
spatial type of expression
--unittrue or -unitcell
windowlength is measured in true length (default)
windowlength is measured in number of cell lengths
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 most often occurring cell value within its window is determined and 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. At a cell, if two or more values occur the same (largest) number of times in its window, the windowlength at that cell is increased with two times the celllength of expression. The most often occurring cell value in this enlarged window is assigned to the cell under consideration. If this enlarged window still does not result in one most often occurring value, the window is progressively enlarged (with steps of two times the cellsize of expression) until a most often occuring cell value is found in the window.
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 which is not the centre cell of the window the value on expression in that cell is included in the computation of the average of the cell values in the window.
This operation belongs to the group of Neigbourhood operators; window operators
- • pcrcalcbindingResult2 = Result2.map;Expr = Expr.map;WinLen2 = WinLen2.map;initialreport Result2 = windowmajority( Expr, WinLen2);• pythonExpr = readmap(“Expr.map”)WinLen2 = readmap(“WinLen2.map”)Result2 = windowmajority( Expr, WinLen2)
- • pcrcalcbindingResult1 = Result1.map;Expr = Expr.map;initialreport Result1 = windowmajority( Expr, 6);• pythonExpr = readmap(“Expr.map”)Result1 = windowmajority( Expr, 6)