如何在R中添加海深数据并按区域填充多边形图层

核心思路 #

以ggplot2为核心框架，先绘制连续色阶的海深底图，再叠加离散色阶的多边形区域图层，通过明确区分填充映射类型、调整图层顺序和透明度，解决两类图层的颜色冲突问题，同时避免逐个添加区域的繁琐操作。

针对不同工具的问题解决与实现 #

1. 使用marmap实现（解决「离散值传入连续填充刻度」错误） #

该错误源于海深底图的连续填充与多边形的离散填充共用了fill映射刻度，需为两类图层分别指定独立的颜色映射：

步骤代码

# 加载依赖包
library(ggplot2)
library(sf)
library(marmap)
library(dplyr)

# 读取已处理好的sf区域数据（CRS=4326）
region_polygons <- st_read("your_region_data.shp") # 替换为你的数据路径

# 获取目标区域的海深数据（示例范围：西经75-60，北纬35-45）
bathy <- getNOAA.bathy(lon1 = -75, lon2 = -60, lat1 = 35, lat2 = 45, resolution = 1)
# 转换为ggplot可用的数据框
bathy_df <- fortify.bathy(bathy)

# 绘制地图
ggplot() +
  # 第一步：绘制海深底图（连续填充）
  geom_raster(data = bathy_df, aes(x = x, y = y, fill = z)) +
  # 设置海深的连续色阶（自定义蓝系渐变）
  scale_fill_continuous(name = "海深（m）",
                        low = "lightblue", high = "darkblue",
                        limits = c(-6000, 0), na.value = "white") +
  # 第二步：叠加区域多边形（离散填充）
  geom_sf(data = region_polygons, aes(fill = Region), alpha = 0.3, color = "black") +
  # 设置区域的自定义离散颜色（匹配MAB/SNE/GB/GOM四类）
  scale_fill_manual(name = "海区",
                    values = c(MAB = "#FF9933", SNE = "#33CC33", GB = "#3366FF", GOM = "#FF33CC")) +
  # 强制使用WGS84坐标
  coord_sf(crs = 4326) +
  # 调整主题
  theme_bw() +
  theme(panel.grid = element_blank())

关键说明

• 明确将海深的fill映射绑定到连续色阶scale_fill_continuous，区域的fill绑定到离散色阶scale_fill_manual，避免刻度冲突。
• 通过alpha参数降低区域填充的透明度，保证海深底图可见。

2. 使用ggOceanMaps实现（解决颜色不足/海深被覆盖问题） #

ggOceanMaps的底图默认会覆盖上层图层，需调整图层顺序并明确自定义颜色数量：

步骤代码

# 加载依赖包
library(ggplot2)
library(sf)
library(ggOceanMaps)

# 读取sf区域数据
region_polygons <- st_read("your_region_data.shp")

# 绘制地图
basemap(limits = c(-75, -60, 35, 45), bathymetry = TRUE) + # 加载海深底图
  # 叠加区域多边形（放在底图之后，避免被覆盖）
  geom_sf(data = region_polygons, aes(fill = Region), alpha = 0.4, color = "black") +
  # 自定义区域颜色（必须匹配所有Region类别）
  scale_fill_manual(name = "海区",
                    values = c(MAB = "#F0E442", SNE = "#009E73", GB = "#0072B2", GOM = "#D55E00")) +
  # 关闭底图的默认填充刻度（避免与区域填充冲突）
  guides(fill = guide_legend(override.aes = list(alpha = 1)))

关键说明

• basemap()先加载海深底图，再用geom_sf叠加区域，确保区域在顶层。
• 自定义颜色时必须包含所有Region的取值（MAB/SNE/GB/GOM），避免「颜色数量不足」错误。
• 用guides()调整图例，消除alpha对图例显示的影响。

通用优化建议 #

• 避免逐个添加区域图层：直接使用geom_sf(aes(fill=Region))即可一次性渲染所有区域，无需按Region拆分数据，更易维护。
• 颜色区分度：选择与海深蓝系差异大的区域颜色（如暖色调），提升可读性。
• 坐标一致性：确保海深数据与sf数据的CRS统一为4326，避免偏移问题。

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