在R中使用ggplot和ggmap在地图上标注企鹅位置遇问题

Hey there! 我来帮你排查下用ggplot和ggmap绘制南极洲企鹅分布时，叠加位置点遇到的「数据被移除+报错」问题～结合常见的踩坑点，给你几个排查方向：

1. 确认数据坐标系和地图匹配 #

ggmap默认返回的地图是WGS84（EPSG:4326）坐标系，也就是常规的经纬度范围（经度-180到180，纬度-90到90）。如果你的企鹅位置数据用的是其他坐标系（比如南极专属的EPSG:3031投影），就会出现点显示不出来或被自动过滤的情况。

可以用sf包转成匹配的坐标系：

library(sf)
# 假设你的企鹅数据框是penguin_data，包含lon（经度）和lat（纬度）列
# 先把数据转成sf对象，指定原数据的坐标系（示例为EPSG:3031）
penguin_sf <- st_as_sf(penguin_data, coords = c("lon", "lat"), crs = 3031)
# 转换成和地图一致的WGS84坐标系
penguin_sf <- st_transform(penguin_sf, crs = 4326)
# 叠加到地图上
ggplot() +
  ggmap(你的地图对象) +
  geom_sf(data = penguin_sf, color = "darkorange", size = 2, inherit.aes = FALSE)

2. 检查图层语法和数据类型 #

有时候是语法细节出问题：

• 确保geom_point()里正确指定数据和映射，不要把企鹅数据放在ggplot()的全局data参数里（ggmap的图层会覆盖全局数据），正确写法示例：

# 先获取南极洲地图（用Stamen地图，无需API密钥）
antarctica_map <- get_map(
  location = c(lon = 0, lat = -90), 
  zoom = 2, 
  maptype = "terrain-background"
)
# 叠加企鹅点
ggmap(antarctica_map) +
  geom_point(
    data = penguin_data, 
    aes(x = lon, y = lat), 
    color = "navy", 
    size = 3,
    alpha = 0.7
  )

• 检查企鹅数据的经纬度是不是数值型：用str(penguin_data)查看，如果是字符型，先转成数值：

penguin_data$lon <- as.numeric(penguin_data$lon)
penguin_data$lat <- as.numeric(penguin_data$lat)

3. 清理数据里的NA值 #

如果你的企鹅数据中存在经纬度为NA的行，geom_point()会自动移除这些数据，可能你看到的「数据被移除」就是这个原因。可以先清理NA：

clean_penguin_data <- na.omit(penguin_data[, c("lon", "lat")])
# 用清理后的数据绘图
ggmap(antarctica_map) +
  geom_point(data = clean_penguin_data, aes(x = lon, y = lat))

4. 检查包版本兼容性 #

旧版本的ggmap和ggplot2可能存在图层叠加的bug，建议更新到最新版：

install.packages(c("ggplot2", "ggmap"))

如果还是解决不了，建议把具体的报错信息贴出来，或者用head(penguin_data)展示数据的前几行，这样能更精准定位问题～

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