在R语言Leaflet中设置插值栅格投影及解决wkt函数报错

我碰到过好几次类似的问题——你用plot()能正常渲染插值后的栅格，但一放到Leaflet里就抛出wkt()函数找不到的错误，这大概率是空间依赖包没加载全，或者栅格的CRS（坐标参考系统）没有被正确识别导致的。咱们一步步来解决：

问题根源分析 #

1. 缺失关键依赖包：wkt()函数是sf包提供的核心函数之一，Leaflet调用addRasterImage()时需要它来处理投影转换，如果没加载sf包就会触发这个错误。
2. 栅格对象定义不完整：你的原代码里grid_raster_small变量没有定义就直接传给了interpolate()，这会导致插值后的栅格缺少完整的空间属性，进一步引发投影识别问题。
3. 网格点CRS未提前设置：在创建krig_grid_small后没有立即设置CRS，后续给插值栅格设置CRS时可能出现不兼容的情况。

修正后的完整可运行代码 #

# 加载所有必需的包（尤其不要漏掉sf）
library(leaflet)
library(raster)
library(fields)
library(sp)
library(sf)

# 生成模拟数据
ozone_df <- data.frame(
  Longitude = runif(200, -112.245075*10000, -111.455581*10000)/10000,
  Latitude = runif(200, 40.063614*10000, 40.827281*10000)/10000,
  Ozone = runif(200, 0, 115)
)

# 创建网格刻度与节点
small_grid_x <- seq(-111.455581, -112.245075, length.out = 500)
small_grid_y <- seq(40.063614, 40.827281, length.out = 500)
krig_grid_small <- expand.grid(small_grid_x, small_grid_y)

# 转为空间点并提前设置CRS（关键步骤）
coordinates(krig_grid_small) <- ~ Var1 + Var2
proj4string(krig_grid_small) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")

# 将空间点转为栅格对象（补充原代码缺失的部分）
grid_raster_small <- rasterFromXYZ(as.data.frame(krig_grid_small))

# 克里金插值
krig_fit_small <- fields::Krig(ozone_df[1:2], ozone_df$Ozone)
ozone_krig_small <- raster::interpolate(grid_raster_small, krig_fit_small)

# 为插值后的栅格设置CRS
crs(ozone_krig_small) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")

# 用plot验证
plot(ozone_krig_small)

# 在Leaflet中加载栅格
leaflet() %>%
  addTiles() %>%
  addRasterImage(ozone_krig_small, project = TRUE) %>%
  setView(lng = -111.941004, lat = 40.610497, zoom = 10) %>%
  addMiniMap()

关键修正点说明 #

• 强制加载sf包：直接解决wkt()函数找不到的问题；
• 补充grid_raster_small的定义：把空间网格点转为Leaflet能识别的栅格对象；
• 提前给krig_grid_small设置CRS：确保后续插值过程中空间属性的一致性，避免投影转换时出现异常。

运行这段代码后，Leaflet应该能正常加载你的臭氧插值栅格了。

内容的提问来源于stack exchange，提问作者bart cubrich