使用R+ggplot结合st_crop裁剪地图未达预期效果的技术求助

首先修正你代码中的一个小错误：st_transform 在管道中的调用方式不对，应该直接传入CRS参数，不需要重复指定数据对象：

lapply(c("rnaturalearthdata", "sf", "ggplot2","dplyr"), require, character.only = TRUE)

# 修正st_transform的调用
world <- ne_coastline(scale = "small", returnclass = "sf") %>%
  st_transform(crs = 4326)
world_cropped <- st_crop(world, xmin = -30, xmax = 90, ymin = 30, ymax = 82)

问题原因 #

st_crop 仅裁剪了地图要素，但ggplot默认会根据裁剪后要素的实际边界自动调整绘图范围，而非你指定的完整方形区域。比如如果裁剪后的海岸线没有延伸到你设定的xmin/xmax/ymin/ymax边界，ggplot就会缩到要素的实际范围，导致无法呈现方形。

替代解决方案（无需用coord_sf设置范围） #

方案1：添加隐形边界框强制绘图范围 #

创建一个与目标范围完全一致的空sf多边形图层，将其添加到绘图中（设置为完全透明），让ggplot识别这个完整的方形边界作为绘图范围：

# 创建目标范围的边界框sf对象
bbox_sf <- st_bbox(c(xmin = -30, xmax = 90, ymin = 30, ymax = 82)) %>%
  st_as_sfc()

# 绘制裁剪后的地图，同时加入隐形边界框
ggplot() +
  geom_sf(data = bbox_sf, fill = NA, color = NA) +  # 隐形边界框，强制范围
  geom_sf(data = world_cropped, fill = NA, color = "black")

方案2：手动构建ggplot的刻度范围 #

利用ggplot的scale_x_continuous和scale_y_continuous手动指定坐标轴范围，同时关闭自动扩展（expand = c(0,0)），确保范围严格匹配你的方形区域：

ggplot() +
  geom_sf(data = world_cropped, fill = NA, color = "black") +
  scale_x_continuous(limits = c(-30, 90), expand = c(0, 0)) +
  scale_y_continuous(limits = c(30, 82), expand = c(0, 0))

注意：这种方法需要确保CRS是地理坐标系（如4326），因为在投影坐标系下，经纬度刻度的设置逻辑不同。

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