extract2: Get zonal statistics from raster

In kongdd/extract2: What the Package Does (One Line, Title Case)

Description

Get zonal statistics from raster

Usage

## S4 method for signature 'Raster'
extract2(r, geoms, fun = sf_mean, weight = TRUE, ...)

## S4 method for signature 'RasterStack'
extract2(r, geoms, fun = sf_mean, weight = TRUE, ...)

## S4 method for signature 'list'
extract2(r, geoms, fun = sf_mean, weight = TRUE, ...)

## S4 method for signature 'character'
extract2(r, geoms, fun = sf_mean, weight = TRUE, ...)

Arguments

r	Raster, RasterBrick or character (raster paths) object
geoms	WKB object, e.g. sf::st_as_binary(sf::st_geometry(basins), EWKB=TRUE)
fun	one of sf_mean(), sf_median(), sf_sum(), sf_sd(), sf_var()
weight	boolean. If true, fraction = fraction * area
...	other parameters to fun

See Also

exactextractr::exact_extract(), sf_func()

Examples

## Not run: 
library(raster)
st <- sf::read_sf("data-raw/shp/Continents.shp")
shp <- rgdal::readOGR("data-raw/shp/Continents.shp")
files <- "data-raw/PML2_yearly_static2017-01-01.tif"
b <- raster::brick(files[1])
geoms = sf::st_as_binary(sf::st_geometry(st), EWKB = TRUE)
r1 = extract2(files, st)

geoms <- overlap(b, st)
r1 = extract2(files, geoms)
r2 = extract2(files, geoms, weight = TRUE)

r3 = extract2(files, geoms)
r4 = extract2(files, shp)

r = brick(files[1]) %>% readAll()
extract2(r, geoms, fun = sf_median)
# extract2(r, geoms, fun = sf_max)
# extract2(r, geoms, fun = sf_min)
extract2(r, geoms, fun = sf_sd)
extract2(r, geoms, fun = sf_sum)
extract2(r, geoms, fun = sf_sum, weight = TRUE)
extract2(r, geoms, fun = sf_var, weight = TRUE)


## End(Not run)

kongdd/extract2 documentation built on May 18, 2021, 5:50 p.m.