sieve: Sieve raster data
In spatialEco: Spatial Analysis and Modelling Utilities

sieve

R Documentation

Sieve raster data

Description

Removes contiguous cells < specified query area

Usage

sieve(x, a, unit = c("m", "km", "ha"))

Arguments

`x`	An integer terra SpatRaster
`a`	Query area to remove
`unit`	The unit to use for area query options are c("m", "km", "ha")

Details

A sieve can be used to establish a minimal mapping unit where contiguous cells < specified query area are set to NA. These NA values can then be filled using focal (majority, median, mean)

Value

A terra SpatRaster with cells < a set to NA

Author(s)

Jeffrey S. Evans <jeffrey_evans@tnc.org>

Examples


library(terra)
elev <- rast(system.file("extdata/elev.tif", package="spatialEco"))
m <- matrix(c(100,200,1,200,300,2,300,400,3,400,
       500,4, 500,600,5), ncol=3, byrow=TRUE)
  x <- classify(elev, m)

# Sieve to a MMU of 60km
sv <- spatialEco::sieve(x, a = 60, unit = "km")
  plot(c(x, sv))

spatialEco documentation built on Nov. 18, 2023, 1:13 a.m.