Рассчитайте общее количество точек выборки в каждой ячейке пространственной сетки.
Я хотел бы сделать сетку и рассчитать общее количество точек выборки в каждой ячейке сетки. Я создал случайно сгенерированные данные и сетку и попытался вычислить количество записей в ячейках сетки, используя как sf, так и растровые пакеты, используя предыдущие аналогичные вопросы SO, но безуспешно. Я также изучил функцию извлечения. Я новичок в пространственном анализе.
library(sf)
library(raster)
library(tidyverse)
library(mapview)
library(mapedit)
#Trial with sf package
# load some spatial data. Administrative Boundary
#https://stackoverflow.com/questions/41787313/how-to-create-a-grid-of- spatial-points
aut <- getData('GADM', country = 'aut', level = 0)
aut <- st_as_sf(aut)
#Try with polygons
grid <- aut %>%
st_make_grid(cellsize = 0.5, what = "polygons") %>%
st_intersection(aut)
#fake data
lat<-runif(1000, 46.5, 48.5)
lon<-runif(1000, 13,16)
pos<-data.frame(lat,lon)
ggplot() +
geom_sf(data = aut) +
geom_sf(data = grid)+
geom_point(data=pos, aes(lon, lat))
#how to count number of records within each cell?
########################################
#Trial with raster package
#https://stackoverflow.com/questions/32889531/r-how-can-i-count-how- many-points-are-in-each-cell-of-my-grid
r<-raster(xmn=13, ymn=46.5, xmx=16, ymx=48.5, res=0.5)
r[] <- 0
#How do I use the pos data here
xy <- spsample(as(extent(r), 'SpatialPolygons'), 100, 'random')
tab <- table(cellFromXY(r, xy))
r[as.numeric(names(tab))] <- tab
plot(r)
points(xy, pch=20)
d <- data.frame(coordinates(r), count=r[])
Я хочу получить таблицу с количеством точек отбора проб.