dem_crop: Crop DEM around points

In habtools: Tools and Metrics for 3D Surfaces and Objects

Crop DEM around points

Description

A function for sampling a DEM by cropping squares of a given size around xy coordinates.

Usage

dem_crop(data, x0, y0, L, plot = FALSE)

Arguments

data	A DEM in RasterLayer format.
x0	A value or vector of central x coordinate(s).
y0	A value or vector of central y coordinate(s).
L	Size of squares to cropped from the DEM.
plot	Logical. Plot the DEM and the cropped sections?

Value

A cropped RasterLayer or list of RasterLayers.

Examples

# around one point
dem_cropped <- dem_crop(horseshoe, -468, 1266, L = 2)
raster::plot(dem_cropped)
points(-468, 1266)

# around multiple points
points <- data.frame(x = c(-467, -465, -466), y = c(1270, 1265, 1268))
dem_list <- dem_crop(horseshoe, points$x, points$y, L = 1, plot = TRUE)

# plot the first element
raster::plot(dem_list[[1]])

habtools documentation built on May 29, 2024, 7:30 a.m.