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

Question

我正在尝试通过传单来绘制加利福尼亚县的地图。 我没有收到任何错误，但是地图未产生正确的结果。

我在https://data.ca.gov/dataset/ca-geographic-boundaries找到了县线的.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`

该程序将生成一条贯穿我整个地图的奇异线，而不是县线。 我已附上输出图像。

失败的地图

Answer 1

运行您的代码时，出现以下警告：

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