我在美国有大约500,000个R 迁徙鸟类的发生数据.
我试图在这些点上叠加网格,然后计算每个网格中出现的次数.计算完计数后,我想将它们引用到网格单元ID.
在R中,我使用该over()函数来获取范围映射中的点,这是一个shapefile.
#Read in occurrence data
data=read.csv("data.csv", header=TRUE)
coordinates(data)=c("LONGITUDE","LATITUDE")
#Get shapefile of the species' range map
range=readOGR(".",layer="data")
proj4string(data)=proj4string(range)
#Get points within the range map
inside.range=!is.na(over(data,as(range,"SpatialPolygons")))
上面的工作正如我所希望的那样,但没有解决我当前的问题:如何处理类型的点SpatialPointsDataFrame,以及栅格的网格.您是否建议使用上面指出的相同方法对栅格网格进行多边形化？或者另一个过程会更有效率？