R addPolygons函数

Question

我正在尝试将一些多边形添加到地图中，但我不确定为什么我的代码无效。我以前在其他.shp文件中使用过相同的代码，它运行正常。我从这里获得了空间文件：http://library.queensu.ca/geo/health_regions。它位于Ontario-＆gt; Local Health Integrated Networks - ＆gt; 2006 - ＆gt;边界文件。

  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是没有投影的。所以我设置了投影并转换成了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")