在R中绘制切萨皮克湾地图无输出结果的技术求助

咱们一步步拆解你的代码为啥没输出预期的切萨皮克湾地图，然后把问题解决掉：

核心问题分析 #

1. 数据源完全不对 #

你用ne_countries()获取的是全球国家边界数据，这个数据集里根本没有切萨皮克湾的水域信息。你设置的坐标范围只是框出了美国东海岸的陆地区域，自然看不到海湾的轮廓。

2. 图像输出触发问题 #

你把ggplot对象赋值给了plot变量，但如果是在R脚本中运行代码，没有主动调用print(plot)或者直接输出这个对象的话，R不会自动显示图像（RStudio交互式窗口可能会自动显示，但脚本模式下必须手动触发）。

修正后的完整可运行代码 #

下面是调整后的代码，能正确显示包含陆地和水域的切萨皮克湾区域：

library("ggplot2")
library("sf")
library("rnaturalearth")
library("rnaturalearthdata")
library("ggspatial")

# 获取美国各州边界数据（覆盖切萨皮克湾所在区域）
us_states <- ne_states(country = "united states of america", returnclass = "sf")
# 获取全球海洋数据（用于显示海湾水域）
oceans <- ne_download(scale = "medium", type = "ocean", returnclass = "sf")

# 绘制地图
plot <- ggplot() +
  # 先画海洋图层（浅蓝色底色）
  geom_sf(data = oceans, fill = "lightblue") +
  # 再画美国各州陆地图层（深绿色）
  geom_sf(data = us_states, fill = "darkgreen") +
  # 锁定切萨皮克湾的坐标范围，关闭多余空白扩展
  coord_sf(xlim = c(-77.5, -75), ylim = c(37, 40), expand = FALSE) +
  # 添加坐标轴标签和标题
  xlab("Longitude") +
  ylab("Latitude") +
  ggtitle("Chesapeake Bay") +
  # 添加海湾文本标注
  annotate(geom = "text", x = -76.1, y = 37.8, label = "Chesapeake Bay",
           color = "brown", size = 3, angle = 90) +
  # 添加指北针
  annotation_north_arrow(location = "tl", pad_x = unit(0.5, "cm"),
                         pad_y = unit(1, "cm"), height = unit(1,"cm"),
                         width = unit(0.5,"cm")) +
  # 设置主题样式
  theme(panel.background = element_rect(fill = "lightblue"),
        panel.grid.major = element_line(linetype = "dashed", size = 0.2))

# 脚本模式下必须手动触发输出，交互式环境可省略
print(plot)

关键调整说明 #

• 替换数据源：用ne_states()获取美国陆地边界，ne_download()获取海洋数据，这样才能呈现出海湾的水域轮廓；
• 调整图层顺序：先画海洋再画陆地，确保陆地在水域上方正常显示；
• 关闭坐标扩展：把expand = TRUE改为FALSE，避免坐标范围外出现多余空白；
• 明确触发输出：脚本中通过print(plot)强制生成图像，解决无输出问题。

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