向已使用 coord_sf 缩放的 ggplot 地图添加比例尺？

Question

我使用 sf 包和 ggplot2 生成了一个地图：

library(ggplot2)
library(sf)
library(rnaturalearth)
state_prov <- rnaturalearth::ne_states(c("united states of america", "canada"), returnclass="sf")
x <- ggplot(data=state_prov) + 
geom_sf()+
coord_sf(xlim=c(-170, -95), ylim=c(40, 75)) 
print(x)

它在 Rstudio 中生成以下地图：

太棒了，但我需要添加一个比例尺。当我尝试使用 ggsn 修改代码时，我根本看不到比例尺。

library(ggplot2)
library(sf)
library(rnaturalearth)
state_prov <- rnaturalearth::ne_states(c("united states of america", "canada"), returnclass="sf")
x <- ggplot(data=state_prov) + 
geom_sf()+
coord_sf(xlim=c(-170, -95), ylim=c(40, 75)) +
ggsn::scalebar(state_prov, location="topleft", dist = 50, dist_unit = "km", 
                 transform=TRUE, model="WGS84", height=0.1)
print(x) 

我尝试更改高度、距离和位置，但没有成功。当我删除对 coord_sf() 的调用时，我可以看到比例尺缩放得很差，这让我相信 ggsn 无法识别地图正在由 coord_sf() 放大。

我该如何解决？ggsn 似乎不容易修改。我愿意使用另一个包或方法，但我确实需要继续以类似的方式调用 ggplot，因为我有一个基于相同结构的更复杂的地图。

谢谢！

Answer 1

正如您所提到的，如果您注释掉coord_sf代码的一部分，则会显示比例尺。我的猜测是ggsn::scalebar必须topleft从整个state_prov数据集中获取它的位置，并且当您使用coord_sf比例尺进行缩放时，比例尺会被裁剪掉。

编辑：当将比例尺放置在具有此比例的纬度/经度投影的地图上时，请注意极端扭曲：/sf/answers/2896149861/

以下是显示比例尺的几个选项。

选项1

使用ggspatial::annotation_scale而不是ggsnwhich 似乎可以识别 中定义的缩放coord_sf。

ggplot(data=state_prov) + 
  geom_sf()+
  coord_sf(xlim=c(-170, -95), ylim=c(40, 75)) +
  ggspatial::annotation_scale(location = 'tl')

选项2

使用您的原始代码，但state_prov在绘图之前进行裁剪，以便scalebar找到正确的topleft.

state_prov_crop <- st_crop(state_prov, xmin=-170, xmax = -95, ymin = 40, ymax = 75)

ggplot(data=state_prov_crop) + 
  geom_sf()+
  #coord_sf(xlim=c(-170, -95), ylim=c(40, 75)) +
  ggsn::scalebar(state_prov_crop, location="topleft", dist = 50, dist_unit = "km", 
                 transform=TRUE, model="WGS84", height=0.1)