R 如何使用.shp将空间多边形应用于传单贴图

我正试图用传单绘制加州各县的地图。我没有收到任何错误，但地图没有产生正确的结果

我在上找到了县线的

.shp

文件。我使用

rgdal

导入了

.shp

文件，并尝试使用

addPolygons（）

函数映射边界

`   # Importing CA County Lines
      ca_counties <- readOGR("Miller/CA_Counties/CA_Counties_TIGER2016.shp")
    # Mapping County Lines
      map1 <- leaflet() %>% 
              addTiles() %>% 
              addPolygons(data=ca_counties)
      map1`

`#导入CA县线
加乌县%
addPolygons（数据=ca_县）
map1`

程序生成的不是县线，而是一条贯穿整个地图的单数线。我已经附上了输出的图像

---

运行代码时，我收到以下警告：

Warning messages:
1: sf layer is not long-lat data 
2: sf layer has inconsistent datum (+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs).
Need '+proj=longlat +datum=WGS84' 

将多边形投影到WGS84适合我。我使用了

sf

包，因为我发现它比

rgdal

更直接

library(sf)
library(leaflet)

ca_counties <- sf::read_sf('CA_Counties_TIGER2016.shp') %>%
  sf::st_transform('+proj=longlat +datum=WGS84')

map1 <- leaflet() %>% 
  addTiles() %>% 
  addPolygons(data=ca_counties)
map1

库（sf）
图书馆（单张）
加乌县%
sf:：st_变换（'+proj=longlat+datum=WGS84'）
map1%
addTiles（）%>%
addPolygons（数据=ca_县）
map1