R语言Leaflet中addPolygons()报错：metaData无适用方法求助

核心错误原因 #

你遇到的Error in UseMethod("metaData")报错，根源是**addPolygons中传入了单独的空间几何列（sfc对象）**，Leaflet的sf适配接口需要接收完整的sf数据框，而非单独的几何列。此外代码里还有几处变量引用和语法错误，一起修正如下：

具体修正点 #

1. 移除addPolygons的data参数：你已经通过管道%>%把merged_df_rep传给了leaflet，不需要再单独指定data = merged_df_rep$geometry，否则会传入孤立的sfc对象触发报错。
2. 修正fillColor的变量引用：原代码里的merged_df_us_1$n_response是不存在的变量，替换为你的目标列merged_df_rep$values（或直接用~values，因为管道已传入数据框）。
3. 修复setView的语法：原代码的参数格式错误，直接传入经纬度数值即可，不需要波浪线。
4. 优化弹窗内容：用htmlEscape处理州名，避免特殊字符导致的弹窗异常，同时用公式接口更简洁。

修正后的完整代码 #

library(ggmap)
library(tidygeocoder)
library(tigris)
library(sf)
library(htmltools)
library(labelled)
library(leaflet) 
library(tidyverse)

# 构建自定义数据
NAME <- c("West Virginia", "Florida", "Illinois", "Minnesota","Maryland", "Rhode Island", "Idaho", "New Hampshire","North Carolina",       "Vermont", "Connecticut", "Delaware","New Mexico", "California", "New Jersey", "Wisconsin", "Oregon", "Nebraska","Pennsylvania","Washington",           "Louisiana", "Georgia", "Alabama" , "Utah","Ohio","Texas", "Colorado", "South Carolina","Oklahoma" , "Tennessee", "Wyoming", "Hawaii", "North Dakota", "Kentucky",  "Maine" , "New York", "Nevada" ,  "Alaska" , "Michigan", "Arkansas","Mississippi","Missouri", "Montana" ,"Kansas", "Indiana", "South Dakota","Massachusetts", "Virginia", "District of Columbia", "Iowa","Arizona")
values <- c(1000:1050)
df_rep <- data.frame(NAME, values)

# 获取并清理州空间数据
map_df <- states()
map_df_clean <- map_df[!(map_df$NAME %in% c("United States Virgin Islands", "Commonwealth of the Northern Mariana Islands", "Guam", "Puerto Rico", "American Samoa")),]

# 合并数据并转换坐标系
merged_df_rep <- left_join(map_df_clean, df_rep, by = "NAME") %>%
  st_transform("+proj=longlat +datum=WGS84")

# 创建配色方案
pal <- colorNumeric("Greens", domain = merged_df_rep$values)

# 构建弹窗（优化版本）
popup_sb <- ~paste0("Total values for ", htmlEscape(NAME), ": ", as.character(values))

# 绘制Leaflet地图
merged_df_rep %>% 
  leaflet() %>% 
  addProviderTiles("CartoDB.Positron") %>% 
  setView(lng = -98.483330, lat = 38.712046, zoom = 4) %>% 
  addPolygons(
    fillColor = ~pal(values),
    fillOpacity = 0.7,
    weight = 0.2, 
    smoothFactor = 0.2,
    popup = popup_sb
  ) %>% 
  addLegend(
    pal = pal, 
    values = ~values,
    position = "bottomright",
    title = "Responses"
  )

额外优化建议 #

• 清理州数据时用%in%替代多个|，代码更简洁易维护。
• 合并数据和转换坐标系用管道串联，减少中间变量。
• 弹窗用公式接口直接引用列名，避免提前生成向量可能导致的行列不匹配问题。

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