嘿，我来帮你搞定用经纬度坐标提取城市名称的事儿！先看你的数据情况——lon和lat居然都是因子（Factor）类型，这可得先处理好类型问题，不然地理编码工具根本认不出这些数值。下面分步骤来操作：

步骤1：将因子类型的经纬度转换为数值型 #

直接用as.numeric()转换因子会返回因子水平的索引值，不是实际的经纬度数值，所以得先转成字符再转数值：

# 转换经纬度类型：因子 → 字符 → 数值
coords$lon <- as.numeric(as.character(coords$lon))
coords$lat <- as.numeric(as.character(coords$lat))

步骤2：逆地理编码获取城市名称 #

这里推荐两个R语言常用的地理编码工具包，按需选择：

方法1：使用tidygeocoder（无需API密钥，推荐） #

这个包基于OpenStreetMap的Nominatim服务，免费且无需额外配置，适合大多数场景：

# 安装并加载包
install.packages("tidygeocoder")
library(tidygeocoder)

# 执行逆地理编码，提取城市和国家信息
coords_with_city <- coords %>%
  reverse_geocode(
    lat = lat, 
    long = lon, 
    method = "osm", 
    address_components = c("city", "country")
  )

运行后，结果框里的city列就是对应的城市名称，country列会返回所属国家，方便你做后续筛选。

方法2：使用ggmap（需Google API密钥） #

如果你需要更精准的地理编码结果，可以用Google的服务，但需要先申请Google Maps API密钥：

# 安装并加载包
install.packages("ggmap")
library(ggmap)

# 注册你的Google API密钥（替换成自己的密钥）
register_google(key = "your_google_api_key")

# 执行逆地理编码
coords_with_city <- coords %>%
  # 拼接经纬度为字符串格式
  mutate(location = paste(lat, lon, sep = ",")) %>%
  # 批量调用地理编码接口
  mutate(geocode_result = purrr::map(location, geocode, output = "more")) %>%
  # 展开嵌套结果
  tidyr::unnest(geocode_result) %>%
  # 保留需要的列（locality对应城市名称）
  select(lon, lat, locality)

注意：Google的API有调用额度限制，超出免费额度会产生费用，使用前要确认好规则。

小提示 #

• 用Nominatim服务时，不要短时间内发送大量请求（官方建议每秒不超过1次），不然会被限流。
• 部分偏远地区的经纬度可能无法返回城市信息，这是地理编码服务的覆盖范围问题，可以尝试调整查询参数或者换用其他服务。

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