sieve {spatialEco}

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 <[e-mail]>

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))