래스터 셀 내 shapefile의 카운트 포인트

Question

R을 사용하여 래스터 파일의 각 셀 내에있는 셰이프 파일에서 포인트 수를 계산할 수 있습니까?

아이디어는 래스터 셀의 원래 값을 유지하면서 래스터의 각 셀에 대해 하나의 행과 해당 셀 내에있는 점 수의 열을 포함하는 데이터 프레임을 검색하는 것입니다.

Answer 1

이 제목에 질문에 대한 답에 StackOverflow에

library(rgdal)   # this package to read and manipulate shapefiles 
library(raster)   # this package for rasters 
shp <- readOGR("my_shape_file.shp")  # read in your points from the shape file 
ras <- raster("my_raster_file")   # read in your raster 
shp <- spTransform(shp, projection(ras)) # make sure that the two spatial objects share the same coordinate system 
cells <- cellFromXY(ras,shp)  # find which cells the points are in 
table(cells)      # and make a table of number of occurences in each cell 

Answer 2

library(raster) 
# example raster and points data 
r <- raster(ncols=10, nrows=5) 
n <- 100 
x <- runif(n) * 360 - 180 
y <- runif(n) * 180 - 90 
xy <- cbind(x, y) 

rp <- rasterize(xy, r, fun=function(x,...)length(x)) 
as.data.frame(rp)