VeloxRaster_extract_points: Extract Values Given Points

In velox: Fast Raster Manipulation and Extraction

Description

Given a set of points, returns all raster values of the cells with which they intersect.

Arguments

sp	A SpatialPoints* object or a sf* POINT object.

Value

A numeric matrix. One row per element in sp, one column per band in the VeloxRaster.

Examples

## Make VeloxRaster with two bands
set.seed(0)
mat1 <- matrix(rnorm(100), 10, 10)
mat2 <- matrix(rnorm(100), 10, 10)
vx <- velox(list(mat1, mat2), extent=c(0,1,0,1), res=c(0.1,0.1),
            crs="+proj=longlat +datum=WGS84 +no_defs")
## Make SpatialPoints
library(sp)
library(rgeos)
coord <- cbind(runif(10), runif(10))
spoint <- SpatialPoints(coords=coord)
## Extract
vx$extract_points(sp=spoint)

velox documentation built on May 2, 2019, 6:01 a.m.