r - 如何通过首先栅格化多边形来计算多多边形的物种丰富度？

Question

我有一个 4 个门的形状文件（每个类 4 个多边形）。每个都包含多个物种范围。意味着我在这 4 个多边形中有近 200 个物种的数据。我的数据看起来像这样（下图）。我从 IUCN 获得了这些数据，并将其剪辑到我的研究区域。数据集的链接是https://www.iucnredlist.org/resources/spatial-data-download

我想计算每个网格单元的物种丰富度（1 公里乘 1 公里），在 R 中对它们进行栅格化。表示每个网格单元中有多少种物种，如下图所示。

如果组合丰富度是不可能的，像爬行动物一样的单一丰富度也可以为我工作等等。我想在 R 中进行分析。

Answer 1

library(raster)
library(sf)

poly <- st_read("directory.shp")

# Creating a raster 
## Write the rows column of polygonn and xmn, xmx, ymn, ymx according to your country data    

r <- raster(nrows=, ncols=, xmn=, xmx=, ymn= , ymx=, 
        crs = "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0", 
        resolution = 0.008333333, vals=NULL)

a <- rasterize(poly, 
     r,
     field = "binomial",
     fun = function (x, ...) length(unique(na.omit(x))))