有没有一种方法可以通过 NA 空间进行栅格渐变？旨在模拟 ArcGIS GWR

Question

我正在尝试模仿您从 ArcGIS 地理加权回归中获得的栅格样式输出，而不是简单地绘制与坐标点对齐的栅格单元格中的系数/值，分配的值离开所有中间单元格 NA，它（至少看起来）用渐变填充它们。有关我的目标，请参阅 ArcGIS 屏幕截图。谢谢！

这是一个代表

library(tidyverse)
library(sf)
library(stars)

df <- quakes %>%
  select(depth, long, lat) %>%
  st_as_sf(coords = c("long", "lat"),
           crs = 4326)

df_rast <- st_rasterize(df)

ggplot() +
  geom_stars(data = df_rast)

ggplot() +
  geom_stars(data = df_rast,
             na.action = na.omit)