RmSmallCellChunks: Remove Small Cell Chunks

In inlmisc: Miscellaneous Functions for the USGS INL Project Office

Description

Remove small cell chunks from a raster layer, where a cell chunk is defined as a group of connected cells with non-missing values. The cell chunk with the largest surface area is preserved and all others removed.

Usage

RmSmallCellChunks(r)

Arguments

r	'RasterLayer'. Raster grid layer with cell values.

Value

An object of class 'RasterLayer' giving r with cell values in the smaller cell chunks set to NA.

Author(s)

J.C. Fisher, U.S. Geological Survey, Idaho Water Science Center

Examples

set.seed(2)
r <- raster::raster(ncols = 10, nrows = 10)
r[] <- round(runif(raster::ncell(r)) * 0.7)
r <- raster::clump(r)
r <- raster::ratify(r)
PlotMap(r)

r_new <- RmSmallCellChunks(r)
PlotMap(r_new)

graphics.off()

inlmisc documentation built on Jan. 25, 2022, 1:14 a.m.