R 添加多边形函数

Question

我正在尝试向地图添加一些多边形，但不确定为什么我的代码无法正常工作。我之前在其他 .shp 文件上使用过相同的代码，并且工作正常。我从这里获得了空间文件：http://library.queensu.ca/geo/health_regions。它位于 Ontario->Local Health Integrated Networks -> 2006 -> Boundary Files 下。

  lhin<- readOGR("lhin_divisions/HR035B08.shp",layer = "HR035B08")
  leaflet() %>% addTiles()%>%setView(lng=-81,lat=45,zoom=6)
  %>%addPolygons(data=lhin,stroke=TRUE, weight=1, col='red')

我觉得我错过了一些简单的东西，但我想不通是什么？

生成的地图应具有此处显示的边界：http://library.queensu.ca/madgic/free/health_regions/2005/images/lhin_ont.pdf

Answer 1

这项工作。你可以改变瓷砖的颜色，我只是用了黑色和白色，因为我不知道我在找什么。 shapefile 在 utm 中，没有投影。所以我设置投影并转换成longlat。

proj4string(lhin) <- CRS("+proj=utm +ellps=WGS84 +datum=WGS84")
lhin <- spTransform(lhin, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
leaflet() %>% addProviderTiles("Stamen.Toner") %>% addPolygons(data=lhin,stroke=TRUE, weight=1, color = ~"red")

R 添加多边形函数

R addPolygons function

r

leaflet

sp