rmSingle {landscapeR} | R Documentation |
Patch creation algorithm can occasionally leave single cells scattered within patches. This function reduces the "salt-pepper" effect, identifying or correcting those cells.
rmSingle(rst, rm = TRUE)
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. |
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.
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)