Выборка случайных точек в R с ограничениями как растра, так и шейп-файла

Я пытаюсь выбрать случайные точки в соответствии с двумя условиями:

в не-NA ячейках растра
внутри полигонов шейп-файла

Я протестировал функцию terra::spatSample(), но не могу выполнить второе условие.

Вот воспроизводимый пример:

library(terra)
set.seed(1)
f <- system.file("ex/elev.tif", package = "terra")
f <- terra::rast(f)

v <- system.file("ex/lux.shp", package = "terra")
v <- vect(v)
v <- subset(v, v$NAME_2 != "Mersch")

test <- terra::crop(f, v, mask = TRUE)

sp <- terra::spatSample(test, 1000, method = "random", replace=FALSE, na.rm=TRUE, as.df=TRUE, as.points=TRUE, xy=TRUE, warn=TRUE, weights=NULL)

plot(test)
plot(v, add = TRUE)
plot(sp, add = TRUE)

r gis raster shapefile terra

19.06.2024 07:00

Стоит ли изучать PHP в 2026-2027 годах?

Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...

Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией

В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.

Приемы CSS-макетирования - floats и Flexbox

Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...

Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest

В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...

Концепция локализации и ее применение в приложениях React ⚡️

Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...

Пользовательский скаляр GraphQL

Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...

Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Один из вариантов — сначала преобразовать f в точку SpatVector, а затем intersect() результирующие точки с v. Таким образом, будут возвращены только точки, пересекающиеся с v. Затем вы можете sample() свои баллы. Обратите внимание, что использование set.seed() в этом репрексе означает, что sample() не будет «случайным».

library(terra)
library(tidyterra)
library(ggplot2)

set.seed(1)
f <- system.file("ex/elev.tif", package = "terra")
f <- rast(f)

# Convert f to points SpatVector and remove NAs
p <- as.points(f, na.rm = TRUE)

v <- system.file("ex/lux.shp", package = "terra")
v <- vect(v)
v <- subset(v, v$NAME_2 != "Mersch")

# Return p if inside v
p1 <- intersect(p, v)

# Return difference between p and p1 (for illustrative purposes)
p2 <- erase(p, p1)

# Sample p1
sp <- sample(p1, 1000)

# Plot
ggplot() +
  geom_spatraster(data = f) +
  geom_spatvector(data = v, alpha = 0.6) +
  geom_spatvector(data = p2, colour = "#920000") +
  geom_spatvector(data = sp) +
  coord_sf(xlim = c(5.8, 5.9),
           ylim = c(49.5, 49.6)) +
  scale_fill_continuous(na.value = "transparent") +
  scale_x_continuous(breaks = seq(5.8, 5.9, length.out = 3)) +
  theme(panel.background = element_blank())

@L Тайрон Отлично! Большое спасибо за Вашу помощь

— 19.06.2024 19:11

Исследовал relate "within" как индекс против intersect, но этот подход более простой.

— 19.06.2024 19:36

19.06.2024 08:59