rmSingle: Remove single tones from patches

In dariomasante/landscapeR: Categorical Landscape Simulation Facility

Description

Patch creation algorithm can occasionally leave single cells scattered within patches. This function reduces the "salt-pepper" effect, identifying or correcting those cells.

Usage

rmSingle(rst, rm = TRUE)

Arguments

rst	input raster landscape.
rm	logical, if TRUE returns the raster without single tones cells, if FALSE a vector of numbers identifying the index of single tones cells.

Value

A raster without single tones cells. If rm=FALSE, it returns a vector of numbers identifying the index of single tones cells. The value assigned to single tone cells is picked from one of the four neighbouring cells, selected at random.

Examples

library(raster)
m = matrix(0, 33, 33)
r = raster(m, xmn=0, xmx=10, ymn=0, ymx=10)
patchSize = 500

## Make a patch and introduce a single tone cell
r = makePatch(r, patchSize, spt=578, rast=TRUE)
r[578] = 0
plot(r)

## Now remove it
plot( rmSingle(r) )

## Single tones can be identified but not removed:
rmSingle(r, rm = FALSE)

dariomasante/landscapeR documentation built on June 5, 2021, 10:47 a.m.