使用plot()可见的Raster在Leaflet中无法显示的问题排查

针对你遇到的栅格在Leaflet中仅显示图例、不显示栅格内容的问题，可以尝试以下几个排查和解决步骤：

1. 明确配色函数的数值范围 #

你当前使用colorNumeric("Blues", domain=NULL)，domain=NULL可能因栅格中的NA值或特殊值导致配色范围识别错误。建议手动指定栅格的有效数值范围：

# 提取栅格的非NA数值作为配色范围
pal <- colorNumeric("Blues", domain = values(r_max, na.rm = TRUE))

# 重新构建地图代码
leaflet() %>% 
  addTiles() %>% 
  setView(lng = 15.6, lat = 28.1, zoom = 3.9) %>%
  addRasterImage(r_max, opacity=0.8, colors=pal) %>%
  addLegend(pal = pal, values = values(r_max, na.rm = TRUE), title = "Max", opacity=1)

2. 检查栅格数据的NA值占比 #

如果栅格中存在大量NA值，可能导致加载后看不到有效内容。可以先检查NA值数量：

sum(is.na(values(r_max)))

如果NA值占比过高，需要确认数据下载或预处理是否存在问题，或者调整地图视角到有有效数据的区域。

3. 手动转换栅格到Web墨卡托投影 #

Leaflet默认使用EPSG:3857（Web墨卡托）投影，虽然addRasterImage会自动转换投影，但手动转换可以避免潜在的适配问题：

library(raster)
# 转换投影
r_max_3857 <- projectRaster(r_max, crs = CRS("+init=EPSG:3857"))

# 使用转换后的栅格加载地图
leaflet() %>% 
  addTiles() %>% 
  setView(lng = 15.6, lat = 28.1, zoom = 3.9) %>%
  addRasterImage(r_max_3857, opacity=0.8, colors=pal) %>%
  addLegend(pal = pal, values = values(r_max, na.rm = TRUE), title = "Max", opacity=1)

4. 简化代码排查问题 #

先去掉自定义配色和图例，用最基础的方式加载栅格，确认是否能显示：

leaflet() %>% 
  addTiles() %>% 
  addRasterImage(r_max, opacity=0.8)

如果这样能显示，说明问题出在配色函数或图例的设置上；如果还是不显示，需要检查栅格文件本身的完整性。

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