为何st_intersects在if条件未满足时仍被执行四次？

核心原因 #

你遇到的问题本质是函数内部逻辑执行顺序错误，或者条件判断未处理大小写匹配，导致本该跳过的st_intersects调用被强制执行：

• 坐标匹配代码提前执行：如果你的函数把st_intersects相关的坐标匹配逻辑写在了if-else分支之外，无论输入的州名是否匹配缩写/全称，每一次函数调用都会先执行坐标匹配代码，自然触发四次警告。
• 大小写不匹配导致条件失效：state.abb和state.name是固定大小写（前者大写、后者首字母大写），如果输入的州名是小写（比如"tx"），直接用state %in% state.abb会返回FALSE，从而进入坐标匹配分支。

修正后的函数代码 #

将坐标匹配逻辑严格限制在最终的else分支内，并统一处理大小写问题，同时消除平面坐标警告：

library(tidyverse)
library(sf)
library(maps)

state_name_corrector <- function(state, lon, lat) {
  # 统一转换为大写，避免大小写匹配误差
  state_upper <- toupper(state)
  
  # 1. 匹配州缩写
  if (state_upper %in% state.abb) {
    return(state_upper)
  }
  
  # 2. 匹配州全称
  state_names_upper <- toupper(state.name)
  if (state_upper %in% state_names_upper) {
    return(state.abb[match(state_upper, state_names_upper)])
  }
  
  # 3. 仅当以上都不匹配时，执行坐标匹配
  # 构造经纬度点并转换为投影坐标系（避免平面假设警告）
  point <- st_sfc(st_point(c(lon, lat)), crs = 4326) %>%
    st_transform(crs = "EPSG:5070") # 转换为Albers等面积投影
  
  # 加载并转换美国州界数据到同一投影
  state_sf <- maps::map("state", plot = FALSE, fill = TRUE) %>%
    st_as_sf() %>%
    st_set_crs(4326) %>%
    st_transform(crs = "EPSG:5070")
  
  # 执行空间匹配
  intersect_idx <- which(st_intersects(point, state_sf, sparse = FALSE))
  matched_state <- str_to_title(state_sf$ID[intersect_idx])
  
  return(state.abb[match(matched_state, state.name)])
}

测试验证 #

用你的测试数据调用函数，只有无法匹配州名的条目才会执行坐标匹配：

# 测试数据：2个缩写、1个全称、1个未知
test_data <- tibble(
  state = c("CA", "NY", "Texas", "unknown"),
  lon = c(-122.4194, -73.9903, -97.7431, -80.1918),
  lat = c(37.7749, 40.7061, 30.2672, 25.7617)
)

# 使用mapply调用
results <- mapply(
  state_name_corrector,
  test_data$state,
  test_data$lon,
  test_data$lat
)

print(results)
# 输出：CA NY TX FL

关键优化点 #

• 惰性执行坐标匹配：只有在前两个条件都不满足时，才会执行空间匹配代码，避免不必要的计算和警告。
• 统一大小写处理：通过toupper()消除输入大小写带来的匹配误差。
• 消除平面警告：将地理坐标系（WGS84）转换为投影坐标系（EPSG:5070），让st_intersects在平面坐标系下执行，符合其设计逻辑。

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