如何使用 R 提取重叠栅格单元所在的多边形？

如何解决如何使用 R 提取重叠栅格单元所在的多边形？

我有一个高分辨率的温度和降水数据栅格层，我想将其覆盖在美国各州的 shapefile 上，并找出栅格层中的每个网格单元出现在哪个州。所以，如果网格单元与北达科他州、南达科他州和明尼苏达州的边界重叠，我需要能够提取该信息。

状态shapefile可以从https://www2.census.gov/geo/tiger/TIGER2020/STATE/下载，光栅文件可以从https://www.northwestknowledge.net/metdata/data/permanent/1980/下载。否则，我一直在尝试使用 sp::over，但这不起作用，并尝试使用 rgeos::

library(raster)
library(rgeos)
library(rgdal)

states <- readOGR('/tl_2020_us_state.shp')
states <- subset(states,!NAME %in% c('Commonwealth of the northern Mariana Islands','Hawaii','Alaska','United States Virgin Islands','Puerto Rico','American Samoa','Guam'))
states <- spTransform(states,"+proj=longlat +datum=wgs84 +no_defs")

burning_index <- raster('/permanent_gridmet_19800101.nc')
states_over_raster <- sp::over(burning_index,states)  # does not work
states_over_raster <- rgeos::overGeomGeom(burning_index,states)  # does not work

似乎 sp::over 仅适用于点，但事实并非如此。我一直在寻找使用栅格数据执行此操作的方法，但找不到任何内容。使用 rgeos 包中的 overGeomGeom 也会失败。对此的任何帮助都会很棒。一如既往，谢谢你。