library(terra)
f <- system.file("extdata/cyl.gpkg", package = "tidyterra")
v <- terra::vect(f)
# Add NAs
v <- v %>% mutate(iso2 = ifelse(cpro <= "09", NA, cpro))
# Init
plot(v, col = "red")
# Mask with lyr.1
v %>%
drop_na(iso2) %>%
plot(col = "red")
# SpatRaster method
# \donttest{
r <- rast(
crs = "EPSG:3857",
extent = c(0, 10, 0, 10),
nlyr = 3,
resolution = c(2.5, 2.5)
)
terra::values(r) <- seq_len(ncell(r) * nlyr(r))
# Add NAs
r[r > 13 & r < 22 | r > 31 & r < 45] <- NA
# Init
plot(r, nc = 3)
# Mask with lyr.1
r %>%
drop_na(lyr.1) %>%
plot(nc = 3)
# Mask with lyr.2
r %>%
drop_na(lyr.2) %>%
plot(nc = 3)
# Mask with lyr.3
r %>%
drop_na(lyr.3) %>%
plot(nc = 3)
# Auto-mask all layers
r %>%
drop_na() %>%
plot(nc = 3)
# }
Run the code above in your browser using DataLab