如何投影 sf 对象以使用精确提取器 R 中的多边形从栅格中提取单元格数据

Question

我需要提取县级栅格单元数据。

dd <- exact_extract(tmax,counties,c('mean','max'))

我正在使用精确的提取器包，它不断向我发出此警告：

Warning message:
In .exact_extract(x, y, fun = fun, ..., include_xy = include_xy,  :
Polygons transformed to raster CRS (EPSG:NA)

我的县文件的 CRS 为“+proj=longlat +ellps=WGS84 +no_defs”。我相信这也称为 EPSG 4326。我什至使用以下代码转换了县多边形文件，但仍然收到此错误。

new_sf_counties <- st_transform(sf_counties, crs = 4326)

这是我需要担心的警告吗？如果是这样，我该如何修复它？县形状文件可以从 https://www2.census.gov/geo/tiger/TIGER2019/COUNTY/ 下载。谢谢你。

Answer 1

对我来说，运行

crs(tmax)<-crs(counties)

以确保栅格数据的 crs 与运行之前的多边形的 crs 相同

dd <- exact_extract(tmax,counties,c('mean','max'))

使警告消息消失。