ggplot绘制水深Shapefile按MIN_ELEV渐变着色仅显单一颜色，如何解决？

问题描述 #

现有一个sf格式的水深Shapefile，数据结构如下：

str(Ab_bath)
Classes ‘sf’ and 'data.frame':  126 obs. of  13 variables:
 $ CLRRGB  : chr  "2844D9FF" "2651DCFF" "215DDEFF" "2266E0FF" ...
 $ COLOR   : int  85 91 91 91 91 55 98 62 62 62 ...
 $ LAYER   : int  30 30 30 30 30 30 30 30 30 30 ...
 $ L_LIMIT : int  0 0 0 0 0 0 0 0 0 0 ...
 $ MAJ_CAT : chr  "Nautical" "Nautical" "Nautical" "Nautical" ...
 $ MAX_ELEV: num  -14.8 -14.5 -14.2 -14 -13.8 ...
 $ MIN_CAT : chr  "Depth Area" "Depth Area" "Depth Area" "Depth Area" ...
 $ MIN_ELEV: num  -15 -14.8 -14.5 -14.2 -14 ...
 $ MODE    : chr  "D_MR" "D_MR" "D_MR" "D_MR" ...
 $ VALUE   : chr  "14.75" "14.5" "14.25" "14" ...
 $ VALUE2  : chr  "15" "14.75" "14.5" "14.25" ...
 $ U_LIMIT : int  8 8 8 8 8 8 8 8 8 8 ...
 $ geometry:sfc_MULTIPOLYGON of length 126; first list element: List of 2
  ..$ :List of 2
  .. ..$ : num [1:3358, 1:2] 0.887 0.887 0.887 0.887 0.887 ...
  .. ..$ : num [1:2678, 1:2] 0.888 0.888 0.888 0.888 0.888 ...
  ..$ :List of 2
  .. ..$ : num [1:1342, 1:2] 0.886 0.886 0.886 0.886 0.886 ...
  .. ..$ : num [1:1049, 1:2] 0.886 0.886 0.886 0.886 0.886 ...
  ..- attr(*, "class")= chr [1:3] "XY" "MULTIPOLYGON" "sfg"
 - attr(*, "sf_column")= chr "geometry"
 - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA NA NA ...
  ..- attr(*, "names")= chr [1:12] "CLRRGB" "COLOR" "LAYER" "L_LIMIT" ...

单一颜色绘图代码可正常运行：

ggplot() +
  geom_sf(data = Ab_bath_utm, fill = "lightblue", color = "blue", alpha = 0.5)+
  coord_sf(datum = 32631) # 设置UTM坐标系

但按MIN_ELEV字段渐变着色时，图形仅显示单一颜色：

ggplot() +
  geom_sf(data = Ab_bath_utm, aes(fill = MIN_ELEV), color = "black", size = 0.1) +
  scale_fill_gradientn(
    colors = c("navy", "blue", "lightblue", "white"),
    name = "Min Elevation (m)"
  ) +
  coord_sf(datum = 32631) +
  theme_minimal()

解决方法 #

1. 显式指定颜色映射范围 #

ggplot自动判断的数值范围可能存在偏差，导致颜色区分度不足。先确认MIN_ELEV的数值范围，再在色标中显式指定：

# 查看MIN_ELEV的数值范围
range(Ab_bath_utm$MIN_ELEV, na.rm = TRUE)

# 修正后的绘图代码
ggplot() +
  geom_sf(data = Ab_bath_utm, aes(fill = MIN_ELEV), color = "black", size = 0.1) +
  scale_fill_gradientn(
    colors = c("navy", "blue", "lightblue", "white"),
    name = "最小海拔 (m)",
    limits = range(Ab_bath_utm$MIN_ELEV, na.rm = TRUE),
    breaks = seq(floor(min(Ab_bath_utm$MIN_ELEV)), ceiling(max(Ab_bath_utm$MIN_ELEV)), by = 2)
  ) +
  coord_sf(datum = 32631) +
  theme_minimal()

2. 使用更稳定的连续色阶 #

推荐使用viridis色系，它对色盲友好且颜色区分度更高，同时通过direction = -1让更深的水深（更小的MIN_ELEV值）对应更深的颜色：

ggplot() +
  geom_sf(data = Ab_bath_utm, aes(fill = MIN_ELEV), color = "black", size = 0.1) +
  scale_fill_viridis_c(
    name = "最小海拔 (m)",
    option = "viridis",
    direction = -1
  ) +
  coord_sf(datum = 32631) +
  theme_minimal()

关键说明 #

• 确认MIN_ELEV为数值类型（从数据结构看是num，无需转换）
• 显式指定色标范围可避免ggplot自动范围判断的偏差
• 水深数据中，数值越小（越负）代表深度越大，调整色标方向可让颜色逻辑更直观

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