如何在 R 中的 map 上插入刻度

Question

我想知道是否可以生成等比例的附加图像，以红色突出显示

library(terra)
library(geodata)
brz <- geodata::gadm("Brazil", level=2, path=".")
mun <- brz[brz$NAME_2 == "Castro", ]
plot(mun, main="", axes=FALSE)

例子：

Answer 1

您可以使用tm_scale_bar中的tmap将比例尺添加到 map。 此外，您可能需要sf将mun从SpatVector转换为sf （tmap 可用的格式）。

library(tmap)
library(sf)

# Transform mun to sf
tm_shape(sf::st_as_sf(mun)) +
  # Draw as polygon with only borders
  tm_borders() +
  # Add scale bar
  tm_scale_bar(position = c("LEFT", "BOTTOM"))

Answer 2

另一种选择是使用库ggspatial和tidyterra 。

library(ggspatial)
library(tidyterra)

ggplot() +
  geom_spatvector(data = mun)  +
  annotation_scale() + # add scale
  annotation_north_arrow(location = "br", which_north = "true") + # add arrow
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    panel.border = element_blank(),
    axis.text.y=element_blank(),  #remove y axis labels
    axis.ticks.y=element_blank(),  #remove y axis ticks
    axis.text.x=element_blank(),  #remove x axis labels
    axis.ticks.x=element_blank()  #remove x axis ticks
  )