R语言中map.where函数无法从地理坐标提取国家名的问题求解

问题根源 #

maps包自带的world地图数据库精度较低，部分区域的边界匹配逻辑不完善，导致你的巴西境内坐标无法被正确识别，返回NA值。

可行解决方案 #

方案1：升级maps的地图精度 #

通过mapdata包获取高精度世界地图数据库，提升匹配成功率：

# 安装依赖包
install.packages("mapdata")
library(mapdata)

# 用高精度数据库匹配国家
data$Country <- map.where(database = "worldHires", data$Longitude, data$Latitude)
# 提取结果中的国家名称（部分结果会包含「地区:国家」格式）
data$Country <- sapply(strsplit(data$Country, ":"), function(x) x[1])

方案2：用sf+rnaturalearth实现精准匹配（推荐） #

这两个包提供的官方地理数据精度更高，匹配逻辑更可靠：

# 安装所需包
install.packages(c("sf", "rnaturalearth"))
library(sf)
library(rnaturalearth)

# 获取标准世界国家边界数据
world_boundaries <- ne_countries(returnclass = "sf")

# 将你的坐标转换为空间对象（采用WGS84坐标系，EPSG:4326）
coords_sf <- st_as_sf(data, coords = c("Longitude", "Latitude"), crs = 4326)

# 空间关联匹配国家信息
matched_result <- st_join(coords_sf, world_boundaries, join = st_within)
# 提取国家全称列到原数据集
data$Country <- matched_result$name_long

方案3：在线反向地理编码（需API密钥） #

如果需要绝对精准的结果，可以调用Google地图的反向地理编码API：

install.packages("ggmap")
library(ggmap)

# 需先在Google Cloud平台申请API密钥并完成注册
# register_google(key = "你的Google Maps API密钥")

# 批量解析坐标获取国家名
data$Country <- apply(data, 1, function(row) {
  geo_res <- revgeocode(c(row["Longitude"], row["Latitude"]), output = "more")
  geo_res$country
})

测试你的坐标：所有给定坐标均位于巴西（Brazil），使用上述任一方案都能正确返回国家名称。

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