Leaflet绘图报错derivePolygons：地图着色展示Average_bleaching数据遇阻

这个derivePolygons报错的核心原因是数据类型与addPolygons的要求不匹配：addPolygons需要输入多边形结构数据（包含闭合的顶点经纬度集合），而你使用的点数据框或SpatialPointsDataFrame是单点数据，无法生成多边形，因此触发错误。以下是针对性的解决方法：

1. 若展示的是点数据的漂白值 #

如果你要标记的是采样点的平均漂白值，应该使用点渲染函数而非多边形函数，示例代码如下：

library(leaflet)

# 假设你的数据框包含lon（经度）、lat（纬度）、Average_bleaching列
leaflet(your_data) %>%
  addTiles() %>%
  addCircles(
    lng = ~lon,
    lat = ~lat,
    # 按漂白值设置颜色渐变
    color = ~colorNumeric("YlOrRd", domain = Average_bleaching)(Average_bleaching),
    radius = 10000, # 根据需求调整点大小
    popup = ~paste("平均漂白值:", round(Average_bleaching, 2))
  )

2. 若要做区域着色（按行政/海域展示平均漂白值） #

如果是要给某类区域（如海域分区、国家）着色展示平均漂白值，必须确保数据是多边形格式，步骤如下：

• 准备对应区域的多边形空间数据（如SpatialPolygonsDataFrame或sf对象）
• 将Average_bleaching字段关联到多边形数据中（确保每个多边形对应一个漂白值）
• 使用sf对象的示例代码：

library(leaflet)
library(sf)

# 假设polygons_sf是包含区域多边形和Average_bleaching的sf对象，坐标系为WGS84（EPSG:4326）
leaflet(polygons_sf) %>%
  addTiles() %>%
  addPolygons(
    fillColor = ~colorNumeric("YlOrRd", domain = Average_bleaching)(Average_bleaching),
    fillOpacity = 0.7,
    color = "#FFFFFF", # 多边形边界颜色
    weight = 1, # 边界粗细
    popup = ~paste("区域:", region_name, "<br>平均漂白值:", round(Average_bleaching, 2))
  )

3. 关键排查步骤 #

• 清理数据中的缺失值：用na.omit(your_data)移除含缺失经纬度或漂白值的行
• 确认空间数据的坐标系：leaflet默认要求WGS84（EPSG:4326），若坐标系不符，用st_transform(your_sf_object, 4326)转换

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